SVMU (Structural Variants from MUmmer) attempts to identify comprehensive sequence variants via the alignment of two contiguous genome assemblies. Its goal is to combine the strengths of different aligners to annotate duplicates, large indels, inversions, small indels, and SNPs from whole genome alignments.
SVMU is being re-developed by the Chakraborty Lab at Texas A&M. We are incorporating new features (e.g., comparing multiple genomes, a VCF output). We expect the new version to be available sometime in Fall 2024. Feel free to continue using the existing version until then and send an email to mahul@tamu.edu if you encounter an issue. SVMU currently works with MUMmer (both v3 and v4 works).
NOTE: If you have an idea or suggestion (including collaboration ideas), write to the above email or Trevor (tdmillar@tamu.edu). However, if you are interested in the svmu versions used in the 2018 Nature Genetics (A4) and 2019 Nature Communication (DSPR) papers, see below:
If you publish results obtained with this pipeline, please cite SVMU as described here https://www.nature.com/articles/s41467-019-12884-1. The version used in the paper is available through commits prior to March 6,2018.
Download and compile the programs -
make
Obtain the mummer and lastz alignments:
nucmer --threads n --prefix sam2ref ref.fasta sample.fasta
lastz ref.fasta[multiple] sample.fasta[multiple] --chain --format=general:name1,strand1,start1,end1,name2,strand2,start2,end2 > sam_lastz.txt
(LASTZ output that svmu reads should have only six columns as mentioned in the lastz command above) (for relatively small genomes, --maxmatch can also be used for nucmer)
Run svmu:
svmu sam2ref.mm.delta ref.fasta sample.fasta snp_mode sam_lastz.txt prefix
snp_mode should be 'h' or 'l'. h = report SNPs; l = no SNPs. currently, this option is turned off [will be activated in the near future].
prefix = a prefix that will be added to the output files.
sv.prefix.txt = A tab-delimited file that summarizes structural mutations (indels, CNVs, inversions) in the sample genome with respect to the reference genome.
small.prefix.txt: A tab-delimited file containing SNPs and small indels that occur within syntenic blocks (or MUMs).
cnv_all.prefix.txt: A tab-delimited file with all the reference genomic regions that are present in higher copy numbers (>1) in the sample genome. Those with "trans" in their names mean either it is a transposable element or non-TE copies of a gene in different chromosomes.
cm.prefix.txt: A bed file with the reference genomic regions that are syntenic between the two genomes.
Finally, If you are using SVMU for your research, please keep in mind that SVMU has not been extensively tested on genomes bigger than Drosophila. So, there is no guarantee that it will work well with other genomes. Currently, it requires ~2.5G memory for the D. melanogaster genome.
KNOWN BUGS/PLANNED FUTURE IMPROVEMENTS: