bcftools icon indicating copy to clipboard operation
bcftools copied to clipboard

INDEL fields in the mpileup/call/norm output report raw reads

Open habibdanesh opened this issue 3 years ago • 1 comments

For the cases of INDEL, I understand that INFO=IDV/IMF fields report the number/fraction of "raw" reads, as explained in the VCF header. But I was expecting the FORMAT=DP/AD/ADF/ADR to be based on the "high-quality" reads--is there a specific reason that they are not?

habibdanesh avatar May 15 '21 05:05 habibdanesh

The FORMAT=DP/AD/ADF/ADR tags are supposed to be high-quality bases, as described in the header. Your question suggests that you are observing something else. In such case we'll need mode details and a small test case to reproduce the problem.

pd3 avatar May 19 '21 12:05 pd3