Num | 1 | ID | task.callpeak_macs2.macs2_rep1.line_71.id_10 | Name | macs2 rep1 | Thread | thread_Root | PID | 22365 | | OK | true | Exit Code | 0 | Retries | | State | FINISHED | Dep. | OK | Cpus | | Mem | | | Start | 2017-01-03 14:25:23 | End | 2017-01-03 16:28:48 | Elapsed | 02:03:24 | Timeout | 00:00:-1 | Wall Timeout | 100 days | | Input files | /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/align/rep1/SRR1370904_1.nodup.tagAlign.gz
/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/HughesInputData/input15MReadsNSCLess1.05.50MSubsample.tagAlign.gz
| Output files | /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/peak/macs2/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.narrowPeak.gz
/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/peak/macs2/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.gappedPeak.gz
/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/peak/macs2/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.broadPeak.gz
/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/signal/macs2/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.fc.signal.bw
/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/signal/macs2/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.pval.signal.bw
| Dependencies | | |
# SYS command. line 73
if [[ -f $(which conda) && $(conda env list | grep aquas_chipseq | wc -l) != "0" ]]; then source 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); renice -n 10 $$
# SYS command. line 74
export LC_COLLATE=C
# SYS command. line 79
macs2 callpeak -t /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/align/rep1/SRR1370904_1.nodup.tagAlign.gz -c /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/HughesInputData/input15MReadsNSCLess1.05.50MSubsample.tagAlign.gz -f BED -n /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/peak/macs2/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign -g hs -p 0.01 --nomodel --shift 0 --extsize 170 --keep-dup all -B --SPMR
# SYS command. line 82
sort -k 8gr,8gr "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/peak/macs2/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign"_peaks.narrowPeak | awk 'BEGIN{OFS="\t"}{$4="Peak_"NR ; print $0}' | gzip -nc > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/peak/macs2/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.narrowPeak.gz
# SYS command. line 85
rm -f "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/peak/macs2/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign"_peaks.xls \
"/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/peak/macs2/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign"_peaks.narrowPeak \
"/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/peak/macs2/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign"_summits.bed
# SYS command. line 92
macs2 callpeak -t /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/align/rep1/SRR1370904_1.nodup.tagAlign.gz -c /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/HughesInputData/input15MReadsNSCLess1.05.50MSubsample.tagAlign.gz -f BED -n /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/peak/macs2/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign -g hs -p 0.01 --broad --nomodel --shift 0 --extsize 170 --keep-dup all
# SYS command. line 95
sort -k 8gr,8gr "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/peak/macs2/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign"_peaks.broadPeak | awk 'BEGIN{OFS="\t"}{$4="Peak_"NR ; print $0}' | gzip -nc > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/peak/macs2/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.broadPeak.gz
# SYS command. line 96
sort -k 14gr,14gr "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/peak/macs2/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign"_peaks.gappedPeak | awk 'BEGIN{OFS="\t"}{$4="Peak_"NR ; print $0}' | gzip -nc > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/peak/macs2/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.gappedPeak.gz
# SYS command. line 99
rm -f "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/peak/macs2/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign"_peaks.xls \
"/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/peak/macs2/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign"_peaks.broadPeak \
"/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/peak/macs2/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign"_peaks.gappedPeak \
"/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/peak/macs2/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign"_summits.bed
# SYS command. line 110
if [[ true == "false" ]]; then \
rm -f "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/peak/macs2/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign"_treat_pileup.bdg "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/peak/macs2/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign"_control_lambda.bdg; \
TASKTIME=$[$(date +%s)-${STARTTIME}]; if [ ${TASKTIME} -lt 60 ]; then echo "Waiting for $[60-${TASKTIME}] seconds."; sleep $[60-${TASKTIME}]; fi; \
exit; \
fi
# SYS command. line 121
macs2 bdgcmp -t "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/peak/macs2/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign"_treat_pileup.bdg -c "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/peak/macs2/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign"_control_lambda.bdg --outdir /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/peak/macs2/rep1 -o "SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign"_FE.bdg -m FE
# SYS command. line 124
slopBed -i "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/peak/macs2/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign"_FE.bdg -g /mnt/data/bds_pipeline_genome_data/hg38/hg38.chrom.sizes -b 0 | awk '{if ($3 != -1) print $0}' | bedClip stdin /mnt/data/bds_pipeline_genome_data/hg38/hg38.chrom.sizes /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/peak/macs2/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.fc.signal.bedgraph
# SYS command. line 125
rm -f "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/peak/macs2/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign"_FE.bdg
# SYS command. line 128
sort -k1,1 -k2,2n /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/peak/macs2/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.fc.signal.bedgraph > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/peak/macs2/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.fc.signal.srt.bedgraph
# SYS command. line 129
bedGraphToBigWig /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/peak/macs2/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.fc.signal.srt.bedgraph /mnt/data/bds_pipeline_genome_data/hg38/hg38.chrom.sizes /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/signal/macs2/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.fc.signal.bw
# SYS command. line 130
rm -f /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/peak/macs2/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.fc.signal.bedgraph /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/peak/macs2/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.fc.signal.srt.bedgraph
# SYS command. line 136
chipReads=$(zcat /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/align/rep1/SRR1370904_1.nodup.tagAlign.gz | wc -l | awk '{printf "%f", $1/1000000}')
# SYS command. line 138
controlReads=$(zcat /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/HughesInputData/input15MReadsNSCLess1.05.50MSubsample.tagAlign.gz | wc -l | awk '{printf "%f", $1/1000000}'); sval=$(echo "${chipReads} ${controlReads}" | awk '$1>$2{printf "%f",$2} $1<=$2{printf "%f",$1}')
# SYS command. line 140
macs2 bdgcmp -t "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/peak/macs2/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign"_treat_pileup.bdg -c "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/peak/macs2/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign"_control_lambda.bdg --outdir /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/peak/macs2/rep1 -o "SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign"_ppois.bdg -m ppois -S "${sval}"
# SYS command. line 143
slopBed -i "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/peak/macs2/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign"_ppois.bdg -g /mnt/data/bds_pipeline_genome_data/hg38/hg38.chrom.sizes -b 0 | awk '{if ($3 != -1) print $0}' | bedClip stdin /mnt/data/bds_pipeline_genome_data/hg38/hg38.chrom.sizes /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/peak/macs2/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.pval.signal.bedgraph
# SYS command. line 144
rm -rf "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/peak/macs2/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign"_ppois.bdg
# SYS command. line 147
sort -k1,1 -k2,2n /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/peak/macs2/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.pval.signal.bedgraph > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/peak/macs2/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.pval.signal.srt.bedgraph
# SYS command. line 148
bedGraphToBigWig /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/peak/macs2/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.pval.signal.srt.bedgraph /mnt/data/bds_pipeline_genome_data/hg38/hg38.chrom.sizes /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/signal/macs2/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.pval.signal.bw
# SYS command. line 149
rm -f /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/peak/macs2/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.pval.signal.bedgraph /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/peak/macs2/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign.pval.signal.srt.bedgraph
# SYS command. line 151
rm -f "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/peak/macs2/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign"_treat_pileup.bdg "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/peak/macs2/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign"_control_lambda.bdg
# SYS command. line 153
TASKTIME=$[$(date +%s)-${STARTTIME}]; if [ ${TASKTIME} -lt 60 ]; then echo "Waiting for $[60-${TASKTIME}] seconds."; sleep $[60-${TASKTIME}]; fi
--------------------Stdout--------------------
32470 (process ID) old priority 0, new priority 10
--------------------Stderr--------------------
INFO @ Tue, 03 Jan 2017 15:29:32:
# Command line: callpeak -t /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/align/rep1/SRR1370904_1.nodup.tagAlign.gz -c /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/HughesInputData/input15MReadsNSCLess1.05.50MSubsample.tagAlign.gz -f BED -n /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/peak/macs2/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign -g hs -p 0.01 --nomodel --shift 0 --extsize 170 --keep-dup all -B --SPMR
# ARGUMENTS LIST:
# name = /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/peak/macs2/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign
# format = BED
# ChIP-seq file = ['/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/align/rep1/SRR1370904_1.nodup.tagAlign.gz']
# control file = ['/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/HughesInputData/input15MReadsNSCLess1.05.50MSubsample.tagAlign.gz']
# effective genome size = 2.70e+09
# band width = 300
# model fold = [5, 50]
# pvalue cutoff = 1.00e-02
# 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 off
# MACS will save fragment pileup signal per million reads
INFO @ Tue, 03 Jan 2017 15:29:32: #1 read tag files...
INFO @ Tue, 03 Jan 2017 15:29:32: #1 read treatment tags...
INFO @ Tue, 03 Jan 2017 15:29:34: 1000000
INFO @ Tue, 03 Jan 2017 15:29:36: 2000000
INFO @ Tue, 03 Jan 2017 15:29:38: 3000000
INFO @ Tue, 03 Jan 2017 15:29:39: 4000000
INFO @ Tue, 03 Jan 2017 15:29:41: 5000000
INFO @ Tue, 03 Jan 2017 15:29:43: 6000000
INFO @ Tue, 03 Jan 2017 15:29:44: 7000000
INFO @ Tue, 03 Jan 2017 15:29:46: 8000000
INFO @ Tue, 03 Jan 2017 15:29:47: 9000000
INFO @ Tue, 03 Jan 2017 15:29:49: 10000000
INFO @ Tue, 03 Jan 2017 15:29:50: 11000000
INFO @ Tue, 03 Jan 2017 15:29:52: 12000000
INFO @ Tue, 03 Jan 2017 15:29:53: 13000000
INFO @ Tue, 03 Jan 2017 15:29:55: 14000000
INFO @ Tue, 03 Jan 2017 15:29:56: 15000000
INFO @ Tue, 03 Jan 2017 15:29:58: 16000000
INFO @ Tue, 03 Jan 2017 15:29:59: 17000000
INFO @ Tue, 03 Jan 2017 15:30:01: 18000000
INFO @ Tue, 03 Jan 2017 15:30:03: 19000000
INFO @ Tue, 03 Jan 2017 15:30:04: 20000000
INFO @ Tue, 03 Jan 2017 15:30:06: #1.2 read input tags...
INFO @ Tue, 03 Jan 2017 15:30:08: 1000000
INFO @ Tue, 03 Jan 2017 15:30:09: 2000000
INFO @ Tue, 03 Jan 2017 15:30:11: 3000000
INFO @ Tue, 03 Jan 2017 15:30:13: 4000000
INFO @ Tue, 03 Jan 2017 15:30:14: 5000000
INFO @ Tue, 03 Jan 2017 15:30:16: 6000000
INFO @ Tue, 03 Jan 2017 15:30:18: 7000000
INFO @ Tue, 03 Jan 2017 15:30:20: 8000000
INFO @ Tue, 03 Jan 2017 15:30:21: 9000000
INFO @ Tue, 03 Jan 2017 15:30:23: 10000000
INFO @ Tue, 03 Jan 2017 15:30:24: 11000000
INFO @ Tue, 03 Jan 2017 15:30:26: 12000000
INFO @ Tue, 03 Jan 2017 15:30:27: 13000000
INFO @ Tue, 03 Jan 2017 15:30:29: 14000000
INFO @ Tue, 03 Jan 2017 15:30:31: 15000000
INFO @ Tue, 03 Jan 2017 15:30:33: 16000000
INFO @ Tue, 03 Jan 2017 15:30:34: 17000000
INFO @ Tue, 03 Jan 2017 15:30:36: 18000000
INFO @ Tue, 03 Jan 2017 15:30:38: 19000000
INFO @ Tue, 03 Jan 2017 15:30:40: 20000000
INFO @ Tue, 03 Jan 2017 15:30:41: 21000000
INFO @ Tue, 03 Jan 2017 15:30:43: 22000000
INFO @ Tue, 03 Jan 2017 15:30:44: 23000000
INFO @ Tue, 03 Jan 2017 15:30:46: 24000000
INFO @ Tue, 03 Jan 2017 15:30:47: 25000000
INFO @ Tue, 03 Jan 2017 15:30:49: 26000000
INFO @ Tue, 03 Jan 2017 15:30:51: 27000000
INFO @ Tue, 03 Jan 2017 15:30:53: 28000000
INFO @ Tue, 03 Jan 2017 15:30:54: 29000000
INFO @ Tue, 03 Jan 2017 15:30:56: 30000000
INFO @ Tue, 03 Jan 2017 15:30:58: 31000000
INFO @ Tue, 03 Jan 2017 15:30:59: 32000000
INFO @ Tue, 03 Jan 2017 15:31:01: 33000000
INFO @ Tue, 03 Jan 2017 15:31:03: 34000000
INFO @ Tue, 03 Jan 2017 15:31:05: 35000000
INFO @ Tue, 03 Jan 2017 15:31:06: 36000000
INFO @ Tue, 03 Jan 2017 15:31:08: 37000000
INFO @ Tue, 03 Jan 2017 15:31:10: 38000000
INFO @ Tue, 03 Jan 2017 15:31:11: 39000000
INFO @ Tue, 03 Jan 2017 15:31:13: 40000000
INFO @ Tue, 03 Jan 2017 15:31:15: 41000000
INFO @ Tue, 03 Jan 2017 15:31:17: 42000000
INFO @ Tue, 03 Jan 2017 15:31:18: 43000000
INFO @ Tue, 03 Jan 2017 15:31:20: 44000000
INFO @ Tue, 03 Jan 2017 15:31:22: 45000000
INFO @ Tue, 03 Jan 2017 15:31:23: 46000000
INFO @ Tue, 03 Jan 2017 15:31:25: 47000000
INFO @ Tue, 03 Jan 2017 15:31:27: 48000000
INFO @ Tue, 03 Jan 2017 15:31:28: 49000000
INFO @ Tue, 03 Jan 2017 15:31:30: 50000000
INFO @ Tue, 03 Jan 2017 15:31:32: #1 tag size is determined as 51 bps
INFO @ Tue, 03 Jan 2017 15:31:32: #1 tag size = 51
INFO @ Tue, 03 Jan 2017 15:31:32: #1 total tags in treatment: 20752479
INFO @ Tue, 03 Jan 2017 15:31:32: #1 total tags in control: 50000000
INFO @ Tue, 03 Jan 2017 15:31:32: #1 finished!
INFO @ Tue, 03 Jan 2017 15:31:32: #2 Build Peak Model...
INFO @ Tue, 03 Jan 2017 15:31:32: #2 Skipped...
INFO @ Tue, 03 Jan 2017 15:31:32: #2 Use 170 as fragment length
INFO @ Tue, 03 Jan 2017 15:31:32: #3 Call peaks...
INFO @ Tue, 03 Jan 2017 15:31:32: #3 Call peaks with given -log10pvalue cutoff: 2.00000 ...
INFO @ Tue, 03 Jan 2017 15:31:32: #3 Pre-compute pvalue-qvalue table...
INFO @ Tue, 03 Jan 2017 15:34:23: #3 In the peak calling step, the following will be performed simultaneously:
INFO @ Tue, 03 Jan 2017 15:34:23: #3 Write bedGraph files for treatment pileup (after scaling if necessary)... /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/peak/macs2/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign_treat_pileup.bdg
INFO @ Tue, 03 Jan 2017 15:34:23: #3 Write bedGraph files for control lambda (after scaling if necessary)... /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/peak/macs2/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign_control_lambda.bdg
INFO @ Tue, 03 Jan 2017 15:34:23: #3 --SPMR is requested, so pileup will be normalized by sequencing depth in million reads.
INFO @ Tue, 03 Jan 2017 15:34:23: #3 Call peaks for each chromosome...
INFO @ Tue, 03 Jan 2017 15:39:35: #4 Write output xls file... /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/peak/macs2/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign_peaks.xls
INFO @ Tue, 03 Jan 2017 15:39:36: #4 Write peak in narrowPeak format file... /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/peak/macs2/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign_peaks.narrowPeak
INFO @ Tue, 03 Jan 2017 15:39:37: #4 Write summits bed file... /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/peak/macs2/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign_summits.bed
INFO @ Tue, 03 Jan 2017 15:39:37: Done!
INFO @ Tue, 03 Jan 2017 15:39:40:
# Command line: callpeak -t /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/align/rep1/SRR1370904_1.nodup.tagAlign.gz -c /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/HughesInputData/input15MReadsNSCLess1.05.50MSubsample.tagAlign.gz -f BED -n /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/peak/macs2/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign -g hs -p 0.01 --broad --nomodel --shift 0 --extsize 170 --keep-dup all
# ARGUMENTS LIST:
# name = /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/peak/macs2/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign
# format = BED
# ChIP-seq file = ['/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/align/rep1/SRR1370904_1.nodup.tagAlign.gz']
# control file = ['/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/HughesInputData/input15MReadsNSCLess1.05.50MSubsample.tagAlign.gz']
# effective genome size = 2.70e+09
# band width = 300
# model fold = [5, 50]
# pvalue cutoff for narrow/strong regions = 1.00e-02
# pvalue cutoff for broad/weak regions = 1.00e-01
# qvalue will not be calculated and reported as -1 in the final output.
# Larger dataset will be scaled towards smaller dataset.
# Range for calculating regional lambda is: 1000 bps and 10000 bps
# Broad region calling is on
INFO @ Tue, 03 Jan 2017 15:39:40: #1 read tag files...
INFO @ Tue, 03 Jan 2017 15:39:40: #1 read treatment tags...
INFO @ Tue, 03 Jan 2017 15:39:42: 1000000
INFO @ Tue, 03 Jan 2017 15:39:43: 2000000
INFO @ Tue, 03 Jan 2017 15:39:45: 3000000
INFO @ Tue, 03 Jan 2017 15:39:46: 4000000
INFO @ Tue, 03 Jan 2017 15:39:48: 5000000
INFO @ Tue, 03 Jan 2017 15:39:49: 6000000
INFO @ Tue, 03 Jan 2017 15:39:51: 7000000
INFO @ Tue, 03 Jan 2017 15:39:52: 8000000
INFO @ Tue, 03 Jan 2017 15:39:54: 9000000
INFO @ Tue, 03 Jan 2017 15:39:55: 10000000
INFO @ Tue, 03 Jan 2017 15:39:57: 11000000
INFO @ Tue, 03 Jan 2017 15:39:58: 12000000
INFO @ Tue, 03 Jan 2017 15:40:00: 13000000
INFO @ Tue, 03 Jan 2017 15:40:01: 14000000
INFO @ Tue, 03 Jan 2017 15:40:03: 15000000
INFO @ Tue, 03 Jan 2017 15:40:04: 16000000
INFO @ Tue, 03 Jan 2017 15:40:06: 17000000
INFO @ Tue, 03 Jan 2017 15:40:07: 18000000
INFO @ Tue, 03 Jan 2017 15:40:09: 19000000
INFO @ Tue, 03 Jan 2017 15:40:10: 20000000
INFO @ Tue, 03 Jan 2017 15:40:12: #1.2 read input tags...
INFO @ Tue, 03 Jan 2017 15:40:13: 1000000
INFO @ Tue, 03 Jan 2017 15:40:15: 2000000
INFO @ Tue, 03 Jan 2017 15:40:17: 3000000
INFO @ Tue, 03 Jan 2017 15:40:18: 4000000
INFO @ Tue, 03 Jan 2017 15:40:20: 5000000
INFO @ Tue, 03 Jan 2017 15:40:21: 6000000
INFO @ Tue, 03 Jan 2017 15:40:23: 7000000
INFO @ Tue, 03 Jan 2017 15:40:25: 8000000
INFO @ Tue, 03 Jan 2017 15:40:26: 9000000
INFO @ Tue, 03 Jan 2017 15:40:28: 10000000
INFO @ Tue, 03 Jan 2017 15:40:29: 11000000
INFO @ Tue, 03 Jan 2017 15:40:31: 12000000
INFO @ Tue, 03 Jan 2017 15:40:32: 13000000
INFO @ Tue, 03 Jan 2017 15:40:34: 14000000
INFO @ Tue, 03 Jan 2017 15:40:35: 15000000
INFO @ Tue, 03 Jan 2017 15:40:37: 16000000
INFO @ Tue, 03 Jan 2017 15:40:38: 17000000
INFO @ Tue, 03 Jan 2017 15:40:40: 18000000
INFO @ Tue, 03 Jan 2017 15:40:41: 19000000
INFO @ Tue, 03 Jan 2017 15:40:43: 20000000
INFO @ Tue, 03 Jan 2017 15:40:44: 21000000
INFO @ Tue, 03 Jan 2017 15:40:46: 22000000
INFO @ Tue, 03 Jan 2017 15:40:47: 23000000
INFO @ Tue, 03 Jan 2017 15:40:49: 24000000
INFO @ Tue, 03 Jan 2017 15:40:50: 25000000
INFO @ Tue, 03 Jan 2017 15:40:52: 26000000
INFO @ Tue, 03 Jan 2017 15:40:54: 27000000
INFO @ Tue, 03 Jan 2017 15:40:55: 28000000
INFO @ Tue, 03 Jan 2017 15:40:57: 29000000
INFO @ Tue, 03 Jan 2017 15:40:58: 30000000
INFO @ Tue, 03 Jan 2017 15:41:00: 31000000
INFO @ Tue, 03 Jan 2017 15:41:01: 32000000
INFO @ Tue, 03 Jan 2017 15:41:03: 33000000
INFO @ Tue, 03 Jan 2017 15:41:04: 34000000
INFO @ Tue, 03 Jan 2017 15:41:06: 35000000
INFO @ Tue, 03 Jan 2017 15:41:07: 36000000
INFO @ Tue, 03 Jan 2017 15:41:09: 37000000
INFO @ Tue, 03 Jan 2017 15:41:10: 38000000
INFO @ Tue, 03 Jan 2017 15:41:12: 39000000
INFO @ Tue, 03 Jan 2017 15:41:13: 40000000
INFO @ Tue, 03 Jan 2017 15:41:15: 41000000
INFO @ Tue, 03 Jan 2017 15:41:16: 42000000
INFO @ Tue, 03 Jan 2017 15:41:18: 43000000
INFO @ Tue, 03 Jan 2017 15:41:19: 44000000
INFO @ Tue, 03 Jan 2017 15:41:21: 45000000
INFO @ Tue, 03 Jan 2017 15:41:23: 46000000
INFO @ Tue, 03 Jan 2017 15:41:24: 47000000
INFO @ Tue, 03 Jan 2017 15:41:26: 48000000
INFO @ Tue, 03 Jan 2017 15:41:27: 49000000
INFO @ Tue, 03 Jan 2017 15:41:29: 50000000
INFO @ Tue, 03 Jan 2017 15:41:30: #1 tag size is determined as 51 bps
INFO @ Tue, 03 Jan 2017 15:41:30: #1 tag size = 51
INFO @ Tue, 03 Jan 2017 15:41:30: #1 total tags in treatment: 20752479
INFO @ Tue, 03 Jan 2017 15:41:30: #1 total tags in control: 50000000
INFO @ Tue, 03 Jan 2017 15:41:30: #1 finished!
INFO @ Tue, 03 Jan 2017 15:41:30: #2 Build Peak Model...
INFO @ Tue, 03 Jan 2017 15:41:30: #2 Skipped...
INFO @ Tue, 03 Jan 2017 15:41:30: #2 Use 170 as fragment length
INFO @ Tue, 03 Jan 2017 15:41:30: #3 Call peaks...
INFO @ Tue, 03 Jan 2017 15:41:30: #3 Call broad peaks with given level1 -log10pvalue cutoff and level2: 2.00000, 1.00000...
INFO @ Tue, 03 Jan 2017 15:41:30: #3 Pre-compute pvalue-qvalue table...
INFO @ Tue, 03 Jan 2017 15:44:53: #3 Call peaks for each chromosome...
INFO @ Tue, 03 Jan 2017 15:47:06: #4 Write output xls file... /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/peak/macs2/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign_peaks.xls
INFO @ Tue, 03 Jan 2017 15:47:12: #4 Write broad peak in broadPeak format file... /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/peak/macs2/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign_peaks.broadPeak
INFO @ Tue, 03 Jan 2017 15:47:20: #4 Write broad peak in bed12/gappedPeak format file... /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/peak/macs2/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign_peaks.gappedPeak
INFO @ Tue, 03 Jan 2017 15:47:24: Done!
INFO @ Tue, 03 Jan 2017 15:47:41: Read and build treatment bedGraph...
INFO @ Tue, 03 Jan 2017 15:48:25: Read and build control bedGraph...
INFO @ Tue, 03 Jan 2017 15:50:09: Build scoreTrackII...
INFO @ Tue, 03 Jan 2017 15:50:55: Calculate scores comparing treatment and control by 'FE'...
INFO @ Tue, 03 Jan 2017 15:56:17: Write bedGraph of scores...
INFO @ Tue, 03 Jan 2017 15:59:55: Finished 'FE'! Please check '/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/peak/macs2/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign_FE.bdg'!
INFO @ Tue, 03 Jan 2017 16:07:54: Read and build treatment bedGraph...
INFO @ Tue, 03 Jan 2017 16:08:35: Read and build control bedGraph...
INFO @ Tue, 03 Jan 2017 16:10:19: Build scoreTrackII...
INFO @ Tue, 03 Jan 2017 16:11:13: Values in your input bedGraph files will be multiplied by 20.752479 ...
INFO @ Tue, 03 Jan 2017 16:16:45: Calculate scores comparing treatment and control by 'ppois'...
INFO @ Tue, 03 Jan 2017 16:17:26: Write bedGraph of scores...
INFO @ Tue, 03 Jan 2017 16:20:33: Finished 'ppois'! Please check '/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/Najafabadi2015Data/hg38Mapping/ZNF684/out/peak/macs2/rep1/SRR1370904_1.nodup.tagAlign_x_input15MReadsNSCLess1.05.50MSubsample.tagAlign_ppois.bdg'!
|