[Bioperl-l] Remote Blast and Report

Scott Markel SMarkel at accelrys.com
Thu Feb 26 02:04:29 UTC 2009


O'car,

There's a polling mechanism you need to use.  See the example in the
Bio::Tools::Run::RemoteBlast module.  Start looking around line 60. 

Scott

Scott Markel, Ph.D.
Principal Bioinformatics Architect  email:  smarkel at accelrys.com
Accelrys (SciTegic R&D)             mobile: +1 858 205 3653
10188 Telesis Court, Suite 100      voice:  +1 858 799 5603
San Diego, CA 92121                 fax:    +1 858 799 5222
USA                                 web:    http://www.accelrys.com

http://www.linkedin.com/in/smarkel
Vice President, Board of Directors:
    International Society for Computational Biology
Co-chair: ISCB Publications Committee
Associate Editor: PLoS Computational Biology
Editorial Board: Briefings in Bioinformatics


> -----Original Message-----
> From: bioperl-l-bounces at lists.open-bio.org [mailto:bioperl-l-
> bounces at lists.open-bio.org] On Behalf Of Ocar Campos
> Sent: Wednesday, 25 February 2009 4:04 PM
> To: Bioperl Mailing List.
> Subject: [Bioperl-l] Remote Blast and Report
> 
> Hello:
> 
>    I'm working in a script to remote blast a file with some sequences, I
> already got the part of sending the query to blast, but I do not get the
> idea of how retrieve a txt report, I mean, like the one you get by running
> a
> blast via web and you can read in a plane text editor.
> 
> This is what I've done so far:
> 
> 
> use Bio::Tools::Run::RemoteBlast;
> use Bio::SearchIO;
> 
> $prog = "tblastx";
> $db = "nr";
> $e_val = "1e-10";
> $remoteBlast = Bio::Tools::Run::RemoteBlast->new(-prog => $prog,
>                          -data => $db,
>                          -expect => $e_val
>                          -readmethod => 'Blast');
> 
> #I select the file to make que query and do the blast.
> $infile = 'file.input.fasta';
> $r = $remoteBlast->submit_blast($infile);
> 
> #this should be the report i get.
> $outfile = 'got.output';
> 
> further than this I've tried some things but none of them work, anybody
> who
> could give an idea of how retrieving the plane text reports please?
> 
> Cheers.
> 
> O'car
> _______________________________________________
> Bioperl-l mailing list
> Bioperl-l at lists.open-bio.org
> http://lists.open-bio.org/mailman/listinfo/bioperl-l




More information about the Bioperl-l mailing list