Why VarScan cannot call variants which appear in 1000genome
0
0
Entering edit mode
6.3 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 • 1.6k views
ADD COMMENT
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
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
0
Entering edit mode

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

ADD REPLY

Login before adding your answer.

Traffic: 2768 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6