[Bioperl-l] bitscore and rawscore for blast reports
Razi Khaja
razi at genet.sickkids.on.ca
Wed Dec 15 15:19:06 EST 2004
Im am trying to parse blast reports using the Bio::SearchIO system.
#!/usr/bin/perl
use strict;
use Bio::SearchIO;
my( $blastFile ) = @ARGV;
my $searchio = new Bio::SearchIO( -file=>$blastFile, -format=>'blast' );
while( my $result = $searchio->next_result() ){
while( my $hit = $result->next_hit() ){
while( my $hsp = $hit->next_hsp() ) {
More information about the Bioperl-l
mailing list