This is a test version of Biostars. For the public version, visit https://www.biostars.org.
XMFA to MAF

Hi everyone :)

i want to convert my XMFA file (output of progressiveMauve) into MAF.

Has anyone had any luck converting xmfa into maf? And how did you do it?

Thanks :)

L

alignment xmfa maf
• 7,057 views • link
written by Lynda lyn • Lynda lyn

thank you Sej Modha :) i'll try it !

• 0 views • link
written by Lynda lyn • Lynda lyn

Please use ADD REPLY/ADD COMMENT when responding to existing posts to keep threads logically organized.

If the solution suggested by @Sej works, please come back and accept the answer (check mark) to indicate that.

• 0 views • link
written by GenoMax GenoMax

2 answers

Thank you Seja Modha for your help :) i found this c++ script to convert xmfa to maf:

https://github.com/Wyss/mauve-py/blob/master/mauve/src/mauveAligner/src/xmfa2maf.cpp

• 0 views • link
written by Lynda lyn • Lynda lyn

Lynda Ait-Zai : I moved this to an answer since it does what you had originally asked for. You can accept this answer (use the green check mark) to provide closure for this thread.

• 0 views • link
written by GenoMax GenoMax

There are no C++ scripts, only C++ programs.

• 0 views • link
written by Raygozak Raygozak

In past, I have used following script to convert xmfa to fasta. https://github.com/lskatz/lskScripts/blob/master/convertAlignment.pl

• 0 views • link
written by Sej Modha Sej Modha

thank you Sej Modha :) i'll try it !

• 0 views • link
written by Lynda lyn • Lynda lyn

Hi Sej Modha , i tried the perl script. its works to have a fasta format but not to get the maf format.

Here is the output of my command: 

$ ./transform_fasta.pl -i threeway.xmfa -o /home/lynda/Documents/threeway.maf -f maf

Converting threeway.xmfa (xmfa) to file /home/lynda/Documents/threeway.maf
------------- EXCEPTION: Bio::Root::NotImplemented -------------
MSG: Abstract method "Bio::AlignIO::maf::write_aln" is not implemented by package Bio::AlignIO::maf.
This is not your fault - author of Bio::AlignIO::maf should be blamed!
STACK: Error::throw
STACK: Bio::Root::Root::throw /usr/local/share/perl/5.22.1/Bio/Root/Root.pm:447
STACK: Bio::Root::RootI::throw_not_implemented /usr/local/share/perl/5.22.1/Bio/Root/RootI.pm:741
STACK: Bio::AlignIO::maf::write_aln /usr/local/share/perl/5.22.1/Bio/AlignIO/maf.pm:189
STACK: main::convertAlnToFile ./transform_fasta.pl:99
STACK: main::convertAln ./transform_fasta.pl:41
STACK: main::main ./transform_fasta.pl:28
STACK: ./transform_fasta.pl:15
----------------------------------------------------------------

Thank you for the help !

• 0 views • link
written by Lynda lyn • Lynda lyn

MFA stands for multiple fasta format so you should provide -f fasta parameter in the command.

• 0 views • link
written by Sej Modha Sej Modha

When i add -f fasta, I get a maf file as output, but by opening the file It doesn't look like a maf file at all, it's a fasta format!. Maf contains more information than fasta format. Bio::AlignIO reads many formats but does not write in every format (the same is true for Bio::SeqIO). AlignIO currently supports output in these formats:

fasta
mase
selex
clustalw
msf
phylip
po
stockholm
XMFA
metafasta
• 0 views • link
written by Lynda lyn • Lynda lyn

I am sorry but I have not written those scripts so I won't able to help much but you could try other built-in format available in Mauve and parse them to achieve what you are trying to do.

Supported alignment output formats are: phylip, clustal, msf, nexus, mega, codon
• 0 views • link
written by Sej Modha Sej Modha

Log in to answer this question.

More posts like this