get_reads.py 文件源码

def largest_overlap(items):
    positions_i = collections.defaultdict(int)
    positions_j = collections.defaultdict(int)
    for disc in items:
        start,end = signi(disc)
        for pos in range(start,end):
            positions_i[pos] += 1
        start,end = signj(disc)
        for pos in range(start,end):
            positions_j[pos] += 1
    i = 0
    overlap_i = 0
    for pos,overlap in iteritems(positions_i):
        if overlap > overlap_i:
            overlap_i = overlap
            i = pos
    j = 0
    overlap_j = 0
    for pos,overlap in iteritems(positions_j):
        if overlap > overlap_j:
            overlap_j = overlap
            j = pos
    return i,j,overlap_i,overlap_j