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

Genotype quality #131

Open
MikeWLloyd opened this issue Dec 14, 2023 · 2 comments
Open

Genotype quality #131

MikeWLloyd opened this issue Dec 14, 2023 · 2 comments

Comments

@MikeWLloyd
Copy link

Can you clarify the following: indel_header.addFormatField("GQ", "1","String", "Genotype quality (currently not supported. Always 0)");

There are non-zero GQ values in the final VCFs. I am confused by the header statement.

@cmnbroad
Copy link

Also, the 4.x versions of the vcf spec all state that GQ is an integer. Having said that, according to github, the particular line referenced above is 8 years old, and its entirely possible that the early 4.x specs incorrectly stated GQ was string (and there is STILL an incorrect example in the 4.3 spec that uses String as well).

@MikeWLloyd
Copy link
Author

For what it is worth, 'PL' is also out of spec.

##FORMAT=<ID=PL,Number=.,Type=Float,Description="Normalized likelihood of the current genotype">

I believe this should be 'Integer' as well.

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

2 participants