[Bioperl-l] Bio::Index::BlastTable->fetch_report

Tristan Lefebure tristan.lefebure at gmail.com
Tue Apr 6 19:45:23 UTC 2010


> Any other way to detect
> and skip sequences that have no hits would also work....

replying to myself: what about having a is_a_query() methods in this
module?

That would make something like:


my $p = 'Sag4|SAL_2304';
my $inx = Bio::Index::BlastTable->new(-filename => $blastdb,
                                   -write_flag => 0);
if($inx->is_a_query($p)) {
       my $blast_result = $inx->fetch_report($p);
}


????



More information about the Bioperl-l mailing list