phylip to fasta

[ac91]

Hi,

Can anyone provide me a script in python or perl that converts phylip interleaved into fasta? I already found some scripts but it seems that they all work with phylip sequential and not with interleaved.

Thanks a lot!!

[Brian Bushnell]

What a coincidence, I have a script which does just that in the BBMap package!

phylip2fasta.sh in=file.phylip out=file.fasta

[ac91]

Thanks so much!

I already saw that there are others useful scripts in bbmap.

[ac91]

What a coincidence, I have a script which does just that in the BBMap package!

phylip2fasta.sh in=file.phylip out=file.fasta

Sorry Brian, but seems that your script does not work well with my files. Can you tell me why?

For example:
file.phy
3 120
seq1 ATGTGTTCAGCTCTAGTCGACAAGTATTTCAACTACCAGGAACCTCCGAATCGACAACTT
seq2 ATGTGTTCGGCTCTTGTAGAGAGACATTTCAGCTTCATAGCGTCTTCG------------
seq3 ATGTGTTCAGCTCTAGTCGACAAATATTTTAACTACCAGGAACCTCCTGATCGACAACAC

CAACTTTCGTTTGAAGACTCCAAAATGGCCTACGACAACGGGTCGGACAGCCAGGGCTCG
------AGGTTTGAAGACTCGAGAATGGCTTACGACAACGGGTCGGACAGCCAGGGTTCG
CAGCTTCAATTTGAGGCCACCAAAATGGCCTACGACAACGGCTCAGATAGCCAAGGCTCC

file.fas
>seq1
ATGTGTTCAGCTCTAGTCGACAAGTATTTCAACTACCAGGAACCTCCGAATCGACAACTT
>seq2
ATGTGTTCGGCTCTTGTAGAGAGACATTTCAGCTTCATAGCGTCTTCG
>seq3
ATGTGTTCAGCTCTAGTCGACAAATATTTTAACTACCAGGAACCTCCTGATCGACAACAC
>CAACTTTCGTTTGAAGACTCCAAAATGGCCTACGACAACGGGTCGGACAGCCAGGGCTCG

>------AGGTTTGAAGACTCGAGAATGGCTTACGACAACGGGTCGGACAGCCAGGGTTCG

>CAGCTTCAATTTGAGGCCACCAAAATGGCCTACGACAACGGCTCAGATAGCCAAGGCTCC

[Brian Bushnell]

I only had a chance to try it on one file, and it worked fine in that case, but that's not a very robust test... can you send me the file you're using?

What you have posted might be enough to replicate the problem, but to ensure it works correctly in entirety, the full file is probably essential.

[dgscofield]

If you haven't already found a solution, I've got a simple script that can interconvert between any of the formats supported in BioPerl's Bio::AlignIO, which you'd need to have installed. For Phylip to Fasta:

convertAlignment.pl -if phylip -of fasta < file.phy > file.fa

It's available in my bioinfo/scripts repository.

[ac91]

I only had a chance to try it on one file, and it worked fine in that case, but that's not a very robust test... can you send me the file you're using?

What you have posted might be enough to replicate the problem, but to ensure it works correctly in entirety, the full file is probably essential.

Sorry for the late reply.
You can find my files attached, but I changed the fasta and phylip to .txt because is impossible to attach files with that extensions.