Why VarScan cannot call variants which appear in 1000genome

0

Entering edit mode

8.8 years ago

sxl919 ▴ 10

HI,

I am using VarScan for variant calling and I use the bam files from 1000genome as control for my project. I run the bam file through VarScan using default setting and then annotate the vcf using Annovar. However, when I try to compare the vcf generated from VarScan and the vcf from 1000 genome, there are some variants which VarScan hasn't called. For example, for SNP rs10000057 is called in the 1000 genome vcf files but it didn't get any calls in the VarScan generated vcf files. Has anyone experience the same situation as I do? Can anyone explain why?

Thank you

VarScan 1000genome • 2.1k views

ADD COMMENT • link updated 16 months ago by Ram 43k • written 8.8 years ago by sxl919 ▴ 10

0

Entering edit mode

What do you mean by VarScan didn't call ? Are they truly present in your bam file (IGV lookup) ?

ADD REPLY • link 8.8 years ago by poisonAlien ★ 3.2k

0

Entering edit mode

Thanks for your comment. I did look at the IGV and the variant is present. However, I was doing variant call individually so I do think that the variant will show up if I do a joint calling.

ADD REPLY • link 8.8 years ago by sxl919 ▴ 10

0

Entering edit mode

hopefully joint call will help (going to know whether VarScan call that position tomorrow)

ADD REPLY • link 8.8 years ago by sxl919 ▴ 10

Login before adding your answer.