Question: PISCES DP is greater than AD
0
gravatar for rosefroehlich
5 months ago by
rosefroehlich0 wrote:

Hallo,

I am using Pisces for variant calling in my TST170 data and have a question regarding the AD and DP values in the vcf. Why is the DP value greater than the sum of the two given AD values for the Alt and Reference Allele? Reading the following on Github made me think that the DP-AD(Alt)=AD(Ref): "Coverage at a particular location is calculated as the sum of all called reads which pass filters (have sufficient mapping quality and base call quality). Coverage is given as DP in the vcf file. Deleted bases have counted towards coverage since version 3.5.5 . SNV count uses the same filters as the coverage count filters. Indel/MNV count is the sum of all reads passing filters where we were able to determine an indel/MNV. The estimate reference count is the coverage minus the variant count. (So, for an indel or MNV, the "reference count" is really the total number of passing reads that did not have the given variant). Variant frequency is calculated by dividing the variant count by the coverage at the given location*.“ https://github.com/Illumina/Pisces/wiki

I would be very thankful for you help.

Best regards, Rose

ADD COMMENTlink written 5 months ago by rosefroehlich0
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 957 users visited in the last hour