Task | Execution | Time | Dependencies | Task program, Errors, StdOut / StdErr |
Num | 1 | ID | task.callpeak_macs2.macs2_g_b_pooled_rep.line_155.id_10 | Name | macs2 g/b pooled_rep | Thread | thread_Root | PID | 46965 | | OK | true | Exit Code | 0 | Retries | | State | FINISHED | Dep. | OK | Cpus | | Mem | | | Start | 2017-04-13 18:37:48 | End | 2017-04-13 19:05:00 | Elapsed | 00:27:12 | Timeout | 00:00:-1 | Wall Timeout | 100 days | | Input files | /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/align/pooled_rep/ENCFF558PHY_ENCFF645WYX.tagAlign.gz
/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/InputData/ss_50M_2000_GRCh38.nodup.tagAlign.gz
| Output files | /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pooled_rep/ENCFF558PHY_ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.gappedPeak.gz
/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pooled_rep/ENCFF558PHY_ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.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/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/align/pooled_rep/ENCFF558PHY_ENCFF645WYX.tagAlign.gz -c /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/InputData/ss_50M_2000_GRCh38.nodup.tagAlign.gz -f BED -n /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pooled_rep/ENCFF558PHY_ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign -g hs -p 0.01 --broad --nomodel --shift 0 --extsize 213 --keep-dup all
# SYS command. line 166
sort -k 8gr,8gr "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pooled_rep/ENCFF558PHY_ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign"_peaks.broadPeak | awk 'BEGIN{OFS="\t"}{$4="Peak_"NR ; print $0}' | gzip -nc > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pooled_rep/ENCFF558PHY_ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.broadPeak.gz
# SYS command. line 167
sort -k 14gr,14gr "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pooled_rep/ENCFF558PHY_ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign"_peaks.gappedPeak | awk 'BEGIN{OFS="\t"}{$4="Peak_"NR ; print $0}' | gzip -nc > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pooled_rep/ENCFF558PHY_ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.gappedPeak.gz
# SYS command. line 170
rm -f "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pooled_rep/ENCFF558PHY_ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign"_peaks.xls \
"/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pooled_rep/ENCFF558PHY_ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign"_peaks.broadPeak \
"/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pooled_rep/ENCFF558PHY_ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign"_peaks.gappedPeak \
"/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pooled_rep/ENCFF558PHY_ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.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--------------------
14014 (process ID) old priority 0, new priority 19
--------------------Stderr--------------------
INFO @ Thu, 13 Apr 2017 18:38:02:
# Command line: callpeak -t /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/align/pooled_rep/ENCFF558PHY_ENCFF645WYX.tagAlign.gz -c /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/InputData/ss_50M_2000_GRCh38.nodup.tagAlign.gz -f BED -n /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pooled_rep/ENCFF558PHY_ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign -g hs -p 0.01 --broad --nomodel --shift 0 --extsize 213 --keep-dup all
# ARGUMENTS LIST:
# name = /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pooled_rep/ENCFF558PHY_ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign
# format = BED
# ChIP-seq file = ['/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/align/pooled_rep/ENCFF558PHY_ENCFF645WYX.tagAlign.gz']
# control file = ['/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/InputData/ss_50M_2000_GRCh38.nodup.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:02: #1 read tag files...
INFO @ Thu, 13 Apr 2017 18:38:02: #1 read treatment tags...
INFO @ Thu, 13 Apr 2017 18:38:04: 1000000
INFO @ Thu, 13 Apr 2017 18:38:06: 2000000
INFO @ Thu, 13 Apr 2017 18:38:08: 3000000
INFO @ Thu, 13 Apr 2017 18:38:10: 4000000
INFO @ Thu, 13 Apr 2017 18:38:12: 5000000
INFO @ Thu, 13 Apr 2017 18:38:14: 6000000
INFO @ Thu, 13 Apr 2017 18:38:16: 7000000
INFO @ Thu, 13 Apr 2017 18:38:18: 8000000
INFO @ Thu, 13 Apr 2017 18:38:20: 9000000
INFO @ Thu, 13 Apr 2017 18:38:22: 10000000
INFO @ Thu, 13 Apr 2017 18:38:24: 11000000
INFO @ Thu, 13 Apr 2017 18:38:26: 12000000
INFO @ Thu, 13 Apr 2017 18:38:28: 13000000
INFO @ Thu, 13 Apr 2017 18:38:30: 14000000
INFO @ Thu, 13 Apr 2017 18:38:32: 15000000
INFO @ Thu, 13 Apr 2017 18:38:34: 16000000
INFO @ Thu, 13 Apr 2017 18:38:36: 17000000
INFO @ Thu, 13 Apr 2017 18:38:38: 18000000
INFO @ Thu, 13 Apr 2017 18:38:40: 19000000
INFO @ Thu, 13 Apr 2017 18:38:42: 20000000
INFO @ Thu, 13 Apr 2017 18:38:44: 21000000
INFO @ Thu, 13 Apr 2017 18:38:46: 22000000
INFO @ Thu, 13 Apr 2017 18:38:48: 23000000
INFO @ Thu, 13 Apr 2017 18:38:50: 24000000
INFO @ Thu, 13 Apr 2017 18:38:52: 25000000
INFO @ Thu, 13 Apr 2017 18:38:54: 26000000
INFO @ Thu, 13 Apr 2017 18:38:56: 27000000
INFO @ Thu, 13 Apr 2017 18:38:58: 28000000
INFO @ Thu, 13 Apr 2017 18:39:00: 29000000
INFO @ Thu, 13 Apr 2017 18:39:02: 30000000
INFO @ Thu, 13 Apr 2017 18:39:04: 31000000
INFO @ Thu, 13 Apr 2017 18:39:06: 32000000
INFO @ Thu, 13 Apr 2017 18:39:08: 33000000
INFO @ Thu, 13 Apr 2017 18:39:10: 34000000
INFO @ Thu, 13 Apr 2017 18:39:12: 35000000
INFO @ Thu, 13 Apr 2017 18:39:14: 36000000
INFO @ Thu, 13 Apr 2017 18:39:16: 37000000
INFO @ Thu, 13 Apr 2017 18:39:18: 38000000
INFO @ Thu, 13 Apr 2017 18:39:20: 39000000
INFO @ Thu, 13 Apr 2017 18:39:22: 40000000
INFO @ Thu, 13 Apr 2017 18:39:23: 41000000
INFO @ Thu, 13 Apr 2017 18:39:25: 42000000
INFO @ Thu, 13 Apr 2017 18:39:27: 43000000
INFO @ Thu, 13 Apr 2017 18:39:29: 44000000
INFO @ Thu, 13 Apr 2017 18:39:31: 45000000
INFO @ Thu, 13 Apr 2017 18:39:33: 46000000
INFO @ Thu, 13 Apr 2017 18:39:35: 47000000
INFO @ Thu, 13 Apr 2017 18:39:37: 48000000
INFO @ Thu, 13 Apr 2017 18:39:39: 49000000
INFO @ Thu, 13 Apr 2017 18:39:41: 50000000
INFO @ Thu, 13 Apr 2017 18:39:43: 51000000
INFO @ Thu, 13 Apr 2017 18:39:45: 52000000
INFO @ Thu, 13 Apr 2017 18:39:47: 53000000
INFO @ Thu, 13 Apr 2017 18:39:49: 54000000
INFO @ Thu, 13 Apr 2017 18:39:51: 55000000
INFO @ Thu, 13 Apr 2017 18:39:53: 56000000
INFO @ Thu, 13 Apr 2017 18:39:55: 57000000
INFO @ Thu, 13 Apr 2017 18:39:57: 58000000
INFO @ Thu, 13 Apr 2017 18:39:59: 59000000
INFO @ Thu, 13 Apr 2017 18:40:01: 60000000
INFO @ Thu, 13 Apr 2017 18:40:03: 61000000
INFO @ Thu, 13 Apr 2017 18:40:05: 62000000
INFO @ Thu, 13 Apr 2017 18:40:07: 63000000
INFO @ Thu, 13 Apr 2017 18:40:08: 64000000
INFO @ Thu, 13 Apr 2017 18:40:10: 65000000
INFO @ Thu, 13 Apr 2017 18:40:12: 66000000
INFO @ Thu, 13 Apr 2017 18:40:14: 67000000
INFO @ Thu, 13 Apr 2017 18:40:16: 68000000
INFO @ Thu, 13 Apr 2017 18:40:18: 69000000
INFO @ Thu, 13 Apr 2017 18:40:20: 70000000
INFO @ Thu, 13 Apr 2017 18:40:22: 71000000
INFO @ Thu, 13 Apr 2017 18:40:24: 72000000
INFO @ Thu, 13 Apr 2017 18:40:26: 73000000
INFO @ Thu, 13 Apr 2017 18:40:28: 74000000
INFO @ Thu, 13 Apr 2017 18:40:30: 75000000
INFO @ Thu, 13 Apr 2017 18:40:32: 76000000
INFO @ Thu, 13 Apr 2017 18:40:34: 77000000
INFO @ Thu, 13 Apr 2017 18:40:36: 78000000
INFO @ Thu, 13 Apr 2017 18:40:38: 79000000
INFO @ Thu, 13 Apr 2017 18:40:40: 80000000
INFO @ Thu, 13 Apr 2017 18:40:42: 81000000
INFO @ Thu, 13 Apr 2017 18:40:44: 82000000
INFO @ Thu, 13 Apr 2017 18:40:46: 83000000
INFO @ Thu, 13 Apr 2017 18:40:47: 84000000
INFO @ Thu, 13 Apr 2017 18:40:49: 85000000
INFO @ Thu, 13 Apr 2017 18:40:51: 86000000
INFO @ Thu, 13 Apr 2017 18:40:53: 87000000
INFO @ Thu, 13 Apr 2017 18:40:55: 88000000
INFO @ Thu, 13 Apr 2017 18:40:57: 89000000
INFO @ Thu, 13 Apr 2017 18:40:59: 90000000
INFO @ Thu, 13 Apr 2017 18:41:01: 91000000
INFO @ Thu, 13 Apr 2017 18:41:03: 92000000
INFO @ Thu, 13 Apr 2017 18:41:05: 93000000
INFO @ Thu, 13 Apr 2017 18:41:07: 94000000
INFO @ Thu, 13 Apr 2017 18:41:09: 95000000
INFO @ Thu, 13 Apr 2017 18:41:11: 96000000
INFO @ Thu, 13 Apr 2017 18:41:13: 97000000
INFO @ Thu, 13 Apr 2017 18:41:15: 98000000
INFO @ Thu, 13 Apr 2017 18:41:17: 99000000
INFO @ Thu, 13 Apr 2017 18:41:19: 100000000
INFO @ Thu, 13 Apr 2017 18:41:21: 101000000
INFO @ Thu, 13 Apr 2017 18:41:23: 102000000
INFO @ Thu, 13 Apr 2017 18:41:25: 103000000
INFO @ Thu, 13 Apr 2017 18:41:27: 104000000
INFO @ Thu, 13 Apr 2017 18:41:29: 105000000
INFO @ Thu, 13 Apr 2017 18:41:30: 106000000
INFO @ Thu, 13 Apr 2017 18:41:32: 107000000
INFO @ Thu, 13 Apr 2017 18:41:34: 108000000
INFO @ Thu, 13 Apr 2017 18:41:45: #1.2 read input tags...
INFO @ Thu, 13 Apr 2017 18:41:47: 1000000
INFO @ Thu, 13 Apr 2017 18:41:48: 2000000
INFO @ Thu, 13 Apr 2017 18:41:50: 3000000
INFO @ Thu, 13 Apr 2017 18:41:52: 4000000
INFO @ Thu, 13 Apr 2017 18:41:54: 5000000
INFO @ Thu, 13 Apr 2017 18:41:56: 6000000
INFO @ Thu, 13 Apr 2017 18:41:58: 7000000
INFO @ Thu, 13 Apr 2017 18:41:59: 8000000
INFO @ Thu, 13 Apr 2017 18:42:01: 9000000
INFO @ Thu, 13 Apr 2017 18:42:03: 10000000
INFO @ Thu, 13 Apr 2017 18:42:05: 11000000
INFO @ Thu, 13 Apr 2017 18:42:07: 12000000
INFO @ Thu, 13 Apr 2017 18:42:09: 13000000
INFO @ Thu, 13 Apr 2017 18:42:10: 14000000
INFO @ Thu, 13 Apr 2017 18:42:12: 15000000
INFO @ Thu, 13 Apr 2017 18:42:14: 16000000
INFO @ Thu, 13 Apr 2017 18:42:16: 17000000
INFO @ Thu, 13 Apr 2017 18:42:18: 18000000
INFO @ Thu, 13 Apr 2017 18:42:20: 19000000
INFO @ Thu, 13 Apr 2017 18:42:21: 20000000
INFO @ Thu, 13 Apr 2017 18:42:23: 21000000
INFO @ Thu, 13 Apr 2017 18:42:25: 22000000
INFO @ Thu, 13 Apr 2017 18:42:27: 23000000
INFO @ Thu, 13 Apr 2017 18:42:29: 24000000
INFO @ Thu, 13 Apr 2017 18:42:31: 25000000
INFO @ Thu, 13 Apr 2017 18:42:33: 26000000
INFO @ Thu, 13 Apr 2017 18:42:35: 27000000
INFO @ Thu, 13 Apr 2017 18:42:36: 28000000
INFO @ Thu, 13 Apr 2017 18:42:38: 29000000
INFO @ Thu, 13 Apr 2017 18:42:40: 30000000
INFO @ Thu, 13 Apr 2017 18:42:42: 31000000
INFO @ Thu, 13 Apr 2017 18:42:44: 32000000
INFO @ Thu, 13 Apr 2017 18:42:46: 33000000
INFO @ Thu, 13 Apr 2017 18:42:47: 34000000
INFO @ Thu, 13 Apr 2017 18:42:49: 35000000
INFO @ Thu, 13 Apr 2017 18:42:51: 36000000
INFO @ Thu, 13 Apr 2017 18:42:53: 37000000
INFO @ Thu, 13 Apr 2017 18:42:55: 38000000
INFO @ Thu, 13 Apr 2017 18:42:57: 39000000
INFO @ Thu, 13 Apr 2017 18:42:58: 40000000
INFO @ Thu, 13 Apr 2017 18:43:00: 41000000
INFO @ Thu, 13 Apr 2017 18:43:02: 42000000
INFO @ Thu, 13 Apr 2017 18:43:04: 43000000
INFO @ Thu, 13 Apr 2017 18:43:06: 44000000
INFO @ Thu, 13 Apr 2017 18:43:08: 45000000
INFO @ Thu, 13 Apr 2017 18:43:09: 46000000
INFO @ Thu, 13 Apr 2017 18:43:11: 47000000
INFO @ Thu, 13 Apr 2017 18:43:13: 48000000
INFO @ Thu, 13 Apr 2017 18:43:15: 49000000
INFO @ Thu, 13 Apr 2017 18:43:17: 50000000
INFO @ Thu, 13 Apr 2017 18:43:19: 51000000
INFO @ Thu, 13 Apr 2017 18:43:20: 52000000
INFO @ Thu, 13 Apr 2017 18:43:22: 53000000
INFO @ Thu, 13 Apr 2017 18:43:24: 54000000
INFO @ Thu, 13 Apr 2017 18:43:26: 55000000
INFO @ Thu, 13 Apr 2017 18:43:28: 56000000
INFO @ Thu, 13 Apr 2017 18:43:30: 57000000
INFO @ Thu, 13 Apr 2017 18:43:31: 58000000
INFO @ Thu, 13 Apr 2017 18:43:33: 59000000
INFO @ Thu, 13 Apr 2017 18:43:35: 60000000
INFO @ Thu, 13 Apr 2017 18:43:37: 61000000
INFO @ Thu, 13 Apr 2017 18:43:39: 62000000
INFO @ Thu, 13 Apr 2017 18:43:41: 63000000
INFO @ Thu, 13 Apr 2017 18:43:43: 64000000
INFO @ Thu, 13 Apr 2017 18:43:45: 65000000
INFO @ Thu, 13 Apr 2017 18:43:47: 66000000
INFO @ Thu, 13 Apr 2017 18:43:48: 67000000
INFO @ Thu, 13 Apr 2017 18:43:50: 68000000
INFO @ Thu, 13 Apr 2017 18:43:52: 69000000
INFO @ Thu, 13 Apr 2017 18:43:54: 70000000
INFO @ Thu, 13 Apr 2017 18:43:56: 71000000
INFO @ Thu, 13 Apr 2017 18:43:58: 72000000
INFO @ Thu, 13 Apr 2017 18:43:59: 73000000
INFO @ Thu, 13 Apr 2017 18:44:01: 74000000
INFO @ Thu, 13 Apr 2017 18:44:04: 75000000
INFO @ Thu, 13 Apr 2017 18:44:05: 76000000
INFO @ Thu, 13 Apr 2017 18:44:07: 77000000
INFO @ Thu, 13 Apr 2017 18:44:09: 78000000
INFO @ Thu, 13 Apr 2017 18:44:11: 79000000
INFO @ Thu, 13 Apr 2017 18:44:13: 80000000
INFO @ Thu, 13 Apr 2017 18:44:15: 81000000
INFO @ Thu, 13 Apr 2017 18:44:16: 82000000
INFO @ Thu, 13 Apr 2017 18:44:18: 83000000
INFO @ Thu, 13 Apr 2017 18:44:20: 84000000
INFO @ Thu, 13 Apr 2017 18:44:22: 85000000
INFO @ Thu, 13 Apr 2017 18:44:24: 86000000
INFO @ Thu, 13 Apr 2017 18:44:26: 87000000
INFO @ Thu, 13 Apr 2017 18:44:27: 88000000
INFO @ Thu, 13 Apr 2017 18:44:29: 89000000
INFO @ Thu, 13 Apr 2017 18:44:31: 90000000
INFO @ Thu, 13 Apr 2017 18:44:33: 91000000
INFO @ Thu, 13 Apr 2017 18:44:35: 92000000
INFO @ Thu, 13 Apr 2017 18:44:37: 93000000
INFO @ Thu, 13 Apr 2017 18:44:38: 94000000
INFO @ Thu, 13 Apr 2017 18:44:40: 95000000
INFO @ Thu, 13 Apr 2017 18:44:42: 96000000
INFO @ Thu, 13 Apr 2017 18:44:44: 97000000
INFO @ Thu, 13 Apr 2017 18:44:46: 98000000
INFO @ Thu, 13 Apr 2017 18:44:47: 99000000
INFO @ Thu, 13 Apr 2017 18:44:49: #1 tag size is determined as 101 bps
INFO @ Thu, 13 Apr 2017 18:44:49: #1 tag size = 101
INFO @ Thu, 13 Apr 2017 18:44:49: #1 total tags in treatment: 108422852
INFO @ Thu, 13 Apr 2017 18:44:49: #1 total tags in control: 99493262
INFO @ Thu, 13 Apr 2017 18:44:49: #1 finished!
INFO @ Thu, 13 Apr 2017 18:44:49: #2 Build Peak Model...
INFO @ Thu, 13 Apr 2017 18:44:49: #2 Skipped...
INFO @ Thu, 13 Apr 2017 18:44:49: #2 Use 213 as fragment length
INFO @ Thu, 13 Apr 2017 18:44:49: #3 Call peaks...
INFO @ Thu, 13 Apr 2017 18:44:49: #3 Call broad peaks with given level1 -log10pvalue cutoff and level2: 2.00000, 1.00000...
INFO @ Thu, 13 Apr 2017 18:44:49: #3 Pre-compute pvalue-qvalue table...
INFO @ Thu, 13 Apr 2017 18:56:21: #3 Call peaks for each chromosome...
INFO @ Thu, 13 Apr 2017 19:04:52: #4 Write output xls file... /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pooled_rep/ENCFF558PHY_ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign_peaks.xls
INFO @ Thu, 13 Apr 2017 19:04:54: #4 Write broad peak in broadPeak format file... /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pooled_rep/ENCFF558PHY_ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign_peaks.broadPeak
INFO @ Thu, 13 Apr 2017 19:04:55: #4 Write broad peak in bed12/gappedPeak format file... /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pooled_rep/ENCFF558PHY_ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign_peaks.gappedPeak
INFO @ Thu, 13 Apr 2017 19:04:55: Done!
|
Num | 2 | ID | task.callpeak_macs2.macs2_g_b_ppr1.line_155.id_11 | Name | macs2 g/b ppr1 | Thread | thread_Root | PID | 46968 | | OK | true | Exit Code | 0 | Retries | | State | FINISHED | Dep. | OK | Cpus | | Mem | | | Start | 2017-04-13 18:37:50 | End | 2017-04-13 18:51:29 | Elapsed | 00:13:38 | Timeout | 00:00:-1 | Wall Timeout | 100 days | | Input files | /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/align/pooled_pseudo_reps/ppr1/ENCFF558PHY.pr1_ENCFF645WYX.pr1.tagAlign.gz
/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/InputData/ss_50M_2000_GRCh38.nodup.tagAlign.gz
| Output files | /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pooled_pseudo_reps/ppr1/ENCFF558PHY.pr1_ENCFF645WYX.pr1.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.gappedPeak.gz
/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pooled_pseudo_reps/ppr1/ENCFF558PHY.pr1_ENCFF645WYX.pr1.tagAlign_x_ss_50M_2000_GRCh38.nodup.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/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/align/pooled_pseudo_reps/ppr1/ENCFF558PHY.pr1_ENCFF645WYX.pr1.tagAlign.gz -c /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/InputData/ss_50M_2000_GRCh38.nodup.tagAlign.gz -f BED -n /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pooled_pseudo_reps/ppr1/ENCFF558PHY.pr1_ENCFF645WYX.pr1.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign -g hs -p 0.01 --broad --nomodel --shift 0 --extsize 213 --keep-dup all
# SYS command. line 166
sort -k 8gr,8gr "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pooled_pseudo_reps/ppr1/ENCFF558PHY.pr1_ENCFF645WYX.pr1.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign"_peaks.broadPeak | awk 'BEGIN{OFS="\t"}{$4="Peak_"NR ; print $0}' | gzip -nc > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pooled_pseudo_reps/ppr1/ENCFF558PHY.pr1_ENCFF645WYX.pr1.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.broadPeak.gz
# SYS command. line 167
sort -k 14gr,14gr "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pooled_pseudo_reps/ppr1/ENCFF558PHY.pr1_ENCFF645WYX.pr1.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign"_peaks.gappedPeak | awk 'BEGIN{OFS="\t"}{$4="Peak_"NR ; print $0}' | gzip -nc > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pooled_pseudo_reps/ppr1/ENCFF558PHY.pr1_ENCFF645WYX.pr1.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.gappedPeak.gz
# SYS command. line 170
rm -f "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pooled_pseudo_reps/ppr1/ENCFF558PHY.pr1_ENCFF645WYX.pr1.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign"_peaks.xls \
"/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pooled_pseudo_reps/ppr1/ENCFF558PHY.pr1_ENCFF645WYX.pr1.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign"_peaks.broadPeak \
"/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pooled_pseudo_reps/ppr1/ENCFF558PHY.pr1_ENCFF645WYX.pr1.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign"_peaks.gappedPeak \
"/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pooled_pseudo_reps/ppr1/ENCFF558PHY.pr1_ENCFF645WYX.pr1.tagAlign_x_ss_50M_2000_GRCh38.nodup.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--------------------
47349 (process ID) old priority 0, new priority 19
--------------------Stderr--------------------
INFO @ Thu, 13 Apr 2017 18:38:01:
# Command line: callpeak -t /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/align/pooled_pseudo_reps/ppr1/ENCFF558PHY.pr1_ENCFF645WYX.pr1.tagAlign.gz -c /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/InputData/ss_50M_2000_GRCh38.nodup.tagAlign.gz -f BED -n /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pooled_pseudo_reps/ppr1/ENCFF558PHY.pr1_ENCFF645WYX.pr1.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign -g hs -p 0.01 --broad --nomodel --shift 0 --extsize 213 --keep-dup all
# ARGUMENTS LIST:
# name = /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pooled_pseudo_reps/ppr1/ENCFF558PHY.pr1_ENCFF645WYX.pr1.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign
# format = BED
# ChIP-seq file = ['/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/align/pooled_pseudo_reps/ppr1/ENCFF558PHY.pr1_ENCFF645WYX.pr1.tagAlign.gz']
# control file = ['/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/InputData/ss_50M_2000_GRCh38.nodup.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:01: #1 read tag files...
INFO @ Thu, 13 Apr 2017 18:38:01: #1 read treatment tags...
INFO @ Thu, 13 Apr 2017 18:38:03: 1000000
INFO @ Thu, 13 Apr 2017 18:38:05: 2000000
INFO @ Thu, 13 Apr 2017 18:38:07: 3000000
INFO @ Thu, 13 Apr 2017 18:38:08: 4000000
INFO @ Thu, 13 Apr 2017 18:38:10: 5000000
INFO @ Thu, 13 Apr 2017 18:38:11: 6000000
INFO @ Thu, 13 Apr 2017 18:38:13: 7000000
INFO @ Thu, 13 Apr 2017 18:38:14: 8000000
INFO @ Thu, 13 Apr 2017 18:38:16: 9000000
INFO @ Thu, 13 Apr 2017 18:38:17: 10000000
INFO @ Thu, 13 Apr 2017 18:38:19: 11000000
INFO @ Thu, 13 Apr 2017 18:38:20: 12000000
INFO @ Thu, 13 Apr 2017 18:38:22: 13000000
INFO @ Thu, 13 Apr 2017 18:38:24: 14000000
INFO @ Thu, 13 Apr 2017 18:38:25: 15000000
INFO @ Thu, 13 Apr 2017 18:38:27: 16000000
INFO @ Thu, 13 Apr 2017 18:38:28: 17000000
INFO @ Thu, 13 Apr 2017 18:38:30: 18000000
INFO @ Thu, 13 Apr 2017 18:38:31: 19000000
INFO @ Thu, 13 Apr 2017 18:38:33: 20000000
INFO @ Thu, 13 Apr 2017 18:38:34: 21000000
INFO @ Thu, 13 Apr 2017 18:38:36: 22000000
INFO @ Thu, 13 Apr 2017 18:38:37: 23000000
INFO @ Thu, 13 Apr 2017 18:38:39: 24000000
INFO @ Thu, 13 Apr 2017 18:38:41: 25000000
INFO @ Thu, 13 Apr 2017 18:38:42: 26000000
INFO @ Thu, 13 Apr 2017 18:38:44: 27000000
INFO @ Thu, 13 Apr 2017 18:38:45: 28000000
INFO @ Thu, 13 Apr 2017 18:38:47: 29000000
INFO @ Thu, 13 Apr 2017 18:38:48: 30000000
INFO @ Thu, 13 Apr 2017 18:38:50: 31000000
INFO @ Thu, 13 Apr 2017 18:38:51: 32000000
INFO @ Thu, 13 Apr 2017 18:38:53: 33000000
INFO @ Thu, 13 Apr 2017 18:38:54: 34000000
INFO @ Thu, 13 Apr 2017 18:38:56: 35000000
INFO @ Thu, 13 Apr 2017 18:38:57: 36000000
INFO @ Thu, 13 Apr 2017 18:38:59: 37000000
INFO @ Thu, 13 Apr 2017 18:39:01: 38000000
INFO @ Thu, 13 Apr 2017 18:39:02: 39000000
INFO @ Thu, 13 Apr 2017 18:39:04: 40000000
INFO @ Thu, 13 Apr 2017 18:39:05: 41000000
INFO @ Thu, 13 Apr 2017 18:39:07: 42000000
INFO @ Thu, 13 Apr 2017 18:39:08: 43000000
INFO @ Thu, 13 Apr 2017 18:39:10: 44000000
INFO @ Thu, 13 Apr 2017 18:39:11: 45000000
INFO @ Thu, 13 Apr 2017 18:39:13: 46000000
INFO @ Thu, 13 Apr 2017 18:39:14: 47000000
INFO @ Thu, 13 Apr 2017 18:39:16: 48000000
INFO @ Thu, 13 Apr 2017 18:39:17: 49000000
INFO @ Thu, 13 Apr 2017 18:39:19: 50000000
INFO @ Thu, 13 Apr 2017 18:39:21: 51000000
INFO @ Thu, 13 Apr 2017 18:39:22: 52000000
INFO @ Thu, 13 Apr 2017 18:39:24: 53000000
INFO @ Thu, 13 Apr 2017 18:39:25: 54000000
INFO @ Thu, 13 Apr 2017 18:39:29: #1.2 read input tags...
INFO @ Thu, 13 Apr 2017 18:39:31: 1000000
INFO @ Thu, 13 Apr 2017 18:39:32: 2000000
INFO @ Thu, 13 Apr 2017 18:39:33: 3000000
INFO @ Thu, 13 Apr 2017 18:39:35: 4000000
INFO @ Thu, 13 Apr 2017 18:39:36: 5000000
INFO @ Thu, 13 Apr 2017 18:39:38: 6000000
INFO @ Thu, 13 Apr 2017 18:39:39: 7000000
INFO @ Thu, 13 Apr 2017 18:39:40: 8000000
INFO @ Thu, 13 Apr 2017 18:39:42: 9000000
INFO @ Thu, 13 Apr 2017 18:39:43: 10000000
INFO @ Thu, 13 Apr 2017 18:39:45: 11000000
INFO @ Thu, 13 Apr 2017 18:39:46: 12000000
INFO @ Thu, 13 Apr 2017 18:39:48: 13000000
INFO @ Thu, 13 Apr 2017 18:39:49: 14000000
INFO @ Thu, 13 Apr 2017 18:39:50: 15000000
INFO @ Thu, 13 Apr 2017 18:39:52: 16000000
INFO @ Thu, 13 Apr 2017 18:39:53: 17000000
INFO @ Thu, 13 Apr 2017 18:39:55: 18000000
INFO @ Thu, 13 Apr 2017 18:39:56: 19000000
INFO @ Thu, 13 Apr 2017 18:39:58: 20000000
INFO @ Thu, 13 Apr 2017 18:39:59: 21000000
INFO @ Thu, 13 Apr 2017 18:40:00: 22000000
INFO @ Thu, 13 Apr 2017 18:40:02: 23000000
INFO @ Thu, 13 Apr 2017 18:40:03: 24000000
INFO @ Thu, 13 Apr 2017 18:40:05: 25000000
INFO @ Thu, 13 Apr 2017 18:40:06: 26000000
INFO @ Thu, 13 Apr 2017 18:40:08: 27000000
INFO @ Thu, 13 Apr 2017 18:40:09: 28000000
INFO @ Thu, 13 Apr 2017 18:40:10: 29000000
INFO @ Thu, 13 Apr 2017 18:40:12: 30000000
INFO @ Thu, 13 Apr 2017 18:40:13: 31000000
INFO @ Thu, 13 Apr 2017 18:40:15: 32000000
INFO @ Thu, 13 Apr 2017 18:40:16: 33000000
INFO @ Thu, 13 Apr 2017 18:40:18: 34000000
INFO @ Thu, 13 Apr 2017 18:40:19: 35000000
INFO @ Thu, 13 Apr 2017 18:40:20: 36000000
INFO @ Thu, 13 Apr 2017 18:40:22: 37000000
INFO @ Thu, 13 Apr 2017 18:40:23: 38000000
INFO @ Thu, 13 Apr 2017 18:40:25: 39000000
INFO @ Thu, 13 Apr 2017 18:40:26: 40000000
INFO @ Thu, 13 Apr 2017 18:40:28: 41000000
INFO @ Thu, 13 Apr 2017 18:40:29: 42000000
INFO @ Thu, 13 Apr 2017 18:40:30: 43000000
INFO @ Thu, 13 Apr 2017 18:40:32: 44000000
INFO @ Thu, 13 Apr 2017 18:40:33: 45000000
INFO @ Thu, 13 Apr 2017 18:40:35: 46000000
INFO @ Thu, 13 Apr 2017 18:40:36: 47000000
INFO @ Thu, 13 Apr 2017 18:40:38: 48000000
INFO @ Thu, 13 Apr 2017 18:40:39: 49000000
INFO @ Thu, 13 Apr 2017 18:40:40: 50000000
INFO @ Thu, 13 Apr 2017 18:40:42: 51000000
INFO @ Thu, 13 Apr 2017 18:40:43: 52000000
INFO @ Thu, 13 Apr 2017 18:40:45: 53000000
INFO @ Thu, 13 Apr 2017 18:40:46: 54000000
INFO @ Thu, 13 Apr 2017 18:40:48: 55000000
INFO @ Thu, 13 Apr 2017 18:40:49: 56000000
INFO @ Thu, 13 Apr 2017 18:40:50: 57000000
INFO @ Thu, 13 Apr 2017 18:40:52: 58000000
INFO @ Thu, 13 Apr 2017 18:40:53: 59000000
INFO @ Thu, 13 Apr 2017 18:40:55: 60000000
INFO @ Thu, 13 Apr 2017 18:40:56: 61000000
INFO @ Thu, 13 Apr 2017 18:40:58: 62000000
INFO @ Thu, 13 Apr 2017 18:40:59: 63000000
INFO @ Thu, 13 Apr 2017 18:41:00: 64000000
INFO @ Thu, 13 Apr 2017 18:41:02: 65000000
INFO @ Thu, 13 Apr 2017 18:41:03: 66000000
INFO @ Thu, 13 Apr 2017 18:41:05: 67000000
INFO @ Thu, 13 Apr 2017 18:41:06: 68000000
INFO @ Thu, 13 Apr 2017 18:41:08: 69000000
INFO @ Thu, 13 Apr 2017 18:41:09: 70000000
INFO @ Thu, 13 Apr 2017 18:41:11: 71000000
INFO @ Thu, 13 Apr 2017 18:41:12: 72000000
INFO @ Thu, 13 Apr 2017 18:41:13: 73000000
INFO @ Thu, 13 Apr 2017 18:41:15: 74000000
INFO @ Thu, 13 Apr 2017 18:41:16: 75000000
INFO @ Thu, 13 Apr 2017 18:41:18: 76000000
INFO @ Thu, 13 Apr 2017 18:41:19: 77000000
INFO @ Thu, 13 Apr 2017 18:41:21: 78000000
INFO @ Thu, 13 Apr 2017 18:41:22: 79000000
INFO @ Thu, 13 Apr 2017 18:41:23: 80000000
INFO @ Thu, 13 Apr 2017 18:41:25: 81000000
INFO @ Thu, 13 Apr 2017 18:41:26: 82000000
INFO @ Thu, 13 Apr 2017 18:41:28: 83000000
INFO @ Thu, 13 Apr 2017 18:41:29: 84000000
INFO @ Thu, 13 Apr 2017 18:41:31: 85000000
INFO @ Thu, 13 Apr 2017 18:41:32: 86000000
INFO @ Thu, 13 Apr 2017 18:41:33: 87000000
INFO @ Thu, 13 Apr 2017 18:41:35: 88000000
INFO @ Thu, 13 Apr 2017 18:41:36: 89000000
INFO @ Thu, 13 Apr 2017 18:41:38: 90000000
INFO @ Thu, 13 Apr 2017 18:41:39: 91000000
INFO @ Thu, 13 Apr 2017 18:41:41: 92000000
INFO @ Thu, 13 Apr 2017 18:41:42: 93000000
INFO @ Thu, 13 Apr 2017 18:41:43: 94000000
INFO @ Thu, 13 Apr 2017 18:41:45: 95000000
INFO @ Thu, 13 Apr 2017 18:41:46: 96000000
INFO @ Thu, 13 Apr 2017 18:41:48: 97000000
INFO @ Thu, 13 Apr 2017 18:41:49: 98000000
INFO @ Thu, 13 Apr 2017 18:41:50: 99000000
INFO @ Thu, 13 Apr 2017 18:41:52: #1 tag size is determined as 101 bps
INFO @ Thu, 13 Apr 2017 18:41:52: #1 tag size = 101
INFO @ Thu, 13 Apr 2017 18:41:52: #1 total tags in treatment: 54211426
INFO @ Thu, 13 Apr 2017 18:41:52: #1 total tags in control: 99493262
INFO @ Thu, 13 Apr 2017 18:41:52: #1 finished!
INFO @ Thu, 13 Apr 2017 18:41:52: #2 Build Peak Model...
INFO @ Thu, 13 Apr 2017 18:41:52: #2 Skipped...
INFO @ Thu, 13 Apr 2017 18:41:52: #2 Use 213 as fragment length
INFO @ Thu, 13 Apr 2017 18:41:52: #3 Call peaks...
INFO @ Thu, 13 Apr 2017 18:41:52: #3 Call broad peaks with given level1 -log10pvalue cutoff and level2: 2.00000, 1.00000...
INFO @ Thu, 13 Apr 2017 18:41:52: #3 Pre-compute pvalue-qvalue table...
INFO @ Thu, 13 Apr 2017 18:47:19: #3 Call peaks for each chromosome...
INFO @ Thu, 13 Apr 2017 18:51:22: #4 Write output xls file... /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pooled_pseudo_reps/ppr1/ENCFF558PHY.pr1_ENCFF645WYX.pr1.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign_peaks.xls
INFO @ Thu, 13 Apr 2017 18:51:23: #4 Write broad peak in broadPeak format file... /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pooled_pseudo_reps/ppr1/ENCFF558PHY.pr1_ENCFF645WYX.pr1.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign_peaks.broadPeak
INFO @ Thu, 13 Apr 2017 18:51:24: #4 Write broad peak in bed12/gappedPeak format file... /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pooled_pseudo_reps/ppr1/ENCFF558PHY.pr1_ENCFF645WYX.pr1.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign_peaks.gappedPeak
INFO @ Thu, 13 Apr 2017 18:51:24: Done!
|
Num | 3 | ID | task.callpeak_macs2.macs2_g_b_ppr2.line_155.id_12 | Name | macs2 g/b ppr2 | Thread | thread_Root | PID | 46969 | | OK | true | Exit Code | 0 | Retries | | State | FINISHED | Dep. | OK | Cpus | | Mem | | | Start | 2017-04-13 18:37:51 | End | 2017-04-13 18:51:37 | Elapsed | 00:13:45 | Timeout | 00:00:-1 | Wall Timeout | 100 days | | Input files | /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/align/pooled_pseudo_reps/ppr2/ENCFF558PHY.pr2_ENCFF645WYX.pr2.tagAlign.gz
/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/InputData/ss_50M_2000_GRCh38.nodup.tagAlign.gz
| Output files | /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pooled_pseudo_reps/ppr2/ENCFF558PHY.pr2_ENCFF645WYX.pr2.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.gappedPeak.gz
/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pooled_pseudo_reps/ppr2/ENCFF558PHY.pr2_ENCFF645WYX.pr2.tagAlign_x_ss_50M_2000_GRCh38.nodup.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/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/align/pooled_pseudo_reps/ppr2/ENCFF558PHY.pr2_ENCFF645WYX.pr2.tagAlign.gz -c /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/InputData/ss_50M_2000_GRCh38.nodup.tagAlign.gz -f BED -n /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pooled_pseudo_reps/ppr2/ENCFF558PHY.pr2_ENCFF645WYX.pr2.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign -g hs -p 0.01 --broad --nomodel --shift 0 --extsize 213 --keep-dup all
# SYS command. line 166
sort -k 8gr,8gr "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pooled_pseudo_reps/ppr2/ENCFF558PHY.pr2_ENCFF645WYX.pr2.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign"_peaks.broadPeak | awk 'BEGIN{OFS="\t"}{$4="Peak_"NR ; print $0}' | gzip -nc > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pooled_pseudo_reps/ppr2/ENCFF558PHY.pr2_ENCFF645WYX.pr2.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.broadPeak.gz
# SYS command. line 167
sort -k 14gr,14gr "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pooled_pseudo_reps/ppr2/ENCFF558PHY.pr2_ENCFF645WYX.pr2.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign"_peaks.gappedPeak | awk 'BEGIN{OFS="\t"}{$4="Peak_"NR ; print $0}' | gzip -nc > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pooled_pseudo_reps/ppr2/ENCFF558PHY.pr2_ENCFF645WYX.pr2.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.gappedPeak.gz
# SYS command. line 170
rm -f "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pooled_pseudo_reps/ppr2/ENCFF558PHY.pr2_ENCFF645WYX.pr2.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign"_peaks.xls \
"/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pooled_pseudo_reps/ppr2/ENCFF558PHY.pr2_ENCFF645WYX.pr2.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign"_peaks.broadPeak \
"/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pooled_pseudo_reps/ppr2/ENCFF558PHY.pr2_ENCFF645WYX.pr2.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign"_peaks.gappedPeak \
"/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pooled_pseudo_reps/ppr2/ENCFF558PHY.pr2_ENCFF645WYX.pr2.tagAlign_x_ss_50M_2000_GRCh38.nodup.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--------------------
48754 (process ID) old priority 0, new priority 19
--------------------Stderr--------------------
INFO @ Thu, 13 Apr 2017 18:38:16:
# Command line: callpeak -t /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/align/pooled_pseudo_reps/ppr2/ENCFF558PHY.pr2_ENCFF645WYX.pr2.tagAlign.gz -c /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/InputData/ss_50M_2000_GRCh38.nodup.tagAlign.gz -f BED -n /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pooled_pseudo_reps/ppr2/ENCFF558PHY.pr2_ENCFF645WYX.pr2.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign -g hs -p 0.01 --broad --nomodel --shift 0 --extsize 213 --keep-dup all
# ARGUMENTS LIST:
# name = /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pooled_pseudo_reps/ppr2/ENCFF558PHY.pr2_ENCFF645WYX.pr2.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign
# format = BED
# ChIP-seq file = ['/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/align/pooled_pseudo_reps/ppr2/ENCFF558PHY.pr2_ENCFF645WYX.pr2.tagAlign.gz']
# control file = ['/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/InputData/ss_50M_2000_GRCh38.nodup.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:16: #1 read tag files...
INFO @ Thu, 13 Apr 2017 18:38:16: #1 read treatment tags...
INFO @ Thu, 13 Apr 2017 18:38:17: 1000000
INFO @ Thu, 13 Apr 2017 18:38:19: 2000000
INFO @ Thu, 13 Apr 2017 18:38:20: 3000000
INFO @ Thu, 13 Apr 2017 18:38:22: 4000000
INFO @ Thu, 13 Apr 2017 18:38:23: 5000000
INFO @ Thu, 13 Apr 2017 18:38:25: 6000000
INFO @ Thu, 13 Apr 2017 18:38:26: 7000000
INFO @ Thu, 13 Apr 2017 18:38:28: 8000000
INFO @ Thu, 13 Apr 2017 18:38:29: 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:45: 19000000
INFO @ Thu, 13 Apr 2017 18:38:46: 20000000
INFO @ Thu, 13 Apr 2017 18:38:48: 21000000
INFO @ Thu, 13 Apr 2017 18:38:49: 22000000
INFO @ Thu, 13 Apr 2017 18:38:51: 23000000
INFO @ Thu, 13 Apr 2017 18:38:52: 24000000
INFO @ Thu, 13 Apr 2017 18:38:54: 25000000
INFO @ Thu, 13 Apr 2017 18:38:55: 26000000
INFO @ Thu, 13 Apr 2017 18:38:57: 27000000
INFO @ Thu, 13 Apr 2017 18:38:58: 28000000
INFO @ Thu, 13 Apr 2017 18:39:00: 29000000
INFO @ Thu, 13 Apr 2017 18:39:01: 30000000
INFO @ Thu, 13 Apr 2017 18:39:03: 31000000
INFO @ Thu, 13 Apr 2017 18:39:04: 32000000
INFO @ Thu, 13 Apr 2017 18:39:06: 33000000
INFO @ Thu, 13 Apr 2017 18:39:08: 34000000
INFO @ Thu, 13 Apr 2017 18:39:09: 35000000
INFO @ Thu, 13 Apr 2017 18:39:11: 36000000
INFO @ Thu, 13 Apr 2017 18:39:12: 37000000
INFO @ Thu, 13 Apr 2017 18:39:14: 38000000
INFO @ Thu, 13 Apr 2017 18:39:15: 39000000
INFO @ Thu, 13 Apr 2017 18:39:17: 40000000
INFO @ Thu, 13 Apr 2017 18:39:18: 41000000
INFO @ Thu, 13 Apr 2017 18:39:20: 42000000
INFO @ Thu, 13 Apr 2017 18:39:21: 43000000
INFO @ Thu, 13 Apr 2017 18:39:23: 44000000
INFO @ Thu, 13 Apr 2017 18:39:24: 45000000
INFO @ Thu, 13 Apr 2017 18:39:26: 46000000
INFO @ Thu, 13 Apr 2017 18:39:27: 47000000
INFO @ Thu, 13 Apr 2017 18:39:29: 48000000
INFO @ Thu, 13 Apr 2017 18:39:30: 49000000
INFO @ Thu, 13 Apr 2017 18:39:32: 50000000
INFO @ Thu, 13 Apr 2017 18:39:33: 51000000
INFO @ Thu, 13 Apr 2017 18:39:35: 52000000
INFO @ Thu, 13 Apr 2017 18:39:37: 53000000
INFO @ Thu, 13 Apr 2017 18:39:38: 54000000
INFO @ Thu, 13 Apr 2017 18:39:42: #1.2 read input tags...
INFO @ Thu, 13 Apr 2017 18:39:43: 1000000
INFO @ Thu, 13 Apr 2017 18:39:45: 2000000
INFO @ Thu, 13 Apr 2017 18:39:46: 3000000
INFO @ Thu, 13 Apr 2017 18:39:47: 4000000
INFO @ Thu, 13 Apr 2017 18:39:49: 5000000
INFO @ Thu, 13 Apr 2017 18:39:50: 6000000
INFO @ Thu, 13 Apr 2017 18:39:52: 7000000
INFO @ Thu, 13 Apr 2017 18:39:53: 8000000
INFO @ Thu, 13 Apr 2017 18:39:54: 9000000
INFO @ Thu, 13 Apr 2017 18:39:56: 10000000
INFO @ Thu, 13 Apr 2017 18:39:57: 11000000
INFO @ Thu, 13 Apr 2017 18:39:59: 12000000
INFO @ Thu, 13 Apr 2017 18:40:00: 13000000
INFO @ Thu, 13 Apr 2017 18:40:01: 14000000
INFO @ Thu, 13 Apr 2017 18:40:03: 15000000
INFO @ Thu, 13 Apr 2017 18:40:04: 16000000
INFO @ Thu, 13 Apr 2017 18:40:05: 17000000
INFO @ Thu, 13 Apr 2017 18:40:07: 18000000
INFO @ Thu, 13 Apr 2017 18:40:08: 19000000
INFO @ Thu, 13 Apr 2017 18:40:10: 20000000
INFO @ Thu, 13 Apr 2017 18:40:11: 21000000
INFO @ Thu, 13 Apr 2017 18:40:12: 22000000
INFO @ Thu, 13 Apr 2017 18:40:14: 23000000
INFO @ Thu, 13 Apr 2017 18:40:15: 24000000
INFO @ Thu, 13 Apr 2017 18:40:17: 25000000
INFO @ Thu, 13 Apr 2017 18:40:18: 26000000
INFO @ Thu, 13 Apr 2017 18:40:19: 27000000
INFO @ Thu, 13 Apr 2017 18:40:21: 28000000
INFO @ Thu, 13 Apr 2017 18:40:22: 29000000
INFO @ Thu, 13 Apr 2017 18:40:23: 30000000
INFO @ Thu, 13 Apr 2017 18:40:25: 31000000
INFO @ Thu, 13 Apr 2017 18:40:26: 32000000
INFO @ Thu, 13 Apr 2017 18:40:28: 33000000
INFO @ Thu, 13 Apr 2017 18:40:29: 34000000
INFO @ Thu, 13 Apr 2017 18:40:30: 35000000
INFO @ Thu, 13 Apr 2017 18:40:32: 36000000
INFO @ Thu, 13 Apr 2017 18:40:33: 37000000
INFO @ Thu, 13 Apr 2017 18:40:35: 38000000
INFO @ Thu, 13 Apr 2017 18:40:36: 39000000
INFO @ Thu, 13 Apr 2017 18:40:37: 40000000
INFO @ Thu, 13 Apr 2017 18:40:39: 41000000
INFO @ Thu, 13 Apr 2017 18:40:40: 42000000
INFO @ Thu, 13 Apr 2017 18:40:41: 43000000
INFO @ Thu, 13 Apr 2017 18:40:43: 44000000
INFO @ Thu, 13 Apr 2017 18:40:44: 45000000
INFO @ Thu, 13 Apr 2017 18:40:46: 46000000
INFO @ Thu, 13 Apr 2017 18:40:47: 47000000
INFO @ Thu, 13 Apr 2017 18:40:48: 48000000
INFO @ Thu, 13 Apr 2017 18:40:50: 49000000
INFO @ Thu, 13 Apr 2017 18:40:51: 50000000
INFO @ Thu, 13 Apr 2017 18:40:53: 51000000
INFO @ Thu, 13 Apr 2017 18:40:54: 52000000
INFO @ Thu, 13 Apr 2017 18:40:55: 53000000
INFO @ Thu, 13 Apr 2017 18:40:57: 54000000
INFO @ Thu, 13 Apr 2017 18:40:58: 55000000
INFO @ Thu, 13 Apr 2017 18:41:00: 56000000
INFO @ Thu, 13 Apr 2017 18:41:01: 57000000
INFO @ Thu, 13 Apr 2017 18:41:02: 58000000
INFO @ Thu, 13 Apr 2017 18:41:04: 59000000
INFO @ Thu, 13 Apr 2017 18:41:05: 60000000
INFO @ Thu, 13 Apr 2017 18:41:06: 61000000
INFO @ Thu, 13 Apr 2017 18:41:08: 62000000
INFO @ Thu, 13 Apr 2017 18:41:09: 63000000
INFO @ Thu, 13 Apr 2017 18:41:11: 64000000
INFO @ Thu, 13 Apr 2017 18:41:12: 65000000
INFO @ Thu, 13 Apr 2017 18:41:13: 66000000
INFO @ Thu, 13 Apr 2017 18:41:15: 67000000
INFO @ Thu, 13 Apr 2017 18:41:16: 68000000
INFO @ Thu, 13 Apr 2017 18:41:18: 69000000
INFO @ Thu, 13 Apr 2017 18:41:19: 70000000
INFO @ Thu, 13 Apr 2017 18:41:20: 71000000
INFO @ Thu, 13 Apr 2017 18:41:22: 72000000
INFO @ Thu, 13 Apr 2017 18:41:23: 73000000
INFO @ Thu, 13 Apr 2017 18:41:25: 74000000
INFO @ Thu, 13 Apr 2017 18:41:26: 75000000
INFO @ Thu, 13 Apr 2017 18:41:27: 76000000
INFO @ Thu, 13 Apr 2017 18:41:29: 77000000
INFO @ Thu, 13 Apr 2017 18:41:30: 78000000
INFO @ Thu, 13 Apr 2017 18:41:31: 79000000
INFO @ Thu, 13 Apr 2017 18:41:33: 80000000
INFO @ Thu, 13 Apr 2017 18:41:34: 81000000
INFO @ Thu, 13 Apr 2017 18:41:36: 82000000
INFO @ Thu, 13 Apr 2017 18:41:37: 83000000
INFO @ Thu, 13 Apr 2017 18:41:38: 84000000
INFO @ Thu, 13 Apr 2017 18:41:40: 85000000
INFO @ Thu, 13 Apr 2017 18:41:41: 86000000
INFO @ Thu, 13 Apr 2017 18:41:43: 87000000
INFO @ Thu, 13 Apr 2017 18:41:44: 88000000
INFO @ Thu, 13 Apr 2017 18:41:45: 89000000
INFO @ Thu, 13 Apr 2017 18:41:47: 90000000
INFO @ Thu, 13 Apr 2017 18:41:48: 91000000
INFO @ Thu, 13 Apr 2017 18:41:50: 92000000
INFO @ Thu, 13 Apr 2017 18:41:51: 93000000
INFO @ Thu, 13 Apr 2017 18:41:52: 94000000
INFO @ Thu, 13 Apr 2017 18:41:54: 95000000
INFO @ Thu, 13 Apr 2017 18:41:55: 96000000
INFO @ Thu, 13 Apr 2017 18:41:56: 97000000
INFO @ Thu, 13 Apr 2017 18:41:58: 98000000
INFO @ Thu, 13 Apr 2017 18:41:59: 99000000
INFO @ Thu, 13 Apr 2017 18:42:01: #1 tag size is determined as 101 bps
INFO @ Thu, 13 Apr 2017 18:42:01: #1 tag size = 101
INFO @ Thu, 13 Apr 2017 18:42:01: #1 total tags in treatment: 54211426
INFO @ Thu, 13 Apr 2017 18:42:01: #1 total tags in control: 99493262
INFO @ Thu, 13 Apr 2017 18:42:01: #1 finished!
INFO @ Thu, 13 Apr 2017 18:42:01: #2 Build Peak Model...
INFO @ Thu, 13 Apr 2017 18:42:01: #2 Skipped...
INFO @ Thu, 13 Apr 2017 18:42:01: #2 Use 213 as fragment length
INFO @ Thu, 13 Apr 2017 18:42:01: #3 Call peaks...
INFO @ Thu, 13 Apr 2017 18:42:01: #3 Call broad peaks with given level1 -log10pvalue cutoff and level2: 2.00000, 1.00000...
INFO @ Thu, 13 Apr 2017 18:42:01: #3 Pre-compute pvalue-qvalue table...
INFO @ Thu, 13 Apr 2017 18:47:26: #3 Call peaks for each chromosome...
INFO @ Thu, 13 Apr 2017 18:51:29: #4 Write output xls file... /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pooled_pseudo_reps/ppr2/ENCFF558PHY.pr2_ENCFF645WYX.pr2.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign_peaks.xls
INFO @ Thu, 13 Apr 2017 18:51:30: #4 Write broad peak in broadPeak format file... /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pooled_pseudo_reps/ppr2/ENCFF558PHY.pr2_ENCFF645WYX.pr2.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign_peaks.broadPeak
INFO @ Thu, 13 Apr 2017 18:51:31: #4 Write broad peak in bed12/gappedPeak format file... /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pooled_pseudo_reps/ppr2/ENCFF558PHY.pr2_ENCFF645WYX.pr2.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign_peaks.gappedPeak
INFO @ Thu, 13 Apr 2017 18:51:32: Done!
|
Num | 4 | ID | task.callpeak_macs2.macs2_g_b_rep1.line_155.id_13 | Name | macs2 g/b rep1 | Thread | thread_Root | PID | 46975 | | OK | true | Exit Code | 0 | Retries | | State | FINISHED | Dep. | OK | Cpus | | Mem | | | Start | 2017-04-13 18:37:54 | End | 2017-04-13 18:52:14 | Elapsed | 00:14:20 | Timeout | 00:00:-1 | Wall Timeout | 100 days | | Input files | /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/align/rep1/ENCFF558PHY.tagAlign.gz
/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/InputData/ss_50M_2000_GRCh38.nodup.tagAlign.gz
| Output files | /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/rep1/ENCFF558PHY.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.gappedPeak.gz
/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/rep1/ENCFF558PHY.tagAlign_x_ss_50M_2000_GRCh38.nodup.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/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/align/rep1/ENCFF558PHY.tagAlign.gz -c /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/InputData/ss_50M_2000_GRCh38.nodup.tagAlign.gz -f BED -n /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/rep1/ENCFF558PHY.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign -g hs -p 0.01 --broad --nomodel --shift 0 --extsize 215 --keep-dup all
# SYS command. line 166
sort -k 8gr,8gr "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/rep1/ENCFF558PHY.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign"_peaks.broadPeak | awk 'BEGIN{OFS="\t"}{$4="Peak_"NR ; print $0}' | gzip -nc > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/rep1/ENCFF558PHY.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.broadPeak.gz
# SYS command. line 167
sort -k 14gr,14gr "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/rep1/ENCFF558PHY.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign"_peaks.gappedPeak | awk 'BEGIN{OFS="\t"}{$4="Peak_"NR ; print $0}' | gzip -nc > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/rep1/ENCFF558PHY.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.gappedPeak.gz
# SYS command. line 170
rm -f "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/rep1/ENCFF558PHY.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign"_peaks.xls \
"/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/rep1/ENCFF558PHY.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign"_peaks.broadPeak \
"/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/rep1/ENCFF558PHY.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign"_peaks.gappedPeak \
"/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/rep1/ENCFF558PHY.tagAlign_x_ss_50M_2000_GRCh38.nodup.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--------------------
48756 (process ID) old priority 0, new priority 19
--------------------Stderr--------------------
INFO @ Thu, 13 Apr 2017 18:38:16:
# Command line: callpeak -t /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/align/rep1/ENCFF558PHY.tagAlign.gz -c /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/InputData/ss_50M_2000_GRCh38.nodup.tagAlign.gz -f BED -n /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/rep1/ENCFF558PHY.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign -g hs -p 0.01 --broad --nomodel --shift 0 --extsize 215 --keep-dup all
# ARGUMENTS LIST:
# name = /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/rep1/ENCFF558PHY.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign
# format = BED
# ChIP-seq file = ['/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/align/rep1/ENCFF558PHY.tagAlign.gz']
# control file = ['/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/InputData/ss_50M_2000_GRCh38.nodup.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:16: #1 read tag files...
INFO @ Thu, 13 Apr 2017 18:38:16: #1 read treatment tags...
INFO @ Thu, 13 Apr 2017 18:38:17: 1000000
INFO @ Thu, 13 Apr 2017 18:38:19: 2000000
INFO @ Thu, 13 Apr 2017 18:38:21: 3000000
INFO @ Thu, 13 Apr 2017 18:38:22: 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:29: 8000000
INFO @ Thu, 13 Apr 2017 18:38:30: 9000000
INFO @ Thu, 13 Apr 2017 18:38:32: 10000000
INFO @ Thu, 13 Apr 2017 18:38:34: 11000000
INFO @ Thu, 13 Apr 2017 18:38:35: 12000000
INFO @ Thu, 13 Apr 2017 18:38:37: 13000000
INFO @ Thu, 13 Apr 2017 18:38:38: 14000000
INFO @ Thu, 13 Apr 2017 18:38:40: 15000000
INFO @ Thu, 13 Apr 2017 18:38:41: 16000000
INFO @ Thu, 13 Apr 2017 18:38:43: 17000000
INFO @ Thu, 13 Apr 2017 18:38:44: 18000000
INFO @ Thu, 13 Apr 2017 18:38:46: 19000000
INFO @ Thu, 13 Apr 2017 18:38:47: 20000000
INFO @ Thu, 13 Apr 2017 18:38:49: 21000000
INFO @ Thu, 13 Apr 2017 18:38:50: 22000000
INFO @ Thu, 13 Apr 2017 18:38:52: 23000000
INFO @ Thu, 13 Apr 2017 18:38:53: 24000000
INFO @ Thu, 13 Apr 2017 18:38:55: 25000000
INFO @ Thu, 13 Apr 2017 18:38:56: 26000000
INFO @ Thu, 13 Apr 2017 18:38:58: 27000000
INFO @ Thu, 13 Apr 2017 18:38:59: 28000000
INFO @ Thu, 13 Apr 2017 18:39:01: 29000000
INFO @ Thu, 13 Apr 2017 18:39:02: 30000000
INFO @ Thu, 13 Apr 2017 18:39:04: 31000000
INFO @ Thu, 13 Apr 2017 18:39:05: 32000000
INFO @ Thu, 13 Apr 2017 18:39:07: 33000000
INFO @ Thu, 13 Apr 2017 18:39:08: 34000000
INFO @ Thu, 13 Apr 2017 18:39:10: 35000000
INFO @ Thu, 13 Apr 2017 18:39:12: 36000000
INFO @ Thu, 13 Apr 2017 18:39:14: 37000000
INFO @ Thu, 13 Apr 2017 18:39:15: 38000000
INFO @ Thu, 13 Apr 2017 18:39:17: 39000000
INFO @ Thu, 13 Apr 2017 18:39:18: 40000000
INFO @ Thu, 13 Apr 2017 18:39:20: 41000000
INFO @ Thu, 13 Apr 2017 18:39:21: 42000000
INFO @ Thu, 13 Apr 2017 18:39:23: 43000000
INFO @ Thu, 13 Apr 2017 18:39:25: 44000000
INFO @ Thu, 13 Apr 2017 18:39:26: 45000000
INFO @ Thu, 13 Apr 2017 18:39:28: 46000000
INFO @ Thu, 13 Apr 2017 18:39:30: 47000000
INFO @ Thu, 13 Apr 2017 18:39:31: 48000000
INFO @ Thu, 13 Apr 2017 18:39:33: 49000000
INFO @ Thu, 13 Apr 2017 18:39:34: 50000000
INFO @ Thu, 13 Apr 2017 18:39:36: 51000000
INFO @ Thu, 13 Apr 2017 18:39:38: 52000000
INFO @ Thu, 13 Apr 2017 18:39:39: 53000000
INFO @ Thu, 13 Apr 2017 18:39:41: 54000000
INFO @ Thu, 13 Apr 2017 18:39:45: #1.2 read input tags...
INFO @ Thu, 13 Apr 2017 18:39:47: 1000000
INFO @ Thu, 13 Apr 2017 18:39:48: 2000000
INFO @ Thu, 13 Apr 2017 18:39:50: 3000000
INFO @ Thu, 13 Apr 2017 18:39:51: 4000000
INFO @ Thu, 13 Apr 2017 18:39:52: 5000000
INFO @ Thu, 13 Apr 2017 18:39:54: 6000000
INFO @ Thu, 13 Apr 2017 18:39:55: 7000000
INFO @ Thu, 13 Apr 2017 18:39:57: 8000000
INFO @ Thu, 13 Apr 2017 18:39:58: 9000000
INFO @ Thu, 13 Apr 2017 18:39:59: 10000000
INFO @ Thu, 13 Apr 2017 18:40:01: 11000000
INFO @ Thu, 13 Apr 2017 18:40:02: 12000000
INFO @ Thu, 13 Apr 2017 18:40:03: 13000000
INFO @ Thu, 13 Apr 2017 18:40:05: 14000000
INFO @ Thu, 13 Apr 2017 18:40:06: 15000000
INFO @ Thu, 13 Apr 2017 18:40:08: 16000000
INFO @ Thu, 13 Apr 2017 18:40:09: 17000000
INFO @ Thu, 13 Apr 2017 18:40:10: 18000000
INFO @ Thu, 13 Apr 2017 18:40:12: 19000000
INFO @ Thu, 13 Apr 2017 18:40:13: 20000000
INFO @ Thu, 13 Apr 2017 18:40:15: 21000000
INFO @ Thu, 13 Apr 2017 18:40:16: 22000000
INFO @ Thu, 13 Apr 2017 18:40:18: 23000000
INFO @ Thu, 13 Apr 2017 18:40:19: 24000000
INFO @ Thu, 13 Apr 2017 18:40:21: 25000000
INFO @ Thu, 13 Apr 2017 18:40:22: 26000000
INFO @ Thu, 13 Apr 2017 18:40:23: 27000000
INFO @ Thu, 13 Apr 2017 18:40:25: 28000000
INFO @ Thu, 13 Apr 2017 18:40:26: 29000000
INFO @ Thu, 13 Apr 2017 18:40:27: 30000000
INFO @ Thu, 13 Apr 2017 18:40:29: 31000000
INFO @ Thu, 13 Apr 2017 18:40:30: 32000000
INFO @ Thu, 13 Apr 2017 18:40:32: 33000000
INFO @ Thu, 13 Apr 2017 18:40:33: 34000000
INFO @ Thu, 13 Apr 2017 18:40:35: 35000000
INFO @ Thu, 13 Apr 2017 18:40:36: 36000000
INFO @ Thu, 13 Apr 2017 18:40:37: 37000000
INFO @ Thu, 13 Apr 2017 18:40:39: 38000000
INFO @ Thu, 13 Apr 2017 18:40:40: 39000000
INFO @ Thu, 13 Apr 2017 18:40:42: 40000000
INFO @ Thu, 13 Apr 2017 18:40:43: 41000000
INFO @ Thu, 13 Apr 2017 18:40:45: 42000000
INFO @ Thu, 13 Apr 2017 18:40:46: 43000000
INFO @ Thu, 13 Apr 2017 18:40:47: 44000000
INFO @ Thu, 13 Apr 2017 18:40:49: 45000000
INFO @ Thu, 13 Apr 2017 18:40:51: 46000000
INFO @ Thu, 13 Apr 2017 18:40:52: 47000000
INFO @ Thu, 13 Apr 2017 18:40:53: 48000000
INFO @ Thu, 13 Apr 2017 18:40:55: 49000000
INFO @ Thu, 13 Apr 2017 18:40:56: 50000000
INFO @ Thu, 13 Apr 2017 18:40:57: 51000000
INFO @ Thu, 13 Apr 2017 18:40:59: 52000000
INFO @ Thu, 13 Apr 2017 18:41:00: 53000000
INFO @ Thu, 13 Apr 2017 18:41:02: 54000000
INFO @ Thu, 13 Apr 2017 18:41:03: 55000000
INFO @ Thu, 13 Apr 2017 18:41:05: 56000000
INFO @ Thu, 13 Apr 2017 18:41:07: 57000000
INFO @ Thu, 13 Apr 2017 18:41:08: 58000000
INFO @ Thu, 13 Apr 2017 18:41:10: 59000000
INFO @ Thu, 13 Apr 2017 18:41:11: 60000000
INFO @ Thu, 13 Apr 2017 18:41:12: 61000000
INFO @ Thu, 13 Apr 2017 18:41:14: 62000000
INFO @ Thu, 13 Apr 2017 18:41:15: 63000000
INFO @ Thu, 13 Apr 2017 18:41:17: 64000000
INFO @ Thu, 13 Apr 2017 18:41:18: 65000000
INFO @ Thu, 13 Apr 2017 18:41:19: 66000000
INFO @ Thu, 13 Apr 2017 18:41:21: 67000000
INFO @ Thu, 13 Apr 2017 18:41:22: 68000000
INFO @ Thu, 13 Apr 2017 18:41:24: 69000000
INFO @ Thu, 13 Apr 2017 18:41:25: 70000000
INFO @ Thu, 13 Apr 2017 18:41:26: 71000000
INFO @ Thu, 13 Apr 2017 18:41:28: 72000000
INFO @ Thu, 13 Apr 2017 18:41:29: 73000000
INFO @ Thu, 13 Apr 2017 18:41:30: 74000000
INFO @ Thu, 13 Apr 2017 18:41:32: 75000000
INFO @ Thu, 13 Apr 2017 18:41:33: 76000000
INFO @ Thu, 13 Apr 2017 18:41:35: 77000000
INFO @ Thu, 13 Apr 2017 18:41:36: 78000000
INFO @ Thu, 13 Apr 2017 18:41:38: 79000000
INFO @ Thu, 13 Apr 2017 18:41:39: 80000000
INFO @ Thu, 13 Apr 2017 18:41:40: 81000000
INFO @ Thu, 13 Apr 2017 18:41:42: 82000000
INFO @ Thu, 13 Apr 2017 18:41:43: 83000000
INFO @ Thu, 13 Apr 2017 18:41:45: 84000000
INFO @ Thu, 13 Apr 2017 18:41:46: 85000000
INFO @ Thu, 13 Apr 2017 18:41:47: 86000000
INFO @ Thu, 13 Apr 2017 18:41:49: 87000000
INFO @ Thu, 13 Apr 2017 18:41:50: 88000000
INFO @ Thu, 13 Apr 2017 18:41:52: 89000000
INFO @ Thu, 13 Apr 2017 18:41:53: 90000000
INFO @ Thu, 13 Apr 2017 18:41:54: 91000000
INFO @ Thu, 13 Apr 2017 18:41:56: 92000000
INFO @ Thu, 13 Apr 2017 18:41:57: 93000000
INFO @ Thu, 13 Apr 2017 18:41:58: 94000000
INFO @ Thu, 13 Apr 2017 18:42:00: 95000000
INFO @ Thu, 13 Apr 2017 18:42:01: 96000000
INFO @ Thu, 13 Apr 2017 18:42:03: 97000000
INFO @ Thu, 13 Apr 2017 18:42:04: 98000000
INFO @ Thu, 13 Apr 2017 18:42:05: 99000000
INFO @ Thu, 13 Apr 2017 18:42:07: #1 tag size is determined as 101 bps
INFO @ Thu, 13 Apr 2017 18:42:07: #1 tag size = 101
INFO @ Thu, 13 Apr 2017 18:42:07: #1 total tags in treatment: 54689796
INFO @ Thu, 13 Apr 2017 18:42:07: #1 total tags in control: 99493262
INFO @ Thu, 13 Apr 2017 18:42:07: #1 finished!
INFO @ Thu, 13 Apr 2017 18:42:07: #2 Build Peak Model...
INFO @ Thu, 13 Apr 2017 18:42:07: #2 Skipped...
INFO @ Thu, 13 Apr 2017 18:42:07: #2 Use 215 as fragment length
INFO @ Thu, 13 Apr 2017 18:42:07: #3 Call peaks...
INFO @ Thu, 13 Apr 2017 18:42:07: #3 Call broad peaks with given level1 -log10pvalue cutoff and level2: 2.00000, 1.00000...
INFO @ Thu, 13 Apr 2017 18:42:07: #3 Pre-compute pvalue-qvalue table...
INFO @ Thu, 13 Apr 2017 18:47:49: #3 Call peaks for each chromosome...
INFO @ Thu, 13 Apr 2017 18:52:07: #4 Write output xls file... /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/rep1/ENCFF558PHY.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign_peaks.xls
INFO @ Thu, 13 Apr 2017 18:52:08: #4 Write broad peak in broadPeak format file... /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/rep1/ENCFF558PHY.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign_peaks.broadPeak
INFO @ Thu, 13 Apr 2017 18:52:09: #4 Write broad peak in bed12/gappedPeak format file... /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/rep1/ENCFF558PHY.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign_peaks.gappedPeak
INFO @ Thu, 13 Apr 2017 18:52:10: Done!
|
Num | 5 | ID | task.callpeak_macs2.macs2_g_b_rep1_pr1.line_155.id_14 | Name | macs2 g/b rep1-pr1 | Thread | thread_Root | PID | 46981 | | OK | true | Exit Code | 0 | Retries | | State | FINISHED | Dep. | OK | Cpus | | Mem | | | Start | 2017-04-13 18:37:57 | End | 2017-04-13 18:55:09 | Elapsed | 00:17:12 | Timeout | 00:00:-1 | Wall Timeout | 100 days | | Input files | /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/align/pseudo_reps/rep1/pr1/ENCFF558PHY.pr1.tagAlign.gz
/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/InputData/ss_50M_2000_GRCh38.nodup.tagAlign.gz
| Output files | /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep1/pr1/ENCFF558PHY.pr1.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.gappedPeak.gz
/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep1/pr1/ENCFF558PHY.pr1.tagAlign_x_ss_50M_2000_GRCh38.nodup.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/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/align/pseudo_reps/rep1/pr1/ENCFF558PHY.pr1.tagAlign.gz -c /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/InputData/ss_50M_2000_GRCh38.nodup.tagAlign.gz -f BED -n /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep1/pr1/ENCFF558PHY.pr1.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign -g hs -p 0.01 --broad --nomodel --shift 0 --extsize 215 --keep-dup all
# SYS command. line 166
sort -k 8gr,8gr "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep1/pr1/ENCFF558PHY.pr1.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign"_peaks.broadPeak | awk 'BEGIN{OFS="\t"}{$4="Peak_"NR ; print $0}' | gzip -nc > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep1/pr1/ENCFF558PHY.pr1.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.broadPeak.gz
# SYS command. line 167
sort -k 14gr,14gr "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep1/pr1/ENCFF558PHY.pr1.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign"_peaks.gappedPeak | awk 'BEGIN{OFS="\t"}{$4="Peak_"NR ; print $0}' | gzip -nc > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep1/pr1/ENCFF558PHY.pr1.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.gappedPeak.gz
# SYS command. line 170
rm -f "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep1/pr1/ENCFF558PHY.pr1.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign"_peaks.xls \
"/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep1/pr1/ENCFF558PHY.pr1.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign"_peaks.broadPeak \
"/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep1/pr1/ENCFF558PHY.pr1.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign"_peaks.gappedPeak \
"/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep1/pr1/ENCFF558PHY.pr1.tagAlign_x_ss_50M_2000_GRCh38.nodup.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--------------------
14457 (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/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/align/pseudo_reps/rep1/pr1/ENCFF558PHY.pr1.tagAlign.gz -c /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/InputData/ss_50M_2000_GRCh38.nodup.tagAlign.gz -f BED -n /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep1/pr1/ENCFF558PHY.pr1.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign -g hs -p 0.01 --broad --nomodel --shift 0 --extsize 215 --keep-dup all
# ARGUMENTS LIST:
# name = /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep1/pr1/ENCFF558PHY.pr1.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign
# format = BED
# ChIP-seq file = ['/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/align/pseudo_reps/rep1/pr1/ENCFF558PHY.pr1.tagAlign.gz']
# control file = ['/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/InputData/ss_50M_2000_GRCh38.nodup.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:19: 1000000
INFO @ Thu, 13 Apr 2017 18:38:21: 2000000
INFO @ Thu, 13 Apr 2017 18:38:23: 3000000
INFO @ Thu, 13 Apr 2017 18:38:25: 4000000
INFO @ Thu, 13 Apr 2017 18:38:27: 5000000
INFO @ Thu, 13 Apr 2017 18:38:29: 6000000
INFO @ Thu, 13 Apr 2017 18:38:31: 7000000
INFO @ Thu, 13 Apr 2017 18:38:33: 8000000
INFO @ Thu, 13 Apr 2017 18:38:35: 9000000
INFO @ Thu, 13 Apr 2017 18:38:37: 10000000
INFO @ Thu, 13 Apr 2017 18:38:39: 11000000
INFO @ Thu, 13 Apr 2017 18:38:41: 12000000
INFO @ Thu, 13 Apr 2017 18:38:43: 13000000
INFO @ Thu, 13 Apr 2017 18:38:45: 14000000
INFO @ Thu, 13 Apr 2017 18:38:47: 15000000
INFO @ Thu, 13 Apr 2017 18:38:49: 16000000
INFO @ Thu, 13 Apr 2017 18:38:51: 17000000
INFO @ Thu, 13 Apr 2017 18:38:53: 18000000
INFO @ Thu, 13 Apr 2017 18:38:55: 19000000
INFO @ Thu, 13 Apr 2017 18:38:57: 20000000
INFO @ Thu, 13 Apr 2017 18:38:59: 21000000
INFO @ Thu, 13 Apr 2017 18:39:01: 22000000
INFO @ Thu, 13 Apr 2017 18:39:03: 23000000
INFO @ Thu, 13 Apr 2017 18:39:05: 24000000
INFO @ Thu, 13 Apr 2017 18:39:07: 25000000
INFO @ Thu, 13 Apr 2017 18:39:09: 26000000
INFO @ Thu, 13 Apr 2017 18:39:11: 27000000
INFO @ Thu, 13 Apr 2017 18:39:14: #1.2 read input tags...
INFO @ Thu, 13 Apr 2017 18:39:16: 1000000
INFO @ Thu, 13 Apr 2017 18:39:18: 2000000
INFO @ Thu, 13 Apr 2017 18:39:20: 3000000
INFO @ Thu, 13 Apr 2017 18:39:22: 4000000
INFO @ Thu, 13 Apr 2017 18:39:24: 5000000
INFO @ Thu, 13 Apr 2017 18:39:26: 6000000
INFO @ Thu, 13 Apr 2017 18:39:28: 7000000
INFO @ Thu, 13 Apr 2017 18:39:29: 8000000
INFO @ Thu, 13 Apr 2017 18:39:31: 9000000
INFO @ Thu, 13 Apr 2017 18:39:33: 10000000
INFO @ Thu, 13 Apr 2017 18:39:35: 11000000
INFO @ Thu, 13 Apr 2017 18:39:37: 12000000
INFO @ Thu, 13 Apr 2017 18:39:39: 13000000
INFO @ Thu, 13 Apr 2017 18:39:41: 14000000
INFO @ Thu, 13 Apr 2017 18:39:43: 15000000
INFO @ Thu, 13 Apr 2017 18:39:45: 16000000
INFO @ Thu, 13 Apr 2017 18:39:46: 17000000
INFO @ Thu, 13 Apr 2017 18:39:48: 18000000
INFO @ Thu, 13 Apr 2017 18:39:50: 19000000
INFO @ Thu, 13 Apr 2017 18:39:52: 20000000
INFO @ Thu, 13 Apr 2017 18:39:54: 21000000
INFO @ Thu, 13 Apr 2017 18:39:56: 22000000
INFO @ Thu, 13 Apr 2017 18:39:58: 23000000
INFO @ Thu, 13 Apr 2017 18:40:00: 24000000
INFO @ Thu, 13 Apr 2017 18:40:02: 25000000
INFO @ Thu, 13 Apr 2017 18:40:03: 26000000
INFO @ Thu, 13 Apr 2017 18:40:05: 27000000
INFO @ Thu, 13 Apr 2017 18:40:07: 28000000
INFO @ Thu, 13 Apr 2017 18:40:09: 29000000
INFO @ Thu, 13 Apr 2017 18:40:11: 30000000
INFO @ Thu, 13 Apr 2017 18:40:13: 31000000
INFO @ Thu, 13 Apr 2017 18:40:15: 32000000
INFO @ Thu, 13 Apr 2017 18:40:17: 33000000
INFO @ Thu, 13 Apr 2017 18:40:19: 34000000
INFO @ Thu, 13 Apr 2017 18:40:20: 35000000
INFO @ Thu, 13 Apr 2017 18:40:22: 36000000
INFO @ Thu, 13 Apr 2017 18:40:24: 37000000
INFO @ Thu, 13 Apr 2017 18:40:26: 38000000
INFO @ Thu, 13 Apr 2017 18:40:28: 39000000
INFO @ Thu, 13 Apr 2017 18:40:30: 40000000
INFO @ Thu, 13 Apr 2017 18:40:32: 41000000
INFO @ Thu, 13 Apr 2017 18:40:34: 42000000
INFO @ Thu, 13 Apr 2017 18:40:35: 43000000
INFO @ Thu, 13 Apr 2017 18:40:37: 44000000
INFO @ Thu, 13 Apr 2017 18:40:39: 45000000
INFO @ Thu, 13 Apr 2017 18:40:41: 46000000
INFO @ Thu, 13 Apr 2017 18:40:43: 47000000
INFO @ Thu, 13 Apr 2017 18:40:45: 48000000
INFO @ Thu, 13 Apr 2017 18:40:47: 49000000
INFO @ Thu, 13 Apr 2017 18:40:48: 50000000
INFO @ Thu, 13 Apr 2017 18:40:50: 51000000
INFO @ Thu, 13 Apr 2017 18:40:52: 52000000
INFO @ Thu, 13 Apr 2017 18:40:54: 53000000
INFO @ Thu, 13 Apr 2017 18:40:56: 54000000
INFO @ Thu, 13 Apr 2017 18:40:58: 55000000
INFO @ Thu, 13 Apr 2017 18:41:00: 56000000
INFO @ Thu, 13 Apr 2017 18:41:02: 57000000
INFO @ Thu, 13 Apr 2017 18:41:04: 58000000
INFO @ Thu, 13 Apr 2017 18:41:06: 59000000
INFO @ Thu, 13 Apr 2017 18:41:08: 60000000
INFO @ Thu, 13 Apr 2017 18:41:09: 61000000
INFO @ Thu, 13 Apr 2017 18:41:11: 62000000
INFO @ Thu, 13 Apr 2017 18:41:13: 63000000
INFO @ Thu, 13 Apr 2017 18:41:15: 64000000
INFO @ Thu, 13 Apr 2017 18:41:17: 65000000
INFO @ Thu, 13 Apr 2017 18:41:19: 66000000
INFO @ Thu, 13 Apr 2017 18:41:21: 67000000
INFO @ Thu, 13 Apr 2017 18:41:23: 68000000
INFO @ Thu, 13 Apr 2017 18:41:25: 69000000
INFO @ Thu, 13 Apr 2017 18:41:27: 70000000
INFO @ Thu, 13 Apr 2017 18:41:29: 71000000
INFO @ Thu, 13 Apr 2017 18:41:30: 72000000
INFO @ Thu, 13 Apr 2017 18:41:32: 73000000
INFO @ Thu, 13 Apr 2017 18:41:34: 74000000
INFO @ Thu, 13 Apr 2017 18:41:36: 75000000
INFO @ Thu, 13 Apr 2017 18:41:38: 76000000
INFO @ Thu, 13 Apr 2017 18:41:40: 77000000
INFO @ Thu, 13 Apr 2017 18:41:42: 78000000
INFO @ Thu, 13 Apr 2017 18:41:44: 79000000
INFO @ Thu, 13 Apr 2017 18:41:46: 80000000
INFO @ Thu, 13 Apr 2017 18:41:47: 81000000
INFO @ Thu, 13 Apr 2017 18:41:49: 82000000
INFO @ Thu, 13 Apr 2017 18:41:51: 83000000
INFO @ Thu, 13 Apr 2017 18:41:53: 84000000
INFO @ Thu, 13 Apr 2017 18:41:55: 85000000
INFO @ Thu, 13 Apr 2017 18:41:57: 86000000
INFO @ Thu, 13 Apr 2017 18:41:58: 87000000
INFO @ Thu, 13 Apr 2017 18:42:00: 88000000
INFO @ Thu, 13 Apr 2017 18:42:02: 89000000
INFO @ Thu, 13 Apr 2017 18:42:04: 90000000
INFO @ Thu, 13 Apr 2017 18:42:06: 91000000
INFO @ Thu, 13 Apr 2017 18:42:08: 92000000
INFO @ Thu, 13 Apr 2017 18:42:09: 93000000
INFO @ Thu, 13 Apr 2017 18:42:11: 94000000
INFO @ Thu, 13 Apr 2017 18:42:13: 95000000
INFO @ Thu, 13 Apr 2017 18:42:15: 96000000
INFO @ Thu, 13 Apr 2017 18:42:17: 97000000
INFO @ Thu, 13 Apr 2017 18:42:19: 98000000
INFO @ Thu, 13 Apr 2017 18:42:21: 99000000
INFO @ Thu, 13 Apr 2017 18:42:23: #1 tag size is determined as 101 bps
INFO @ Thu, 13 Apr 2017 18:42:23: #1 tag size = 101
INFO @ Thu, 13 Apr 2017 18:42:23: #1 total tags in treatment: 27344898
INFO @ Thu, 13 Apr 2017 18:42:23: #1 total tags in control: 99493262
INFO @ Thu, 13 Apr 2017 18:42:23: #1 finished!
INFO @ Thu, 13 Apr 2017 18:42:23: #2 Build Peak Model...
INFO @ Thu, 13 Apr 2017 18:42:23: #2 Skipped...
INFO @ Thu, 13 Apr 2017 18:42:23: #2 Use 215 as fragment length
INFO @ Thu, 13 Apr 2017 18:42:23: #3 Call peaks...
INFO @ Thu, 13 Apr 2017 18:42:23: #3 Call broad peaks with given level1 -log10pvalue cutoff and level2: 2.00000, 1.00000...
INFO @ Thu, 13 Apr 2017 18:42:23: #3 Pre-compute pvalue-qvalue table...
INFO @ Thu, 13 Apr 2017 18:49:10: #3 Call peaks for each chromosome...
INFO @ Thu, 13 Apr 2017 18:54:50: #4 Write output xls file... /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep1/pr1/ENCFF558PHY.pr1.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign_peaks.xls
INFO @ Thu, 13 Apr 2017 18:54:53: #4 Write broad peak in broadPeak format file... /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep1/pr1/ENCFF558PHY.pr1.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign_peaks.broadPeak
INFO @ Thu, 13 Apr 2017 18:54:58: #4 Write broad peak in bed12/gappedPeak format file... /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep1/pr1/ENCFF558PHY.pr1.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign_peaks.gappedPeak
INFO @ Thu, 13 Apr 2017 18:55:00: Done!
|
Num | 6 | ID | task.callpeak_macs2.macs2_g_b_rep1_pr2.line_155.id_15 | Name | macs2 g/b rep1-pr2 | Thread | thread_Root | PID | 47016 | | OK | true | Exit Code | 0 | Retries | | State | FINISHED | Dep. | OK | Cpus | | Mem | | | Start | 2017-04-13 18:51:33 | End | 2017-04-13 19:03:56 | Elapsed | 00:12:23 | Timeout | 00:00:-1 | Wall Timeout | 100 days | | Input files | /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/align/pseudo_reps/rep1/pr2/ENCFF558PHY.pr2.tagAlign.gz
/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/InputData/ss_50M_2000_GRCh38.nodup.tagAlign.gz
| Output files | /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep1/pr2/ENCFF558PHY.pr2.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.gappedPeak.gz
/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep1/pr2/ENCFF558PHY.pr2.tagAlign_x_ss_50M_2000_GRCh38.nodup.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/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/align/pseudo_reps/rep1/pr2/ENCFF558PHY.pr2.tagAlign.gz -c /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/InputData/ss_50M_2000_GRCh38.nodup.tagAlign.gz -f BED -n /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep1/pr2/ENCFF558PHY.pr2.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign -g hs -p 0.01 --broad --nomodel --shift 0 --extsize 215 --keep-dup all
# SYS command. line 166
sort -k 8gr,8gr "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep1/pr2/ENCFF558PHY.pr2.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign"_peaks.broadPeak | awk 'BEGIN{OFS="\t"}{$4="Peak_"NR ; print $0}' | gzip -nc > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep1/pr2/ENCFF558PHY.pr2.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.broadPeak.gz
# SYS command. line 167
sort -k 14gr,14gr "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep1/pr2/ENCFF558PHY.pr2.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign"_peaks.gappedPeak | awk 'BEGIN{OFS="\t"}{$4="Peak_"NR ; print $0}' | gzip -nc > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep1/pr2/ENCFF558PHY.pr2.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.gappedPeak.gz
# SYS command. line 170
rm -f "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep1/pr2/ENCFF558PHY.pr2.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign"_peaks.xls \
"/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep1/pr2/ENCFF558PHY.pr2.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign"_peaks.broadPeak \
"/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep1/pr2/ENCFF558PHY.pr2.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign"_peaks.gappedPeak \
"/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep1/pr2/ENCFF558PHY.pr2.tagAlign_x_ss_50M_2000_GRCh38.nodup.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--------------------
24782 (process ID) old priority 0, new priority 19
--------------------Stderr--------------------
INFO @ Thu, 13 Apr 2017 18:51:46:
# Command line: callpeak -t /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/align/pseudo_reps/rep1/pr2/ENCFF558PHY.pr2.tagAlign.gz -c /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/InputData/ss_50M_2000_GRCh38.nodup.tagAlign.gz -f BED -n /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep1/pr2/ENCFF558PHY.pr2.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign -g hs -p 0.01 --broad --nomodel --shift 0 --extsize 215 --keep-dup all
# ARGUMENTS LIST:
# name = /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep1/pr2/ENCFF558PHY.pr2.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign
# format = BED
# ChIP-seq file = ['/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/align/pseudo_reps/rep1/pr2/ENCFF558PHY.pr2.tagAlign.gz']
# control file = ['/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/InputData/ss_50M_2000_GRCh38.nodup.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:51:46: #1 read tag files...
INFO @ Thu, 13 Apr 2017 18:51:46: #1 read treatment tags...
INFO @ Thu, 13 Apr 2017 18:51:48: 1000000
INFO @ Thu, 13 Apr 2017 18:51:50: 2000000
INFO @ Thu, 13 Apr 2017 18:51:52: 3000000
INFO @ Thu, 13 Apr 2017 18:51:54: 4000000
INFO @ Thu, 13 Apr 2017 18:51:56: 5000000
INFO @ Thu, 13 Apr 2017 18:51:58: 6000000
INFO @ Thu, 13 Apr 2017 18:51:59: 7000000
INFO @ Thu, 13 Apr 2017 18:52:01: 8000000
INFO @ Thu, 13 Apr 2017 18:52:02: 9000000
INFO @ Thu, 13 Apr 2017 18:52:04: 10000000
INFO @ Thu, 13 Apr 2017 18:52:05: 11000000
INFO @ Thu, 13 Apr 2017 18:52:07: 12000000
INFO @ Thu, 13 Apr 2017 18:52:08: 13000000
INFO @ Thu, 13 Apr 2017 18:52:10: 14000000
INFO @ Thu, 13 Apr 2017 18:52:12: 15000000
INFO @ Thu, 13 Apr 2017 18:52:13: 16000000
INFO @ Thu, 13 Apr 2017 18:52:15: 17000000
INFO @ Thu, 13 Apr 2017 18:52:16: 18000000
INFO @ Thu, 13 Apr 2017 18:52:18: 19000000
INFO @ Thu, 13 Apr 2017 18:52:19: 20000000
INFO @ Thu, 13 Apr 2017 18:52:21: 21000000
INFO @ Thu, 13 Apr 2017 18:52:22: 22000000
INFO @ Thu, 13 Apr 2017 18:52:24: 23000000
INFO @ Thu, 13 Apr 2017 18:52:26: 24000000
INFO @ Thu, 13 Apr 2017 18:52:28: 25000000
INFO @ Thu, 13 Apr 2017 18:52:29: 26000000
INFO @ Thu, 13 Apr 2017 18:52:31: 27000000
INFO @ Thu, 13 Apr 2017 18:52:33: #1.2 read input tags...
INFO @ Thu, 13 Apr 2017 18:52:34: 1000000
INFO @ Thu, 13 Apr 2017 18:52:36: 2000000
INFO @ Thu, 13 Apr 2017 18:52:37: 3000000
INFO @ Thu, 13 Apr 2017 18:52:39: 4000000
INFO @ Thu, 13 Apr 2017 18:52:40: 5000000
INFO @ Thu, 13 Apr 2017 18:52:41: 6000000
INFO @ Thu, 13 Apr 2017 18:52:43: 7000000
INFO @ Thu, 13 Apr 2017 18:52:44: 8000000
INFO @ Thu, 13 Apr 2017 18:52:46: 9000000
INFO @ Thu, 13 Apr 2017 18:52:47: 10000000
INFO @ Thu, 13 Apr 2017 18:52:48: 11000000
INFO @ Thu, 13 Apr 2017 18:52:50: 12000000
INFO @ Thu, 13 Apr 2017 18:52:51: 13000000
INFO @ Thu, 13 Apr 2017 18:52:53: 14000000
INFO @ Thu, 13 Apr 2017 18:52:54: 15000000
INFO @ Thu, 13 Apr 2017 18:52:55: 16000000
INFO @ Thu, 13 Apr 2017 18:52:57: 17000000
INFO @ Thu, 13 Apr 2017 18:52:58: 18000000
INFO @ Thu, 13 Apr 2017 18:52:59: 19000000
INFO @ Thu, 13 Apr 2017 18:53:01: 20000000
INFO @ Thu, 13 Apr 2017 18:53:02: 21000000
INFO @ Thu, 13 Apr 2017 18:53:04: 22000000
INFO @ Thu, 13 Apr 2017 18:53:05: 23000000
INFO @ Thu, 13 Apr 2017 18:53:06: 24000000
INFO @ Thu, 13 Apr 2017 18:53:08: 25000000
INFO @ Thu, 13 Apr 2017 18:53:09: 26000000
INFO @ Thu, 13 Apr 2017 18:53:11: 27000000
INFO @ Thu, 13 Apr 2017 18:53:12: 28000000
INFO @ Thu, 13 Apr 2017 18:53:13: 29000000
INFO @ Thu, 13 Apr 2017 18:53:15: 30000000
INFO @ Thu, 13 Apr 2017 18:53:16: 31000000
INFO @ Thu, 13 Apr 2017 18:53:18: 32000000
INFO @ Thu, 13 Apr 2017 18:53:19: 33000000
INFO @ Thu, 13 Apr 2017 18:53:20: 34000000
INFO @ Thu, 13 Apr 2017 18:53:22: 35000000
INFO @ Thu, 13 Apr 2017 18:53:23: 36000000
INFO @ Thu, 13 Apr 2017 18:53:24: 37000000
INFO @ Thu, 13 Apr 2017 18:53:26: 38000000
INFO @ Thu, 13 Apr 2017 18:53:27: 39000000
INFO @ Thu, 13 Apr 2017 18:53:29: 40000000
INFO @ Thu, 13 Apr 2017 18:53:30: 41000000
INFO @ Thu, 13 Apr 2017 18:53:31: 42000000
INFO @ Thu, 13 Apr 2017 18:53:33: 43000000
INFO @ Thu, 13 Apr 2017 18:53:34: 44000000
INFO @ Thu, 13 Apr 2017 18:53:35: 45000000
INFO @ Thu, 13 Apr 2017 18:53:37: 46000000
INFO @ Thu, 13 Apr 2017 18:53:38: 47000000
INFO @ Thu, 13 Apr 2017 18:53:40: 48000000
INFO @ Thu, 13 Apr 2017 18:53:41: 49000000
INFO @ Thu, 13 Apr 2017 18:53:42: 50000000
INFO @ Thu, 13 Apr 2017 18:53:44: 51000000
INFO @ Thu, 13 Apr 2017 18:53:45: 52000000
INFO @ Thu, 13 Apr 2017 18:53:46: 53000000
INFO @ Thu, 13 Apr 2017 18:53:48: 54000000
INFO @ Thu, 13 Apr 2017 18:53:49: 55000000
INFO @ Thu, 13 Apr 2017 18:53:51: 56000000
INFO @ Thu, 13 Apr 2017 18:53:52: 57000000
INFO @ Thu, 13 Apr 2017 18:53:53: 58000000
INFO @ Thu, 13 Apr 2017 18:53:55: 59000000
INFO @ Thu, 13 Apr 2017 18:53:56: 60000000
INFO @ Thu, 13 Apr 2017 18:53:58: 61000000
INFO @ Thu, 13 Apr 2017 18:53:59: 62000000
INFO @ Thu, 13 Apr 2017 18:54:00: 63000000
INFO @ Thu, 13 Apr 2017 18:54:02: 64000000
INFO @ Thu, 13 Apr 2017 18:54:03: 65000000
INFO @ Thu, 13 Apr 2017 18:54:05: 66000000
INFO @ Thu, 13 Apr 2017 18:54:06: 67000000
INFO @ Thu, 13 Apr 2017 18:54:07: 68000000
INFO @ Thu, 13 Apr 2017 18:54:09: 69000000
INFO @ Thu, 13 Apr 2017 18:54:10: 70000000
INFO @ Thu, 13 Apr 2017 18:54:11: 71000000
INFO @ Thu, 13 Apr 2017 18:54:13: 72000000
INFO @ Thu, 13 Apr 2017 18:54:14: 73000000
INFO @ Thu, 13 Apr 2017 18:54:16: 74000000
INFO @ Thu, 13 Apr 2017 18:54:17: 75000000
INFO @ Thu, 13 Apr 2017 18:54:18: 76000000
INFO @ Thu, 13 Apr 2017 18:54:20: 77000000
INFO @ Thu, 13 Apr 2017 18:54:21: 78000000
INFO @ Thu, 13 Apr 2017 18:54:23: 79000000
INFO @ Thu, 13 Apr 2017 18:54:24: 80000000
INFO @ Thu, 13 Apr 2017 18:54:25: 81000000
INFO @ Thu, 13 Apr 2017 18:54:27: 82000000
INFO @ Thu, 13 Apr 2017 18:54:28: 83000000
INFO @ Thu, 13 Apr 2017 18:54:29: 84000000
INFO @ Thu, 13 Apr 2017 18:54:31: 85000000
INFO @ Thu, 13 Apr 2017 18:54:32: 86000000
INFO @ Thu, 13 Apr 2017 18:54:34: 87000000
INFO @ Thu, 13 Apr 2017 18:54:35: 88000000
INFO @ Thu, 13 Apr 2017 18:54:36: 89000000
INFO @ Thu, 13 Apr 2017 18:54:38: 90000000
INFO @ Thu, 13 Apr 2017 18:54:39: 91000000
INFO @ Thu, 13 Apr 2017 18:54:41: 92000000
INFO @ Thu, 13 Apr 2017 18:54:42: 93000000
INFO @ Thu, 13 Apr 2017 18:54:43: 94000000
INFO @ Thu, 13 Apr 2017 18:54:45: 95000000
INFO @ Thu, 13 Apr 2017 18:54:46: 96000000
INFO @ Thu, 13 Apr 2017 18:54:47: 97000000
INFO @ Thu, 13 Apr 2017 18:54:49: 98000000
INFO @ Thu, 13 Apr 2017 18:54:50: 99000000
INFO @ Thu, 13 Apr 2017 18:54:52: #1 tag size is determined as 101 bps
INFO @ Thu, 13 Apr 2017 18:54:52: #1 tag size = 101
INFO @ Thu, 13 Apr 2017 18:54:52: #1 total tags in treatment: 27344898
INFO @ Thu, 13 Apr 2017 18:54:52: #1 total tags in control: 99493262
INFO @ Thu, 13 Apr 2017 18:54:52: #1 finished!
INFO @ Thu, 13 Apr 2017 18:54:52: #2 Build Peak Model...
INFO @ Thu, 13 Apr 2017 18:54:52: #2 Skipped...
INFO @ Thu, 13 Apr 2017 18:54:52: #2 Use 215 as fragment length
INFO @ Thu, 13 Apr 2017 18:54:52: #3 Call peaks...
INFO @ Thu, 13 Apr 2017 18:54:52: #3 Call broad peaks with given level1 -log10pvalue cutoff and level2: 2.00000, 1.00000...
INFO @ Thu, 13 Apr 2017 18:54:52: #3 Pre-compute pvalue-qvalue table...
INFO @ Thu, 13 Apr 2017 18:59:48: #3 Call peaks for each chromosome...
INFO @ Thu, 13 Apr 2017 19:03:43: #4 Write output xls file... /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep1/pr2/ENCFF558PHY.pr2.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign_peaks.xls
INFO @ Thu, 13 Apr 2017 19:03:45: #4 Write broad peak in broadPeak format file... /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep1/pr2/ENCFF558PHY.pr2.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign_peaks.broadPeak
INFO @ Thu, 13 Apr 2017 19:03:47: #4 Write broad peak in bed12/gappedPeak format file... /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep1/pr2/ENCFF558PHY.pr2.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign_peaks.gappedPeak
INFO @ Thu, 13 Apr 2017 19:03:48: Done!
|
Num | 7 | ID | task.callpeak_macs2.macs2_g_b_rep2.line_155.id_16 | Name | macs2 g/b rep2 | Thread | thread_Root | PID | 47017 | | OK | true | Exit Code | 0 | Retries | | State | FINISHED | Dep. | OK | Cpus | | Mem | | | Start | 2017-04-13 18:51:41 | End | 2017-04-13 19:14:01 | Elapsed | 00:22:20 | Timeout | 00:00:-1 | Wall Timeout | 100 days | | Input files | /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/align/rep2/ENCFF645WYX.tagAlign.gz
/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/InputData/ss_50M_2000_GRCh38.nodup.tagAlign.gz
| Output files | /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/rep2/ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.gappedPeak.gz
/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/rep2/ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.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/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/align/rep2/ENCFF645WYX.tagAlign.gz -c /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/InputData/ss_50M_2000_GRCh38.nodup.tagAlign.gz -f BED -n /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/rep2/ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign -g hs -p 0.01 --broad --nomodel --shift 0 --extsize 210 --keep-dup all
# SYS command. line 166
sort -k 8gr,8gr "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/rep2/ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign"_peaks.broadPeak | awk 'BEGIN{OFS="\t"}{$4="Peak_"NR ; print $0}' | gzip -nc > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/rep2/ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.broadPeak.gz
# SYS command. line 167
sort -k 14gr,14gr "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/rep2/ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign"_peaks.gappedPeak | awk 'BEGIN{OFS="\t"}{$4="Peak_"NR ; print $0}' | gzip -nc > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/rep2/ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.gappedPeak.gz
# SYS command. line 170
rm -f "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/rep2/ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign"_peaks.xls \
"/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/rep2/ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign"_peaks.broadPeak \
"/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/rep2/ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign"_peaks.gappedPeak \
"/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/rep2/ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.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--------------------
19594 (process ID) old priority 0, new priority 19
--------------------Stderr--------------------
INFO @ Thu, 13 Apr 2017 18:52:02:
# Command line: callpeak -t /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/align/rep2/ENCFF645WYX.tagAlign.gz -c /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/InputData/ss_50M_2000_GRCh38.nodup.tagAlign.gz -f BED -n /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/rep2/ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign -g hs -p 0.01 --broad --nomodel --shift 0 --extsize 210 --keep-dup all
# ARGUMENTS LIST:
# name = /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/rep2/ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign
# format = BED
# ChIP-seq file = ['/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/align/rep2/ENCFF645WYX.tagAlign.gz']
# control file = ['/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/InputData/ss_50M_2000_GRCh38.nodup.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:52:02: #1 read tag files...
INFO @ Thu, 13 Apr 2017 18:52:02: #1 read treatment tags...
INFO @ Thu, 13 Apr 2017 18:52:04: 1000000
INFO @ Thu, 13 Apr 2017 18:52:06: 2000000
INFO @ Thu, 13 Apr 2017 18:52:08: 3000000
INFO @ Thu, 13 Apr 2017 18:52:10: 4000000
INFO @ Thu, 13 Apr 2017 18:52:12: 5000000
INFO @ Thu, 13 Apr 2017 18:52:14: 6000000
INFO @ Thu, 13 Apr 2017 18:52:16: 7000000
INFO @ Thu, 13 Apr 2017 18:52:18: 8000000
INFO @ Thu, 13 Apr 2017 18:52:20: 9000000
INFO @ Thu, 13 Apr 2017 18:52:22: 10000000
INFO @ Thu, 13 Apr 2017 18:52:24: 11000000
INFO @ Thu, 13 Apr 2017 18:52:26: 12000000
INFO @ Thu, 13 Apr 2017 18:52:28: 13000000
INFO @ Thu, 13 Apr 2017 18:52:30: 14000000
INFO @ Thu, 13 Apr 2017 18:52:32: 15000000
INFO @ Thu, 13 Apr 2017 18:52:34: 16000000
INFO @ Thu, 13 Apr 2017 18:52:36: 17000000
INFO @ Thu, 13 Apr 2017 18:52:38: 18000000
INFO @ Thu, 13 Apr 2017 18:52:40: 19000000
INFO @ Thu, 13 Apr 2017 18:52:42: 20000000
INFO @ Thu, 13 Apr 2017 18:52:44: 21000000
INFO @ Thu, 13 Apr 2017 18:52:46: 22000000
INFO @ Thu, 13 Apr 2017 18:52:48: 23000000
INFO @ Thu, 13 Apr 2017 18:52:50: 24000000
INFO @ Thu, 13 Apr 2017 18:52:52: 25000000
INFO @ Thu, 13 Apr 2017 18:52:54: 26000000
INFO @ Thu, 13 Apr 2017 18:52:56: 27000000
INFO @ Thu, 13 Apr 2017 18:52:58: 28000000
INFO @ Thu, 13 Apr 2017 18:53:00: 29000000
INFO @ Thu, 13 Apr 2017 18:53:02: 30000000
INFO @ Thu, 13 Apr 2017 18:53:04: 31000000
INFO @ Thu, 13 Apr 2017 18:53:06: 32000000
INFO @ Thu, 13 Apr 2017 18:53:08: 33000000
INFO @ Thu, 13 Apr 2017 18:53:10: 34000000
INFO @ Thu, 13 Apr 2017 18:53:12: 35000000
INFO @ Thu, 13 Apr 2017 18:53:14: 36000000
INFO @ Thu, 13 Apr 2017 18:53:16: 37000000
INFO @ Thu, 13 Apr 2017 18:53:18: 38000000
INFO @ Thu, 13 Apr 2017 18:53:20: 39000000
INFO @ Thu, 13 Apr 2017 18:53:22: 40000000
INFO @ Thu, 13 Apr 2017 18:53:24: 41000000
INFO @ Thu, 13 Apr 2017 18:53:26: 42000000
INFO @ Thu, 13 Apr 2017 18:53:28: 43000000
INFO @ Thu, 13 Apr 2017 18:53:30: 44000000
INFO @ Thu, 13 Apr 2017 18:53:32: 45000000
INFO @ Thu, 13 Apr 2017 18:53:34: 46000000
INFO @ Thu, 13 Apr 2017 18:53:36: 47000000
INFO @ Thu, 13 Apr 2017 18:53:38: 48000000
INFO @ Thu, 13 Apr 2017 18:53:40: 49000000
INFO @ Thu, 13 Apr 2017 18:53:42: 50000000
INFO @ Thu, 13 Apr 2017 18:53:43: 51000000
INFO @ Thu, 13 Apr 2017 18:53:45: 52000000
INFO @ Thu, 13 Apr 2017 18:53:47: 53000000
INFO @ Thu, 13 Apr 2017 18:53:53: #1.2 read input tags...
INFO @ Thu, 13 Apr 2017 18:53:55: 1000000
INFO @ Thu, 13 Apr 2017 18:53:57: 2000000
INFO @ Thu, 13 Apr 2017 18:53:59: 3000000
INFO @ Thu, 13 Apr 2017 18:54:01: 4000000
INFO @ Thu, 13 Apr 2017 18:54:03: 5000000
INFO @ Thu, 13 Apr 2017 18:54:05: 6000000
INFO @ Thu, 13 Apr 2017 18:54:06: 7000000
INFO @ Thu, 13 Apr 2017 18:54:08: 8000000
INFO @ Thu, 13 Apr 2017 18:54:10: 9000000
INFO @ Thu, 13 Apr 2017 18:54:12: 10000000
INFO @ Thu, 13 Apr 2017 18:54:14: 11000000
INFO @ Thu, 13 Apr 2017 18:54:15: 12000000
INFO @ Thu, 13 Apr 2017 18:54:17: 13000000
INFO @ Thu, 13 Apr 2017 18:54:19: 14000000
INFO @ Thu, 13 Apr 2017 18:54:21: 15000000
INFO @ Thu, 13 Apr 2017 18:54:23: 16000000
INFO @ Thu, 13 Apr 2017 18:54:25: 17000000
INFO @ Thu, 13 Apr 2017 18:54:27: 18000000
INFO @ Thu, 13 Apr 2017 18:54:29: 19000000
INFO @ Thu, 13 Apr 2017 18:54:31: 20000000
INFO @ Thu, 13 Apr 2017 18:54:33: 21000000
INFO @ Thu, 13 Apr 2017 18:54:34: 22000000
INFO @ Thu, 13 Apr 2017 18:54:36: 23000000
INFO @ Thu, 13 Apr 2017 18:54:38: 24000000
INFO @ Thu, 13 Apr 2017 18:54:40: 25000000
INFO @ Thu, 13 Apr 2017 18:54:42: 26000000
INFO @ Thu, 13 Apr 2017 18:54:44: 27000000
INFO @ Thu, 13 Apr 2017 18:54:45: 28000000
INFO @ Thu, 13 Apr 2017 18:54:47: 29000000
INFO @ Thu, 13 Apr 2017 18:54:49: 30000000
INFO @ Thu, 13 Apr 2017 18:54:51: 31000000
INFO @ Thu, 13 Apr 2017 18:54:53: 32000000
INFO @ Thu, 13 Apr 2017 18:54:55: 33000000
INFO @ Thu, 13 Apr 2017 18:54:56: 34000000
INFO @ Thu, 13 Apr 2017 18:54:58: 35000000
INFO @ Thu, 13 Apr 2017 18:55:00: 36000000
INFO @ Thu, 13 Apr 2017 18:55:02: 37000000
INFO @ Thu, 13 Apr 2017 18:55:04: 38000000
INFO @ Thu, 13 Apr 2017 18:55:06: 39000000
INFO @ Thu, 13 Apr 2017 18:55:07: 40000000
INFO @ Thu, 13 Apr 2017 18:55:09: 41000000
INFO @ Thu, 13 Apr 2017 18:55:11: 42000000
INFO @ Thu, 13 Apr 2017 18:55:13: 43000000
INFO @ Thu, 13 Apr 2017 18:55:15: 44000000
INFO @ Thu, 13 Apr 2017 18:55:17: 45000000
INFO @ Thu, 13 Apr 2017 18:55:19: 46000000
INFO @ Thu, 13 Apr 2017 18:55:20: 47000000
INFO @ Thu, 13 Apr 2017 18:55:22: 48000000
INFO @ Thu, 13 Apr 2017 18:55:24: 49000000
INFO @ Thu, 13 Apr 2017 18:55:26: 50000000
INFO @ Thu, 13 Apr 2017 18:55:28: 51000000
INFO @ Thu, 13 Apr 2017 18:55:30: 52000000
INFO @ Thu, 13 Apr 2017 18:55:32: 53000000
INFO @ Thu, 13 Apr 2017 18:55:33: 54000000
INFO @ Thu, 13 Apr 2017 18:55:35: 55000000
INFO @ Thu, 13 Apr 2017 18:55:37: 56000000
INFO @ Thu, 13 Apr 2017 18:55:39: 57000000
INFO @ Thu, 13 Apr 2017 18:55:41: 58000000
INFO @ Thu, 13 Apr 2017 18:55:43: 59000000
INFO @ Thu, 13 Apr 2017 18:55:45: 60000000
INFO @ Thu, 13 Apr 2017 18:55:46: 61000000
INFO @ Thu, 13 Apr 2017 18:55:48: 62000000
INFO @ Thu, 13 Apr 2017 18:55:50: 63000000
INFO @ Thu, 13 Apr 2017 18:55:52: 64000000
INFO @ Thu, 13 Apr 2017 18:55:54: 65000000
INFO @ Thu, 13 Apr 2017 18:55:56: 66000000
INFO @ Thu, 13 Apr 2017 18:55:57: 67000000
INFO @ Thu, 13 Apr 2017 18:55:59: 68000000
INFO @ Thu, 13 Apr 2017 18:56:01: 69000000
INFO @ Thu, 13 Apr 2017 18:56:03: 70000000
INFO @ Thu, 13 Apr 2017 18:56:05: 71000000
INFO @ Thu, 13 Apr 2017 18:56:07: 72000000
INFO @ Thu, 13 Apr 2017 18:56:08: 73000000
INFO @ Thu, 13 Apr 2017 18:56:10: 74000000
INFO @ Thu, 13 Apr 2017 18:56:12: 75000000
INFO @ Thu, 13 Apr 2017 18:56:14: 76000000
INFO @ Thu, 13 Apr 2017 18:56:16: 77000000
INFO @ Thu, 13 Apr 2017 18:56:18: 78000000
INFO @ Thu, 13 Apr 2017 18:56:19: 79000000
INFO @ Thu, 13 Apr 2017 18:56:21: 80000000
INFO @ Thu, 13 Apr 2017 18:56:23: 81000000
INFO @ Thu, 13 Apr 2017 18:56:25: 82000000
INFO @ Thu, 13 Apr 2017 18:56:27: 83000000
INFO @ Thu, 13 Apr 2017 18:56:29: 84000000
INFO @ Thu, 13 Apr 2017 18:56:30: 85000000
INFO @ Thu, 13 Apr 2017 18:56:32: 86000000
INFO @ Thu, 13 Apr 2017 18:56:34: 87000000
INFO @ Thu, 13 Apr 2017 18:56:36: 88000000
INFO @ Thu, 13 Apr 2017 18:56:38: 89000000
INFO @ Thu, 13 Apr 2017 18:56:40: 90000000
INFO @ Thu, 13 Apr 2017 18:56:41: 91000000
INFO @ Thu, 13 Apr 2017 18:56:43: 92000000
INFO @ Thu, 13 Apr 2017 18:56:45: 93000000
INFO @ Thu, 13 Apr 2017 18:56:47: 94000000
INFO @ Thu, 13 Apr 2017 18:56:49: 95000000
INFO @ Thu, 13 Apr 2017 18:56:51: 96000000
INFO @ Thu, 13 Apr 2017 18:56:52: 97000000
INFO @ Thu, 13 Apr 2017 18:56:54: 98000000
INFO @ Thu, 13 Apr 2017 18:56:56: 99000000
INFO @ Thu, 13 Apr 2017 18:56:58: #1 tag size is determined as 101 bps
INFO @ Thu, 13 Apr 2017 18:56:58: #1 tag size = 101
INFO @ Thu, 13 Apr 2017 18:56:58: #1 total tags in treatment: 53733056
INFO @ Thu, 13 Apr 2017 18:56:58: #1 total tags in control: 99493262
INFO @ Thu, 13 Apr 2017 18:56:58: #1 finished!
INFO @ Thu, 13 Apr 2017 18:56:58: #2 Build Peak Model...
INFO @ Thu, 13 Apr 2017 18:56:58: #2 Skipped...
INFO @ Thu, 13 Apr 2017 18:56:58: #2 Use 210 as fragment length
INFO @ Thu, 13 Apr 2017 18:56:58: #3 Call peaks...
INFO @ Thu, 13 Apr 2017 18:56:58: #3 Call broad peaks with given level1 -log10pvalue cutoff and level2: 2.00000, 1.00000...
INFO @ Thu, 13 Apr 2017 18:56:58: #3 Pre-compute pvalue-qvalue table...
INFO @ Thu, 13 Apr 2017 19:06:38: #3 Call peaks for each chromosome...
INFO @ Thu, 13 Apr 2017 19:13:51: #4 Write output xls file... /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/rep2/ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign_peaks.xls
INFO @ Thu, 13 Apr 2017 19:13:53: #4 Write broad peak in broadPeak format file... /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/rep2/ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign_peaks.broadPeak
INFO @ Thu, 13 Apr 2017 19:13:54: #4 Write broad peak in bed12/gappedPeak format file... /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/rep2/ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign_peaks.gappedPeak
INFO @ Thu, 13 Apr 2017 19:13:56: Done!
|
Num | 8 | ID | task.callpeak_macs2.macs2_g_b_rep2_pr1.line_155.id_17 | Name | macs2 g/b rep2-pr1 | Thread | thread_Root | PID | 47018 | | OK | true | Exit Code | 0 | Retries | | State | FINISHED | Dep. | OK | Cpus | | Mem | | | Start | 2017-04-13 18:52:19 | End | 2017-04-13 19:04:34 | Elapsed | 00:12:15 | Timeout | 00:00:-1 | Wall Timeout | 100 days | | Input files | /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/align/pseudo_reps/rep2/pr1/ENCFF645WYX.pr1.tagAlign.gz
/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/InputData/ss_50M_2000_GRCh38.nodup.tagAlign.gz
| Output files | /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep2/pr1/ENCFF645WYX.pr1.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.gappedPeak.gz
/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep2/pr1/ENCFF645WYX.pr1.tagAlign_x_ss_50M_2000_GRCh38.nodup.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/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/align/pseudo_reps/rep2/pr1/ENCFF645WYX.pr1.tagAlign.gz -c /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/InputData/ss_50M_2000_GRCh38.nodup.tagAlign.gz -f BED -n /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep2/pr1/ENCFF645WYX.pr1.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign -g hs -p 0.01 --broad --nomodel --shift 0 --extsize 210 --keep-dup all
# SYS command. line 166
sort -k 8gr,8gr "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep2/pr1/ENCFF645WYX.pr1.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign"_peaks.broadPeak | awk 'BEGIN{OFS="\t"}{$4="Peak_"NR ; print $0}' | gzip -nc > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep2/pr1/ENCFF645WYX.pr1.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.broadPeak.gz
# SYS command. line 167
sort -k 14gr,14gr "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep2/pr1/ENCFF645WYX.pr1.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign"_peaks.gappedPeak | awk 'BEGIN{OFS="\t"}{$4="Peak_"NR ; print $0}' | gzip -nc > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep2/pr1/ENCFF645WYX.pr1.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.gappedPeak.gz
# SYS command. line 170
rm -f "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep2/pr1/ENCFF645WYX.pr1.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign"_peaks.xls \
"/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep2/pr1/ENCFF645WYX.pr1.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign"_peaks.broadPeak \
"/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep2/pr1/ENCFF645WYX.pr1.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign"_peaks.gappedPeak \
"/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep2/pr1/ENCFF645WYX.pr1.tagAlign_x_ss_50M_2000_GRCh38.nodup.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--------------------
26604 (process ID) old priority 0, new priority 19
--------------------Stderr--------------------
INFO @ Thu, 13 Apr 2017 18:52:30:
# Command line: callpeak -t /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/align/pseudo_reps/rep2/pr1/ENCFF645WYX.pr1.tagAlign.gz -c /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/InputData/ss_50M_2000_GRCh38.nodup.tagAlign.gz -f BED -n /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep2/pr1/ENCFF645WYX.pr1.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign -g hs -p 0.01 --broad --nomodel --shift 0 --extsize 210 --keep-dup all
# ARGUMENTS LIST:
# name = /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep2/pr1/ENCFF645WYX.pr1.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign
# format = BED
# ChIP-seq file = ['/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/align/pseudo_reps/rep2/pr1/ENCFF645WYX.pr1.tagAlign.gz']
# control file = ['/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/InputData/ss_50M_2000_GRCh38.nodup.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:52:30: #1 read tag files...
INFO @ Thu, 13 Apr 2017 18:52:30: #1 read treatment tags...
INFO @ Thu, 13 Apr 2017 18:52:32: 1000000
INFO @ Thu, 13 Apr 2017 18:52:34: 2000000
INFO @ Thu, 13 Apr 2017 18:52:35: 3000000
INFO @ Thu, 13 Apr 2017 18:52:37: 4000000
INFO @ Thu, 13 Apr 2017 18:52:38: 5000000
INFO @ Thu, 13 Apr 2017 18:52:40: 6000000
INFO @ Thu, 13 Apr 2017 18:52:41: 7000000
INFO @ Thu, 13 Apr 2017 18:52:43: 8000000
INFO @ Thu, 13 Apr 2017 18:52:44: 9000000
INFO @ Thu, 13 Apr 2017 18:52:46: 10000000
INFO @ Thu, 13 Apr 2017 18:52:47: 11000000
INFO @ Thu, 13 Apr 2017 18:52:49: 12000000
INFO @ Thu, 13 Apr 2017 18:52:50: 13000000
INFO @ Thu, 13 Apr 2017 18:52:52: 14000000
INFO @ Thu, 13 Apr 2017 18:52:53: 15000000
INFO @ Thu, 13 Apr 2017 18:52:55: 16000000
INFO @ Thu, 13 Apr 2017 18:52:56: 17000000
INFO @ Thu, 13 Apr 2017 18:52:58: 18000000
INFO @ Thu, 13 Apr 2017 18:52:59: 19000000
INFO @ Thu, 13 Apr 2017 18:53:01: 20000000
INFO @ Thu, 13 Apr 2017 18:53:02: 21000000
INFO @ Thu, 13 Apr 2017 18:53:04: 22000000
INFO @ Thu, 13 Apr 2017 18:53:05: 23000000
INFO @ Thu, 13 Apr 2017 18:53:07: 24000000
INFO @ Thu, 13 Apr 2017 18:53:08: 25000000
INFO @ Thu, 13 Apr 2017 18:53:10: 26000000
INFO @ Thu, 13 Apr 2017 18:53:13: #1.2 read input tags...
INFO @ Thu, 13 Apr 2017 18:53:14: 1000000
INFO @ Thu, 13 Apr 2017 18:53:16: 2000000
INFO @ Thu, 13 Apr 2017 18:53:17: 3000000
INFO @ Thu, 13 Apr 2017 18:53:19: 4000000
INFO @ Thu, 13 Apr 2017 18:53:20: 5000000
INFO @ Thu, 13 Apr 2017 18:53:21: 6000000
INFO @ Thu, 13 Apr 2017 18:53:23: 7000000
INFO @ Thu, 13 Apr 2017 18:53:24: 8000000
INFO @ Thu, 13 Apr 2017 18:53:26: 9000000
INFO @ Thu, 13 Apr 2017 18:53:27: 10000000
INFO @ Thu, 13 Apr 2017 18:53:29: 11000000
INFO @ Thu, 13 Apr 2017 18:53:30: 12000000
INFO @ Thu, 13 Apr 2017 18:53:31: 13000000
INFO @ Thu, 13 Apr 2017 18:53:33: 14000000
INFO @ Thu, 13 Apr 2017 18:53:34: 15000000
INFO @ Thu, 13 Apr 2017 18:53:36: 16000000
INFO @ Thu, 13 Apr 2017 18:53:37: 17000000
INFO @ Thu, 13 Apr 2017 18:53:38: 18000000
INFO @ Thu, 13 Apr 2017 18:53:40: 19000000
INFO @ Thu, 13 Apr 2017 18:53:41: 20000000
INFO @ Thu, 13 Apr 2017 18:53:42: 21000000
INFO @ Thu, 13 Apr 2017 18:53:44: 22000000
INFO @ Thu, 13 Apr 2017 18:53:45: 23000000
INFO @ Thu, 13 Apr 2017 18:53:47: 24000000
INFO @ Thu, 13 Apr 2017 18:53:48: 25000000
INFO @ Thu, 13 Apr 2017 18:53:49: 26000000
INFO @ Thu, 13 Apr 2017 18:53:51: 27000000
INFO @ Thu, 13 Apr 2017 18:53:52: 28000000
INFO @ Thu, 13 Apr 2017 18:53:54: 29000000
INFO @ Thu, 13 Apr 2017 18:53:55: 30000000
INFO @ Thu, 13 Apr 2017 18:53:56: 31000000
INFO @ Thu, 13 Apr 2017 18:53:58: 32000000
INFO @ Thu, 13 Apr 2017 18:53:59: 33000000
INFO @ Thu, 13 Apr 2017 18:54:01: 34000000
INFO @ Thu, 13 Apr 2017 18:54:02: 35000000
INFO @ Thu, 13 Apr 2017 18:54:03: 36000000
INFO @ Thu, 13 Apr 2017 18:54:05: 37000000
INFO @ Thu, 13 Apr 2017 18:54:06: 38000000
INFO @ Thu, 13 Apr 2017 18:54:08: 39000000
INFO @ Thu, 13 Apr 2017 18:54:09: 40000000
INFO @ Thu, 13 Apr 2017 18:54:11: 41000000
INFO @ Thu, 13 Apr 2017 18:54:12: 42000000
INFO @ Thu, 13 Apr 2017 18:54:14: 43000000
INFO @ Thu, 13 Apr 2017 18:54:15: 44000000
INFO @ Thu, 13 Apr 2017 18:54:17: 45000000
INFO @ Thu, 13 Apr 2017 18:54:18: 46000000
INFO @ Thu, 13 Apr 2017 18:54:20: 47000000
INFO @ Thu, 13 Apr 2017 18:54:21: 48000000
INFO @ Thu, 13 Apr 2017 18:54:23: 49000000
INFO @ Thu, 13 Apr 2017 18:54:24: 50000000
INFO @ Thu, 13 Apr 2017 18:54:25: 51000000
INFO @ Thu, 13 Apr 2017 18:54:27: 52000000
INFO @ Thu, 13 Apr 2017 18:54:28: 53000000
INFO @ Thu, 13 Apr 2017 18:54:30: 54000000
INFO @ Thu, 13 Apr 2017 18:54:31: 55000000
INFO @ Thu, 13 Apr 2017 18:54:32: 56000000
INFO @ Thu, 13 Apr 2017 18:54:34: 57000000
INFO @ Thu, 13 Apr 2017 18:54:35: 58000000
INFO @ Thu, 13 Apr 2017 18:54:37: 59000000
INFO @ Thu, 13 Apr 2017 18:54:38: 60000000
INFO @ Thu, 13 Apr 2017 18:54:39: 61000000
INFO @ Thu, 13 Apr 2017 18:54:41: 62000000
INFO @ Thu, 13 Apr 2017 18:54:42: 63000000
INFO @ Thu, 13 Apr 2017 18:54:44: 64000000
INFO @ Thu, 13 Apr 2017 18:54:45: 65000000
INFO @ Thu, 13 Apr 2017 18:54:46: 66000000
INFO @ Thu, 13 Apr 2017 18:54:48: 67000000
INFO @ Thu, 13 Apr 2017 18:54:49: 68000000
INFO @ Thu, 13 Apr 2017 18:54:51: 69000000
INFO @ Thu, 13 Apr 2017 18:54:52: 70000000
INFO @ Thu, 13 Apr 2017 18:54:53: 71000000
INFO @ Thu, 13 Apr 2017 18:54:55: 72000000
INFO @ Thu, 13 Apr 2017 18:54:57: 73000000
INFO @ Thu, 13 Apr 2017 18:54:58: 74000000
INFO @ Thu, 13 Apr 2017 18:54:59: 75000000
INFO @ Thu, 13 Apr 2017 18:55:01: 76000000
INFO @ Thu, 13 Apr 2017 18:55:02: 77000000
INFO @ Thu, 13 Apr 2017 18:55:04: 78000000
INFO @ Thu, 13 Apr 2017 18:55:05: 79000000
INFO @ Thu, 13 Apr 2017 18:55:07: 80000000
INFO @ Thu, 13 Apr 2017 18:55:08: 81000000
INFO @ Thu, 13 Apr 2017 18:55:10: 82000000
INFO @ Thu, 13 Apr 2017 18:55:11: 83000000
INFO @ Thu, 13 Apr 2017 18:55:13: 84000000
INFO @ Thu, 13 Apr 2017 18:55:14: 85000000
INFO @ Thu, 13 Apr 2017 18:55:16: 86000000
INFO @ Thu, 13 Apr 2017 18:55:17: 87000000
INFO @ Thu, 13 Apr 2017 18:55:19: 88000000
INFO @ Thu, 13 Apr 2017 18:55:21: 89000000
INFO @ Thu, 13 Apr 2017 18:55:22: 90000000
INFO @ Thu, 13 Apr 2017 18:55:23: 91000000
INFO @ Thu, 13 Apr 2017 18:55:25: 92000000
INFO @ Thu, 13 Apr 2017 18:55:26: 93000000
INFO @ Thu, 13 Apr 2017 18:55:28: 94000000
INFO @ Thu, 13 Apr 2017 18:55:29: 95000000
INFO @ Thu, 13 Apr 2017 18:55:31: 96000000
INFO @ Thu, 13 Apr 2017 18:55:32: 97000000
INFO @ Thu, 13 Apr 2017 18:55:34: 98000000
INFO @ Thu, 13 Apr 2017 18:55:35: 99000000
INFO @ Thu, 13 Apr 2017 18:55:37: #1 tag size is determined as 101 bps
INFO @ Thu, 13 Apr 2017 18:55:37: #1 tag size = 101
INFO @ Thu, 13 Apr 2017 18:55:37: #1 total tags in treatment: 26866528
INFO @ Thu, 13 Apr 2017 18:55:37: #1 total tags in control: 99493262
INFO @ Thu, 13 Apr 2017 18:55:37: #1 finished!
INFO @ Thu, 13 Apr 2017 18:55:37: #2 Build Peak Model...
INFO @ Thu, 13 Apr 2017 18:55:37: #2 Skipped...
INFO @ Thu, 13 Apr 2017 18:55:37: #2 Use 210 as fragment length
INFO @ Thu, 13 Apr 2017 18:55:37: #3 Call peaks...
INFO @ Thu, 13 Apr 2017 18:55:37: #3 Call broad peaks with given level1 -log10pvalue cutoff and level2: 2.00000, 1.00000...
INFO @ Thu, 13 Apr 2017 18:55:37: #3 Pre-compute pvalue-qvalue table...
INFO @ Thu, 13 Apr 2017 19:00:28: #3 Call peaks for each chromosome...
INFO @ Thu, 13 Apr 2017 19:04:22: #4 Write output xls file... /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep2/pr1/ENCFF645WYX.pr1.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign_peaks.xls
INFO @ Thu, 13 Apr 2017 19:04:24: #4 Write broad peak in broadPeak format file... /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep2/pr1/ENCFF645WYX.pr1.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign_peaks.broadPeak
INFO @ Thu, 13 Apr 2017 19:04:26: #4 Write broad peak in bed12/gappedPeak format file... /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep2/pr1/ENCFF645WYX.pr1.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign_peaks.gappedPeak
INFO @ Thu, 13 Apr 2017 19:04:28: Done!
|
Num | 9 | ID | task.callpeak_macs2.macs2_g_b_rep2_pr2.line_155.id_18 | Name | macs2 g/b rep2-pr2 | Thread | thread_Root | PID | 47024 | | OK | true | Exit Code | 0 | Retries | | State | FINISHED | Dep. | OK | Cpus | | Mem | | | Start | 2017-04-13 18:55:12 | End | 2017-04-13 19:07:33 | Elapsed | 00:12:20 | Timeout | 00:00:-1 | Wall Timeout | 100 days | | Input files | /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/align/pseudo_reps/rep2/pr2/ENCFF645WYX.pr2.tagAlign.gz
/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/InputData/ss_50M_2000_GRCh38.nodup.tagAlign.gz
| Output files | /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep2/pr2/ENCFF645WYX.pr2.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.gappedPeak.gz
/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep2/pr2/ENCFF645WYX.pr2.tagAlign_x_ss_50M_2000_GRCh38.nodup.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/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/align/pseudo_reps/rep2/pr2/ENCFF645WYX.pr2.tagAlign.gz -c /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/InputData/ss_50M_2000_GRCh38.nodup.tagAlign.gz -f BED -n /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep2/pr2/ENCFF645WYX.pr2.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign -g hs -p 0.01 --broad --nomodel --shift 0 --extsize 210 --keep-dup all
# SYS command. line 166
sort -k 8gr,8gr "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep2/pr2/ENCFF645WYX.pr2.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign"_peaks.broadPeak | awk 'BEGIN{OFS="\t"}{$4="Peak_"NR ; print $0}' | gzip -nc > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep2/pr2/ENCFF645WYX.pr2.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.broadPeak.gz
# SYS command. line 167
sort -k 14gr,14gr "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep2/pr2/ENCFF645WYX.pr2.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign"_peaks.gappedPeak | awk 'BEGIN{OFS="\t"}{$4="Peak_"NR ; print $0}' | gzip -nc > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep2/pr2/ENCFF645WYX.pr2.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.gappedPeak.gz
# SYS command. line 170
rm -f "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep2/pr2/ENCFF645WYX.pr2.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign"_peaks.xls \
"/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep2/pr2/ENCFF645WYX.pr2.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign"_peaks.broadPeak \
"/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep2/pr2/ENCFF645WYX.pr2.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign"_peaks.gappedPeak \
"/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep2/pr2/ENCFF645WYX.pr2.tagAlign_x_ss_50M_2000_GRCh38.nodup.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--------------------
34102 (process ID) old priority 0, new priority 19
--------------------Stderr--------------------
INFO @ Thu, 13 Apr 2017 18:55:32:
# Command line: callpeak -t /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/align/pseudo_reps/rep2/pr2/ENCFF645WYX.pr2.tagAlign.gz -c /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/InputData/ss_50M_2000_GRCh38.nodup.tagAlign.gz -f BED -n /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep2/pr2/ENCFF645WYX.pr2.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign -g hs -p 0.01 --broad --nomodel --shift 0 --extsize 210 --keep-dup all
# ARGUMENTS LIST:
# name = /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep2/pr2/ENCFF645WYX.pr2.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign
# format = BED
# ChIP-seq file = ['/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/align/pseudo_reps/rep2/pr2/ENCFF645WYX.pr2.tagAlign.gz']
# control file = ['/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/InputData/ss_50M_2000_GRCh38.nodup.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:55:32: #1 read tag files...
INFO @ Thu, 13 Apr 2017 18:55:32: #1 read treatment tags...
INFO @ Thu, 13 Apr 2017 18:55:34: 1000000
INFO @ Thu, 13 Apr 2017 18:55:35: 2000000
INFO @ Thu, 13 Apr 2017 18:55:37: 3000000
INFO @ Thu, 13 Apr 2017 18:55:38: 4000000
INFO @ Thu, 13 Apr 2017 18:55:40: 5000000
INFO @ Thu, 13 Apr 2017 18:55:41: 6000000
INFO @ Thu, 13 Apr 2017 18:55:43: 7000000
INFO @ Thu, 13 Apr 2017 18:55:44: 8000000
INFO @ Thu, 13 Apr 2017 18:55:46: 9000000
INFO @ Thu, 13 Apr 2017 18:55:47: 10000000
INFO @ Thu, 13 Apr 2017 18:55:49: 11000000
INFO @ Thu, 13 Apr 2017 18:55:50: 12000000
INFO @ Thu, 13 Apr 2017 18:55:52: 13000000
INFO @ Thu, 13 Apr 2017 18:55:53: 14000000
INFO @ Thu, 13 Apr 2017 18:55:55: 15000000
INFO @ Thu, 13 Apr 2017 18:55:56: 16000000
INFO @ Thu, 13 Apr 2017 18:55:58: 17000000
INFO @ Thu, 13 Apr 2017 18:56:00: 18000000
INFO @ Thu, 13 Apr 2017 18:56:01: 19000000
INFO @ Thu, 13 Apr 2017 18:56:03: 20000000
INFO @ Thu, 13 Apr 2017 18:56:04: 21000000
INFO @ Thu, 13 Apr 2017 18:56:06: 22000000
INFO @ Thu, 13 Apr 2017 18:56:07: 23000000
INFO @ Thu, 13 Apr 2017 18:56:09: 24000000
INFO @ Thu, 13 Apr 2017 18:56:10: 25000000
INFO @ Thu, 13 Apr 2017 18:56:12: 26000000
INFO @ Thu, 13 Apr 2017 18:56:15: #1.2 read input tags...
INFO @ Thu, 13 Apr 2017 18:56:16: 1000000
INFO @ Thu, 13 Apr 2017 18:56:18: 2000000
INFO @ Thu, 13 Apr 2017 18:56:19: 3000000
INFO @ Thu, 13 Apr 2017 18:56:20: 4000000
INFO @ Thu, 13 Apr 2017 18:56:22: 5000000
INFO @ Thu, 13 Apr 2017 18:56:23: 6000000
INFO @ Thu, 13 Apr 2017 18:56:25: 7000000
INFO @ Thu, 13 Apr 2017 18:56:26: 8000000
INFO @ Thu, 13 Apr 2017 18:56:27: 9000000
INFO @ Thu, 13 Apr 2017 18:56:29: 10000000
INFO @ Thu, 13 Apr 2017 18:56:30: 11000000
INFO @ Thu, 13 Apr 2017 18:56:31: 12000000
INFO @ Thu, 13 Apr 2017 18:56:33: 13000000
INFO @ Thu, 13 Apr 2017 18:56:34: 14000000
INFO @ Thu, 13 Apr 2017 18:56:36: 15000000
INFO @ Thu, 13 Apr 2017 18:56:37: 16000000
INFO @ Thu, 13 Apr 2017 18:56:38: 17000000
INFO @ Thu, 13 Apr 2017 18:56:40: 18000000
INFO @ Thu, 13 Apr 2017 18:56:41: 19000000
INFO @ Thu, 13 Apr 2017 18:56:43: 20000000
INFO @ Thu, 13 Apr 2017 18:56:44: 21000000
INFO @ Thu, 13 Apr 2017 18:56:45: 22000000
INFO @ Thu, 13 Apr 2017 18:56:47: 23000000
INFO @ Thu, 13 Apr 2017 18:56:48: 24000000
INFO @ Thu, 13 Apr 2017 18:56:50: 25000000
INFO @ Thu, 13 Apr 2017 18:56:51: 26000000
INFO @ Thu, 13 Apr 2017 18:56:52: 27000000
INFO @ Thu, 13 Apr 2017 18:56:54: 28000000
INFO @ Thu, 13 Apr 2017 18:56:55: 29000000
INFO @ Thu, 13 Apr 2017 18:56:56: 30000000
INFO @ Thu, 13 Apr 2017 18:56:58: 31000000
INFO @ Thu, 13 Apr 2017 18:56:59: 32000000
INFO @ Thu, 13 Apr 2017 18:57:01: 33000000
INFO @ Thu, 13 Apr 2017 18:57:02: 34000000
INFO @ Thu, 13 Apr 2017 18:57:03: 35000000
INFO @ Thu, 13 Apr 2017 18:57:05: 36000000
INFO @ Thu, 13 Apr 2017 18:57:06: 37000000
INFO @ Thu, 13 Apr 2017 18:57:07: 38000000
INFO @ Thu, 13 Apr 2017 18:57:09: 39000000
INFO @ Thu, 13 Apr 2017 18:57:10: 40000000
INFO @ Thu, 13 Apr 2017 18:57:12: 41000000
INFO @ Thu, 13 Apr 2017 18:57:13: 42000000
INFO @ Thu, 13 Apr 2017 18:57:14: 43000000
INFO @ Thu, 13 Apr 2017 18:57:16: 44000000
INFO @ Thu, 13 Apr 2017 18:57:17: 45000000
INFO @ Thu, 13 Apr 2017 18:57:19: 46000000
INFO @ Thu, 13 Apr 2017 18:57:20: 47000000
INFO @ Thu, 13 Apr 2017 18:57:21: 48000000
INFO @ Thu, 13 Apr 2017 18:57:23: 49000000
INFO @ Thu, 13 Apr 2017 18:57:24: 50000000
INFO @ Thu, 13 Apr 2017 18:57:25: 51000000
INFO @ Thu, 13 Apr 2017 18:57:27: 52000000
INFO @ Thu, 13 Apr 2017 18:57:28: 53000000
INFO @ Thu, 13 Apr 2017 18:57:30: 54000000
INFO @ Thu, 13 Apr 2017 18:57:31: 55000000
INFO @ Thu, 13 Apr 2017 18:57:32: 56000000
INFO @ Thu, 13 Apr 2017 18:57:34: 57000000
INFO @ Thu, 13 Apr 2017 18:57:35: 58000000
INFO @ Thu, 13 Apr 2017 18:57:37: 59000000
INFO @ Thu, 13 Apr 2017 18:57:38: 60000000
INFO @ Thu, 13 Apr 2017 18:57:39: 61000000
INFO @ Thu, 13 Apr 2017 18:57:41: 62000000
INFO @ Thu, 13 Apr 2017 18:57:42: 63000000
INFO @ Thu, 13 Apr 2017 18:57:43: 64000000
INFO @ Thu, 13 Apr 2017 18:57:45: 65000000
INFO @ Thu, 13 Apr 2017 18:57:46: 66000000
INFO @ Thu, 13 Apr 2017 18:57:48: 67000000
INFO @ Thu, 13 Apr 2017 18:57:49: 68000000
INFO @ Thu, 13 Apr 2017 18:57:50: 69000000
INFO @ Thu, 13 Apr 2017 18:57:52: 70000000
INFO @ Thu, 13 Apr 2017 18:57:53: 71000000
INFO @ Thu, 13 Apr 2017 18:57:55: 72000000
INFO @ Thu, 13 Apr 2017 18:57:56: 73000000
INFO @ Thu, 13 Apr 2017 18:57:57: 74000000
INFO @ Thu, 13 Apr 2017 18:57:59: 75000000
INFO @ Thu, 13 Apr 2017 18:58:00: 76000000
INFO @ Thu, 13 Apr 2017 18:58:01: 77000000
INFO @ Thu, 13 Apr 2017 18:58:03: 78000000
INFO @ Thu, 13 Apr 2017 18:58:04: 79000000
INFO @ Thu, 13 Apr 2017 18:58:06: 80000000
INFO @ Thu, 13 Apr 2017 18:58:07: 81000000
INFO @ Thu, 13 Apr 2017 18:58:08: 82000000
INFO @ Thu, 13 Apr 2017 18:58:10: 83000000
INFO @ Thu, 13 Apr 2017 18:58:11: 84000000
INFO @ Thu, 13 Apr 2017 18:58:13: 85000000
INFO @ Thu, 13 Apr 2017 18:58:14: 86000000
INFO @ Thu, 13 Apr 2017 18:58:15: 87000000
INFO @ Thu, 13 Apr 2017 18:58:17: 88000000
INFO @ Thu, 13 Apr 2017 18:58:18: 89000000
INFO @ Thu, 13 Apr 2017 18:58:19: 90000000
INFO @ Thu, 13 Apr 2017 18:58:21: 91000000
INFO @ Thu, 13 Apr 2017 18:58:22: 92000000
INFO @ Thu, 13 Apr 2017 18:58:24: 93000000
INFO @ Thu, 13 Apr 2017 18:58:25: 94000000
INFO @ Thu, 13 Apr 2017 18:58:26: 95000000
INFO @ Thu, 13 Apr 2017 18:58:28: 96000000
INFO @ Thu, 13 Apr 2017 18:58:29: 97000000
INFO @ Thu, 13 Apr 2017 18:58:30: 98000000
INFO @ Thu, 13 Apr 2017 18:58:32: 99000000
INFO @ Thu, 13 Apr 2017 18:58:33: #1 tag size is determined as 101 bps
INFO @ Thu, 13 Apr 2017 18:58:33: #1 tag size = 101
INFO @ Thu, 13 Apr 2017 18:58:33: #1 total tags in treatment: 26866528
INFO @ Thu, 13 Apr 2017 18:58:33: #1 total tags in control: 99493262
INFO @ Thu, 13 Apr 2017 18:58:33: #1 finished!
INFO @ Thu, 13 Apr 2017 18:58:33: #2 Build Peak Model...
INFO @ Thu, 13 Apr 2017 18:58:33: #2 Skipped...
INFO @ Thu, 13 Apr 2017 18:58:33: #2 Use 210 as fragment length
INFO @ Thu, 13 Apr 2017 18:58:33: #3 Call peaks...
INFO @ Thu, 13 Apr 2017 18:58:33: #3 Call broad peaks with given level1 -log10pvalue cutoff and level2: 2.00000, 1.00000...
INFO @ Thu, 13 Apr 2017 18:58:33: #3 Pre-compute pvalue-qvalue table...
INFO @ Thu, 13 Apr 2017 19:03:31: #3 Call peaks for each chromosome...
INFO @ Thu, 13 Apr 2017 19:07:20: #4 Write output xls file... /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep2/pr2/ENCFF645WYX.pr2.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign_peaks.xls
INFO @ Thu, 13 Apr 2017 19:07:23: #4 Write broad peak in broadPeak format file... /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep2/pr2/ENCFF645WYX.pr2.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign_peaks.broadPeak
INFO @ Thu, 13 Apr 2017 19:07:24: #4 Write broad peak in bed12/gappedPeak format file... /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pseudo_reps/rep2/pr2/ENCFF645WYX.pr2.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign_peaks.gappedPeak
INFO @ Thu, 13 Apr 2017 19:07:25: Done!
|
Num | 10 | ID | task.callpeak_naive_overlap.naive_overlap.line_143.id_19 | Name | naive_overlap | Thread | thread_Root | PID | 47031 | | OK | true | Exit Code | 0 | Retries | | State | FINISHED | Dep. | OK | Cpus | | Mem | | | Start | 2017-04-13 19:14:04 | End | 2017-04-13 19:15:15 | Elapsed | 00:01:10 | Timeout | 00:00:-1 | Wall Timeout | 100 days | | Input files | /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pooled_rep/ENCFF558PHY_ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.gappedPeak.gz
/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/rep1/ENCFF558PHY.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.gappedPeak.gz
/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/rep2/ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.gappedPeak.gz
/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pooled_pseudo_reps/ppr1/ENCFF558PHY.pr1_ENCFF645WYX.pr1.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.gappedPeak.gz
/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pooled_pseudo_reps/ppr2/ENCFF558PHY.pr2_ENCFF645WYX.pr2.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.gappedPeak.gz
| Output files | /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/overlap/ENCFF558PHY_ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.naive_overlap.filt.gappedPeak.gz
| Dependencies | | |
# SYS command. line 145
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 148
intersectBed -wo -a <(zcat -f /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pooled_rep/ENCFF558PHY_ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.gappedPeak.gz) -b <(zcat -f /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/rep1/ENCFF558PHY.tagAlign_x_ss_50M_2000_GRCh38.nodup.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/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/rep2/ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.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/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/overlap/ENCFF558PHY_ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.PooledInReps.gappedPeak.gz
# SYS command. line 151
intersectBed -wo -a <(zcat -f /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pooled_rep/ENCFF558PHY_ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.gappedPeak.gz) -b <(zcat -f /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pooled_pseudo_reps/ppr1/ENCFF558PHY.pr1_ENCFF645WYX.pr1.tagAlign_x_ss_50M_2000_GRCh38.nodup.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/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pooled_pseudo_reps/ppr2/ENCFF558PHY.pr2_ENCFF645WYX.pr2.tagAlign_x_ss_50M_2000_GRCh38.nodup.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/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/overlap/ENCFF558PHY_ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.PooledInPsRep1AndPsRep2.gappedPeak.gz
# SYS command. line 154
zcat /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/overlap/ENCFF558PHY_ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.PooledInReps.gappedPeak.gz /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/overlap/ENCFF558PHY_ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.PooledInPsRep1AndPsRep2.gappedPeak.gz | sort | uniq | gzip -nc > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/overlap/ENCFF558PHY_ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.naive_overlap.gappedPeak.gz
# SYS command. line 156
bedtools intersect -v -a <(zcat -f /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/overlap/ENCFF558PHY_ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.naive_overlap.gappedPeak.gz) -b <(zcat -f /mnt/data/bds_pipeline_genome_data/hg38_ENCODE/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/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/overlap/ENCFF558PHY_ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.naive_overlap.filt.gappedPeak.gz
# SYS command. line 158
rm -f /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/overlap/ENCFF558PHY_ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.PooledInReps.gappedPeak.gz /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/overlap/ENCFF558PHY_ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.PooledInPsRep1AndPsRep2.gappedPeak.gz
# SYS command. line 160
TASKTIME=$[$(date +%s)-${STARTTIME}]; if [ ${TASKTIME} -lt 60 ]; then echo "Waiting for $[60-${TASKTIME}] seconds."; sleep $[60-${TASKTIME}]; fi
--------------------Stdout--------------------
25324 (process ID) old priority 0, new priority 19
Waiting for 52 seconds.
|
Num | 11 | ID | task.callpeak_naive_overlap.naive_overlap.line_143.id_20 | Name | naive_overlap | Thread | thread_Root | PID | 47032 | | OK | true | Exit Code | 0 | Retries | | State | FINISHED | Dep. | OK | Cpus | | Mem | | | Start | 2017-04-13 19:14:06 | End | 2017-04-13 19:15:30 | Elapsed | 00:01:23 | Timeout | 00:00:-1 | Wall Timeout | 100 days | | Input files | /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pooled_rep/ENCFF558PHY_ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.broadPeak.gz
/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/rep1/ENCFF558PHY.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.broadPeak.gz
/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/rep2/ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.broadPeak.gz
/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pooled_pseudo_reps/ppr1/ENCFF558PHY.pr1_ENCFF645WYX.pr1.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.broadPeak.gz
/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pooled_pseudo_reps/ppr2/ENCFF558PHY.pr2_ENCFF645WYX.pr2.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.broadPeak.gz
| Output files | /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/overlap/ENCFF558PHY_ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.naive_overlap.filt.broadPeak.gz
| Dependencies | | |
# SYS command. line 145
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 148
intersectBed -wo -a <(zcat -f /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pooled_rep/ENCFF558PHY_ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.broadPeak.gz) -b <(zcat -f /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/rep1/ENCFF558PHY.tagAlign_x_ss_50M_2000_GRCh38.nodup.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/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/rep2/ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.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/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/overlap/ENCFF558PHY_ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.PooledInReps.broadPeak.gz
# SYS command. line 151
intersectBed -wo -a <(zcat -f /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pooled_rep/ENCFF558PHY_ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.broadPeak.gz) -b <(zcat -f /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pooled_pseudo_reps/ppr1/ENCFF558PHY.pr1_ENCFF645WYX.pr1.tagAlign_x_ss_50M_2000_GRCh38.nodup.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/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/pooled_pseudo_reps/ppr2/ENCFF558PHY.pr2_ENCFF645WYX.pr2.tagAlign_x_ss_50M_2000_GRCh38.nodup.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/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/overlap/ENCFF558PHY_ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.PooledInPsRep1AndPsRep2.broadPeak.gz
# SYS command. line 154
zcat /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/overlap/ENCFF558PHY_ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.PooledInReps.broadPeak.gz /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/overlap/ENCFF558PHY_ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.PooledInPsRep1AndPsRep2.broadPeak.gz | sort | uniq | gzip -nc > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/overlap/ENCFF558PHY_ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.naive_overlap.broadPeak.gz
# SYS command. line 156
bedtools intersect -v -a <(zcat -f /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/overlap/ENCFF558PHY_ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.naive_overlap.broadPeak.gz) -b <(zcat -f /mnt/data/bds_pipeline_genome_data/hg38_ENCODE/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/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/overlap/ENCFF558PHY_ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.naive_overlap.filt.broadPeak.gz
# SYS command. line 158
rm -f /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/overlap/ENCFF558PHY_ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.PooledInReps.broadPeak.gz /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/overlap/ENCFF558PHY_ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.PooledInPsRep1AndPsRep2.broadPeak.gz
# SYS command. line 160
TASKTIME=$[$(date +%s)-${STARTTIME}]; if [ ${TASKTIME} -lt 60 ]; then echo "Waiting for $[60-${TASKTIME}] seconds."; sleep $[60-${TASKTIME}]; fi
--------------------Stdout--------------------
26068 (process ID) old priority 0, new priority 19
Waiting for 54 seconds.
|
Num | 12 | ID | task.report.peak2hammock.line_361.id_23 | Name | peak2hammock | Thread | thread_Root | PID | 29607 | | OK | true | Exit Code | 0 | Retries | | State | FINISHED | Dep. | OK | Cpus | | Mem | | | Start | 2017-04-13 19:15:34 | End | 2017-04-13 19:15:45 | Elapsed | 00:00:10 | Timeout | 00:00:-1 | Wall Timeout | 100 days | | Input files | /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/overlap/ENCFF558PHY_ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.naive_overlap.filt.gappedPeak.gz
| Output files | /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/overlap/ENCFF558PHY_ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.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/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/overlap/ENCFF558PHY_ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.naive_overlap.filt.gappedPeak.gz | sed '/^\(chr\)/!d' | sort -k1,1V -k2,2n > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/overlap/ENCFF558PHY_ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.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/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/overlap/ENCFF558PHY_ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.naive_overlap.filt.gappedPeak.tmp /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/overlap/ENCFF558PHY_ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.naive_overlap.filt.gappedPeak.hammock
# SYS command. line 367
rm -f /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/peak/macs2/overlap/ENCFF558PHY_ENCFF645WYX.tagAlign_x_ss_50M_2000_GRCh38.nodup.tagAlign.naive_overlap.filt.gappedPeak.tmp
--------------------Stdout--------------------
29611 (process ID) old priority 0, new priority 19
|
Num | 13 | ID | task.graphviz.report.line_97.id_24 | Name | report | Thread | thread_Root | PID | 30132 | | OK | true | Exit Code | 0 | Retries | | State | FINISHED | Dep. | OK | Cpus | | Mem | | | Start | 2017-04-13 19:15:46 | End | 2017-04-13 19:15:54 | 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/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/report/ZIC2_MACS2_workflow.dot > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/report/ZIC2_MACS2_workflow.svg 2> /dev/null || echo "svg: /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/ZIC2_MACS2/out/report/ZIC2_MACS2_workflow.svg" # to suppress dot font error (exit code=1)
--------------------Stdout--------------------
30142 (process ID) old priority 0, new priority 19
|