Question: How to convert MUMmer output to SAM, GFF, or any other common format
0
gravatar for Michael Dondrup
4.6 years ago by
Bergen, Norway
Michael Dondrup47k wrote:

I am aligning contigs against a reference genome and would like import the output of Nucmer-MUMmer as a track in GBrowse. MUMmer uses its own idiosyncratic output formats (delta format), and to my surprise I was unable to find any working parser for Mummer in BioPerl or BioPython.

I found some old requests, like

but seemingly nothing of this code ever made it. Does anyone know more about it?

parser format mummer • 5.1k views
ADD COMMENTlink modified 4.6 years ago by kloetzl1.1k • written 4.6 years ago by Michael Dondrup47k

Note, http://www.bioperl.org/ and http://biopython.org/ seem to be both down atm, so I will check them once they come back.

ADD REPLYlink written 4.6 years ago by Michael Dondrup47k

This is closely related: Mummer To Viewable Alignment Format (Fasta Or Aln...) but doesn't answer my question.

ADD REPLYlink written 4.6 years ago by Michael Dondrup47k
2
gravatar for kloetzl
4.6 years ago by
kloetzl1.1k
European Union
kloetzl1.1k wrote:

The Debian package for Mummer comes with the delta2maf program which does exactly what it claims to do. That program comes from the Mugsy suite if you are on a different system.

ADD COMMENTlink modified 4.6 years ago • written 4.6 years ago by kloetzl1.1k

Dear koetzl, thank you for your answer. I have found delta2maf as ELF binary in mugsy on sourceforge, but I am not using Debian but CentOS, and the binary Seg faults. I guess I should re-compile it, do you know where the original sources are? I was also searching https://packages.debian.org and it looks like there is no source for delta2maf other than this patch: https://anonscm.debian.org/viewvc/debian-med/trunk/packages/mummer/trunk/debian/patches/addition_from_mugsy.patch?view=markup

ADD REPLYlink written 4.5 years ago by Michael Dondrup47k

Oh source is here: svn://svn.code.sf.net/p/mugsy/code/trunk to compile:

svn checkout svn://svn.code.sf.net/p/mugsy/code/trunk/MUMmer3.20/
cd MUMmer3.20/src/tigr
make delta2maf
g++  -O3 tigrinc.cc -c -o tigrinc.o
g++  -O3 translate.cc -c -o translate.o
g++  -O3 delta.cc -c -o delta.o
g++  -O3 delta2maf.cc tigrinc.o translate.o delta.o -o      /export/home/michaeld/compile/MUMmer3.20/src/tigr/delta2maf; chmod 755 /export/home/michaeld/compile/MUMmer3.20/src/tigr/delta2maf

./delta2maf
Segmentation fault (core dumped)
   :(
ADD REPLYlink modified 4.5 years ago • written 4.5 years ago by Michael Dondrup47k
1

Yeah, the program segfaults when no args are provided. :/ You have to call it with the delta file as first argument. ./delta2maf foo.delta

ADD REPLYlink written 4.5 years ago by kloetzl1.1k

Yep, it is working! I think the maf format is still a bit verbose (for 80k contigs vs. 30k scaffolds, yields a 1.4GB file) to import into GBrowse, because it contains all the aligned sequences, I think I can further convert it to GFF though, and therefore can accept this answer.

ADD REPLYlink modified 4.5 years ago • written 4.5 years ago by Michael Dondrup47k

I have recently been using (a plugin I made) mafviewer in jbrowse for viewing these whole genome alignments :)! https://github.com/cmdcolin/mafviewer

ADD REPLYlink written 3.5 years ago by cmdcolin1.4k
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: 990 users visited in the last hour