Task | Execution | Time | Dependencies | Task program, Errors, StdOut / StdErr |
Num | 1 | ID | task.callpeak_macs2.macs2_g_b_rep1.line_155.id_10 | Name | macs2 g/b rep1 | Thread | thread_Root | PID | 46995 | | OK | true | Exit Code | 0 | Retries | | State | FINISHED | Dep. | OK | Cpus | | Mem | | | Start | 2017-04-13 18:38:04 | End | 2017-04-13 18:45:43 | Elapsed | 00:07:38 | Timeout | 00:00:-1 | Wall Timeout | 100 days | | Input files | /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322/out/align/rep1/SRR1370894_1.nodup.tagAlign.gz
/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/HughesInputData/input15MReadsNSCLess1.05.50MSubsample.tagAlign.gz
| Output files | /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/rep1/SRR1370894_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.gappedPeak.gz
/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/rep1/SRR1370894_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.broadPeak.gz
| Dependencies | | |
# SYS command. line 157
if [[ -f $(which /software/miniconda3/bin/conda) && $(/software/miniconda3/bin/conda env list | grep aquas_chipseq | wc -l) != "0" ]]; then source /software/miniconda3/bin/activate aquas_chipseq; sleep 5; fi; export PATH=/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/src/chipseq_pipeline/.:/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/src/chipseq_pipeline/modules:/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/src/chipseq_pipeline/utils:${PATH}:/bin:/usr/bin:/usr/local/bin:${HOME}/.bds; set -o pipefail; STARTTIME=$(date +%s); if (( $(nice)<19 )); then renice -n 19 $$; fi
# SYS command. line 158
export LC_COLLATE=C
# SYS command. line 163
macs2 callpeak -t /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322/out/align/rep1/SRR1370894_1.nodup.tagAlign.gz -c /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/HughesInputData/input15MReadsNSCLess1.05.50MSubsample.tagAlign.gz -f BED -n /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/rep1/SRR1370894_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign -g hs -p 0.01 --broad --nomodel --shift 0 --extsize 160 --keep-dup all
# SYS command. line 166
sort -k 8gr,8gr "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/rep1/SRR1370894_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign"_peaks.broadPeak | awk 'BEGIN{OFS="\t"}{$4="Peak_"NR ; print $0}' | gzip -nc > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/rep1/SRR1370894_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.broadPeak.gz
# SYS command. line 167
sort -k 14gr,14gr "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/rep1/SRR1370894_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign"_peaks.gappedPeak | awk 'BEGIN{OFS="\t"}{$4="Peak_"NR ; print $0}' | gzip -nc > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/rep1/SRR1370894_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.gappedPeak.gz
# SYS command. line 170
rm -f "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/rep1/SRR1370894_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign"_peaks.xls \
"/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/rep1/SRR1370894_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign"_peaks.broadPeak \
"/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/rep1/SRR1370894_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign"_peaks.gappedPeak \
"/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/rep1/SRR1370894_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign"_summits.bed
# SYS command. line 175
TASKTIME=$[$(date +%s)-${STARTTIME}]; if [ ${TASKTIME} -lt 60 ]; then echo "Waiting for $[60-${TASKTIME}] seconds."; sleep $[60-${TASKTIME}]; fi
--------------------Stdout--------------------
24546 (process ID) old priority 0, new priority 19
--------------------Stderr--------------------
INFO @ Thu, 13 Apr 2017 18:38:17:
# Command line: callpeak -t /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322/out/align/rep1/SRR1370894_1.nodup.tagAlign.gz -c /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/HughesInputData/input15MReadsNSCLess1.05.50MSubsample.tagAlign.gz -f BED -n /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/rep1/SRR1370894_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign -g hs -p 0.01 --broad --nomodel --shift 0 --extsize 160 --keep-dup all
# ARGUMENTS LIST:
# name = /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/rep1/SRR1370894_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign
# format = BED
# ChIP-seq file = ['/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322/out/align/rep1/SRR1370894_1.nodup.tagAlign.gz']
# control file = ['/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/HughesInputData/input15MReadsNSCLess1.05.50MSubsample.tagAlign.gz']
# effective genome size = 2.70e+09
# band width = 300
# model fold = [5, 50]
# pvalue cutoff for narrow/strong regions = 1.00e-02
# pvalue cutoff for broad/weak regions = 1.00e-01
# qvalue will not be calculated and reported as -1 in the final output.
# Larger dataset will be scaled towards smaller dataset.
# Range for calculating regional lambda is: 1000 bps and 10000 bps
# Broad region calling is on
INFO @ Thu, 13 Apr 2017 18:38:17: #1 read tag files...
INFO @ Thu, 13 Apr 2017 18:38:17: #1 read treatment tags...
INFO @ Thu, 13 Apr 2017 18:38:18: 1000000
INFO @ Thu, 13 Apr 2017 18:38:20: 2000000
INFO @ Thu, 13 Apr 2017 18:38:21: 3000000
INFO @ Thu, 13 Apr 2017 18:38:23: 4000000
INFO @ Thu, 13 Apr 2017 18:38:24: 5000000
INFO @ Thu, 13 Apr 2017 18:38:26: 6000000
INFO @ Thu, 13 Apr 2017 18:38:27: 7000000
INFO @ Thu, 13 Apr 2017 18:38:28: 8000000
INFO @ Thu, 13 Apr 2017 18:38:30: 9000000
INFO @ Thu, 13 Apr 2017 18:38:31: 10000000
INFO @ Thu, 13 Apr 2017 18:38:33: 11000000
INFO @ Thu, 13 Apr 2017 18:38:34: 12000000
INFO @ Thu, 13 Apr 2017 18:38:36: 13000000
INFO @ Thu, 13 Apr 2017 18:38:37: 14000000
INFO @ Thu, 13 Apr 2017 18:38:39: 15000000
INFO @ Thu, 13 Apr 2017 18:38:40: 16000000
INFO @ Thu, 13 Apr 2017 18:38:42: 17000000
INFO @ Thu, 13 Apr 2017 18:38:43: 18000000
INFO @ Thu, 13 Apr 2017 18:38:44: 19000000
INFO @ Thu, 13 Apr 2017 18:38:46: 20000000
INFO @ Thu, 13 Apr 2017 18:38:47: #1.2 read input tags...
INFO @ Thu, 13 Apr 2017 18:38:49: 1000000
INFO @ Thu, 13 Apr 2017 18:38:50: 2000000
INFO @ Thu, 13 Apr 2017 18:38:52: 3000000
INFO @ Thu, 13 Apr 2017 18:38:54: 4000000
INFO @ Thu, 13 Apr 2017 18:38:55: 5000000
INFO @ Thu, 13 Apr 2017 18:38:57: 6000000
INFO @ Thu, 13 Apr 2017 18:38:58: 7000000
INFO @ Thu, 13 Apr 2017 18:39:00: 8000000
INFO @ Thu, 13 Apr 2017 18:39:01: 9000000
INFO @ Thu, 13 Apr 2017 18:39:03: 10000000
INFO @ Thu, 13 Apr 2017 18:39:04: 11000000
INFO @ Thu, 13 Apr 2017 18:39:06: 12000000
INFO @ Thu, 13 Apr 2017 18:39:07: 13000000
INFO @ Thu, 13 Apr 2017 18:39:09: 14000000
INFO @ Thu, 13 Apr 2017 18:39:10: 15000000
INFO @ Thu, 13 Apr 2017 18:39:12: 16000000
INFO @ Thu, 13 Apr 2017 18:39:13: 17000000
INFO @ Thu, 13 Apr 2017 18:39:15: 18000000
INFO @ Thu, 13 Apr 2017 18:39:16: 19000000
INFO @ Thu, 13 Apr 2017 18:39:18: 20000000
INFO @ Thu, 13 Apr 2017 18:39:19: 21000000
INFO @ Thu, 13 Apr 2017 18:39:21: 22000000
INFO @ Thu, 13 Apr 2017 18:39:22: 23000000
INFO @ Thu, 13 Apr 2017 18:39:24: 24000000
INFO @ Thu, 13 Apr 2017 18:39:25: 25000000
INFO @ Thu, 13 Apr 2017 18:39:27: 26000000
INFO @ Thu, 13 Apr 2017 18:39:29: 27000000
INFO @ Thu, 13 Apr 2017 18:39:30: 28000000
INFO @ Thu, 13 Apr 2017 18:39:32: 29000000
INFO @ Thu, 13 Apr 2017 18:39:34: 30000000
INFO @ Thu, 13 Apr 2017 18:39:35: 31000000
INFO @ Thu, 13 Apr 2017 18:39:37: 32000000
INFO @ Thu, 13 Apr 2017 18:39:38: 33000000
INFO @ Thu, 13 Apr 2017 18:39:39: 34000000
INFO @ Thu, 13 Apr 2017 18:39:41: 35000000
INFO @ Thu, 13 Apr 2017 18:39:42: 36000000
INFO @ Thu, 13 Apr 2017 18:39:44: 37000000
INFO @ Thu, 13 Apr 2017 18:39:45: 38000000
INFO @ Thu, 13 Apr 2017 18:39:47: 39000000
INFO @ Thu, 13 Apr 2017 18:39:49: 40000000
INFO @ Thu, 13 Apr 2017 18:39:51: 41000000
INFO @ Thu, 13 Apr 2017 18:39:52: 42000000
INFO @ Thu, 13 Apr 2017 18:39:53: 43000000
INFO @ Thu, 13 Apr 2017 18:39:55: 44000000
INFO @ Thu, 13 Apr 2017 18:39:56: 45000000
INFO @ Thu, 13 Apr 2017 18:39:58: 46000000
INFO @ Thu, 13 Apr 2017 18:40:00: 47000000
INFO @ Thu, 13 Apr 2017 18:40:01: 48000000
INFO @ Thu, 13 Apr 2017 18:40:03: 49000000
INFO @ Thu, 13 Apr 2017 18:40:04: 50000000
INFO @ Thu, 13 Apr 2017 18:40:06: #1 tag size is determined as 49 bps
INFO @ Thu, 13 Apr 2017 18:40:06: #1 tag size = 49
INFO @ Thu, 13 Apr 2017 18:40:06: #1 total tags in treatment: 20891566
INFO @ Thu, 13 Apr 2017 18:40:06: #1 total tags in control: 50000000
INFO @ Thu, 13 Apr 2017 18:40:06: #1 finished!
INFO @ Thu, 13 Apr 2017 18:40:06: #2 Build Peak Model...
INFO @ Thu, 13 Apr 2017 18:40:06: #2 Skipped...
INFO @ Thu, 13 Apr 2017 18:40:06: #2 Use 160 as fragment length
INFO @ Thu, 13 Apr 2017 18:40:06: #3 Call peaks...
INFO @ Thu, 13 Apr 2017 18:40:06: #3 Call broad peaks with given level1 -log10pvalue cutoff and level2: 2.00000, 1.00000...
INFO @ Thu, 13 Apr 2017 18:40:06: #3 Pre-compute pvalue-qvalue table...
INFO @ Thu, 13 Apr 2017 18:42:58: #3 Call peaks for each chromosome...
INFO @ Thu, 13 Apr 2017 18:45:14: #4 Write output xls file... /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/rep1/SRR1370894_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign_peaks.xls
INFO @ Thu, 13 Apr 2017 18:45:21: #4 Write broad peak in broadPeak format file... /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/rep1/SRR1370894_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign_peaks.broadPeak
INFO @ Thu, 13 Apr 2017 18:45:27: #4 Write broad peak in bed12/gappedPeak format file... /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/rep1/SRR1370894_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign_peaks.gappedPeak
INFO @ Thu, 13 Apr 2017 18:45:32: Done!
|
Num | 2 | ID | task.callpeak_macs2.macs2_g_b_rep1_pr1.line_155.id_11 | Name | macs2 g/b rep1-pr1 | Thread | thread_Root | PID | 46996 | | OK | true | Exit Code | 0 | Retries | | State | FINISHED | Dep. | OK | Cpus | | Mem | | | Start | 2017-04-13 18:38:07 | End | 2017-04-13 18:45:29 | Elapsed | 00:07:22 | Timeout | 00:00:-1 | Wall Timeout | 100 days | | Input files | /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/align/pseudo_reps/rep1/pr1/SRR1370894_1.nodup.pr1.tagAlign.gz
/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/HughesInputData/input15MReadsNSCLess1.05.50MSubsample.tagAlign.gz
| Output files | /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/pseudo_reps/rep1/pr1/SRR1370894_1.nodup.pr1.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.gappedPeak.gz
/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/pseudo_reps/rep1/pr1/SRR1370894_1.nodup.pr1.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.broadPeak.gz
| Dependencies | | |
# SYS command. line 157
if [[ -f $(which /software/miniconda3/bin/conda) && $(/software/miniconda3/bin/conda env list | grep aquas_chipseq | wc -l) != "0" ]]; then source /software/miniconda3/bin/activate aquas_chipseq; sleep 5; fi; export PATH=/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/src/chipseq_pipeline/.:/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/src/chipseq_pipeline/modules:/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/src/chipseq_pipeline/utils:${PATH}:/bin:/usr/bin:/usr/local/bin:${HOME}/.bds; set -o pipefail; STARTTIME=$(date +%s); if (( $(nice)<19 )); then renice -n 19 $$; fi
# SYS command. line 158
export LC_COLLATE=C
# SYS command. line 163
macs2 callpeak -t /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/align/pseudo_reps/rep1/pr1/SRR1370894_1.nodup.pr1.tagAlign.gz -c /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/HughesInputData/input15MReadsNSCLess1.05.50MSubsample.tagAlign.gz -f BED -n /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/pseudo_reps/rep1/pr1/SRR1370894_1.nodup.pr1.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign -g hs -p 0.01 --broad --nomodel --shift 0 --extsize 160 --keep-dup all
# SYS command. line 166
sort -k 8gr,8gr "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/pseudo_reps/rep1/pr1/SRR1370894_1.nodup.pr1.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign"_peaks.broadPeak | awk 'BEGIN{OFS="\t"}{$4="Peak_"NR ; print $0}' | gzip -nc > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/pseudo_reps/rep1/pr1/SRR1370894_1.nodup.pr1.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.broadPeak.gz
# SYS command. line 167
sort -k 14gr,14gr "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/pseudo_reps/rep1/pr1/SRR1370894_1.nodup.pr1.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign"_peaks.gappedPeak | awk 'BEGIN{OFS="\t"}{$4="Peak_"NR ; print $0}' | gzip -nc > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/pseudo_reps/rep1/pr1/SRR1370894_1.nodup.pr1.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.gappedPeak.gz
# SYS command. line 170
rm -f "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/pseudo_reps/rep1/pr1/SRR1370894_1.nodup.pr1.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign"_peaks.xls \
"/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/pseudo_reps/rep1/pr1/SRR1370894_1.nodup.pr1.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign"_peaks.broadPeak \
"/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/pseudo_reps/rep1/pr1/SRR1370894_1.nodup.pr1.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign"_peaks.gappedPeak \
"/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/pseudo_reps/rep1/pr1/SRR1370894_1.nodup.pr1.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign"_summits.bed
# SYS command. line 175
TASKTIME=$[$(date +%s)-${STARTTIME}]; if [ ${TASKTIME} -lt 60 ]; then echo "Waiting for $[60-${TASKTIME}] seconds."; sleep $[60-${TASKTIME}]; fi
--------------------Stdout--------------------
25169 (process ID) old priority 0, new priority 19
--------------------Stderr--------------------
INFO @ Thu, 13 Apr 2017 18:38:32:
# Command line: callpeak -t /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/align/pseudo_reps/rep1/pr1/SRR1370894_1.nodup.pr1.tagAlign.gz -c /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/HughesInputData/input15MReadsNSCLess1.05.50MSubsample.tagAlign.gz -f BED -n /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/pseudo_reps/rep1/pr1/SRR1370894_1.nodup.pr1.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign -g hs -p 0.01 --broad --nomodel --shift 0 --extsize 160 --keep-dup all
# ARGUMENTS LIST:
# name = /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/pseudo_reps/rep1/pr1/SRR1370894_1.nodup.pr1.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign
# format = BED
# ChIP-seq file = ['/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/align/pseudo_reps/rep1/pr1/SRR1370894_1.nodup.pr1.tagAlign.gz']
# control file = ['/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/HughesInputData/input15MReadsNSCLess1.05.50MSubsample.tagAlign.gz']
# effective genome size = 2.70e+09
# band width = 300
# model fold = [5, 50]
# pvalue cutoff for narrow/strong regions = 1.00e-02
# pvalue cutoff for broad/weak regions = 1.00e-01
# qvalue will not be calculated and reported as -1 in the final output.
# Larger dataset will be scaled towards smaller dataset.
# Range for calculating regional lambda is: 1000 bps and 10000 bps
# Broad region calling is on
INFO @ Thu, 13 Apr 2017 18:38:32: #1 read tag files...
INFO @ Thu, 13 Apr 2017 18:38:32: #1 read treatment tags...
INFO @ Thu, 13 Apr 2017 18:38:34: 1000000
INFO @ Thu, 13 Apr 2017 18:38:35: 2000000
INFO @ Thu, 13 Apr 2017 18:38:37: 3000000
INFO @ Thu, 13 Apr 2017 18:38:39: 4000000
INFO @ Thu, 13 Apr 2017 18:38:40: 5000000
INFO @ Thu, 13 Apr 2017 18:38:42: 6000000
INFO @ Thu, 13 Apr 2017 18:38:43: 7000000
INFO @ Thu, 13 Apr 2017 18:38:45: 8000000
INFO @ Thu, 13 Apr 2017 18:38:47: 9000000
INFO @ Thu, 13 Apr 2017 18:38:48: 10000000
INFO @ Thu, 13 Apr 2017 18:38:50: #1.2 read input tags...
INFO @ Thu, 13 Apr 2017 18:38:51: 1000000
INFO @ Thu, 13 Apr 2017 18:38:53: 2000000
INFO @ Thu, 13 Apr 2017 18:38:54: 3000000
INFO @ Thu, 13 Apr 2017 18:38:56: 4000000
INFO @ Thu, 13 Apr 2017 18:38:57: 5000000
INFO @ Thu, 13 Apr 2017 18:38:59: 6000000
INFO @ Thu, 13 Apr 2017 18:39:00: 7000000
INFO @ Thu, 13 Apr 2017 18:39:02: 8000000
INFO @ Thu, 13 Apr 2017 18:39:03: 9000000
INFO @ Thu, 13 Apr 2017 18:39:05: 10000000
INFO @ Thu, 13 Apr 2017 18:39:06: 11000000
INFO @ Thu, 13 Apr 2017 18:39:08: 12000000
INFO @ Thu, 13 Apr 2017 18:39:09: 13000000
INFO @ Thu, 13 Apr 2017 18:39:11: 14000000
INFO @ Thu, 13 Apr 2017 18:39:12: 15000000
INFO @ Thu, 13 Apr 2017 18:39:14: 16000000
INFO @ Thu, 13 Apr 2017 18:39:15: 17000000
INFO @ Thu, 13 Apr 2017 18:39:17: 18000000
INFO @ Thu, 13 Apr 2017 18:39:18: 19000000
INFO @ Thu, 13 Apr 2017 18:39:20: 20000000
INFO @ Thu, 13 Apr 2017 18:39:21: 21000000
INFO @ Thu, 13 Apr 2017 18:39:23: 22000000
INFO @ Thu, 13 Apr 2017 18:39:24: 23000000
INFO @ Thu, 13 Apr 2017 18:39:26: 24000000
INFO @ Thu, 13 Apr 2017 18:39:27: 25000000
INFO @ Thu, 13 Apr 2017 18:39:29: 26000000
INFO @ Thu, 13 Apr 2017 18:39:31: 27000000
INFO @ Thu, 13 Apr 2017 18:39:32: 28000000
INFO @ Thu, 13 Apr 2017 18:39:33: 29000000
INFO @ Thu, 13 Apr 2017 18:39:35: 30000000
INFO @ Thu, 13 Apr 2017 18:39:36: 31000000
INFO @ Thu, 13 Apr 2017 18:39:38: 32000000
INFO @ Thu, 13 Apr 2017 18:39:40: 33000000
INFO @ Thu, 13 Apr 2017 18:39:41: 34000000
INFO @ Thu, 13 Apr 2017 18:39:43: 35000000
INFO @ Thu, 13 Apr 2017 18:39:44: 36000000
INFO @ Thu, 13 Apr 2017 18:39:46: 37000000
INFO @ Thu, 13 Apr 2017 18:39:47: 38000000
INFO @ Thu, 13 Apr 2017 18:39:49: 39000000
INFO @ Thu, 13 Apr 2017 18:39:50: 40000000
INFO @ Thu, 13 Apr 2017 18:39:52: 41000000
INFO @ Thu, 13 Apr 2017 18:39:53: 42000000
INFO @ Thu, 13 Apr 2017 18:39:55: 43000000
INFO @ Thu, 13 Apr 2017 18:39:56: 44000000
INFO @ Thu, 13 Apr 2017 18:39:58: 45000000
INFO @ Thu, 13 Apr 2017 18:40:00: 46000000
INFO @ Thu, 13 Apr 2017 18:40:01: 47000000
INFO @ Thu, 13 Apr 2017 18:40:03: 48000000
INFO @ Thu, 13 Apr 2017 18:40:04: 49000000
INFO @ Thu, 13 Apr 2017 18:40:06: 50000000
INFO @ Thu, 13 Apr 2017 18:40:07: #1 tag size is determined as 51 bps
INFO @ Thu, 13 Apr 2017 18:40:07: #1 tag size = 51
INFO @ Thu, 13 Apr 2017 18:40:07: #1 total tags in treatment: 10445783
INFO @ Thu, 13 Apr 2017 18:40:07: #1 total tags in control: 50000000
INFO @ Thu, 13 Apr 2017 18:40:07: #1 finished!
INFO @ Thu, 13 Apr 2017 18:40:07: #2 Build Peak Model...
INFO @ Thu, 13 Apr 2017 18:40:07: #2 Skipped...
INFO @ Thu, 13 Apr 2017 18:40:07: #2 Use 160 as fragment length
INFO @ Thu, 13 Apr 2017 18:40:07: #3 Call peaks...
INFO @ Thu, 13 Apr 2017 18:40:07: #3 Call broad peaks with given level1 -log10pvalue cutoff and level2: 2.00000, 1.00000...
INFO @ Thu, 13 Apr 2017 18:40:07: #3 Pre-compute pvalue-qvalue table...
INFO @ Thu, 13 Apr 2017 18:42:42: #3 Call peaks for each chromosome...
INFO @ Thu, 13 Apr 2017 18:44:50: #4 Write output xls file... /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/pseudo_reps/rep1/pr1/SRR1370894_1.nodup.pr1.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign_peaks.xls
INFO @ Thu, 13 Apr 2017 18:44:56: #4 Write broad peak in broadPeak format file... /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/pseudo_reps/rep1/pr1/SRR1370894_1.nodup.pr1.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign_peaks.broadPeak
INFO @ Thu, 13 Apr 2017 18:45:04: #4 Write broad peak in bed12/gappedPeak format file... /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/pseudo_reps/rep1/pr1/SRR1370894_1.nodup.pr1.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign_peaks.gappedPeak
INFO @ Thu, 13 Apr 2017 18:45:13: Done!
|
Num | 3 | ID | task.callpeak_macs2.macs2_g_b_rep1_pr2.line_155.id_12 | Name | macs2 g/b rep1-pr2 | Thread | thread_Root | PID | 46998 | | OK | true | Exit Code | 0 | Retries | | State | FINISHED | Dep. | OK | Cpus | | Mem | | | Start | 2017-04-13 18:38:08 | End | 2017-04-13 18:45:31 | Elapsed | 00:07:23 | Timeout | 00:00:-1 | Wall Timeout | 100 days | | Input files | /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/align/pseudo_reps/rep1/pr2/SRR1370894_1.nodup.pr2.tagAlign.gz
/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/HughesInputData/input15MReadsNSCLess1.05.50MSubsample.tagAlign.gz
| Output files | /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/pseudo_reps/rep1/pr2/SRR1370894_1.nodup.pr2.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.gappedPeak.gz
/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/pseudo_reps/rep1/pr2/SRR1370894_1.nodup.pr2.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.broadPeak.gz
| Dependencies | | |
# SYS command. line 157
if [[ -f $(which /software/miniconda3/bin/conda) && $(/software/miniconda3/bin/conda env list | grep aquas_chipseq | wc -l) != "0" ]]; then source /software/miniconda3/bin/activate aquas_chipseq; sleep 5; fi; export PATH=/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/src/chipseq_pipeline/.:/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/src/chipseq_pipeline/modules:/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/src/chipseq_pipeline/utils:${PATH}:/bin:/usr/bin:/usr/local/bin:${HOME}/.bds; set -o pipefail; STARTTIME=$(date +%s); if (( $(nice)<19 )); then renice -n 19 $$; fi
# SYS command. line 158
export LC_COLLATE=C
# SYS command. line 163
macs2 callpeak -t /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/align/pseudo_reps/rep1/pr2/SRR1370894_1.nodup.pr2.tagAlign.gz -c /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/HughesInputData/input15MReadsNSCLess1.05.50MSubsample.tagAlign.gz -f BED -n /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/pseudo_reps/rep1/pr2/SRR1370894_1.nodup.pr2.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign -g hs -p 0.01 --broad --nomodel --shift 0 --extsize 160 --keep-dup all
# SYS command. line 166
sort -k 8gr,8gr "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/pseudo_reps/rep1/pr2/SRR1370894_1.nodup.pr2.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign"_peaks.broadPeak | awk 'BEGIN{OFS="\t"}{$4="Peak_"NR ; print $0}' | gzip -nc > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/pseudo_reps/rep1/pr2/SRR1370894_1.nodup.pr2.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.broadPeak.gz
# SYS command. line 167
sort -k 14gr,14gr "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/pseudo_reps/rep1/pr2/SRR1370894_1.nodup.pr2.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign"_peaks.gappedPeak | awk 'BEGIN{OFS="\t"}{$4="Peak_"NR ; print $0}' | gzip -nc > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/pseudo_reps/rep1/pr2/SRR1370894_1.nodup.pr2.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.gappedPeak.gz
# SYS command. line 170
rm -f "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/pseudo_reps/rep1/pr2/SRR1370894_1.nodup.pr2.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign"_peaks.xls \
"/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/pseudo_reps/rep1/pr2/SRR1370894_1.nodup.pr2.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign"_peaks.broadPeak \
"/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/pseudo_reps/rep1/pr2/SRR1370894_1.nodup.pr2.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign"_peaks.gappedPeak \
"/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/pseudo_reps/rep1/pr2/SRR1370894_1.nodup.pr2.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign"_summits.bed
# SYS command. line 175
TASKTIME=$[$(date +%s)-${STARTTIME}]; if [ ${TASKTIME} -lt 60 ]; then echo "Waiting for $[60-${TASKTIME}] seconds."; sleep $[60-${TASKTIME}]; fi
--------------------Stdout--------------------
25170 (process ID) old priority 0, new priority 19
--------------------Stderr--------------------
INFO @ Thu, 13 Apr 2017 18:38:32:
# Command line: callpeak -t /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/align/pseudo_reps/rep1/pr2/SRR1370894_1.nodup.pr2.tagAlign.gz -c /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/HughesInputData/input15MReadsNSCLess1.05.50MSubsample.tagAlign.gz -f BED -n /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/pseudo_reps/rep1/pr2/SRR1370894_1.nodup.pr2.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign -g hs -p 0.01 --broad --nomodel --shift 0 --extsize 160 --keep-dup all
# ARGUMENTS LIST:
# name = /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/pseudo_reps/rep1/pr2/SRR1370894_1.nodup.pr2.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign
# format = BED
# ChIP-seq file = ['/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/align/pseudo_reps/rep1/pr2/SRR1370894_1.nodup.pr2.tagAlign.gz']
# control file = ['/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/HughesInputData/input15MReadsNSCLess1.05.50MSubsample.tagAlign.gz']
# effective genome size = 2.70e+09
# band width = 300
# model fold = [5, 50]
# pvalue cutoff for narrow/strong regions = 1.00e-02
# pvalue cutoff for broad/weak regions = 1.00e-01
# qvalue will not be calculated and reported as -1 in the final output.
# Larger dataset will be scaled towards smaller dataset.
# Range for calculating regional lambda is: 1000 bps and 10000 bps
# Broad region calling is on
INFO @ Thu, 13 Apr 2017 18:38:32: #1 read tag files...
INFO @ Thu, 13 Apr 2017 18:38:32: #1 read treatment tags...
INFO @ Thu, 13 Apr 2017 18:38:34: 1000000
INFO @ Thu, 13 Apr 2017 18:38:36: 2000000
INFO @ Thu, 13 Apr 2017 18:38:37: 3000000
INFO @ Thu, 13 Apr 2017 18:38:39: 4000000
INFO @ Thu, 13 Apr 2017 18:38:40: 5000000
INFO @ Thu, 13 Apr 2017 18:38:42: 6000000
INFO @ Thu, 13 Apr 2017 18:38:44: 7000000
INFO @ Thu, 13 Apr 2017 18:38:45: 8000000
INFO @ Thu, 13 Apr 2017 18:38:47: 9000000
INFO @ Thu, 13 Apr 2017 18:38:48: 10000000
INFO @ Thu, 13 Apr 2017 18:38:50: #1.2 read input tags...
INFO @ Thu, 13 Apr 2017 18:38:52: 1000000
INFO @ Thu, 13 Apr 2017 18:38:53: 2000000
INFO @ Thu, 13 Apr 2017 18:38:55: 3000000
INFO @ Thu, 13 Apr 2017 18:38:57: 4000000
INFO @ Thu, 13 Apr 2017 18:38:58: 5000000
INFO @ Thu, 13 Apr 2017 18:39:00: 6000000
INFO @ Thu, 13 Apr 2017 18:39:02: 7000000
INFO @ Thu, 13 Apr 2017 18:39:03: 8000000
INFO @ Thu, 13 Apr 2017 18:39:05: 9000000
INFO @ Thu, 13 Apr 2017 18:39:06: 10000000
INFO @ Thu, 13 Apr 2017 18:39:08: 11000000
INFO @ Thu, 13 Apr 2017 18:39:10: 12000000
INFO @ Thu, 13 Apr 2017 18:39:11: 13000000
INFO @ Thu, 13 Apr 2017 18:39:13: 14000000
INFO @ Thu, 13 Apr 2017 18:39:15: 15000000
INFO @ Thu, 13 Apr 2017 18:39:16: 16000000
INFO @ Thu, 13 Apr 2017 18:39:18: 17000000
INFO @ Thu, 13 Apr 2017 18:39:19: 18000000
INFO @ Thu, 13 Apr 2017 18:39:21: 19000000
INFO @ Thu, 13 Apr 2017 18:39:22: 20000000
INFO @ Thu, 13 Apr 2017 18:39:24: 21000000
INFO @ Thu, 13 Apr 2017 18:39:26: 22000000
INFO @ Thu, 13 Apr 2017 18:39:27: 23000000
INFO @ Thu, 13 Apr 2017 18:39:29: 24000000
INFO @ Thu, 13 Apr 2017 18:39:31: 25000000
INFO @ Thu, 13 Apr 2017 18:39:32: 26000000
INFO @ Thu, 13 Apr 2017 18:39:34: 27000000
INFO @ Thu, 13 Apr 2017 18:39:36: 28000000
INFO @ Thu, 13 Apr 2017 18:39:37: 29000000
INFO @ Thu, 13 Apr 2017 18:39:39: 30000000
INFO @ Thu, 13 Apr 2017 18:39:40: 31000000
INFO @ Thu, 13 Apr 2017 18:39:42: 32000000
INFO @ Thu, 13 Apr 2017 18:39:43: 33000000
INFO @ Thu, 13 Apr 2017 18:39:45: 34000000
INFO @ Thu, 13 Apr 2017 18:39:46: 35000000
INFO @ Thu, 13 Apr 2017 18:39:48: 36000000
INFO @ Thu, 13 Apr 2017 18:39:49: 37000000
INFO @ Thu, 13 Apr 2017 18:39:51: 38000000
INFO @ Thu, 13 Apr 2017 18:39:52: 39000000
INFO @ Thu, 13 Apr 2017 18:39:54: 40000000
INFO @ Thu, 13 Apr 2017 18:39:55: 41000000
INFO @ Thu, 13 Apr 2017 18:39:57: 42000000
INFO @ Thu, 13 Apr 2017 18:39:59: 43000000
INFO @ Thu, 13 Apr 2017 18:40:00: 44000000
INFO @ Thu, 13 Apr 2017 18:40:02: 45000000
INFO @ Thu, 13 Apr 2017 18:40:03: 46000000
INFO @ Thu, 13 Apr 2017 18:40:05: 47000000
INFO @ Thu, 13 Apr 2017 18:40:06: 48000000
INFO @ Thu, 13 Apr 2017 18:40:08: 49000000
INFO @ Thu, 13 Apr 2017 18:40:09: 50000000
INFO @ Thu, 13 Apr 2017 18:40:11: #1 tag size is determined as 49 bps
INFO @ Thu, 13 Apr 2017 18:40:11: #1 tag size = 49
INFO @ Thu, 13 Apr 2017 18:40:11: #1 total tags in treatment: 10445783
INFO @ Thu, 13 Apr 2017 18:40:11: #1 total tags in control: 50000000
INFO @ Thu, 13 Apr 2017 18:40:11: #1 finished!
INFO @ Thu, 13 Apr 2017 18:40:11: #2 Build Peak Model...
INFO @ Thu, 13 Apr 2017 18:40:11: #2 Skipped...
INFO @ Thu, 13 Apr 2017 18:40:11: #2 Use 160 as fragment length
INFO @ Thu, 13 Apr 2017 18:40:11: #3 Call peaks...
INFO @ Thu, 13 Apr 2017 18:40:11: #3 Call broad peaks with given level1 -log10pvalue cutoff and level2: 2.00000, 1.00000...
INFO @ Thu, 13 Apr 2017 18:40:11: #3 Pre-compute pvalue-qvalue table...
INFO @ Thu, 13 Apr 2017 18:42:46: #3 Call peaks for each chromosome...
INFO @ Thu, 13 Apr 2017 18:44:53: #4 Write output xls file... /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/pseudo_reps/rep1/pr2/SRR1370894_1.nodup.pr2.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign_peaks.xls
INFO @ Thu, 13 Apr 2017 18:45:04: #4 Write broad peak in broadPeak format file... /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/pseudo_reps/rep1/pr2/SRR1370894_1.nodup.pr2.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign_peaks.broadPeak
INFO @ Thu, 13 Apr 2017 18:45:12: #4 Write broad peak in bed12/gappedPeak format file... /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/pseudo_reps/rep1/pr2/SRR1370894_1.nodup.pr2.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign_peaks.gappedPeak
INFO @ Thu, 13 Apr 2017 18:45:17: Done!
|
Num | 4 | ID | task.callpeak_naive_overlap.naive_overlap.line_59.id_13 | Name | naive_overlap | Thread | thread_Root | PID | 47001 | | OK | true | Exit Code | 0 | Retries | | State | FINISHED | Dep. | OK | Cpus | | Mem | | | Start | 2017-04-13 18:45:45 | End | 2017-04-13 18:47:00 | Elapsed | 00:01:14 | Timeout | 01:00:00 | Wall Timeout | 100 days | | Input files | /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/rep1/SRR1370894_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.gappedPeak.gz
/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/pseudo_reps/rep1/pr1/SRR1370894_1.nodup.pr1.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.gappedPeak.gz
/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/pseudo_reps/rep1/pr2/SRR1370894_1.nodup.pr2.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.gappedPeak.gz
| Output files | /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/overlap/SRR1370894_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.naive_overlap.filt.gappedPeak.gz
| Dependencies | | |
# SYS command. line 61
if [[ -f $(which /software/miniconda3/bin/conda) && $(/software/miniconda3/bin/conda env list | grep aquas_chipseq | wc -l) != "0" ]]; then source /software/miniconda3/bin/activate aquas_chipseq; sleep 5; fi; export PATH=/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/src/chipseq_pipeline/.:/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/src/chipseq_pipeline/modules:/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/src/chipseq_pipeline/utils:${PATH}:/bin:/usr/bin:/usr/local/bin:${HOME}/.bds; set -o pipefail; STARTTIME=$(date +%s); if (( $(nice)<19 )); then renice -n 19 $$; fi
# SYS command. line 64
intersectBed -wo -a <(zcat -f /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/rep1/SRR1370894_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.gappedPeak.gz) -b <(zcat -f /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/pseudo_reps/rep1/pr1/SRR1370894_1.nodup.pr1.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.gappedPeak.gz) | \
awk 'BEGIN{FS="\t";OFS="\t"} {s1=$3-$2; s2=$18-$17; if (($31/s1 >= 0.5) || ($31/s2 >= 0.5)) {print $0}}' | cut -f 1-15 | sort | uniq | \
intersectBed -wo -a stdin -b <(zcat -f /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/pseudo_reps/rep1/pr2/SRR1370894_1.nodup.pr2.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.gappedPeak.gz) | \
awk 'BEGIN{FS="\t";OFS="\t"} {s1=$3-$2; s2=$18-$17; if (($31/s1 >= 0.5) || ($31/s2 >= 0.5)) {print $0}}' | cut -f 1-15 | sort | uniq | gzip -nc > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/overlap/SRR1370894_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.PooledInRep1AndRep2.gappedPeak.gz
# SYS command. line 70
zcat /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/overlap/SRR1370894_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.PooledInRep1AndRep2.gappedPeak.gz | sort | uniq | gzip -nc > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/overlap/SRR1370894_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.naive_overlap.gappedPeak.gz
# SYS command. line 72
bedtools intersect -v -a <(zcat -f /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/overlap/SRR1370894_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.naive_overlap.gappedPeak.gz) -b <(zcat -f /mnt/data/annotations/by_release/hg20.GRCh38/hg38.blacklist.bed.gz) | grep -P 'chr[\dXY]+[ \t]' | awk 'BEGIN{OFS="\t"} {if ($5>1000) $5=1000; print $0}' | gzip -nc > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/overlap/SRR1370894_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.naive_overlap.filt.gappedPeak.gz
# SYS command. line 74
rm -f /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/overlap/SRR1370894_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.PooledInRep1AndRep2.gappedPeak.gz
# SYS command. line 76
TASKTIME=$[$(date +%s)-${STARTTIME}]; if [ ${TASKTIME} -lt 60 ]; then echo "Waiting for $[60-${TASKTIME}] seconds."; sleep $[60-${TASKTIME}]; fi
--------------------Stdout--------------------
10003 (process ID) old priority 0, new priority 19
Waiting for 40 seconds.
--------------------Stderr--------------------
***** WARNING: File /dev/fd/62 has inconsistent naming convention for record:
GL000219.1 178041 178976 Peak_3795 195 . 178041 178976 0 2 1,1 0,934 9.89096 21.42591 19.53124
***** WARNING: File /dev/fd/62 has inconsistent naming convention for record:
GL000219.1 178041 178976 Peak_3795 195 . 178041 178976 0 2 1,1 0,934 9.89096 21.42591 19.53124
***** WARNING: File /dev/fd/63 has inconsistent naming convention for record:
GL000219.1 93433 93745 Peak_4391 172 . 93433 93745 0 3 1,182,1 0,66,311 9.40746 19.26130 17.20913
***** WARNING: File /dev/fd/63 has inconsistent naming convention for record:
GL000219.1 93433 93745 Peak_4391 172 . 93433 93745 0 3 1,182,1 0,66,311 9.40746 19.26130 17.20913
***** WARNING: File /dev/fd/63 has inconsistent naming convention for record:
GL000195.1 137499 137916 Peak_206243 3 . 137499 137916 0 3 1,165,1 0,251,416 1.19003 1.08146 0.33098
***** WARNING: File /dev/fd/63 has inconsistent naming convention for record:
GL000195.1 137499 137916 Peak_206243 3 . 137499 137916 0 3 1,165,1 0,251,416 1.19003 1.08146 0.33098
|
Num | 5 | ID | task.callpeak_naive_overlap.naive_overlap.line_59.id_14 | Name | naive_overlap | Thread | thread_Root | PID | 47002 | | OK | true | Exit Code | 0 | Retries | | State | FINISHED | Dep. | OK | Cpus | | Mem | | | Start | 2017-04-13 18:45:45 | End | 2017-04-13 18:47:00 | Elapsed | 00:01:14 | Timeout | 01:00:00 | Wall Timeout | 100 days | | Input files | /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/rep1/SRR1370894_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.broadPeak.gz
/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/pseudo_reps/rep1/pr1/SRR1370894_1.nodup.pr1.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.broadPeak.gz
/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/pseudo_reps/rep1/pr2/SRR1370894_1.nodup.pr2.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.broadPeak.gz
| Output files | /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/overlap/SRR1370894_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.naive_overlap.filt.broadPeak.gz
| Dependencies | | |
# SYS command. line 61
if [[ -f $(which /software/miniconda3/bin/conda) && $(/software/miniconda3/bin/conda env list | grep aquas_chipseq | wc -l) != "0" ]]; then source /software/miniconda3/bin/activate aquas_chipseq; sleep 5; fi; export PATH=/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/src/chipseq_pipeline/.:/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/src/chipseq_pipeline/modules:/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/src/chipseq_pipeline/utils:${PATH}:/bin:/usr/bin:/usr/local/bin:${HOME}/.bds; set -o pipefail; STARTTIME=$(date +%s); if (( $(nice)<19 )); then renice -n 19 $$; fi
# SYS command. line 64
intersectBed -wo -a <(zcat -f /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/rep1/SRR1370894_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.broadPeak.gz) -b <(zcat -f /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/pseudo_reps/rep1/pr1/SRR1370894_1.nodup.pr1.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.broadPeak.gz) | \
awk 'BEGIN{FS="\t";OFS="\t"} {s1=$3-$2; s2=$12-$11; if (($19/s1 >= 0.5) || ($19/s2 >= 0.5)) {print $0}}' | cut -f 1-9 | sort | uniq | \
intersectBed -wo -a stdin -b <(zcat -f /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/pseudo_reps/rep1/pr2/SRR1370894_1.nodup.pr2.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.broadPeak.gz) | \
awk 'BEGIN{FS="\t";OFS="\t"} {s1=$3-$2; s2=$12-$11; if (($19/s1 >= 0.5) || ($19/s2 >= 0.5)) {print $0}}' | cut -f 1-9 | sort | uniq | gzip -nc > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/overlap/SRR1370894_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.PooledInRep1AndRep2.broadPeak.gz
# SYS command. line 70
zcat /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/overlap/SRR1370894_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.PooledInRep1AndRep2.broadPeak.gz | sort | uniq | gzip -nc > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/overlap/SRR1370894_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.naive_overlap.broadPeak.gz
# SYS command. line 72
bedtools intersect -v -a <(zcat -f /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/overlap/SRR1370894_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.naive_overlap.broadPeak.gz) -b <(zcat -f /mnt/data/annotations/by_release/hg20.GRCh38/hg38.blacklist.bed.gz) | grep -P 'chr[\dXY]+[ \t]' | awk 'BEGIN{OFS="\t"} {if ($5>1000) $5=1000; print $0}' | gzip -nc > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/overlap/SRR1370894_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.naive_overlap.filt.broadPeak.gz
# SYS command. line 74
rm -f /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/overlap/SRR1370894_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.PooledInRep1AndRep2.broadPeak.gz
# SYS command. line 76
TASKTIME=$[$(date +%s)-${STARTTIME}]; if [ ${TASKTIME} -lt 60 ]; then echo "Waiting for $[60-${TASKTIME}] seconds."; sleep $[60-${TASKTIME}]; fi
--------------------Stdout--------------------
10004 (process ID) old priority 0, new priority 19
Waiting for 44 seconds.
--------------------Stderr--------------------
***** WARNING: File /dev/fd/62 has inconsistent naming convention for record:
GL000219.1 178041 178976 Peak_3795 195 . 9.89096 21.42591 19.53124
***** WARNING: File /dev/fd/62 has inconsistent naming convention for record:
GL000219.1 178041 178976 Peak_3795 195 . 9.89096 21.42591 19.53124
***** WARNING: File /dev/fd/63 has inconsistent naming convention for record:
GL000219.1 93433 93745 Peak_4391 172 . 9.40746 19.26130 17.20913
***** WARNING: File /dev/fd/63 has inconsistent naming convention for record:
GL000219.1 93433 93745 Peak_4391 172 . 9.40746 19.26130 17.20913
***** WARNING: File /dev/fd/63 has inconsistent naming convention for record:
GL000195.1 137499 137916 Peak_206243 3 . 1.19003 1.08146 0.33098
***** WARNING: File /dev/fd/63 has inconsistent naming convention for record:
GL000195.1 137499 137916 Peak_206243 3 . 1.19003 1.08146 0.33098
|
Num | 6 | ID | task.report.peak2hammock.line_361.id_17 | Name | peak2hammock | Thread | thread_Root | PID | 12826 | | OK | true | Exit Code | 0 | Retries | | State | FINISHED | Dep. | OK | Cpus | | Mem | | | Start | 2017-04-13 18:47:02 | End | 2017-04-13 18:47:20 | Elapsed | 00:00:17 | Timeout | 00:00:-1 | Wall Timeout | 100 days | | Input files | /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/overlap/SRR1370894_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.naive_overlap.filt.gappedPeak.gz
| Output files | /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/overlap/SRR1370894_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.naive_overlap.filt.gappedPeak.hammock.gz
| Dependencies | | |
# SYS command. line 363
if [[ -f $(which /software/miniconda3/bin/conda) && $(/software/miniconda3/bin/conda env list | grep aquas_chipseq | wc -l) != "0" ]]; then source /software/miniconda3/bin/activate aquas_chipseq; sleep 5; fi; export PATH=/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/src/chipseq_pipeline/.:/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/src/chipseq_pipeline/modules:/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/src/chipseq_pipeline/utils:${PATH}:/bin:/usr/bin:/usr/local/bin:${HOME}/.bds; set -o pipefail; STARTTIME=$(date +%s); if (( $(nice)<19 )); then renice -n 19 $$; fi
# SYS command. line 364
zcat /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/overlap/SRR1370894_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.naive_overlap.filt.gappedPeak.gz | sed '/^\(chr\)/!d' | sort -k1,1V -k2,2n > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/overlap/SRR1370894_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.naive_overlap.filt.gappedPeak.tmp
# SYS command. line 366
/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/src/chipseq_pipeline/utils/gappedpeak.py /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/overlap/SRR1370894_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.naive_overlap.filt.gappedPeak.tmp /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/overlap/SRR1370894_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.naive_overlap.filt.gappedPeak.hammock
# SYS command. line 367
rm -f /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/peak/macs2/overlap/SRR1370894_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.naive_overlap.filt.gappedPeak.tmp
--------------------Stdout--------------------
12830 (process ID) old priority 0, new priority 19
|
Num | 7 | ID | task.graphviz.report.line_97.id_18 | Name | report | Thread | thread_Root | PID | 13654 | | OK | true | Exit Code | 0 | Retries | | State | FINISHED | Dep. | OK | Cpus | | Mem | | | Start | 2017-04-13 18:47:21 | End | 2017-04-13 18:47:29 | Elapsed | 00:00:08 | Timeout | 00:00:-1 | Wall Timeout | 100 days | | Input files | | Output files | | Dependencies | | |
# SYS command. line 98
if [[ -f $(which /software/miniconda3/bin/conda) && $(/software/miniconda3/bin/conda env list | grep aquas_chipseq | wc -l) != "0" ]]; then source /software/miniconda3/bin/activate aquas_chipseq; sleep 5; fi; export PATH=/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/src/chipseq_pipeline/.:/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/src/chipseq_pipeline/modules:/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/src/chipseq_pipeline/utils:${PATH}:/bin:/usr/bin:/usr/local/bin:${HOME}/.bds; set -o pipefail; STARTTIME=$(date +%s); if (( $(nice)<19 )); then renice -n 19 $$; fi
# SYS command. line 99
dot -Tsvg /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/report/ZNF322_MACS2_workflow.dot > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/report/ZNF322_MACS2_workflow.svg 2> /dev/null || echo "svg: /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF322_MACS2/out/report/ZNF322_MACS2_workflow.svg" # to suppress dot font error (exit code=1)
--------------------Stdout--------------------
13658 (process ID) old priority 0, new priority 19
|