Esempio n. 1
0
def iterator_filter_overlapping_target( psls, options ):

    ninput, noutput, ndiscarded = 0, 0, 0    
    for block in Blat.iterator_target_overlap( psls, options.threshold_merge_distance ):
        l = len(block)
        ninput += l
        if l > 1: 
            ndiscarded += l
        else:
            yield block[0]
            noutput += 1

    E.info( "iterator_filter_overlapping_target: ninput=%i, noutput=%i, ndiscarded=%i" % (ninput, noutput,ndiscarded) )
Esempio n. 2
0
def iterator_filter_overlapping_target(psls, options):

    ninput, noutput, ndiscarded = 0, 0, 0
    for block in Blat.iterator_target_overlap(
            psls, options.threshold_merge_distance):
        l = len(block)
        ninput += l
        if l > 1:
            ndiscarded += l
        else:
            yield block[0]
            noutput += 1

    E.info("iterator_filter_overlapping_target: ninput=%i, noutput=%i, "
           "ndiscarded=%i" % (ninput, noutput, ndiscarded))