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

Tristan Lefebure tristan.lefebure at gmail.com
Tue Apr 6 15:45:23 EDT 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

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