[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