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 requires a method filtersWereApplied #741

Open
heuermh opened this issue Nov 14, 2016 · 5 comments
Open

Genotype requires a method filtersWereApplied #741

heuermh opened this issue Nov 14, 2016 · 5 comments

Comments

@heuermh
Copy link
Contributor

heuermh commented Nov 14, 2016

Where do the assertions in #663 come from?

"A Genotype always has filters. If there are no non-PASS filters, it's PASS."

Per VCF specification:
"FT : sample genotype filter indicating if this genotype was “called” (similar in concept to the FILTER field). Again, use PASS to indicate that all filters have been passed, a semi-colon separated list of codes for filters that fail, or ‘.’ to indicate that filters have not been applied. These values should be described in the meta- information in the same way as FILTERs (String, no white-space or semi-colons permitted)"

Genotype requires a method boolean filtersWereApplied() similar to that in VariantContext to support this distinction.

@fnothaft
Copy link
Contributor

+1, subscribing

@yfarjoun
Copy link
Contributor

I guess the assertion was a way to say that you should always be able to
ask about a genotype whether it is filtered or not.

do you think that there should also be a "wereFiltersApplied()" function
that would return false if the FT tag is missing or "."?

On Mon, Nov 14, 2016 at 4:46 PM, Frank Austin Nothaft <
notifications@github.com> wrote:

+1, subscribing


You are receiving this because you are subscribed to this thread.
Reply to this email directly, view it on GitHub
#741 (comment),
or mute the thread
https://github.com/notifications/unsubscribe-auth/ACnk0rBwoRNUDZDi101YNUbKqu1O6knlks5q-NakgaJpZM4KxqYf
.

@heuermh
Copy link
Contributor Author

heuermh commented Nov 15, 2016

do you think that there should also be a "wereFiltersApplied()" function that would return false if the FT tag is missing or "."?

Yes, that would be my suggestion. Because of the special treatment FT receives as one of PRIMARY_KEYS (https://github.com/samtools/htsjdk/blob/master/src/main/java/htsjdk/variant/variantcontext/Genotype.java#L54) the actual attribute value appears not to be accessible.

@heuermh
Copy link
Contributor Author

heuermh commented Apr 19, 2017

This one has been sitting for a while, any update?

@yfarjoun
Copy link
Contributor

I think that it's a reasonable idea. Genotype filters haven't loomed large in our world...(though this might soon change)...So I'm not sure that we would prioritize submitting the PR...but we will review it if you have a suggestion.

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