-
Notifications
You must be signed in to change notification settings - Fork 0
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
📚 update docs for genotype filtering #57
base: dm-exome-panel-genotyping
Are you sure you want to change the base?
Conversation
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Looks pretty good overall, a couple small questions and comments
docs/GATK_GERMLINE_README.md
Outdated
small-scale experiments, such as targeted gene panels or exome studies with | ||
fewer than 30 exomes." Therefore, VQSR is only activated in this workflow when | ||
the input gVCFs for this workflow come from whole genome sequencing experiments | ||
or when the user provides 30 or more exome gVCFs. |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
If you provide 30+ exomes, target, do you get a joint called VCF, or individual calls with VQSR?
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
They all get jointly processed. 30+ exomes will go to VQSR, 30+ targeted will still be hard filtered.
Targeted seq just is too high depth and too low data to work properly in VQSR.
Hard Filtering is really only constrained by having sufficient depth. In the | ||
case of exome and targeted sequencing, the depths are more than sufficient. Our | ||
current approach for hard filtering mirrors the default approach outlined in | ||
the GATK documentation. However as they point out, "You absolutely SHOULD |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Is there link to the various filtering possibilities/syntax?
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
There are but they are GATK doc links, which means they'll just end up dead sooner or later.
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
I'm debating making archive links to everything, but in the ticket.
bea6e65
to
60db58d
Compare
Description
Updating the docs and pub app formatting for the genotyping workflow.
Part of https://d3b.atlassian.net/browse/BIXU-3800
Type of change
How Has This Been Tested?
Please describe the tests that you ran to verify your changes. Provide instructions so we can reproduce. Please also list any relevant details for your test configuration
Test Configuration:
Checklist: