diff --git a/workflow/scripts/runDeepTools.py b/workflow/scripts/runDeepTools.py index 164b0b7d9d976d0fd5f485f37ddfc1a323490244..81e02b45127096b5a40e45655d8eab15f42575c0 100644 --- a/workflow/scripts/runDeepTools.py +++ b/workflow/scripts/runDeepTools.py @@ -51,10 +51,10 @@ def run_signal(files, labels, genome): bw_commands = [] for f in files: bw_commands.append("bamCoverage -bs 10 -b "+f+" -o "+f.replace("bam","bw")) - #work_pool = Pool(min(len(files), 12)) - #work_pool.map(bam2bw_wrapper, bw_commands) - #work_pool.close() - #work_pool.join() + work_pool = Pool(min(len(files), 12)) + work_pool.map(bam2bw_wrapper, bw_commands) + work_pool.close() + work_pool.join() cm_command = "computeMatrix scale-regions -R "+gene_bed+" -a 3000 -b 3000 --regionBodyLength 5000 --skipZeros -S *.bw -o samples.deeptools_generegionscalematrix.gz" #p = subprocess.Popen(cm_command, shell=True) diff --git a/workflow/scripts/runMemechip.py b/workflow/scripts/runMemechip.py index d0f66f60abd60aa2bb2ca64af270b7e2bf82b97d..4b17dc0d5b65b001c803005a0302db0ec6aa4099 100644 --- a/workflow/scripts/runMemechip.py +++ b/workflow/scripts/runMemechip.py @@ -39,8 +39,10 @@ def main(): def run(infile, genome, limit, output) infile = pybedtools.BedTool(infile) #genome = twobitreader.TwoBitFile(genome) - output = open(args.outfile,"w") + output = open(output+".fa","w") rowcount = 1 + if limit ==-1: + limit = len(infile) for record in infile: while rowcount <=limit: rowcount += 1 @@ -56,7 +58,7 @@ def run(infile, genome, limit, output) SeqIO.write(newfa,output+".fa","fasta") output.close() #Call memechip - meme_command = "meme-chip -oc "+output+" -meme-minw 5 -meme-maxw 15 -meme-nmotifs 10" + meme_command = "meme-chip -oc "+output+"_memechip"+" -meme-minw 5 -meme-maxw 15 -meme-nmotifs 10 "+output+".fa" p = subprocess.Popen(meme_command, shell=True) p.communicate()