MUMmer4

MUMmer4 performs rapid alignment of large DNA sequences to enable whole-genome and assembly comparisons at terabase-scale.

Key Features:

• Enhanced Data Structure: Replaces the 32-bit suffix tree with a 48-bit suffix array, enabling handling of sequence lengths up to 141 terabase pairs (Tbp).
• Parallel Processing: Uses parallel processing to improve alignment speed for large query sets.
• Whole-Genome Alignment: Performs whole-genome sequence alignment for large-scale genomic comparisons.
• Read-to-Reference Alignment: Aligns sequencing reads to reference genomes, with potentially lower sensitivity and accuracy than dedicated read aligners.
• Assembly Comparison: Compares different assemblies of the same genome to identify concordant and discordant regions.
• nucmer Program: Includes the nucmer alignment program for nucleotide-level alignments.

Scientific Applications:

• Whole-genome comparison: Aligning whole genome sequences for comparative genomics and structural variation analysis.
• Read mapping: Mapping sequencing reads to reference genomes for variant detection and assembly validation.
• Assembly-to-assembly comparison: Comparing different assemblies of the same genome to assess assembly differences and improvements.
• Comparative genomics example: Pairwise genome comparisons such as human versus chimpanzee, which have been reported as 98% identical across 96% of their length.

Methodology:

Implements a 48-bit suffix array in place of a 32-bit suffix tree, employs parallel processing for speed, and includes the nucmer alignment program.