[bip] pygr: reload AnnotationDB

Diane Trout diane at caltech.edu
Wed Mar 26 14:34:46 PDT 2008


>> perfect. thanks.
> so to make it painfully obvious, so i dont forget:
> given an AnnotationDB, adb, that was created with BlastDB('/ 
> somewhere/asdf.fa'),
> one should access that BlastDB with
>
>>>> seq_db = adb.seqDict.prefixDict['asdf']

Yes.

If you attach annotations they might show up as 'noname1'

 >>> kg = cnestedlist.NLMSA('db/knownGene.txt','r')
 >>> kg.seqDict.prefixDict.keys()
['hg18', 'noname2', 'noname1']
 >>> kg.seqDict.prefixDict['hg18'].keys()
['chr6_random',
  'chr19_random',
  'chr8_random',
  'chr5_h2_hap1',

diane



More information about the biology-in-python mailing list