[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