Question: SNP calling from multiple bacterial genomes
1
gravatar for biotech
5.5 years ago by
biotech520
United States
biotech520 wrote:

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 • 3.8k views
ADD COMMENTlink modified 2.6 years ago • written 5.5 years ago by biotech520
3
gravatar for Erik Garrison
5.5 years ago by
Erik Garrison2.1k
Somerville, MA
Erik Garrison2.1k wrote:

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 COMMENTlink written 5.5 years ago by Erik Garrison2.1k
1
gravatar for moorem
3.5 years ago by
moorem220
United Kingdom
moorem220 wrote:

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 COMMENTlink written 3.5 years ago by moorem220
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1332 users visited in the last hour