Entering edit mode
8.6 years ago
William
★
5.2k
Is there a tool already out there that can add the QualityOfDepth(QD) annotation to a VCF file created by freebayes?
The general idea of variant QualityByDepth(QD) is explained here and the detailed formula is the best answer (highlighted in green) here.
Freebayes outputs somewhat different fields than GATK. Instead of AC (allele count) per genotype:
##FORMAT=<ID=RO,Number=1,Type=Integer,Description="Reference allele observation count">
##FORMAT=<ID=AO,Number=A,Type=Integer,Description="Alternate allele observation count">
The basic QD formula should be something like.
QD = QUAL / Sum (Sum (RO+AO)) (over just the samples that have a read that contributes to the variant call).