How to convert MUMmer output to SAM, GFF, or any other common format
1
1
Entering edit mode
6.1 years ago

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?

mummer format parser • 6.8k views
ADD COMMENT
0
Entering edit mode

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 REPLY
0
Entering edit mode

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

ADD REPLY
2
Entering edit mode
6.1 years ago
kloetzl ★ 1.1k

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 COMMENT
0
Entering edit mode

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 REPLY
0
Entering edit mode

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 REPLY
1
Entering edit mode

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 REPLY
0
Entering edit mode

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 REPLY
0
Entering edit mode

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

ADD REPLY

Login before adding your answer.

Traffic: 2315 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6