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
====================================================================