Esempio n. 1
0
 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]))
Esempio n. 2
0
 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]))