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

Read Depth in VCF #305

Closed
mproberts99 opened this issue May 9, 2024 · 4 comments
Closed

Read Depth in VCF #305

mproberts99 opened this issue May 9, 2024 · 4 comments
Assignees
Labels
bug Something isn't working

Comments

@mproberts99
Copy link

mproberts99 commented May 9, 2024

Hi,
I found a typo in the output VCF that is causing errors when parsing. The AF parameter has NUMBER=1, so it must have exactly 1 entry, but with the recent update of adding af for all alternate alleles, it has more than one value so it should be NUMBER=G I believe.
##FORMAT=<ID=AF,Number=G,Type=Float,Description="Observed allele frequency in reads, for each ALT allele, in the same order as listed, or the REF allele for a RefCall">
Also, because we use amplicon sequencing, we have to lower the MQ filter to 0 so having the DP automatically display reads with MQ>20 rather than what was specified in the command may be misleading?

@aquaskyline aquaskyline added the bug Something isn't working label May 10, 2024
@zhengzhenxian
Copy link
Collaborator

@mproberts99

Thanks for reporting this! The number of AF has been updated in the header, and we also revised the DP description in the header for clarification. Pls kindly re-pull the latest docker image, thanks!

Zhenxian

@mproberts99
Copy link
Author

Thanks so much!

@mproberts99
Copy link
Author

Ah sorry @zhengzhenxian , it looks like it should actually be NUMBER=R. I got this error:
Incorrect number of FORMAT/AF values at chr5:70070208, cannot merge. The tag is defined as Number=G, but found
2 values and 3 alleles. See also http://samtools.github.io/bcftools/howtos/FAQ.html#incorrect-nfields

@mproberts99 mproberts99 reopened this May 16, 2024
@zhengzhenxian
Copy link
Collaborator

@mproberts99

It seems assigning ‘A’ would be more flexible, I tested with bcftools merge and it worked well. I have updated the docker image and pls kindly re-pull, thanks!

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
bug Something isn't working
Projects
None yet
Development

No branches or pull requests

3 participants