[Bioperl-l] how to get annotations (especially ensembl IDs) for a list of genes?

Sendu Bala bix at sendu.me.uk
Sat Jul 29 07:25:38 UTC 2006


zhihua li wrote:
> Hi all,
> 
> I have a list of like 300 genes (actually their refseq IDs).  Now I 
> wanna get more information (annotations) for each of the genes. 
> Speficially, I want a mapping of the refseq IDs to Ensembl gene IDs.
> 
> I know how to do it through a web page. But I'm wondering if I can also 
> do it via bioperl

One possible way is to use the Ensembl perl API:
http://www.ensembl.org/info/software/core/core_tutorial.html

You'd get a gene or transcript adapator and use
fetch_all_by_external_name() iirc.

I'm aware that not every entrez id can be mapped that way, but perhaps
most if not all refseqs will work.



More information about the Bioperl-l mailing list