[Bioperl-l] retrieving by acc from a local multifasta file
paulo.david at netvisao.pt
Fri Aug 6 07:47:42 EDT 2004
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
( for details check
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 )
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?
More information about the Bioperl-l