Question: How to get complete SNPs list using Mauve
gravatar for Paul
5 months ago by
Paul70 wrote:

I am having multiple contigs for 1000s of genome, and each genome was mapped with the reference genome. I am using command line Mauve to extract all the SNPs from the alignments of contigs with the Reference genome. the command line used was

java org.gel.mauve.analysis.SnpExporter -f aln.xmfa -o aln.snps

And the list of SNPs I get for each genome is like the following:

SNP Ref GenWidePos1     sequence2Contig sequence2PosInContg sequence2GenWidePos2
GT  Ref 55649           Contig_18       130                 3255611
CT  Ref 55553           Contig_76       226                 3255707
AG  Ref 54394           Contig_72       1394                3256875
AG  Ref 51949           Contig_35       3839                3259320
TC  Ref 50557           Cont

But there more SNPs present between the contigs and the reference, but mauve gives me only half of the SNPs and it stops in the middle. Like the last line in the above mentioned example.

mauve snp java assembly genome • 349 views
ADD COMMENTlink modified 5 months ago • written 5 months ago by Paul70
Please log in to add an answer.


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