[Bioperl-l] dna seq for an AA entry

Ewan Birney birney@ebi.ac.uk
Sun, 27 May 2001 10:50:13 +0100 (BST)


On Sat, 26 May 2001, Hilmar Lapp wrote:

> "Saurabh D. Patel" wrote:
> > 
> > Now that my bioperl setup is working, I have another (newbie) question:
> > given a genbank accession number/amino acid sequence, is there any easy
> > way to get corresponding dna sequence?  I can think of a kludge where I
> 
> If there is a crosslink somewhere in the feature table (which I
> think there is in Genbank, but I'm not sure particularly where off
> the top of my head, and swissprot also has crosslinks), pull out
> that piece of annotation (e.g. scan the attached features) and
> retrieve the sequence entry using Bio::DB::GenBank or
> Bio::DB::EMBL (or Bio::DB::Universal :-))

I find this easier to do in the EMBL entries because the DR lines are more
explicit - they are in the $seq->annotation->each_DBLink() inside bioperl.

In genbank you will need to parse the seqfeatures and find the db_xref
tag.



> 
> 	Hilmar
> 
> -- 
> -----------------------------------------------------------------
> Hilmar Lapp                              email: hilmarl@yahoo.com
> GNF, San Diego, Ca. 92122                phone: +1 858 812 1757
> -----------------------------------------------------------------
> _______________________________________________
> Bioperl-l mailing list
> Bioperl-l@bioperl.org
> http://bioperl.org/mailman/listinfo/bioperl-l
> 

-----------------------------------------------------------------
Ewan Birney. Mobile: +44 (0)7970 151230, Work: +44 1223 494420
<birney@ebi.ac.uk>. 
-----------------------------------------------------------------