def test_schema(self):
"Test schema"
sp_hbb1 = testutil.datafile('sp_hbb1')
sp2 = seqdb.BlastDB(sp_hbb1)
sp2.__doc__ = 'another sp'
worldbase.Bio.Seq.sp2 = sp2
sp = worldbase.Bio.Seq.Swissprot.sp42()
m = mapping.Mapping(sourceDB=sp, targetDB=sp2)
m.__doc__ = 'sp -> sp2'
worldbase.Bio.Seq.testmap = m
worldbase.schema.Bio.Seq.testmap = metabase.OneToManyRelation(sp, sp2)
worldbase.commit()
worldbase.clear_cache()
sp3 = seqdb.BlastDB(sp_hbb1)
sp3.__doc__ = 'sp number 3'
worldbase.Bio.Seq.sp3 = sp3
sp2 = worldbase.Bio.Seq.sp2()
m = mapping.Mapping(sourceDB=sp3, targetDB=sp2)
m.__doc__ = 'sp3 -> sp2'
worldbase.Bio.Seq.testmap2 = m
worldbase.schema.Bio.Seq.testmap2 = metabase.OneToManyRelation(
sp3, sp2)
l = worldbase._mdb.resourceCache.keys()
l.sort()
assert l == ['Bio.Seq.sp2', 'Bio.Seq.sp3', 'Bio.Seq.testmap2']
worldbase.commit()
g = worldbase._mdb.writer.storage.graph
expected = set([
'Bio.Annotation.annoDB', 'Bio.Seq.Swissprot.sp42', 'Bio.Seq.sp2',
'Bio.Seq.sp3'
])
found = set(g.keys())
self.EQ(len(expected - found), 0)
def test_schema(self):
"Test schema"
sp_hbb1 = testutil.datafile('sp_hbb1')
sp2 = seqdb.BlastDB(sp_hbb1)
sp2.__doc__ = 'another sp'
pygr.Data.Bio.Seq.sp2 = sp2
sp = pygr.Data.Bio.Seq.Swissprot.sp42()
m = mapping.Mapping(sourceDB=sp, targetDB=sp2)
m.__doc__ = 'sp -> sp2'
pygr.Data.Bio.Seq.testmap = m
pygr.Data.schema.Bio.Seq.testmap = pygr.Data.OneToManyRelation(sp, sp2)
pygr.Data.save()
pygr.Data.clear_cache()
sp3 = seqdb.BlastDB(sp_hbb1)
sp3.__doc__ = 'sp number 3'
pygr.Data.Bio.Seq.sp3 = sp3
sp2 = pygr.Data.Bio.Seq.sp2()
m = mapping.Mapping(sourceDB=sp3, targetDB=sp2)
m.__doc__ = 'sp3 -> sp2'
pygr.Data.Bio.Seq.testmap2 = m
pygr.Data.schema.Bio.Seq.testmap2 = pygr.Data.OneToManyRelation(sp3,
sp2)
# List all cached resources.
l = pygr.Data.getResource.resourceCache.keys()
l.sort()
assert l == ['Bio.Seq.sp2', 'Bio.Seq.sp3', 'Bio.Seq.testmap2']
pygr.Data.save()
g = pygr.Data.getResource.writer.storage.graph
expected = set(['Bio.Annotation.annoDB',
'Bio.Seq.Swissprot.sp42', 'Bio.Seq.sp2', 'Bio.Seq.sp3'])
found = set(g.keys())
self.EQ(len(expected - found), 0)
def populate_swissprot():
"Populate the current worldbase with swissprot data"
# build BlastDB out of the sequences
sp_hbb1 = testutil.datafile('sp_hbb1')
sp = seqdb.BlastDB(sp_hbb1)
sp.__doc__ = 'little swissprot'
worldbase.Bio.Seq.Swissprot.sp42 = sp
# also store a fragment
hbb = sp['HBB1_TORMA']
ival = hbb[10:35]
ival.__doc__ = 'fragment'
worldbase.Bio.Seq.frag = ival
# build a mapping to itself
m = mapping.Mapping(sourceDB=sp, targetDB=sp)
trypsin = sp['PRCA_ANAVA']
m[hbb] = trypsin
m.__doc__ = 'map sp to itself'
worldbase.Bio.Seq.spmap = m
# create an annotation database and bind as exons attribute
worldbase.schema.Bio.Seq.spmap = metabase.OneToManyRelation(
sp, sp, bindAttrs=('buddy', ))
annoDB = seqdb.AnnotationDB({1: ('HBB1_TORMA', 10, 50)},
sp,
sliceAttrDict=dict(id=0, start=1, stop=2))
exon = annoDB[1]
# generate the names where these will be stored
tempdir = testutil.TempDir('exonAnnot')
filename = tempdir.subfile('cnested')
nlmsa = cnestedlist.NLMSA(filename,
'w',
pairwiseMode=True,
bidirectional=False)
nlmsa.addAnnotation(exon)
nlmsa.build()
annoDB.__doc__ = 'a little annotation db'
nlmsa.__doc__ = 'a little map'
worldbase.Bio.Annotation.annoDB = annoDB
worldbase.Bio.Annotation.map = nlmsa
worldbase.schema.Bio.Annotation.map = \
metabase.ManyToManyRelation(sp, annoDB, bindAttrs=('exons', ))