Question: Conversion Of Alignment Files
1
gravatar for joe.loquasto
6.1 years ago by
joe.loquasto40
joe.loquasto40 wrote:

Hi,

I would like to use RAxML to construct a phylogentic tree. I have Mauve alignment files. Can I either convert to phylip? or what is the best way to construct a full genome alignment in this format.

Thank you Joe

alignment phylogenetics • 3.2k views
ADD COMMENTlink modified 4.7 years ago by Chrispin Chaguza230 • written 6.1 years ago by joe.loquasto40

A side note: RAxML is commonly used , it is very fast.  However, it employs homogeneous models.  If your data is not homogenous, you should consider heterogenous models such as CAT employed in PhyloBayes or p4.

 

 

ADD REPLYlink written 4.9 years ago by moranr250
0
gravatar for Josh Herr
6.1 years ago by
Josh Herr5.6k
University of Nebraska
Josh Herr5.6k wrote:

Hi Joe! MAUVE alignment files are in XMFA format and are a little troublesome to convert (in my opinion). When I have done this in the past I used a script written in Perl by Lee Katz. I'm not aware of any other converters for XFMA files, but I wouldn't be surprised if some are out there... Hope that helps.

ADD COMMENTlink written 6.1 years ago by Josh Herr5.6k
0
gravatar for Chrispin Chaguza
4.7 years ago by
University of Liverpool, UK
Chrispin Chaguza230 wrote:

Seems like BioPerl accepts XMFA input file format (http://search.cpan.org/~cjfields/BioPerl-1.6.924/Bio/AlignIO/xmfa.pm). You can use Bio::AlignIO to convert it to other formats e.g. Phylip in order to use it with RAxML.

Full genome alignments are only easy to construct for smaller genomes e.g. viruses but extremely difficult to construct for large genomes e.g. bacteria (or eukaryotic sequences) because it's computationally demanding among other factors. You can read this paper by Croucher et al (methods section) to see if you could use a similar approach like the one they used to construct phylogenetic tree for 600+ bacterial genomes (without doing a whole genome alignment) http://www.nature.com/ng/journal/v45/n6/full/ng.2625.html

 

ADD COMMENTlink written 4.7 years ago by Chrispin Chaguza230
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: 1714 users visited in the last hour