Mugsy: multiple whole genome alignment tool

Jit — Fri, 08 Dec 2017 17:41:14 -0600

Mugsy is a multiple whole genome aligner. Mugsy uses Nucmer for pairwise alignment, a custom graph based segmentation procedure for identifying collinear regions, and the segment-based progressive multiple alignment strategy from Seqan::TCoffee. Mugsy accepts draft genomes in the form of multi-FASTA files and does not require a reference genome.

To cite Mugsy, use:

Angiuoli SV and Salzberg SL. Mugsy: Fast multiple alignment of closely related whole genomes.Bioinformatics 2011 27(3):334-4

Address of the bookmark: http://mugsy.sourceforge.net/

Comment by BioStar

BioStar — Tue, 27 Oct 2020 19:29:52 -0500

Set up the environment

bash
source /path/to/install/mugsyenv.sh

Run mugsy

mugsy --directory /data/output --prefix mygenomes genome1.fasta genome2.fasta genome3.fasta

This example will align three genomes and output a file /data/output/mygenomes.maf. The --directory setting is also used for storing temporary files during the run.

The prefix of each input filename will be used as the genome name in the output files (eg. genome1 from genome1.fsa). Header lines in the FASTA files should not contain ':' or '-' to avoid parsing problems.

View or parse MAF output

The output of Mugsy is Multiple Alignment Format (MAF).
One option for browsing is GMAJ, which provides a reference-based stand-alone viewer for MAF files.

java -jar gmaj.jar mygenomes.maf

BOL: Mugsy: multiple whole genome alignment tool

Mugsy: multiple whole genome alignment tool

Comment by BioStar