[Bioperl-l] get html-output on my prompt
ak2083
ankh.egypt.public at googlemail.com
Sun Jun 24 13:44:49 UTC 2012
Hello,
i try to blast a protein sequence (in sequence.fasta) against the blast-db.
If i parse this
with the command 'perl file.pl', then the programm displayed html-code on my
prompt and the out-file
was not created. What could i do so that the results write into the
out-file?
my $report = $blast_obj->submit_blast('sequence.fasta');
my $x = 0;
print STDERR "waiting..." if ($x > 0);
while (my @rids = $blast_obj->each_rid() && $x == 5)
{
foreach my $rid (@rids)
{
my $rc = $blast_obj->retrieve_blast($rid);
if (!ref($rc))
{
if ($rc < 0)
{ $blast_obj->remove_rid($rid); }
print STDERR "." if ($x > 0);
sleep 5;
}
else
{
my $result = $rc->next_result();
my $filename = $result->query_name(), '\.out';
$blast_obj->save_output($filename);
$blast_obj->remove_rid($rid);
print "\nQuery Name: ", $result->query_name();
while (my $hit = $result->next_hit())
{
next unless ($x > 0);
print "\thit name is ", $hit->name, "\n";
while (my $hsp = $hit->next_hsp())
{
print "\t\tscore is ", $hsp->score(), "\n";
}
}
}
}
}
--
View this message in context: http://old.nabble.com/get-html-output-on-my-prompt-tp34061631p34061631.html
Sent from the Perl - Bioperl-L mailing list archive at Nabble.com.
More information about the Bioperl-l
mailing list