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