SNP calling from multiple bacterial genomes
2
1
Entering edit mode
10.6 years ago
biotech ▴ 570

Hi there,

I'm trying to find a strategy to call core SNPs in a collection of hundreds of bacterial genomes. I've been able to get some of them using kSNP tool Scalable SNP Analyses of 100+ Bacterial or Viral Genomes, but we think that tool is loosing many of them.

I'm thinking of using a reference-based method.

Thanks so much for your help, Best, Bernardo

Citation: Gardner SN, Slezak T (2010) Scalable SNP Analyses of 100+ Bacterial or Viral Genomes. J Forensic Res 1:107. doi:10.4172/2157-7145.1000107

snp bacteria metagenomics • 5.7k views
ADD COMMENT
3
Entering edit mode
10.5 years ago
Erik Garrison ★ 2.4k

You should definitely use multiple methods to generate results. In comparison between methods you'll be more likely to catch issues generated by a particular approach, and the marginal contribution of very different variant detection methods often far outweighs the trouble of running them.

I develop a reference-based detection method, freebayes, which is designed for and has been used in exactly this kind of context.

Running freebayes should be very easy provided you have a reference and have aligned your data against it. Please ping me or the freebayes mailing list if you have any questions.

ADD COMMENT
1
Entering edit mode
8.5 years ago
moorem ▴ 240

There are lots of methods by which to align short reads and call variants, of which GATK performs the best. There can be little concordance between methods at times, particularly with indel calls so it's important.

I've written a short blog post on using GATK with the updated version of haplotype caller which can now be set to ploidy =1.

https://approachedinthelimit.wordpress.com/2015/10/09/variant-calling-with-gatk/

Let me know if this is any use to you!

ADD COMMENT

Login before adding your answer.

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