[Bioperl-l] deprecated blast parse?
Wiepert, Mathieu
Wiepert.Mathieu@mayo.edu
Tue, 26 Feb 2002 12:22:57 -0600
Hi,
I have code (surreptitiously copied from Jason's examples) that is giving me
a deprecation error. I have a refresh from this morning of bioperl-live.
Method subject deprecated: use hit() instead
STACK Bio::SeqFeature::SimilarityPair::subject
/home/mxw02/bioperl_latest/lib/site_perl/5.6.0/Bio/SeqFeature/SimilarityPair
.pm:342
STACK toplevel bpblast.pl:32
453 dbj|D53816.1|D53816
Code snippet is
my @params = ('outfile' => 'seq.out', 'program' => 'tblastn', 'database' =>
'est_human');
my $factory = Bio::Tools::Run::StandAloneBlast->new(@params);
my $str = Bio::SeqIO->new(-file=>'seq.fa' , '-format' => 'Fasta' );
my $input = $str->next_seq();
my $blast_report = $factory->blastall($input);
my $searchio = new Bio::SearchIO(-format => 'blast',
# -result_factory => $blast_report);
-file => 'seq.out');
while( my $result = $searchio->next_result ) {
while( my $hit = $result->next_hit ) {
while( my $hsp = $hit->next_hsp ) {
print $hsp->score . " " . $hsp->subject->seqname .
"\n";
}
}
}
I am sure i did something wrong...
-Mat