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

[Paulo Almeida]

Maria,

If your fasta file is in a format that can be parsed by Bioperl, I think 
you just need to load the sequences into a Bio::Seq object, using 
Bio::SeqIO.

( for details check 
http://search.cpan.org/~birney/bioperl-1.4/Bio/SeqIO.pm )

Then you can get the accesion number of each sequence with 
$seq->accesion_number and only retrieve those that match your list.

( for details check http://search.cpan.org/~birney/bioperl-1.4/Bio/Seq.pm )

Good luck,
Paulo

Maria Persico wrote:

>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?
>