def test_get_val2keys(self):
v2ks = lcommon.get_val2keys(self.dict3)
for v, ks in v2ks.iteritems():
v2ks[v] = sorted(ks)
dict3 = {1: ['sam12'], 2: ['sam12', 'sam23'], 3: ['sam23']}
self.assertEqual(v2ks, dict3)
self.assertRaises(lcommon.InputError, lcommon.get_val2key_1to1, self.dict3)
def run(self):
c2s = pickle.load(gzip.open(self.infile, 'rb'))
filter_c2sams = filter_clones(c2s, self.sams, self.minsam, self.maxsam)
sam2clones = lcommon.get_val2keys(filter_c2sams)
for sam, clones in sam2clones.iteritems():
infile = os.path.join(self.db_dir, sam, self.v)
outfile = os.path.join(self.outdir, sam, self.v)
self.addChildTarget(OverlapClones_Sample(infile, outfile, clones))