Re: [MUMmer-help] whole genome alignment question

[Diogo C. <dio...@gm...>]

Hi Adam,

Thanks for the tip.  Just tried it out and looking at the .report file I can see the average identity between the 2 genomes.  Thanks!


Best,
d.



----------
Diogo Camacho
e: dio...@gm...
p: (617) 945-4383





> On Mar 3, 2015, at 10:30 AM, Adam Phillippy <aph...@gm...> wrote:
> 
> Hi Diogo,
> If you run the 'dnadiff' script that comes with MUMmer it will generate a .report file that will include global similarity metrics based on Nucmer alignments.
> 
> Best,
> -Adam
> 
> 
> On Tue, Mar 3, 2015 at 9:50 AM, Diogo Camacho <dio...@gm...> wrote:
> Hi,
> 
> This is a very naive, simple question, but I can’t seem to locate a good answer for it on the manual.  I have a set of bacterial genomes that I want to align against a different (single) bacterial genome.  The obvious tool to do this seems to be NUCmer and I did a test run of 1 genome against the “reference” genome by doing:
> 
> nucmer -maxmatch -p test-output <REFERENCE GENOME> <GENOME TO BE COMPARED>
> 
> NUCmer runs fine and when I do show-coords I get the table with the alignments and the respective identity between the 2 genomes I am comparing.  But what I would like would be to get a single, global, percent identity score between the 2 genomes, as I basically want to compare the two genomes at a global level.  Is there a simple way to do this using the MUMer suite?
> 
> Many thanks!
> 
> 
> Best,
> d.
> 
> 
> 
> ----------
> Diogo Camacho
> e: dio...@gm...
> p: (617) 945-4383
> 
> 
> 
> 
> 
> 
> ------------------------------------------------------------------------------
> Dive into the World of Parallel Programming The Go Parallel Website, sponsored
> by Intel and developed in partnership with Slashdot Media, is your hub for all
> things parallel software development, from weekly thought leadership blogs to
> news, videos, case studies, tutorials and more. Take a look and join the
> conversation now. http://goparallel.sourceforge.net/
> _______________________________________________
> MUMmer-help mailing list
> MUM...@li...
> https://lists.sourceforge.net/lists/listinfo/mummer-help
>