Bioperl: Processing blast report of multiple input sequences

Kris Boulez krbou@pgsgent.be
Thu, 18 May 2000 16:39:57 +0200


Is a Bio::Tools::Blast object supposed to be able to parse a blast
output file containing the output of blasting multiple sequences (i.e.
the input to blast (BLASTN 2.0.9) was a Fasta file containing multiple
sequences. When doing 

  my $blastobj = Bio::Tools::Blast->new( -file => "$in_file",
                                         -parse => 1,
                                         -signif => '1',
                                          );
          
  foreach my $hit ($blastobj->hits) {
                         [ do somehting with $hit ] 
  }

I only see the hits for the first sequence. Is this intentional and
should I split up the report my self for each input sequence ?

Kris,
=========== Bioperl Project Mailing List Message Footer =======
Project URL: http://bio.perl.org/
For info about how to (un)subscribe, where messages are archived, etc:
http://www.techfak.uni-bielefeld.de/bcd/Perl/Bio/vsns-bcd-perl.html
====================================================================