This is a test version of Biostars. For the public version, visit https://www.biostars.org.
how to add QualityOfDepth(QD) annotation to freebayes vcf

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).

vcf qualitybydepth freebayes

1 answer

I believe GATK Variant Annotator will do that for you.

Log in to answer this question.