Question: Group a set of plots using MUMmer
0
gravatar for Rox
4.0 years ago by
Rox1.3k
France / Toulouse / GeT-Plage
Rox1.3k wrote:

Hello everyone !

This is my first try using MUMmer, I would like to visualize something close to this publication : http://biorxiv.org/content/biorxiv/early/2016/12/02/090969.full.pdf (page 18). I'm trying to align the 6 longest contigs of my assembly to see where it aligns on D.melanogaster genome.

I've already tried something like that using mummer :

#Launching mummer
./mummer -b melanogaster.fasta suzukii_6first_contigs.fasta > mummer_mel_against_top6_canu.mums
#Plot using mummer plot
./mummerplot -x "[0,250000]" -y "[0,250000]" -postscript -p mummer_mel_against_top6_canu mummer_mel_against_top6_canu.mums'

Each time I've tried something, the -x and -y options where required, but I don't know at which contigs from my assembly or which chromosomes to the ref it correspond too...

I've followed the example proposed here : http://mummer.sourceforge.net/examples/#nucmernucmer , but I can't find an example of how to create the same kind of plot I saw in the bioarXiv paper.

Anyone with some experience with mummer can help me understanding how the plots are working ?

Cheers,

Roxane

alignment mummer mummer-plot • 1.5k views
ADD COMMENTlink modified 4.0 years ago by skbrimer650 • written 4.0 years ago by Rox1.3k
0
gravatar for skbrimer
4.0 years ago by
skbrimer650
United States
skbrimer650 wrote:

It looks like it is actually 6 MUMer plots. The paper says they took their 6 largest contigs and mapped it against the 6 D.melanogaster chromosomal arms. I do not know anything about the fruit fly genome so I do not know if you need to make a multi-fasta of all the arms or not, but you would only map one contig to 6 arms. So to recreate the figure you would need to make that plot for each contig and then place neatly together after the fact.

ADD COMMENTlink written 4.0 years ago by skbrimer650
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1678 users visited in the last hour
_