MAUVEALIGNER

NAME

mauveAligner - efficiently constructing multiple genome alignments

SYNOPSIS

mauveAligner [options] <seq1 filename> <sml1 filename> ... <seqN filename> <smlN filename>

DESCRIPTION

The mauveAligner and progressiveMauve alignment algorithms have been implemented as command-line programs included with the downloadable Mauve software. When run from the command-line, these programs provide options not yet available in the graphical interface.

OPTIONS

--output=<file> Output file name.

Prints to screen by default

--mums Find MUMs only, do not attempt to determine locally collinear blocks
(LCBs)
--no-recursion Don’t perform recursive anchor identification (implies
--no-gapped-alignment)
--no-lcb-extension If determining LCBs, don’t attempt to extend the LCBs
--seed-size=<number> Initial seed match size, default is log_2( average
seq. length )
--max-extension-iterations=<number> Limit LCB extensions to this number of
attempts, default is 4
--eliminate-inclusions Eliminate linked inclusions in subset matches.
--weight=<number> Minimum LCB weight in base pairs per sequence
--match-input=<file> Use specified match file instead of searching for
matches
--lcb-match-input Indicates that the match input file contains matches
that have been clustered into LCBs
--lcb-input=<file> Use specified lcb file instead of constructing LCBs
(skips LCB generation)
--scratch-path=<path> For large genomes, use a directory for storage of
temporary data. Should be given two or more times to with different paths.
--id-matrix=<file> Generate LCB stats and write them to the specified file
--island-size=<number> Find islands larger than the given number
--island-output=<file> Output islands the given file (requires
--island-size)
--backbone-size=<number> Find stretches of backbone longer than the given
number of b.p.
--max-backbone-gap=<number> Allow backbone to be interrupted by gaps up to
this length in b.p.
--backbone-output=<file> Output islands the given file (requires
--island-size)
--coverage-output=<file> Output a coverage list to the specified file (-
for stdout)
--repeats Generates a repeat map. Only one sequence can be specified
--output-guide-tree=<file> Write out a guide tree to the designated file
--collinear Assume that input sequences are collinear--they have no
rearrangements

Gapped alignment controls:

--no-gapped-alignment Don’t perform a gapped alignment
--max-gapped-aligner-length=<number> Maximum number of base pairs to
attempt aligning with the gapped aligner
--min-recursive-gap-length=<number> Minimum size of gaps that Mauve will
perform recursive MUM anchoring on (Default is 200)

Signed permutation matrix options:

--permutation-matrix-output=<file> Write out the LCBs as a signed
permutation matrix to the given file
--permutation-matrix-min-weight=<number> A permutation matrix will be
written for every set of LCBs with weight between this value and the value
of --weight

Alignment output options:

--alignment-output-dir=<directory> Outputs a set of alignment files (one
per LCB) to a given directory
--alignment-output-format=<directory> Selects the output format for
--alignment-output-dir
--output-alignment=<file> Write out an XMFA format alignment to the
designated file

Supported alignment output formats are: phylip, clustal, msf, nexus, mega, codon