[Bioperl-l] Parsing ClustalW reports

Rhazes Spell rus1@duke.edu
Sat, 29 Jun 2002 15:15:08 -0400 (EDT)


I am having trouble parsing a ClustalW report.  After creating the stream
I call the next_aln method, but I don't get an alignment object.  I
assume that I making a simple newbie mistake.  If anyone could be of assistance I would
appreciate it.

Here is a simple script and the output generated:

use strict;
use Bio::AlignIO;

my $stream  = Bio::AlignIO->new(-format =>'clustalw',
				-file=>'act1.aln');
my $aln = $stream->next_aln();

foreach my $seq ( $aln->each_seq() ) {
    print $seq->desc();
}

OUTPUT:
perl clustalParse.pl
Can't call method "each_seq" on an undefined value at clustalParse.pl line
21, <GEN0> line 44.



Thanks in advance
Rhazes Spell

**************************************
Duke University
Department of Biomedical Engineering
Duke University Medical Center,
 Center for Genome Technology

rhazes.spell@duke.edu
919.684.2720
***************************************