def batch_scan(points, xdist=20, ydist=20, N=5):
    chr_pair_points = group_hits(points)

    clusters = []
    for chr_pair in sorted(chr_pair_points.keys()):
        points = chr_pair_points[chr_pair]
        clusters.extend(synteny_scan(points, xdist, ydist, N))

    return clusters
runs synteny_scan() per chromosome pair
dataset/ETHPy150Open tanghaibao/jcvi/compara/synteny.py/batch_scan