| Num | 1 | | ID | task.callpeak_macs2.macs2_rep1.line_67.id_10 | | Name | macs2 rep1 | | Thread | thread_Root | | PID | 4558 | | | OK | true | | Exit Code | 0 | | Retries | | | State | FINISHED | | Dep. | OK | | Cpus | | | Mem | | | | Start | 2016-11-07 13:06:12 | | End | 2016-11-07 13:47:53 | | Elapsed | 00:41:40 | | Timeout | 00:00:-1 | | Wall Timeout | 100 days | | | Input files | /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/align/rep1/ZF5_allReps.nodup.tagAlign.gz
/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/InputData/r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign.gz
| | Output files | /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/peak/macs2/rep1/ZF5_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign.narrowPeak.gz
/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/peak/macs2/rep1/ZF5_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign.gappedPeak.gz
/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/signal/macs2/rep1/ZF5_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign.fc.signal.bw
/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/signal/macs2/rep1/ZF5_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign.pval.signal.bw
| | Dependencies | | |
# SYS command. line 69
if [[ -f $(which conda) && $(conda env list | grep aquas_chipseq | wc -l) != "0" ]]; then source activate aquas_chipseq; fi; export PATH=/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/src/TF_chipseq_pipeline/.:/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/src/TF_chipseq_pipeline/modules:/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/src/TF_chipseq_pipeline/utils:${PATH}:/bin:/usr/bin:/usr/local/bin:${HOME}/.bds; set -o pipefail; STARTTIME=$(date +%s); renice -n 10 $$
# SYS command. line 70
export LC_COLLATE=C
# SYS command. line 75
macs2 callpeak -t /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/align/rep1/ZF5_allReps.nodup.tagAlign.gz -c /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/InputData/r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign.gz -f BED -n /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/peak/macs2/rep1/ZF5_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign -g mm -p 0.01 --nomodel --shift 0 --extsize 95 --keep-dup all -B --SPMR
# SYS command. line 78
sort -k 8gr,8gr "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/peak/macs2/rep1/ZF5_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign"_peaks.narrowPeak | awk 'BEGIN{OFS="\t"}{$4="Peak_"NR ; print $0}' | gzip -c > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/peak/macs2/rep1/ZF5_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign.narrowPeak.gz
# SYS command. line 81
rm -f "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/peak/macs2/rep1/ZF5_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign"_peaks.xls \
"/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/peak/macs2/rep1/ZF5_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign"_peaks.narrowPeak \
"/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/peak/macs2/rep1/ZF5_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign"_summits.bed
# SYS command. line 88
macs2 callpeak -t /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/align/rep1/ZF5_allReps.nodup.tagAlign.gz -c /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/InputData/r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign.gz -f BED -n /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/peak/macs2/rep1/ZF5_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign -g mm -p 0.01 --broad --nomodel --shift 0 --extsize 95 --keep-dup all
# SYS command. line 91
sort -k 8gr,8gr "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/peak/macs2/rep1/ZF5_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign"_peaks.broadPeak | awk 'BEGIN{OFS="\t"}{$4="Peak_"NR ; print $0}' | gzip -c > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/peak/macs2/rep1/ZF5_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign.broadPeak.gz
# SYS command. line 92
sort -k 14gr,14gr "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/peak/macs2/rep1/ZF5_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign"_peaks.gappedPeak | awk 'BEGIN{OFS="\t"}{$4="Peak_"NR ; print $0}' | gzip -c > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/peak/macs2/rep1/ZF5_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign.gappedPeak.gz
# SYS command. line 95
rm -f "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/peak/macs2/rep1/ZF5_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign"_peaks.xls \
"/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/peak/macs2/rep1/ZF5_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign"_peaks.broadPeak \
"/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/peak/macs2/rep1/ZF5_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign"_peaks.gappedPeak \
"/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/peak/macs2/rep1/ZF5_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign"_summits.bed
# SYS command. line 100
if [[ true == "false" ]]; then \
rm -f "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/peak/macs2/rep1/ZF5_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign"_treat_pileup.bdg "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/peak/macs2/rep1/ZF5_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.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 111
macs2 bdgcmp -t "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/peak/macs2/rep1/ZF5_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign"_treat_pileup.bdg -c "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/peak/macs2/rep1/ZF5_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign"_control_lambda.bdg --outdir /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/peak/macs2/rep1 -o "ZF5_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign"_FE.bdg -m FE
# SYS command. line 114
slopBed -i "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/peak/macs2/rep1/ZF5_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign"_FE.bdg -g /mnt/data/bds_pipeline_genome_data/mm10/mm10.chrom.sizes -b 0 | awk '{if ($3 != -1) print $0}' | bedClip stdin /mnt/data/bds_pipeline_genome_data/mm10/mm10.chrom.sizes /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/peak/macs2/rep1/ZF5_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign.fc.signal.bedgraph
# SYS command. line 115
rm -f "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/peak/macs2/rep1/ZF5_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign"_FE.bdg
# SYS command. line 118
sort -k1,1 -k2,2n /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/peak/macs2/rep1/ZF5_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign.fc.signal.bedgraph > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/peak/macs2/rep1/ZF5_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign.fc.signal.srt.bedgraph
# SYS command. line 119
bedGraphToBigWig /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/peak/macs2/rep1/ZF5_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign.fc.signal.srt.bedgraph /mnt/data/bds_pipeline_genome_data/mm10/mm10.chrom.sizes /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/signal/macs2/rep1/ZF5_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign.fc.signal.bw
# SYS command. line 120
rm -f /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/peak/macs2/rep1/ZF5_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign.fc.signal.bedgraph /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/peak/macs2/rep1/ZF5_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign.fc.signal.srt.bedgraph
# SYS command. line 126
chipReads=$(zcat /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/align/rep1/ZF5_allReps.nodup.tagAlign.gz | wc -l | awk '{printf "%f", $1/1000000}')
# SYS command. line 128
controlReads=$(zcat /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/InputData/r120621A_l1i7_aB180__wt__bioBirAonly.nodup.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 130
macs2 bdgcmp -t "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/peak/macs2/rep1/ZF5_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign"_treat_pileup.bdg -c "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/peak/macs2/rep1/ZF5_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign"_control_lambda.bdg --outdir /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/peak/macs2/rep1 -o "ZF5_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign"_ppois.bdg -m ppois -S "${sval}"
# SYS command. line 133
slopBed -i "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/peak/macs2/rep1/ZF5_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign"_ppois.bdg -g /mnt/data/bds_pipeline_genome_data/mm10/mm10.chrom.sizes -b 0 | awk '{if ($3 != -1) print $0}' | bedClip stdin /mnt/data/bds_pipeline_genome_data/mm10/mm10.chrom.sizes /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/peak/macs2/rep1/ZF5_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign.pval.signal.bedgraph
# SYS command. line 134
rm -rf "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/peak/macs2/rep1/ZF5_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign"_ppois.bdg
# SYS command. line 137
sort -k1,1 -k2,2n /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/peak/macs2/rep1/ZF5_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign.pval.signal.bedgraph > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/peak/macs2/rep1/ZF5_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign.pval.signal.srt.bedgraph
# SYS command. line 138
bedGraphToBigWig /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/peak/macs2/rep1/ZF5_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign.pval.signal.srt.bedgraph /mnt/data/bds_pipeline_genome_data/mm10/mm10.chrom.sizes /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/signal/macs2/rep1/ZF5_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign.pval.signal.bw
# SYS command. line 139
rm -f /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/peak/macs2/rep1/ZF5_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign.pval.signal.bedgraph /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/peak/macs2/rep1/ZF5_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign.pval.signal.srt.bedgraph
# SYS command. line 141
rm -f "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/peak/macs2/rep1/ZF5_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign"_treat_pileup.bdg "/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/peak/macs2/rep1/ZF5_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign"_control_lambda.bdg
# SYS command. line 143
TASKTIME=$[$(date +%s)-${STARTTIME}]; if [ ${TASKTIME} -lt 60 ]; then echo "Waiting for $[60-${TASKTIME}] seconds."; sleep $[60-${TASKTIME}]; fi
--------------------Stdout--------------------
6169 (process ID) old priority 0, new priority 10
--------------------Stderr--------------------
INFO @ Mon, 07 Nov 2016 13:06:33:
# Command line: callpeak -t /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/align/rep1/ZF5_allReps.nodup.tagAlign.gz -c /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/InputData/r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign.gz -f BED -n /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/peak/macs2/rep1/ZF5_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign -g mm -p 0.01 --nomodel --shift 0 --extsize 95 --keep-dup all -B --SPMR
# ARGUMENTS LIST:
# name = /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/peak/macs2/rep1/ZF5_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign
# format = BED
# ChIP-seq file = ['/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/align/rep1/ZF5_allReps.nodup.tagAlign.gz']
# control file = ['/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/InputData/r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign.gz']
# effective genome size = 1.87e+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 @ Mon, 07 Nov 2016 13:06:33: #1 read tag files...
INFO @ Mon, 07 Nov 2016 13:06:33: #1 read treatment tags...
INFO @ Mon, 07 Nov 2016 13:06:35: 1000000
INFO @ Mon, 07 Nov 2016 13:06:37: 2000000
INFO @ Mon, 07 Nov 2016 13:06:38: 3000000
INFO @ Mon, 07 Nov 2016 13:06:40: 4000000
INFO @ Mon, 07 Nov 2016 13:06:41: 5000000
INFO @ Mon, 07 Nov 2016 13:06:43: 6000000
INFO @ Mon, 07 Nov 2016 13:06:44: 7000000
INFO @ Mon, 07 Nov 2016 13:06:46: 8000000
INFO @ Mon, 07 Nov 2016 13:06:47: 9000000
INFO @ Mon, 07 Nov 2016 13:06:49: 10000000
INFO @ Mon, 07 Nov 2016 13:06:50: 11000000
INFO @ Mon, 07 Nov 2016 13:06:52: 12000000
INFO @ Mon, 07 Nov 2016 13:06:54: 13000000
INFO @ Mon, 07 Nov 2016 13:06:55: 14000000
INFO @ Mon, 07 Nov 2016 13:06:57: 15000000
INFO @ Mon, 07 Nov 2016 13:06:58: 16000000
INFO @ Mon, 07 Nov 2016 13:07:00: 17000000
INFO @ Mon, 07 Nov 2016 13:07:01: 18000000
INFO @ Mon, 07 Nov 2016 13:07:03: 19000000
INFO @ Mon, 07 Nov 2016 13:07:04: 20000000
INFO @ Mon, 07 Nov 2016 13:07:06: 21000000
INFO @ Mon, 07 Nov 2016 13:07:08: 22000000
INFO @ Mon, 07 Nov 2016 13:07:09: 23000000
INFO @ Mon, 07 Nov 2016 13:07:11: 24000000
INFO @ Mon, 07 Nov 2016 13:07:13: 25000000
INFO @ Mon, 07 Nov 2016 13:07:14: 26000000
INFO @ Mon, 07 Nov 2016 13:07:16: 27000000
INFO @ Mon, 07 Nov 2016 13:07:18: 28000000
INFO @ Mon, 07 Nov 2016 13:07:19: 29000000
INFO @ Mon, 07 Nov 2016 13:07:21: 30000000
INFO @ Mon, 07 Nov 2016 13:07:22: 31000000
INFO @ Mon, 07 Nov 2016 13:07:24: 32000000
INFO @ Mon, 07 Nov 2016 13:07:26: 33000000
INFO @ Mon, 07 Nov 2016 13:07:27: 34000000
INFO @ Mon, 07 Nov 2016 13:07:29: 35000000
INFO @ Mon, 07 Nov 2016 13:07:31: 36000000
INFO @ Mon, 07 Nov 2016 13:07:32: 37000000
INFO @ Mon, 07 Nov 2016 13:07:34: 38000000
INFO @ Mon, 07 Nov 2016 13:07:35: 39000000
INFO @ Mon, 07 Nov 2016 13:07:37: 40000000
INFO @ Mon, 07 Nov 2016 13:07:38: 41000000
INFO @ Mon, 07 Nov 2016 13:07:40: 42000000
INFO @ Mon, 07 Nov 2016 13:07:42: 43000000
INFO @ Mon, 07 Nov 2016 13:07:44: 44000000
INFO @ Mon, 07 Nov 2016 13:07:45: 45000000
INFO @ Mon, 07 Nov 2016 13:07:47: 46000000
INFO @ Mon, 07 Nov 2016 13:07:48: 47000000
INFO @ Mon, 07 Nov 2016 13:07:50: 48000000
INFO @ Mon, 07 Nov 2016 13:07:52: 49000000
INFO @ Mon, 07 Nov 2016 13:07:53: 50000000
INFO @ Mon, 07 Nov 2016 13:07:55: 51000000
INFO @ Mon, 07 Nov 2016 13:07:57: 52000000
INFO @ Mon, 07 Nov 2016 13:07:58: 53000000
INFO @ Mon, 07 Nov 2016 13:08:00: 54000000
INFO @ Mon, 07 Nov 2016 13:08:02: 55000000
INFO @ Mon, 07 Nov 2016 13:08:03: 56000000
INFO @ Mon, 07 Nov 2016 13:08:05: 57000000
INFO @ Mon, 07 Nov 2016 13:08:06: 58000000
INFO @ Mon, 07 Nov 2016 13:08:08: 59000000
INFO @ Mon, 07 Nov 2016 13:08:09: #1.2 read input tags...
INFO @ Mon, 07 Nov 2016 13:08:11: 1000000
INFO @ Mon, 07 Nov 2016 13:08:13: 2000000
INFO @ Mon, 07 Nov 2016 13:08:14: 3000000
INFO @ Mon, 07 Nov 2016 13:08:15: #1 tag size is determined as 49 bps
INFO @ Mon, 07 Nov 2016 13:08:15: #1 tag size = 49
INFO @ Mon, 07 Nov 2016 13:08:15: #1 total tags in treatment: 59491042
INFO @ Mon, 07 Nov 2016 13:08:15: #1 total tags in control: 3384182
INFO @ Mon, 07 Nov 2016 13:08:15: #1 finished!
INFO @ Mon, 07 Nov 2016 13:08:15: #2 Build Peak Model...
INFO @ Mon, 07 Nov 2016 13:08:15: #2 Skipped...
INFO @ Mon, 07 Nov 2016 13:08:15: #2 Use 95 as fragment length
INFO @ Mon, 07 Nov 2016 13:08:15: #3 Call peaks...
INFO @ Mon, 07 Nov 2016 13:08:15: #3 Call peaks with given -log10pvalue cutoff: 2.00000 ...
INFO @ Mon, 07 Nov 2016 13:08:15: #3 Pre-compute pvalue-qvalue table...
INFO @ Mon, 07 Nov 2016 13:09:49: #3 In the peak calling step, the following will be performed simultaneously:
INFO @ Mon, 07 Nov 2016 13:09:49: #3 Write bedGraph files for treatment pileup (after scaling if necessary)... /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/peak/macs2/rep1/ZF5_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign_treat_pileup.bdg
INFO @ Mon, 07 Nov 2016 13:09:49: #3 Write bedGraph files for control lambda (after scaling if necessary)... /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/peak/macs2/rep1/ZF5_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign_control_lambda.bdg
INFO @ Mon, 07 Nov 2016 13:09:49: #3 --SPMR is requested, so pileup will be normalized by sequencing depth in million reads.
INFO @ Mon, 07 Nov 2016 13:09:49: #3 Call peaks for each chromosome...
INFO @ Mon, 07 Nov 2016 13:14:36: #4 Write output xls file... /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/peak/macs2/rep1/ZF5_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign_peaks.xls
INFO @ Mon, 07 Nov 2016 13:14:36: #4 Write peak in narrowPeak format file... /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/peak/macs2/rep1/ZF5_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign_peaks.narrowPeak
INFO @ Mon, 07 Nov 2016 13:14:36: #4 Write summits bed file... /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/peak/macs2/rep1/ZF5_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign_summits.bed
INFO @ Mon, 07 Nov 2016 13:14:36: Done!
INFO @ Mon, 07 Nov 2016 13:14:38:
# Command line: callpeak -t /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/align/rep1/ZF5_allReps.nodup.tagAlign.gz -c /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/InputData/r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign.gz -f BED -n /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/peak/macs2/rep1/ZF5_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign -g mm -p 0.01 --broad --nomodel --shift 0 --extsize 95 --keep-dup all
# ARGUMENTS LIST:
# name = /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/peak/macs2/rep1/ZF5_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign
# format = BED
# ChIP-seq file = ['/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/align/rep1/ZF5_allReps.nodup.tagAlign.gz']
# control file = ['/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/InputData/r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign.gz']
# effective genome size = 1.87e+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 @ Mon, 07 Nov 2016 13:14:38: #1 read tag files...
INFO @ Mon, 07 Nov 2016 13:14:38: #1 read treatment tags...
INFO @ Mon, 07 Nov 2016 13:14:40: 1000000
INFO @ Mon, 07 Nov 2016 13:14:41: 2000000
INFO @ Mon, 07 Nov 2016 13:14:43: 3000000
INFO @ Mon, 07 Nov 2016 13:14:45: 4000000
INFO @ Mon, 07 Nov 2016 13:14:47: 5000000
INFO @ Mon, 07 Nov 2016 13:14:48: 6000000
INFO @ Mon, 07 Nov 2016 13:14:50: 7000000
INFO @ Mon, 07 Nov 2016 13:14:51: 8000000
INFO @ Mon, 07 Nov 2016 13:14:53: 9000000
INFO @ Mon, 07 Nov 2016 13:14:55: 10000000
INFO @ Mon, 07 Nov 2016 13:14:56: 11000000
INFO @ Mon, 07 Nov 2016 13:14:58: 12000000
INFO @ Mon, 07 Nov 2016 13:15:00: 13000000
INFO @ Mon, 07 Nov 2016 13:15:01: 14000000
INFO @ Mon, 07 Nov 2016 13:15:03: 15000000
INFO @ Mon, 07 Nov 2016 13:15:05: 16000000
INFO @ Mon, 07 Nov 2016 13:15:06: 17000000
INFO @ Mon, 07 Nov 2016 13:15:08: 18000000
INFO @ Mon, 07 Nov 2016 13:15:10: 19000000
INFO @ Mon, 07 Nov 2016 13:15:11: 20000000
INFO @ Mon, 07 Nov 2016 13:15:13: 21000000
INFO @ Mon, 07 Nov 2016 13:15:15: 22000000
INFO @ Mon, 07 Nov 2016 13:15:16: 23000000
INFO @ Mon, 07 Nov 2016 13:15:18: 24000000
INFO @ Mon, 07 Nov 2016 13:15:20: 25000000
INFO @ Mon, 07 Nov 2016 13:15:21: 26000000
INFO @ Mon, 07 Nov 2016 13:15:23: 27000000
INFO @ Mon, 07 Nov 2016 13:15:25: 28000000
INFO @ Mon, 07 Nov 2016 13:15:26: 29000000
INFO @ Mon, 07 Nov 2016 13:15:28: 30000000
INFO @ Mon, 07 Nov 2016 13:15:29: 31000000
INFO @ Mon, 07 Nov 2016 13:15:31: 32000000
INFO @ Mon, 07 Nov 2016 13:15:33: 33000000
INFO @ Mon, 07 Nov 2016 13:15:35: 34000000
INFO @ Mon, 07 Nov 2016 13:15:36: 35000000
INFO @ Mon, 07 Nov 2016 13:15:38: 36000000
INFO @ Mon, 07 Nov 2016 13:15:40: 37000000
INFO @ Mon, 07 Nov 2016 13:15:41: 38000000
INFO @ Mon, 07 Nov 2016 13:15:43: 39000000
INFO @ Mon, 07 Nov 2016 13:15:45: 40000000
INFO @ Mon, 07 Nov 2016 13:15:46: 41000000
INFO @ Mon, 07 Nov 2016 13:15:48: 42000000
INFO @ Mon, 07 Nov 2016 13:15:50: 43000000
INFO @ Mon, 07 Nov 2016 13:15:51: 44000000
INFO @ Mon, 07 Nov 2016 13:15:53: 45000000
INFO @ Mon, 07 Nov 2016 13:15:55: 46000000
INFO @ Mon, 07 Nov 2016 13:15:56: 47000000
INFO @ Mon, 07 Nov 2016 13:15:58: 48000000
INFO @ Mon, 07 Nov 2016 13:16:00: 49000000
INFO @ Mon, 07 Nov 2016 13:16:01: 50000000
INFO @ Mon, 07 Nov 2016 13:16:03: 51000000
INFO @ Mon, 07 Nov 2016 13:16:05: 52000000
INFO @ Mon, 07 Nov 2016 13:16:06: 53000000
INFO @ Mon, 07 Nov 2016 13:16:08: 54000000
INFO @ Mon, 07 Nov 2016 13:16:10: 55000000
INFO @ Mon, 07 Nov 2016 13:16:11: 56000000
INFO @ Mon, 07 Nov 2016 13:16:13: 57000000
INFO @ Mon, 07 Nov 2016 13:16:14: 58000000
INFO @ Mon, 07 Nov 2016 13:16:16: 59000000
INFO @ Mon, 07 Nov 2016 13:16:18: #1.2 read input tags...
INFO @ Mon, 07 Nov 2016 13:16:19: 1000000
INFO @ Mon, 07 Nov 2016 13:16:21: 2000000
INFO @ Mon, 07 Nov 2016 13:16:23: 3000000
INFO @ Mon, 07 Nov 2016 13:16:23: #1 tag size is determined as 49 bps
INFO @ Mon, 07 Nov 2016 13:16:23: #1 tag size = 49
INFO @ Mon, 07 Nov 2016 13:16:23: #1 total tags in treatment: 59491042
INFO @ Mon, 07 Nov 2016 13:16:23: #1 total tags in control: 3384182
INFO @ Mon, 07 Nov 2016 13:16:23: #1 finished!
INFO @ Mon, 07 Nov 2016 13:16:23: #2 Build Peak Model...
INFO @ Mon, 07 Nov 2016 13:16:23: #2 Skipped...
INFO @ Mon, 07 Nov 2016 13:16:23: #2 Use 95 as fragment length
INFO @ Mon, 07 Nov 2016 13:16:23: #3 Call peaks...
INFO @ Mon, 07 Nov 2016 13:16:23: #3 Call broad peaks with given level1 -log10pvalue cutoff and level2: 2.00000, 1.00000...
INFO @ Mon, 07 Nov 2016 13:16:23: #3 Pre-compute pvalue-qvalue table...
INFO @ Mon, 07 Nov 2016 13:17:34: #3 Call peaks for each chromosome...
INFO @ Mon, 07 Nov 2016 13:18:04: #4 Write output xls file... /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/peak/macs2/rep1/ZF5_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign_peaks.xls
INFO @ Mon, 07 Nov 2016 13:18:04: #4 Write broad peak in broadPeak format file... /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/peak/macs2/rep1/ZF5_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign_peaks.broadPeak
INFO @ Mon, 07 Nov 2016 13:18:04: #4 Write broad peak in bed12/gappedPeak format file... /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/peak/macs2/rep1/ZF5_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign_peaks.gappedPeak
INFO @ Mon, 07 Nov 2016 13:18:05: Done!
INFO @ Mon, 07 Nov 2016 13:18:07: Read and build treatment bedGraph...
INFO @ Mon, 07 Nov 2016 13:20:07: Read and build control bedGraph...
INFO @ Mon, 07 Nov 2016 13:20:20: Build scoreTrackII...
INFO @ Mon, 07 Nov 2016 13:21:02: Calculate scores comparing treatment and control by 'FE'...
INFO @ Mon, 07 Nov 2016 13:26:38: Write bedGraph of scores...
INFO @ Mon, 07 Nov 2016 13:30:17: Finished 'FE'! Please check '/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/peak/macs2/rep1/ZF5_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign_FE.bdg'!
INFO @ Mon, 07 Nov 2016 13:38:08: Read and build treatment bedGraph...
INFO @ Mon, 07 Nov 2016 13:40:05: Read and build control bedGraph...
INFO @ Mon, 07 Nov 2016 13:40:17: Build scoreTrackII...
INFO @ Mon, 07 Nov 2016 13:40:54: Values in your input bedGraph files will be multiplied by 3.384182 ...
INFO @ Mon, 07 Nov 2016 13:45:56: Calculate scores comparing treatment and control by 'ppois'...
INFO @ Mon, 07 Nov 2016 13:46:31: Write bedGraph of scores...
INFO @ Mon, 07 Nov 2016 13:47:02: Finished 'ppois'! Please check '/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF5/out/peak/macs2/rep1/ZF5_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign_ppois.bdg'!
|