[Bioperl-l] bl2seq and next_aln()
Gregory Drake Wilson
gdw1 at cornell.edu
Wed Aug 3 12:00:13 EDT 2005
I am trying to parse a bl2seq file but am only being returned one of the
alignments when there are 2+.
Code:
my @params = (program => 'blastn' , 'outfile' => 'bl2seq.out');
my $factory = Bio::Tools::Run::StandAloneBlast->new(@params);
my $report = $factory->bl2seq($seq1, $seq2);
my $str = Bio::AlignIO->new(-file=> 'bl2seq.out','-format' =>
'bl2seq');
while ( my $aln = $str->next_aln() ) {
print $aln->consensus_iupac()."\n";
}
Opening 'bl2seq.out' shows mutiple alignments, yet this code only returns
the first one in the file. Any thoughts?
Greg
More information about the Bioperl-l
mailing list