[Bioperl-l] UniGene

[Sean Davis]

Badr,

The simplest way is to go to the ftp site for unigene:

ftp://ftp.ncbi.nih.gov/repository/UniGene

Get the file for the organism you are interested that ends in .gb_cid_lid. 
Just choose the file for the organism of interest.  For example, the first 
few lines of Cre.gb_cid_lid look like:

AY171232 3219 -
AY171231 3003 -
AY171230 4370 -
AY171229 2671 -
AY184800 2793 -
AY184799 6486 -
AY184798 206 -
AY184797 3607 -
AY184796 2281 -
AY177787 2380 -
AY212923 4329 -
AB091079 3370 -

The first column contains genbank accessions.  The second contains unigene 
cluster ids.  So, for the first genbank AY171232, the unigene accession is 
Cre.3219 (You have to prepend the 2- or 3- letter organism code).  You can 
read these into a hash that is keyed on the genbank accession with a value 
that is the reference to an array of unigene cluster IDs for each genbank (a 
genbank can be in multiple unigene clusters).

Hope that helps.
Sean

----- Original Message ----- 
From: "badr al-daihani" <aldaihani at hotmail.co.uk>
To: <bioperl-l at bioperl.org>
Sent: Sunday, April 17, 2005 5:04 PM
Subject: [Bioperl-l] UniGene


> Hi folks
>
> would you please tell me how to retrieve the unigene number of a gene 
> (UniGene)
> knowing the GenBankaccession number ?
>
>
> Best regards
>
> Badr
>
> _________________________________________________________________
> It's fast, it's easy and it's free. Get MSN Messenger today! 
> http://www.msn.co.uk/messenger
>
> _______________________________________________
> Bioperl-l mailing list
> Bioperl-l at portal.open-bio.org
> http://portal.open-bio.org/mailman/listinfo/bioperl-l
>