[Bioperl-l] Bplite
Jason Eric Stajich
jason@cgt.mc.duke.edu
Sat, 24 Nov 2001 10:14:21 -0500 (EST)
a) Case counts!
line should read
use Bio::Tools::BPlite
instead of
use Bio::Tools::Bplite
b) You aren't actually calling print anywhere in your script so you won't get
any results. The code in the SYNOPSIS is meant to be an example of how
the API is used, but you will need to build up a script so that it
prints the data in the way you expect to see it.
c) You can either specify the input from STDIN by running the script like:
% perl mybplitescript.pl < atop.txt
or change your script so that it reads
my $report = new Bio::Tools::BPlite(-file => "atop.txt");
Best of luck.
On Sat, 24 Nov 2001, nabula easter wrote:
> Dear list,
> I have a blast report (atop1.txt) and i am trying to
> parse this blast report for sequence name, its high
> E-value. for this i tried Bplite module.
>
> As per bioperl docs, i used this following script, but
> unfortunately i did not get any result. i know there
> is a mistake in this, but i am unable to find that.
> Please help me.
> my blast report: atop1.txt
> Script:
>
> use Bio::Tools::Bplite;
> $report = new Bplite (-fh=>\*STDIN);
> $report->query;
> while (my $sbjct = $report->nextSbjct){
> $sbjct->name;
> while (my $hsp = $sbjct->nextHS) {$hsp->score;}
> }
>
>
> thank you
> NE
>
> __________________________________________________
> Do You Yahoo!?
> Yahoo! GeoCities - quick and easy web site hosting, just $8.95/month.
> http://geocities.yahoo.com/ps/info1
> _______________________________________________
> Bioperl-l mailing list
> Bioperl-l@bioperl.org
> http://bioperl.org/mailman/listinfo/bioperl-l
>
--
Jason Stajich
Duke University
jason@cgt.mc.duke.edu