[Bioperl-l] retrieving by acc from a local multifasta file

Remo Sanges sanges at biogem.it
Fri Aug 6 07:06:45 EDT 2004

On Aug 6, 2004, at 12:50 PM, Maria Persico wrote:

> Hi All,
> This may be a stupid problem but for me it's something difficult:
> I have a list of swissprot accessions(my_acc) and I want to extract 
> from
> uniprot_sprot.fasta only sequences of my list.
> How can do this with bioperl?

Hi Maria,

you could create an index of your fasta file
and then retrieve sequences from it.

Read the tutorial section 'III.1.2 Indexing and accessing local 

and the documentation of Bio::Index::Fasta module



More information about the Bioperl-l mailing list