#!usr/bin!perl -w
#local blast search

#intecration of bio-perl-tools
use Bio::SeqIO;
use Bio::Tools::Run::StandAloneBlast;


#read as-sequence via commandline
my($as) = $ARGV[0];


my $Seq_in = Bio::SeqIO->new (-file => $ARGV[0], 
             -format => 'fasta');
my $query = $Seq_in->next_seq(); 

my $factory = Bio::Tools::Run::StandAloneBlast->new(
               
 'program'  => 'blastp',
               
 'database' => 'tobydb',
               
  _READMETHOD => "Blast"
             );
$factory->outfile('blast.out');
my $blast_report = $factory->blastall($query);
my $result = $blast_report->next_result;
while( my $hit = $result->next_hit()) {
    print "\thit name: ", $hit->name(), 
    " significance: ", $hit->significance(), "\n";
}

exit;

use File::Temp qw(tempfile);
my($fh,$file) = tempfile(UNLINK => 1);
print $fh $ARGV[0];