Пример #1
0
def main():
    os.chdir('./bin')
    dirc = sys.argv[1]
    cfg = load_mirquant_config_file(sys.argv[3])
    name = os.path.basename(dirc.split('./IntermediateFiles/')[0])
    mirquant_output = sample_output_paths(cfg['paths']['output'], name)
    logName = '{}_collectRes.log'.format(os.path.basename(dirc))
    initiate_logging('{}/collect_results_logs/'.format(mirquant_output['log']),
                     logName)
    bt = '{}.results'.format(sys.argv[1])
    res_li = resource_paths(cfg['parameters']['species'], cfg['paths'],
                            cfg['parameters'])
    mir_file, TRNAfile = res_li[1], res_li[5]
    outDir = os.path.dirname(bt)
    filesLib = glob.glob('{}/../*LIB.fa'.format(os.path.dirname(dirc)))[0]
    print_run_info(dirc, bt, cfg['parameters']['species'], outDir, filesLib)
    btWins, EM = result_file_dict(bt, dirc)
    mirs = get_miR_info(mir_file)
    counters = initialize_counters()

    bedFile = []
    for res in sorted(btWins):
        bedFile, counters = mainChunk(res, counters, bedFile, dirc, EM,
                                      TRNAfile, mirs, outDir, mirquant_output)

    write_summary_to_log(counters)
    write_shrimp_results_bed(sys.argv[1], bedFile)
    os.system('rm {}'.format(sys.argv[2]))
Пример #2
0
def main(arg):
    start = time.time()
    os.chdir('./bin')
    cfg = load_mirquant_config_file(arg.conf)
    scfg = load_sys_config_file(arg.conf)
    job = build_job(scfg['job'])
    t_job = build_job(scfg['job_threaded'])

    dr, dr_i, fi_base = set_up_output_folder(arg.sample,
                                             cfg['paths']['output'])
    out_di = sample_output_paths(cfg['paths']['output'], fi_base)
    initiate_logging(out_di['log'], 'chainSubmission.log')
    res_li = resource_paths(cfg['parameters']['species'], cfg['paths'],
                            cfg['parameters'])
    tRNA, tmRNA, BI = res_li[4], res_li[3], res_li[7]
    define_input_varibles(cfg)

    lib = set_lib(dr_i, fi_base)
    MINrna, MAXrna = cutadapt_cmd(arg.sample, lib, cfg['cutadapt'],
                                  out_di['log'], arg.conf)
    logging.info('cutadapt = {}'.format(time.time() - start))
    separate_by_read_length(MINrna, MAXrna, lib, out_di['output'])
    for length in range(MINrna, MAXrna + 1):
        fi = '{}_{}.fq'.format(lib, length)
        bowtie(fi, length, BI, cfg['bowtie'])
    logging.info('bowtie = {}'.format(time.time() - start))
    window_creation(MINrna, MAXrna, lib, BI, tRNA, tmRNA)
    logging.info('window generation = {}'.format(time.time() - start))
    mapping_statistics(arg.sample, lib, dr, fi_base)
    run_shrimp_alignment(MINrna, MAXrna, lib, out_di['log'], out_di['temp'],
                         t_job, arg.conf)
    bt_postProcEM.main('{}_merge.bed'.format(lib), '{}_allGS.bed'.format(lib),
                       out_di['temp'])
    reduce_shrimp_res(out_di['temp'], dr_i, job)
    logging.info('finish = {}'.format(time.time() - start))
Пример #3
0
def main():
    os.chdir('./bin') 
    dirc = sys.argv[1]
    cfg = load_mirquant_config_file(sys.argv[3])
    name = os.path.basename(dirc.split('./IntermediateFiles/')[0])
    mirquant_output = sample_output_paths(cfg['paths']['output'], name)
    logName  = '{}_collectRes.log'.format(os.path.basename(dirc))
    initiate_logging('{}/collect_results_logs/'.format(mirquant_output['log']), logName)
    bt = '{}.results'.format(sys.argv[1])
    res_li = resource_paths(cfg['parameters']['species'], cfg['paths'], cfg['parameters'])
    mir_file, TRNAfile = res_li[1], res_li[5]
    outDir = os.path.dirname(bt)
    filesLib = glob.glob('{}/../*LIB.fa'.format(os.path.dirname(dirc)))[0] 
    print_run_info(dirc, bt, cfg['parameters']['species'], outDir, filesLib)
    btWins, EM = result_file_dict(bt, dirc)
    mirs =  get_miR_info(mir_file)
    counters = initialize_counters()

    bedFile = [] 
    for res in sorted(btWins):
        bedFile, counters = mainChunk(res, counters, bedFile, dirc, EM, TRNAfile, mirs, outDir, mirquant_output)

    write_summary_to_log(counters)
    write_shrimp_results_bed(sys.argv[1], bedFile)
    os.system('rm {}'.format(sys.argv[2])) 
Пример #4
0
def main(arg):
    os.chdir('./bin')
    cfg = load_mirquant_config_file(arg.conf)
    scfg = load_sys_config_file(arg.conf)
    job = build_job(scfg['job'])

    dr, dr_i, fi_base = set_up_output_folder(arg.sample, cfg['paths']['output'])
    out_di = sample_output_paths(cfg['paths']['output'], fi_base)
    initiate_logging(out_di['log'], 'chainSubmission.log')
    res_li = resource_paths(cfg['parameters']['species'], cfg['paths'], cfg['parameters'])
    tRNA, tmRNA, BI = res_li[4], res_li[3], res_li[7]
    define_input_varibles(cfg)

    generate_adapter_files.main(arg.sample, out_di['log'], arg.conf)
    MAXrna = get_maxRNA_length(arg.sample, cfg['cutadapt'])
    lib = set_lib(dr_i, fi_base)
    MINrna = cutadapt_cmd(arg.sample, lib, cfg['cutadapt'])
    separate_by_read_length(MINrna, MAXrna, lib, out_di['output'])
    for length in range(MINrna, MAXrna + 1):
        fi = '{}_{}.fq'.format(lib, length)
        bowtie(fi, length, BI, cfg['bowtie'])
    window_creation(MINrna, MAXrna, lib, BI, tRNA, tmRNA)
    mapping_statistics(arg.sample, lib, dr, fi_base)
    run_shrimp_alignment(MINrna, MAXrna, lib, out_di['log'], out_di['temp'], job, arg.conf)
    bt_postProcEM.main('{}_merge.bed'.format(lib), '{}_allGS.bed'.format(lib), out_di['temp'])
    reduce_shrimp_res(out_di['temp'], dr_i, job)
Пример #5
0
def main(size, lib, base, conf, log_dir):
    cfg = load_mirquant_config_file(conf)
    initiate_logging(log_dir, "SHRiMP_{}.log".format(size))

    reads = "{}{}.noHit".format(base, size)
    working_dir = os.getcwd()

    make_out_dir(lib, size)
    seedList = make_base_seed(size)
    group_Name = shrimp_submission("pypath", cfg["shrimp"], seedList, lib, log_dir, reads)
    shrimp_postProcGS.main(conf, os.getcwd())
Пример #6
0
def main(size, lib, base, conf, log_dir):
    cfg = load_mirquant_config_file(conf)
    initiate_logging(log_dir, 'SHRiMP_{}.log'.format(size))

    reads = '{}{}.noHit'.format(base, size)
    working_dir = os.getcwd()

    make_out_dir(lib, size)
    seedList = make_base_seed(size)
    group_Name = shrimp_submission('pypath', cfg['shrimp'], seedList, lib,
                                   log_dir, reads)
    shrimp_postProcGS.main(conf, os.getcwd())