[bip] pygr: reload AnnotationDB

[Brent Pedersen]

On Wed, Mar 26, 2008 at 12:28 PM, Diane Trout <diane at ghic.org> wrote:
> >
>  > def read():
>  >
>  # My read code is using cnestedlist.NLMSA(path, 'r')
>  >    alignment = cnestedlist.NLMSA('/tmp/test')
>  >
>  # i then use
>  #   genome = alignment.seqDict.prefixDict[genome_name]
>  # instead of
>
> >     seq_db = seqdb.BlastDB('seq.fa')
>  # for my queries
>
>  (I looked at the keys in prefixDict to figure out what the right
>  genome_name was.)
>
>  Diane
>

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']