def test_bedarcsdb_tiles(): """ Retrieve a 1D tile from a 2d file """ filename = op.join("data", "arrowhead_domains_short.txt.multires.db") pos = 2 hgbad.tiles(filename, ["b.2.{}".format(pos)])
def test_bedarcsdb_tiles(): ''' Retrieve a 1D tile from a 2d file ''' filename = op.join('data', 'arrowhead_domains_short.txt.multires.db') pos = 2 hgbad.tiles(filename, ["b.2.{}".format(pos)])
def get_counts(filename, zoom, pos): data = hgbad.tiles(filename, ["b.{}.{}".format(zoom, pos)]) # print('data:', len(data[0][1][pos])) counts = col.defaultdict(int) for i, d in enumerate(data[0][1][pos]): chr1, chr2 = sorted([d['fields'][0], d['fields'][3]]) #print('d:', d) counts[(chr1, chr2)] += 1 return counts