Question: How to parsing mugsy (.maf) alignment output to RaxML input
0
gravatar for jeccy.J
3.3 years ago by
jeccy.J50
Argentinia
jeccy.J50 wrote:

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 • 1.6k views
ADD COMMENTlink modified 2.9 years ago by rafa.rios.5050 • written 3.3 years ago by jeccy.J50
0
gravatar for rafa.rios.50
2.9 years ago by
rafa.rios.5050
Colombia
rafa.rios.5050 wrote:

To do that I used biopython, which has a module that is not in the main distribution yet, the module is MafIO (http://biopython.org/wiki/Multiple_Alignment_Format) and with the module AlignIO you can convert from maf to phylip.

ADD COMMENTlink written 2.9 years ago by rafa.rios.5050
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: 1534 users visited in the last hour