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