What is the best way for bacterial variant calling?
1
0
Entering edit mode
3.4 years ago

I am trying to analyze Illumina sequenced bacterial data from known species to find mutations and I am stuck. Which tool I should use? It seems like half of the tools hugely over-represents the number of variants, even with strict parameters (Breseq, Freebayes) and the other half finds barely any variants (Snippy, VarScan) even with very loose parameters. I read several papers about this and tried a lot of software but still I feel like all the tools belong to either category. What can I do?

mutation prokaryote illumina variant calling • 1.0k views
ADD COMMENT
1
Entering edit mode
3.4 years ago

In general, the approach is to filter the variant calls yourself, rather than expect the variant caller to find just the right ones.

Freebayes does indeed find a lot of variants but filtering them after will allow you to hone in the ones that are most reliable and credible:

https://thegenomefactory.blogspot.com/2018/10/a-unix-one-liner-to-call-bacterial.html

ADD COMMENT

Login before adding your answer.

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