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

Alternate allele length of zero #30

Open
tfwillems opened this issue Mar 18, 2014 · 1 comment
Open

Alternate allele length of zero #30

tfwillems opened this issue Mar 18, 2014 · 1 comment
Labels

Comments

@tfwillems
Copy link
Contributor

When an alternate allele has a length of 0 and there's only 1 alternate allele, the alternate allele field is printed in the VCF as the empty string. While this is technically not incorrect, it causes awk to think there's one less field than every other line in the VCF file. I suspect this will break a lot of subsequent tools, and we may want to consider requiring that alleles have a length >= 1

@tfwillems tfwillems added the bug label Mar 18, 2014
@mgymrek
Copy link
Owner

mgymrek commented Mar 18, 2014

Thanks for point this out. While technically we could have an allele of size 0, I think in almost all such cases it would be a calling error on our part.

Easiest option is to change MIN_ALLELE_SIZE to 1 in ReadContainer.cpp to not allow length 0 alleles.

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

No branches or pull requests

2 participants