Question: the number of the reads covering the snp after HaplotypeCaller
0
gravatar for 543541656
10 weeks ago by
54354165620
54354165620 wrote:

when I used the following command:

(base) fyx@huangji-5885H-V5:~/bm4_mw/bsa$ samtools view s_add_rmdup.bam Chr11:28164766-28164766 | wc -l
43

we can see that the position 28164766 having 43 reads,but after running the following command:

gatk  HaplotypeCaller  \
-R rice.fasta \
-I input.bam \
--intervals chr11
-O output.g.vcf.gz \

wo can see the folling vcf file,the raw DP of the snp locating 28164766 is 52 and the filtered DP is 48,both are larger than 43,I want to know that's why???

Chr11   28164766    .   G   A   1132.64 .   AC=1;AF=0.500;AN=2;BaseQRankSum=-1.144;DP=52;ExcessHet=3.0103;FS=1.177;MLEAC=1;MLEAF=0.500;MQ=57.58;MQRankSum=4.345;QD=23.60;ReadPosRankSum=2.939;SOR=0.501 GT:AD:DP:GQ:PL  0/1:16,32:48:99:1140,0,576
depth snp gatk haplotypecaller • 136 views
ADD COMMENTlink modified 10 weeks ago by Pierre Lindenbaum134k • written 10 weeks ago by 54354165620
1
gravatar for pau.rodriguez.sodupe
10 weeks ago by
pau.rodriguez.sodupe20 wrote:

Hello

HaplotypeCaller performs local realignment for variant calling. That results in an altered/corrected bam (different to the input bam) from which determines Allele Depths...

If you want to see the bam used for calling you should use the --bamout ${OUTPUT_BAMOUT} argument. This bam will only contain the reads overlaping positions called in the output vcf (normally a much smaller bam file).

Best,

Pau

ADD COMMENTlink modified 10 weeks ago • written 10 weeks ago by pau.rodriguez.sodupe20
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: 2251 users visited in the last hour
_