def call1 (cvsfile, write_func, *args ):
    argv = args[0]
    if len(argv) < 5:
        sys.stderr.write("Need 5 extra arguments for 'call', options <loc column> <score column> <cutoff> <min length> <max gap>\ne.g. command: <0> <1> <0.5> <10000> <2000>, means to use the first column as genome coordinations to call enriched regions from the second column, the threshold to call enriched region is 0.5, the minimum length of region is 10k, and the maximum gap to link two nearby regions is 2k.\n")
        sys.exit()
    cor_column = cvsfile.fieldnames[int(argv[0])]
    var_column = cvsfile.fieldnames[int(argv[1])]
    cutoff = float(argv[2])
    min_len = int(argv[3])
    max_gap = int(argv[4])
    wtrack = WigTrackI()
    add_func = wtrack.add_loc
    for l in cvsfile:
        cor = l.setdefault(cor_column,None)
        var = l.setdefault(var_column,None)        
        if cor and var and cor != "NA" and var != "NA":
            (chrom,start,end) = cor.split(".")
            add_func(chrom,int(start),float(var))
    wtrack.span = int(end)-int(start)
    write_func( "# regions called from %s:%s\n" % (argv[1],var_column) )
    
    bpeaks = wtrack.call_peaks(cutoff=cutoff,min_length=min_len,max_gap=max_gap)
    write_func( bpeaks.tobed() )

# ------------------------------------
# Classes
# ------------------------------------

# ------------------------------------
# Main function
# ------------------------------------
Build wheels.
dataset/ETHPy150Open taoliu/taolib/Scripts/ce_histone_matrix.py/call1--crossover with--dataset/ETHPy150Open GeekTrainer/Flask/Work/Trivia - Module 5/env/Lib/site-packages/pip/wheel.py/WheelBuilder.build