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