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]))
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))
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]))
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)
def main(conf):
os.chdir('./bin')
cfg = load_mirquant_config_file(conf)
samples = return_sample_results_directories(cfg['paths']['project'])
for sample in samples:
samp_name = os.path.basename(sample[:-1])
out_di = sample_output_paths(cfg['paths']['output'], samp_name)
run_summary2Tab_clust(cfg['paths'], cfg['parameters']['species'], sample, conf)
move_files_to_out_dir(out_di, sample, samp_name)
write_summary_table(sample)
def main(conf):
os.chdir('./bin')
cfg = load_mirquant_config_file(conf)
samples = return_sample_results_directories(cfg['paths']['project'])
for sample in samples:
print '\nProcessing sample {}...'.format(sample)
samp_name = os.path.basename(sample[:-1])
out_di = sample_output_paths(cfg['paths']['output'], samp_name)
run_summary2Tab_clust(cfg['paths'], cfg['parameters']['species'],
sample, conf)
move_files_to_out_dir(out_di, sample, samp_name)
write_summary_table(sample)
def main(conf, shrimp_dir):
check_input()
cfg = load_mirquant_config_file(conf)
res_li = resource_paths(cfg['parameters']['species'], cfg['paths'], cfg['parameters'])
sample = os.path.basename(shrimp_dir.split('./IntermediateFiles/')[0])
out_di = sample_output_paths(cfg['paths']['output'], sample)
logging.info('\n\n### Processing SHRiMP results ###\n')
mir_fi = res_li[1]
mirList, mirStrand = load_mir_info(mir_fi)
tagCount, hits, maps = load_SHRiMP_res(shrimp_dir)
hits, tags, pCount = get_best_alignments(hits, maps)
write_processed_shrimp_output(hits, tags, mirList, mirStrand, shrimp_dir)
remove_temp_file(os.path.basename(shrimp_dir).split('_')[1], out_di['temp'])
logging.info('Total SHRiMP alignments = {}; proportional count = {}'.format(tagCount, pCount))
def main(conf, shrimp_dir):
print 'Shrimp post-processing'
check_input()
cfg = load_mirquant_config_file(conf)
res_li = resource_paths(cfg['parameters']['species'], cfg['paths'], cfg['parameters'])
sample = os.path.basename(shrimp_dir.split('./IntermediateFiles/')[0])
out_di = sample_output_paths(cfg['paths']['output'], sample)
logging.info('\n\n### Processing SHRiMP results ###\n')
mir_fi = res_li[1]
mirList, mirStrand = load_mir_info(mir_fi)
tagCount, hits, maps = load_SHRiMP_res(shrimp_dir)
hits, tags, pCount = get_best_alignments(hits, maps)
write_processed_shrimp_output(hits, tags, mirList, mirStrand, shrimp_dir)
remove_temp_file(os.path.basename(shrimp_dir).split('_')[1], out_di['temp'])
logging.info('Total SHRiMP alignments = {}; proportional count = {}'.format(tagCount, pCount))
