problem with SNPs and Indels detectiong useng WGS data
1
0
Entering edit mode
6 weeks ago
Sara ▴ 160

we have WGS data and looking for SNPs and Indels globally. we found Indels and SNPs using PCR and also when we load the sequencing data on IGV, we can see those events by eye (of course with low frequency). but when we used the following tools we can not detect them in the resulting VCF file:

lofreq , mutect, strelka and bcftool

do you know what could be the problem?

WGS • 313 views
ADD COMMENT
0
Entering edit mode

you asked many questions on this forum but you validated only one answer. See here

ADD REPLY
0
Entering edit mode

OK thanks for letting me know. do you have any answer for this question?

ADD REPLY
0
Entering edit mode

Disagreement between variant callers is expected, and it should also be expected that these [variant callers] will 'miss' (fail to call) genuine variants. For such a question as this, however, you have to explain better the entire pipeline that you've used (including sample preparation, alignment tool and metrics, etc), and show as much relevant code as possible.

ADD REPLY
0
Entering edit mode
6 weeks ago
EagleEye 6.9k

Hi,

I had this issue while ago for detecting very-rare variants. I was using freebayes and I had to adjust the default parameters (--min-alternate-count 1 --min-alternate-fraction 0.0000001). I am sure when you play with the default parameters you will manage to detect with any tools (you just have to dig in a bit!).

ADD COMMENT

Login before adding your answer.

Traffic: 1199 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