[Bioperl-l] Caching sequences

Sendu Bala bix at sendu.me.uk
Fri Oct 27 08:35:13 EDT 2006


Nathan S. Haigh wrote:
> I have a script that is capable of downloading sequences from GenBank
> based on GI numbers. I retrieve them if fasta format in order to save
> bandwidth, but I'd like to take this one step further and cache the
> sequences in case the user want to rerun the script using some of the
> GI's they used previously.
> 
> Does anyone have any guidance on how best to do this?

You'd probably write the sequences out in some suitable format and 
access them via Bio::Index

Or, I'm sure bioperl-db excels at this kind of thing, but is a little 
more involved if this is only a simple situation.


More information about the Bioperl-l mailing list