Question: Depth of coverage in Bam file does not match depth in VCF file
0
gravatar for jsneaththompson
18 months ago by
jsneaththompson60 wrote:

My lab has a variant calling pipeline where sorted bam files are used as input for Pindel. The pindel output files (D and SI) are then converted to vcfs for annotation.

However, when the sorted bam files are loaded into IGV the loaded coverage file gives one value for depth at a variant's location which is not the same as the value obtained from the vcf files (AD + RD).

I thought that maybe the depth shown in the vcf files might be an average of coverage over the length of the indel, but I can't see that that's the case. Is there an issue here, or are IGV coverage and VCF coverage not directly comparable? Would appreciate any advice; below is the code that generates the first vcf files from the sorted bam files:

for k in `cat samplesPindel.txt`; do
    ls ../sorted_bam/${k}.sorted.bam > ${k}_pindelinput.txt
sed -i 's|\(.\+\)\.sorted.bam|\1\.sorted.bam\t350\t\1|' ${k}_pindelinput.txt
done

#Run Pindel
for i in `cat samplesPindel.txt`; do
    pindel -f ${hg38} -i ${i}_pindelinput.txt -c ALL -o ${i}
#Create vcf files of deletion and insertion out of pindel
    pindel2vcf -p ${i}_D -r ${hg38} -R GRCh38 -d 201312 -G -v ${i}_D.vcf
    pindel2vcf -p ${i}_SI -r ${hg38} -R GRCh38 -d 201312 -G -v ${i}_SI.vcf
done
ADD COMMENTlink modified 18 months ago by Ming Tang2.4k • written 18 months ago by jsneaththompson60

most caller will ignore the duplicate read, failing-qualoty-reads, discordant reads, low MAPQ, etc... so the number of reads could be lower in the VCF/DP field that in the BAM.

ADD REPLYlink written 18 months ago by Pierre Lindenbaum114k

In addition to Pierre's answer, there are also various ways of calculating coverage. For example, deletion events could be classified as covered, or might not be. Secondary/supplementary alignments might be considered covered, or they might be ignored. Variant-callers might trim the ends of reads, or they might not. They can also perform realignment.

So, it's not surprising that the results differ; I wouldn't worry about it.

ADD REPLYlink written 18 months ago by Brian Bushnell16k
0
gravatar for Ming Tang
18 months ago by
Ming Tang2.4k
Houston/MD Anderson Cancer Center
Ming Tang2.4k wrote:

as Pierre says, most callers do some filtering of the reads, so the reads number in the VCF file will be fewer than in the bam file.

ADD COMMENTlink written 18 months ago by Ming Tang2.4k
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: 1287 users visited in the last hour