Creating Phylogenetic tree from MAF format file

[mez]

I hope this is relevant to the discussions here.

I am working on a metagenome project and have genome sequences of several strains of a bacteria species, un-assembled, but in contig reads. I used the mugsy program to align them and I have been trying to find a tool to create a tree out of the MAF format file I generated, but I am having issues because of the nature of the original sequences being in contigs rather than assembled.

I am hoping someone jas worked with unassembled genomes with contig reads and could be able to point me to how I could generate an acceptable multi-alignment format from the data

Thanks as always.

mez.

[themerlin]

I'm a little unclear about your question. If you have contigs, then you are working with an assembly. Do you mean that you are working with overlapping paired-end reads?

[mez]

Sorry, I guess I should say the reads have been assembled into contigs, I was having trouble aligning the genomes of the different strains, but I am making progress on that now. Just trying to figure how best to use the MAF format file generated to create a phylogenetic tree.

[themerlin]

I'm not sure that you can create a tree from a MAF file with overlapping reads. With Mugsy, you usually provide separate files for each genome and then infer a tree from sequence conserved in each block. If you can get your MAF file in this format, you can create a FASTA alignment with Galaxy or with scripts in the bx-python package (https://bitbucket.org/james_taylor/bx-python/wiki/Home).