was hoping to get some advice on working with the AD field for variants in VCFs produced by Pindel and Pindel2VCF.
My lab has a variant calling pipeline that uses GATK's UnifiedGenotyper for calling SNVs and small Indels, as well as using Pindel to call deletions and short insertions. As part of our VCF filtering, we remove variants with an Variant Allele Fraction below a certain threshold. However, it has been brought to my attention by other members of the lab that when Pindel and UnifiedGenotyper call the same variant, they have different AD values and so different VAF values (I'm not sure if it's relevant, but the AD for UnifiedGenotyper is the same as found in IGV).
Normally the difference is small enough that it doesn't affect our filtering, but there have been some cases where a variant called by both branches of the pipeline has an VAF greater than our cutoff threshold in one VCF, but has an VAF lower than threshold in the other VCF. To avoid this issue, which value for AD should I be using for calculating VAF? The AD from UnifiedGenotyper, or the AD from Pindel?
Thanks in advance for any suggestions.
I've also posted this question on SeqAnswers (http://seqanswers.com/forums/showthread.php?p=208171#post208171)