[Bioperl-l] problems parsing interlevaed phylip file with Bio::AlignIO
Hannes Hettling
hannes.hettling at gmail.com
Fri May 23 10:11:19 EDT 2014
Dear all,
I am having trouble reading in an interleaved phylip file. It works with a
smaller example, but not with the file attached to this post (93 taxa,
36338 sites).
Here is my code:
use Bio::AlignIO;
my $in = Bio::AlignIO->new(-file => 'myphylip.phy' ,
-format => 'phylip',
-interleaved => 1,
-idlength => 10);
while ( my $aln = $in->next_aln() ) {
print "Alignment found!\n";
print ref ($aln)."\n";
}
The variable $in is a Bio::AlignIO::phylip object, as expected. The phylip
file seems to be valid, I was able to convert it to fasta (using
http://sequenceconversion.bugaco.com/converter/biology/sequences/phylip_to_fasta.php)
and then I could read the fasta with AlignIO.
Does anyone have an idea what I could me missing?
Many thanks in advance,
Hannes
-------------- next part --------------
A non-text attachment was scrubbed...
Name: myphylip.phy
Type: application/octet-stream
Size: 3383535 bytes
Desc: not available
URL: <http://lists.open-bio.org/pipermail/bioperl-l/attachments/20140523/e76cd446/attachment-0001.obj>
More information about the Bioperl-l
mailing list