This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to parsing mugsy (.maf) alignment output to RaxML input

Hi all,

I was trying to run a phylogeny using RaxML of few bacterial genome. I have generated the alignment file using mugsy and convert .maf to phylip using Bioperl. Now the problem is mugsy alignment generated a multiple block local alignment not a global alignment file. Can anyone help me how can I convert that multi block alignment to single global alignment file. To run mugsy I used the command bellow

$mugsy --directory /data/output --prefix mygenomes genome1.fasta genome2.fasta genome3.fasta ..

Thanks advance.

mugsy raxml phylogeny
• 3,903 views • link
updated by Ram Ram • written by jeccy.J jeccy.J

1 answer

To do that I used biopython, which has a module that is not in the main distribution yet, the module is MafIO and with the module AlignIO you can convert from maf to phylip.

• 0 views • link
updated by Ram Ram • written by rafa.rios.50 • rafa.rios.50

Log in to answer this question.

More posts like this