Question: Mauve alignment vs samtools for SNP calling
0
gravatar for gabriel.jabud
5 months ago by
United States
gabriel.jabud20 wrote:

If I have two assembled genomes (one reference), what is the difference between doing SNP calling using something like Samtools from the raw Illumina reads and doing a Mauve alignment after assembly of each genome? Would Samtools be much more accurate? It seems easier to just do a Mauve alignment and extract the SNPs since I've already done the assemblies but maybe this is oversimplifying it. These are bacterial genomes if that matters.

ADD COMMENTlink modified 5 months ago by jrj.healey13k • written 5 months ago by gabriel.jabud20
2
gravatar for jrj.healey
5 months ago by
jrj.healey13k
United Kingdom
jrj.healey13k wrote:

You should use the reads and samtools if you have them. Part of identifying SNPs etc is determining how well they're supported (i.e. are they real or sequencing errors for instance). If you aligned 2 genomes with mauve, firstly it will be horribly inaccurate as you say, but also you'll have no information about whether that SNP is supported in most or all of the reads covering that position.

ADD COMMENTlink written 5 months ago by jrj.healey13k
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: 1074 users visited in the last hour