pipeline.py 文件源码

def run_ec_on_bin(*args):
    cloneset, mismatch_rate, confidence, max_Q = args
    try:
        gc.disable()
        seqlen = len(next(iter(cloneset)).seq)
        logger.info("Starting QMerge on cloneset: seqlen: %s, \
#clones: %s, #sequences: %s, #bases: %s, #mutation_count: %s"%(seqlen,
len(cloneset), cloneset.sequence_count, cloneset.base_count,
cloneset.mutation_count))

        cloneset = run_qmerge_on_bin(cloneset, mismatch_rate, confidence,
                max_Q)
        logger.info("Starting IMerge on cloneset: seqlen: %s, \
#clones: %s, #sequences: %s, #bases: %s, #mutation_count: %s"%(seqlen,
len(cloneset), cloneset.sequence_count, cloneset.base_count,
cloneset.mutation_count))
        cloneset = run_imerge_on_bin(cloneset, mismatch_rate, confidence)
    finally:
        gc.enable()
    return cloneset, mismatch_rate