Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

recombinant HCV genome not detected as such #35

Open
sschmutz opened this issue Nov 30, 2018 · 1 comment
Open

recombinant HCV genome not detected as such #35

sschmutz opened this issue Nov 30, 2018 · 1 comment
Assignees

Comments

@sschmutz
Copy link
Member

Example 1000421130

despite there's a recombinant sequence (in case of HCV just one, HCV-1b/2k) which is very close to the actual sequence (created with SmaltAlign), MinVar chose a non recombinant sequence as reference:
in minvar.log INFO 2018/11/29 09:13:23 prepare.py: main() 472: Using non recombinant

This then leads to a cns_ambiguous.fasta sequence which contains a lot of Ns
The recombinant form is also not listed in subtype_evidence.csv

@mihuber mihuber assigned mihuber and MaryamZaheri and unassigned mihuber Dec 3, 2018
@MaryamZaheri
Copy link
Member

The logic in MinVar is that if the viral reads are aligned to a non_recombinant reference more than 50% then do not try the recombinant reference sequences.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

3 participants