Tweet
In a recent paper entitled "Exome sequencing as a tool for Mendelian disease gene discovery", the authors write:
"The analysis pipeline consisted of: alignment using the Burrows–Wheeler alignment tool; recalibration; realignment around insertion–deletions and merging with the Genome Analysis Toolkit (GATK)91; and removal of duplicates with PICARD. Variants were called using the following parameters: quality score > 50, allele balance ratio < 0.75; homopolymer run > 3; and quality by depth < 8."
I have three questions:
1. Once the removal of duplicates step/PICARD step is done, you will be left with a bam (or sam) file. So what subsequent steps are taken in order to find the four numbers that they're filtering on?
2. If I run things through the GATK pipeline up through the UnifiedGenotyper step, I end up with a vcf file. This site describes elements in that file:
Is "allele balance ratio" the same as "##INFO=<ID=AF,Number=.,Type=Float,Description="Allele Frequency, for each ALT allele, in the same order as listed">"?
3. The only mention I see of "quality by depth" is here:
##INFO=<ID=QD,Number=1,Type=Float,Description="Variant Confidence/Quality by Depth">
I assume that "Variant Confidence" is the same as "QUAL", in which case I could calculate "Quality by Depth", but this seems wrong because they want LOW "quality by depth" scores. Do they instead mean this QD ratio when they say "quality by depth"?
Thank you.
Eric
"The analysis pipeline consisted of: alignment using the Burrows–Wheeler alignment tool; recalibration; realignment around insertion–deletions and merging with the Genome Analysis Toolkit (GATK)91; and removal of duplicates with PICARD. Variants were called using the following parameters: quality score > 50, allele balance ratio < 0.75; homopolymer run > 3; and quality by depth < 8."
I have three questions:
1. Once the removal of duplicates step/PICARD step is done, you will be left with a bam (or sam) file. So what subsequent steps are taken in order to find the four numbers that they're filtering on?
2. If I run things through the GATK pipeline up through the UnifiedGenotyper step, I end up with a vcf file. This site describes elements in that file:
Is "allele balance ratio" the same as "##INFO=<ID=AF,Number=.,Type=Float,Description="Allele Frequency, for each ALT allele, in the same order as listed">"?
3. The only mention I see of "quality by depth" is here:
##INFO=<ID=QD,Number=1,Type=Float,Description="Variant Confidence/Quality by Depth">
I assume that "Variant Confidence" is the same as "QUAL", in which case I could calculate "Quality by Depth", but this seems wrong because they want LOW "quality by depth" scores. Do they instead mean this QD ratio when they say "quality by depth"?
Thank you.
Eric
Comment