[Bioperl-l] retrieven ids

Luba Pardo lubapardo at gmail.com
Thu Mar 1 07:48:27 UTC 2007


Thank you very much.
L. Pardo

On 28/02/07, Dave Messina <dmessina at wustl.edu> wrote:
>
> Whenever I'm unsure of how to do something, I first look to see if one of
> the  HOWTOs on bioperl.org covers it. In this case, the Features HOWTO has
> example code which I think will do what you want.
>
> Genbank records typically have the coding sequence of a protein as a
> feature, so I would do something like:
>
> - use the RefSeq protein IDs to query Entrez and get back the Genbank
> records.
>
> - read the Features HOWTO to refresh my memory on the syntax for grabbing
> features.
>
> That HOWTO is at:
> http://www.bioperl.org/wiki/HOWTO:Feature-Annotation
>
> - whip up a little script to loop through the Genbank records one at a
> time with SeqIO and pull out the cDNA sequence features.
>
>
> Dave
>
>
>



More information about the Bioperl-l mailing list