[Bioperl-l] SearchIO: write/read database

Dave Messina David.Messina at sbc.su.se
Tue May 20 17:10:26 UTC 2008


Echoing what Sendu said, the most straightforward way to read output from a
program without writing to a file first is to open a pipe from the program
as a filehandle.

e.g.
my @args = ('/usr/bin/blastp', '/path/to/my/blastdb',
'/path/to/my/blastquery');
open (my $blast_fh, '-|', @args) or die "couldn't open blast stream";

Now you can use that filehandle to read in data and manipulate just as if
you had read it in from a file (except of course since it's a stream you
can't rewind).


Dave



More information about the Bioperl-l mailing list