I am having trouble with Mummer and Dgenies to provide me with a nice continuous dot plot. Basically, I have a bunch of de novo assembled genomes from Flye where the number of contigs is essentially to chromosome level. However, some of the entire contigs are inverted, compared to one another or to the reference and I'm not sure how to visualize these on a dotplot as one continuous diagonal ranging from chromosome 1-12. It seems to start the dotplot at an arbitrarily chosen chromosome for each plot rather than being arranged. I'm also not sure how to 'invert' the inverted chromosomes/contigs to correctly fit the dot plot for better visualization. I've also used RagTag. Even if I can't get the contings into the correct orientation, I would like to at least order the dot plot from chromosomes 1-12.
Can anyone give me some advice as to how to tackle this?
•
written
by
Drew •
4 answers
With dgenies you can left click on any alignment square of the dotplot, this will show a local menu, and reverse it using the "reverse" item.
I'm also not sure how to 'invert' the inverted chromosomes/contigs to correctly fit the dot plot for better visualization.
Just reverse complement the sequences that are in opposite orientation to the reference.
Even if I can't get the contings into the correct orientation, I would like to at least order the dot plot from chromosomes 1-12.
There's a Sort contigs button in D-genies. That should sort the contigs in the same order as the reference. Check also, if the reference itself is sorted properly.
this blog post has a method for flipping reverse orientation contigs. it is it's own plotting tool instead of dgenies but perhaps could be leveraged http://jmonlong.github.io/Hippocamplus/2017/09/19/mummerplots-with-ggplot2/
You might be interested in my pipeline (https://nf-co.re/pairgenomealign), which runs last-dotplot in a mode that attempts to flip sequences to diagonalise the plot.
Log in to answer this question.
-
Genome Assembly QC from BAM files
written by SomeOne •
-
Genome Assembly Verification, Chimeric contigs
written by SomeOne •
-
Canu genome assembly output
written by SomeOne •
-
Obtain satellite DNA monomer from dot plot
written by kirillkirilenko •
-
Plasmid or Chromosome?
written by Jaykumar •
-
Group a set of plots using MUMmer
written by Rox
-
Speeding up NUCmer for de novo contig to reference alignment
written by lkw222 •
-
Mummer alignment of two chromosomes
written by wjidea •
-
Tutorial: Mini-tutorial: How to fix too-long sequence names in SyMAP
written by Philipp Bayer
-
How To Compare Complete Genomes With Mummer Or Other Tools
written by Layla •
