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

Where can I find "Phred-scaled quality score" (QUAL)? #1506

Closed
majkiw opened this issue Apr 27, 2017 · 2 comments
Closed

Where can I find "Phred-scaled quality score" (QUAL)? #1506

majkiw opened this issue Apr 27, 2017 · 2 comments
Milestone

Comments

@majkiw
Copy link
Contributor

@majkiw majkiw commented Apr 27, 2017

Hey,
Following instructions in https://github.com/bigdatagenomics/adam/wiki/FAQ-(Frequently-Asked-Questions) and http://bdgenomics.org/mail/ I created a question in Google Groups https://groups.google.com/forum/#!topic/adam-developers/WNu5a8wIwz8
Since I don't think anyone looks there anymore it I will ask it here again.

When reading VCF I cannot find "Phred-scaled quality score" (QUAL) and even htsjdkVC.getPhredScaledQual doesn't seem to be accessed anywhere.
I also don't see the values when printing the whole Genotype object.
Is the QUAL field omitted on purpose or am I missing something obvious?

Thank you,
Michał Wysocki

@majkiw majkiw changed the title Where can I find "Phred-scaled quality score" (QUAL) Where can I find "Phred-scaled quality score" (QUAL)? Apr 27, 2017
@heuermh
Copy link
Member

@heuermh heuermh commented Apr 27, 2017

Sorry for not seeing your post on the mailing list!

In the VCF specification, QUAL is defined

6. QUAL - quality: Phred-scaled quality score for the assertion made
in ALT. i.e. −10log10 prob(call in ALT is wrong). If ALT is ‘.’ (no variant) then
this is −10log10 prob(variant), and if ALT is not ‘.’ this is −10log10 prob(no variant).
If unknown, the missing value should be specified. (Float)

We use a split-allelic model, where each ALT at a site is a separate Variant. There is no way to distribute this per-site quality score over all the separate variants so we have to drop it.

@majkiw
Copy link
Contributor Author

@majkiw majkiw commented Apr 28, 2017

Thank you for the quick response.
I read the definition before but I didn't realize split-allelic approach breaks it.

@majkiw majkiw closed this Apr 28, 2017
@heuermh heuermh modified the milestone: 0.23.0 Jul 22, 2017
@heuermh heuermh added this to Completed in Release 0.23.0 Jan 4, 2018
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
Linked pull requests

Successfully merging a pull request may close this issue.

None yet
2 participants
You can’t perform that action at this time.