[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