def __addNibs(self, paths):
"add nib sequences to object"
# /cluster/data/hg17/nib/chrX.nib chrX 154824264
lines = procOps.callProcLines(["nibSize"] + paths)
for line in lines:
row = line.split("\t")
self[row[1]] = GenomeSeq(self, row[0], row[1], int(row[2]))
def __addTwoBit(self, path):
"add twoBit sequences to object"
self.genomeDb = path[0]
# chrX 154824264
lines = procOps.callProcLines(["twoBitInfo", path, "stdout"])
for line in lines:
row = line.split("\t")
self[row[0]] = GenomeSeq(self, path, row[0], int(row[1]))