[Bioperl-l] Appending blast results for multi sequence data set.

Charles Hauser chauser@duke.edu
21 Oct 2002 10:56:45 -0400


Hi,

Trying to blast multiple sequences and collect all results in a single
file. 
If I set @params w/ 'outfile' => '>>../20021010/BLAST/Xbn.out'
			OR	=> ">>../20021010/BLAST/Xbn.out"

 blastall crashes:

MSG: blastall call crashed: 256 /usr/bin/blastall -p  blastn  -d  ../20021010/BLAST/X  -i  /tmp/4OjYf9fQp6  -e  1e-20  -o  >>../20021010/BLAST/Xbn.out 




my @params = ('program' => 'blastn','database' => 'X','outfile' => '../20021010/BLAST/Xbn.out','_READMETHOD' => 'Blast','e' => '1e-20');
my $factory = Bio::Tools::Run::StandAloneBlast->new(@params);
my $Seq_in = Bio::SeqIO->new (-file => $ARGV[0], -format => 'fasta');

while ( my $query = $Seq_in->next_seq() ) {
    my $blast_report = $factory->blastall($query);
    my $result = $blast_report->next_result;
    while( my $hit = $result->next_hit()) { 
	parse report....
    }
}


Suggestions?

Charles