[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