MUM&Co is a simple bash script that uses Whole Genome Alignment information provided by MUMmer (v4) to detect variants.
VERSION >= 3 UPDATE
Only uses MUMmer4 now and contains a thread count option
Contains a VCF output file with all calls currently being imprecise
Contains another output file containing the calls alongside the respective DNA impacted
This new step requires samtools installation
Now calls the reverse of tandem duplications, tandem contractions (>50bp)
MUM&Co is able to detect:
Deletions, insertions, tandem duplications and tandem contractions (>=50bp & <=150kb)
Inversions (>=1kb) and translocations (>=10kb)
MUM&Co requires installation of MUMmer4 and samtools.
MUM&Co will look for the MUMmer toolkit's and samtools scripts path using 'which xxxxx'.
An error warning will print and the script will stop if these paths cannot be found
This path can be editted directly in the script if required.
In order to help with downstream analysis:
Renaming and re-orientation of the query genome contigs to correspond to their reference counterparts
Tools such as RaGOO and Ragout can do this alongside scaffolding of contigs (this is not currently recommended for short-read based assemblies)
Options:
-r or --reference_genome path to reference genome
-q or --query_genome path to query genome
-g or --genome_size size of genome
-o or --output output prefix (default: mumandco)
-t or --threads thread number (default: 1)
-ml or --minlen minimum length of alignments (default: 50bp)
-b or --blast adds the blast option to identify is insertions or deletions look repetitive or novel
Test run script:
bash mumandco_v*.sh -r ./yeast.tidy.fa -q ./yeast_tidy_DEL100.fa -g 12500000 -o DEL100_test -t 2 -b
OUTPUT FOLDER:
Folder with alignments used for SV detection
Txt file with summary of SVs detected
TSV file with all the detected SVs
TSV file with all detected SVs plus the DNA associated with the event (all from reference except insertions)
VCF file with all calls currently being imprecise
TSV NOTES:
The last column in the TSV file contains notes.
'complicated' : multiple calls within the same region; generally overlapping insertions and deletions
'double' : several calls at the same coordinates; generally tandem duplications or contractions with multiple copy changes
']chrX:xxxxxx]' : a VCF inspired notation for the association of the translocation fragments with the other fragments
e.g. for chr1 with its right border at 250000bp assocaited with chr2 at 100000bp;
the note would be as follows for chr 1: ']chr2:100000]' and for chr2 : '[chr1:250000['
As such, each translocation fragment as called as an event, is now a breakend-like call and will be duplicated if both borders are involved in translocations
VCF TRA EVENT:
The later notation for the TSV file is currently being added to the alt column in the VCF for 'TRA' events.
Currently it is not a called a breakend site (contains no nucleotide at edge) but can be interpreted similarly
Note:
MUMmer4 is now required due to the hard wired thread option not available during alignment with MUMmer3
The blast option (-b /--blast) using BLAST to search for insertion and deletion events in the reference/query in order to label them as either mobile or novel events.
Reference:
Samuel O’Donnell, Gilles Fischer, MUM&Co: accurate detection of all SV types through whole-genome alignment, Bioinformatics, Volume 36, Issue 10, 15 May 2020, Pages 3242–3243, https://doi.org/10.1093/bioinformatics/btaa115