In Vcf file, generate by samtools mpileup, there is an AD field (Allelic depth). This field is suppose to replace the previous DP4 field (Number of high-quality ref-forward, ref-reverse, alt-forward and alt-reverse bases).

But when I use -t DP4, AD in mpileup, I have a different number of values in the 2 fields:

    GT:PL:DP:DP4:AD 0/1:27,0,107:5:4,0,1,0:4,1,0

And I don't understand why one data is missing in the AD field... Is it not supposed to be like the DP4, so:

> high-quality ref-forward, ref-reverse, alt-forward and alt-reverse bases

?

After that I am not sure: what there are 3 values in the AD field ? And what are the differences between AD and DP4?

Thanks for your help.

Rob
I think DP4 shows ref and all alt reads, but AD further separates the different alt alleles (a single record can have multiple alt alleles).

If you check the AD and DP4 descriptions in the VCF header, it should explain both fields.

I actually don't have AD in my samtools 1.3 VCF. They may have taken it out or it might depend on some other setting. Thus, I can't really confirm.