BigDataScript report: chipseq.bds

Script file /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/src/TF_chipseq_pipeline/chipseq.bds
Program ID chipseq.bds.20161109_175132_676
Start time 2016-11-09 17:51:32
Run time 03:29:39.202
Tasks executed 3
Tasks failed 2
Tasks failed names 
markdup_bam_picard ctl1
bwa_aln rep1
Arguments* [-title, SRR1370912_1, -nth, 5, -macs2_for_pooled_rep_only, -species, hg38, -url_base, http://mitra.stanford.edu/kundaje/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370912_1/out, -fastq, /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/RawReadsFiles/SRR1370912_1.fastq.gz, -ctl_bam, /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/InputData/ss_50M_2000_GRCh38.bam, -final_stage, xcor, $]
System* sge
Cpus* -1
Exit value 1
* Values in global scope when program finished execution.

Timeline

Task Graph

Parallel (threads) details

Thread ID Parent Tasks
thread_Root None
thread_41 thread_Root 
 chipseq.bds.20161109_175132_676_parallel_41/task.align_bwa.bwa_aln_rep1.line_108.id_10
thread_42 thread_Root 
 chipseq.bds.20161109_175132_676_parallel_42/task.postalign_bam.dedup_bam_1_ctl1.line_95.id_11
chipseq.bds.20161109_175132_676_parallel_42/task.postalign_bam.markdup_bam_picard_ctl1.line_347.id_12

Task details

Task Execution Time Dependencies Task program, Errors, StdOut / StdErr
Num 1
ID task.align_bwa.bwa_aln_rep1.line_108.id_10
Name bwa_aln rep1
Thread thread_41
PID 6000
OK false
Exit Code 3
Retries
State RUNNING
Dep. WAIT
Cpus
Mem
Start 2016-11-09 17:51:36
End 2016-11-09 17:51:36
Elapsed 00:00:00
Timeout 00:00:-1
Wall Timeout 100 days
Input files /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/RawReadsFiles/SRR1370912_1.fastq.gz
Output files /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370912_1/out/align/rep1/SRR1370912_1.sai
Dependencies 
 
# SYS command. line 110

 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 113

 bwa aln -q 5 -l 32 -k 2 -t 1 /mnt/data/bds_pipeline_genome_data/hg38/bwa_index/GRCh38.p3.genome.fa /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/RawReadsFiles/SRR1370912_1.fastq.gz > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370912_1/out/align/rep1/SRR1370912_1.sai

# SYS command. line 115

 TASKTIME=$[$(date +%s)-${STARTTIME}]; if [ ${TASKTIME} -lt 60 ]; then echo "Waiting for $[60-${TASKTIME}] seconds."; sleep $[60-${TASKTIME}]; fi
 
   
--------------------Stdout--------------------
11867 (process ID) old priority 0, new priority 10

 
--------------------Stderr--------------------
[bwa_aln] 17bp reads: max_diff = 2
[bwa_aln] 38bp reads: max_diff = 3
[bwa_aln] 64bp reads: max_diff = 4
[bwa_aln] 93bp reads: max_diff = 5
[bwa_aln] 124bp reads: max_diff = 6
[bwa_aln] 157bp reads: max_diff = 7
[bwa_aln] 190bp reads: max_diff = 8
[bwa_aln] 225bp reads: max_diff = 9
[bwa_read_seq] 0.3% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 57.01 sec
[bwa_aln_core] write to the disk... 0.04 sec
[bwa_aln_core] 262144 sequences have been processed.
[bwa_read_seq] 0.3% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 51.72 sec
[bwa_aln_core] write to the disk... 0.03 sec
[bwa_aln_core] 524288 sequences have been processed.
[bwa_read_seq] 0.3% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 50.76 sec
[bwa_aln_core] write to the disk... 0.03 sec
[bwa_aln_core] 786432 sequences have been processed.
[bwa_read_seq] 0.3% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 50.58 sec
[bwa_aln_core] write to the disk... 0.04 sec
[bwa_aln_core] 1048576 sequences have been processed.
[bwa_read_seq] 0.4% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 50.50 sec
[bwa_aln_core] write to the disk... 0.03 sec
[bwa_aln_core] 1310720 sequences have been processed.
[bwa_read_seq] 1.1% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 52.09 sec
[bwa_aln_core] write to the disk... 0.03 sec
[bwa_aln_core] 1572864 sequences have been processed.
[bwa_read_seq] 0.7% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 50.70 sec
[bwa_aln_core] write to the disk... 0.04 sec
[bwa_aln_core] 1835008 sequences have been processed.
[bwa_read_seq] 0.4% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 50.12 sec
[bwa_aln_core] write to the disk... 0.04 sec
[bwa_aln_core] 2097152 sequences have been processed.
[bwa_read_seq] 0.3% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 50.25 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 2359296 sequences have been processed.
[bwa_read_seq] 0.3% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 49.71 sec
[bwa_aln_core] write to the disk... 0.03 sec
[bwa_aln_core] 2621440 sequences have been processed.
[bwa_read_seq] 0.3% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 49.50 sec
[bwa_aln_core] write to the disk... 0.03 sec
[bwa_aln_core] 2883584 sequences have been processed.
[bwa_read_seq] 0.3% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 49.83 sec
[bwa_aln_core] write to the disk... 0.03 sec
[bwa_aln_core] 3145728 sequences have been processed.
[bwa_read_seq] 0.4% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 50.13 sec
[bwa_aln_core] write to the disk... 0.03 sec
[bwa_aln_core] 3407872 sequences have been processed.
[bwa_read_seq] 0.4% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 50.38 sec
[bwa_aln_core] write to the disk... 0.04 sec
[bwa_aln_core] 3670016 sequences have been processed.
[bwa_read_seq] 0.5% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 50.33 sec
[bwa_aln_core] write to the disk... 0.03 sec
[bwa_aln_core] 3932160 sequences have been processed.
[bwa_read_seq] 0.3% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 50.84 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 4194304 sequences have been processed.
[bwa_read_seq] 0.3% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 50.08 sec
[bwa_aln_core] write to the disk... 0.03 sec
[bwa_aln_core] 4456448 sequences have been processed.
[bwa_read_seq] 0.3% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 49.72 sec
[bwa_aln_core] write to the disk... 0.03 sec
[bwa_aln_core] 4718592 sequences have been processed.
[bwa_read_seq] 0.3% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 49.31 sec
[bwa_aln_core] write to the disk... 0.03 sec
[bwa_aln_core] 4980736 sequences have been processed.
[bwa_read_seq] 0.4% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 49.84 sec
[bwa_aln_core] write to the disk... 0.03 sec
[bwa_aln_core] 5242880 sequences have been processed.
[bwa_read_seq] 0.8% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 51.16 sec
[bwa_aln_core] write to the disk... 0.04 sec
[bwa_aln_core] 5505024 sequences have been processed.
[bwa_read_seq] 0.7% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 50.95 sec
[bwa_aln_core] write to the disk... 0.03 sec
[bwa_aln_core] 5767168 sequences have been processed.
[bwa_read_seq] 0.5% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 51.12 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 6029312 sequences have been processed.
[bwa_read_seq] 0.4% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 51.17 sec
[bwa_aln_core] write to the disk... 0.03 sec
[bwa_aln_core] 6291456 sequences have been processed.
[bwa_read_seq] 0.4% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 50.13 sec
[bwa_aln_core] write to the disk... 0.03 sec
[bwa_aln_core] 6553600 sequences have been processed.
[bwa_read_seq] 0.4% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 59.45 sec
[bwa_aln_core] write to the disk... 0.04 sec
[bwa_aln_core] 6815744 sequences have been processed.
[bwa_read_seq] 0.5% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 61.72 sec
[bwa_aln_core] write to the disk... 0.04 sec
[bwa_aln_core] 7077888 sequences have been processed.
[bwa_read_seq] 0.7% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 81.58 sec
[bwa_aln_core] write to the disk... 0.04 sec
[bwa_aln_core] 7340032 sequences have been processed.
[bwa_read_seq] 1.9% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 115.64 sec
[bwa_aln_core] write to the disk... 0.03 sec
[bwa_aln_core] 7602176 sequences have been processed.
[bwa_read_seq] 1.1% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 95.32 sec
[bwa_aln_core] write to the disk... 0.04 sec
[bwa_aln_core] 7864320 sequences have been processed.
[bwa_read_seq] 0.4% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 51.76 sec
[bwa_aln_core] write to the disk... 0.03 sec
[bwa_aln_core] 8126464 sequences have been processed.
[bwa_read_seq] 0.4% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 51.21 sec
[bwa_aln_core] write to the disk... 0.04 sec
[bwa_aln_core] 8388608 sequences have been processed.
[bwa_read_seq] 0.4% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 50.57 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 8650752 sequences have been processed.
[bwa_read_seq] 0.4% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 52.41 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 8912896 sequences have been processed.
[bwa_read_seq] 0.4% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 52.23 sec
[bwa_aln_core] write to the disk... 0.03 sec
[bwa_aln_core] 9175040 sequences have been processed.
[bwa_read_seq] 0.5% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 58.27 sec
[bwa_aln_core] write to the disk... 0.04 sec
[bwa_aln_core] 9437184 sequences have been processed.
[bwa_read_seq] 0.5% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 68.53 sec
[bwa_aln_core] write to the disk... 0.04 sec
[bwa_aln_core] 9699328 sequences have been processed.
[bwa_read_seq] 0.5% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 65.76 sec
[bwa_aln_core] write to the disk... 0.04 sec
[bwa_aln_core] 9961472 sequences have been processed.
[bwa_read_seq] 0.4% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 66.96 sec
[bwa_aln_core] write to the disk... 0.03 sec
[bwa_aln_core] 10223616 sequences have been processed.
[bwa_read_seq] 0.4% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 55.21 sec
[bwa_aln_core] write to the disk... 0.03 sec
[bwa_aln_core] 10485760 sequences have been processed.
[bwa_read_seq] 0.4% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 55.28 sec
[bwa_aln_core] write to the disk... 0.04 sec
[bwa_aln_core] 10747904 sequences have been processed.
[bwa_read_seq] 0.5% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 60.01 sec
[bwa_aln_core] write to the disk... 0.04 sec
[bwa_aln_core] 11010048 sequences have been processed.
[bwa_read_seq] 0.8% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 74.80 sec
[bwa_aln_core] write to the disk... 0.04 sec
[bwa_aln_core] 11272192 sequences have been processed.
[bwa_read_seq] 1.1% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 100.01 sec
[bwa_aln_core] write to the disk... 0.04 sec
[bwa_aln_core] 11534336 sequences have been processed.
[bwa_read_seq] 0.8% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 79.28 sec
[bwa_aln_core] write to the disk... 0.04 sec
[bwa_aln_core] 11796480 sequences have been processed.
[bwa_read_seq] 0.3% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 51.03 sec
[bwa_aln_core] write to the disk... 0.03 sec
[bwa_aln_core] 12058624 sequences have been processed.
[bwa_read_seq] 0.3% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 50.23 sec
[bwa_aln_core] write to the disk... 0.03 sec
[bwa_aln_core] 12320768 sequences have been processed.
[bwa_read_seq] 0.3% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 50.01 sec
[bwa_aln_core] write to the disk... 0.03 sec
[bwa_aln_core] 12582912 sequences have been processed.
[bwa_read_seq] 0.3% bases are trimmed.
[bwa_aln_core] calculate SA coordinate...
Num 2
ID task.postalign_bam.dedup_bam_1_ctl1.line_95.id_11
Name dedup_bam_1 ctl1
Thread thread_42
PID 6003
OK true
Exit Code 0
Retries
State FINISHED
Dep. OK
Cpus 4
Mem
Start 2016-11-09 17:51:36
End 2016-11-09 18:48:55
Elapsed 00:57:18
Timeout 00:00:-1
Wall Timeout 100 days
Input files /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/InputData/ss_50M_2000_GRCh38.bam
Output files /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370912_1/out/align/ctl1/ss_50M_2000_GRCh38.filt.bam
Dependencies 
 
# SYS command. line 97

 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 105

 if [[ 0 > 0 ]]; then \
				sambamba sort -t 4 /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/InputData/ss_50M_2000_GRCh38.bam -n -o /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370912_1/out/align/ctl1/ss_50M_2000_GRCh38.qnmsrt.bam; \
				samtools view -h /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370912_1/out/align/ctl1/ss_50M_2000_GRCh38.qnmsrt.bam | $(which assign_multimappers.py) -k 0 | \
				samtools view -F 1804 -Su /dev/stdin | \
				sambamba sort -t 4 /dev/stdin -o /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370912_1/out/align/ctl1/ss_50M_2000_GRCh38.filt.bam; \
				rm -f /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370912_1/out/align/ctl1/ss_50M_2000_GRCh38.qnmsrt.bam; \
			    else \
			    	samtools view -F 1804 -q 30 -u /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/InputData/ss_50M_2000_GRCh38.bam | \
			    	sambamba sort -t 4 /dev/stdin -o /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370912_1/out/align/ctl1/ss_50M_2000_GRCh38.filt.bam; \
			    fi

# SYS command. line 116

 TASKTIME=$[$(date +%s)-${STARTTIME}]; if [ ${TASKTIME} -lt 60 ]; then echo "Waiting for $[60-${TASKTIME}] seconds."; sleep $[60-${TASKTIME}]; fi
 
   
--------------------Stdout--------------------
26087 (process ID) old priority 0, new priority 10
Num 3
ID task.postalign_bam.markdup_bam_picard_ctl1.line_347.id_12
Name markdup_bam_picard ctl1
Thread thread_42
PID 6275
OK false
Exit Code 3
Retries
State KILLED
Dep. ERROR
Cpus 1
Mem
Start 2016-11-09 18:48:56
End 2016-11-09 20:17:08
Elapsed 01:28:12
Timeout 00:00:-1
Wall Timeout 100 days
Input files /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370912_1/out/align/ctl1/ss_50M_2000_GRCh38.filt.bam
Output files /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370912_1/out/align/ctl1/ss_50M_2000_GRCh38.dupmark.bam /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370912_1/out/qc/ctl1/ss_50M_2000_GRCh38.dup.qc
Dependencies 
 
# SYS command. line 349

 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 353

 export _JAVA_OPTIONS="-Xms256M -Xmx4G -XX:ParallelGCThreads=1"

# SYS command. line 354

 export MAX_JAVA_MEM="12G"

# SYS command. line 359

 if [ -f "${PICARDROOT}/picard.jar" ]; then \
		      MARKDUP="${PICARDROOT}/picard.jar MarkDuplicates"; \
		    elif [ -f "${PICARDROOT}/MarkDuplicates.jar" ]; then \
		      MARKDUP="${PICARDROOT}/MarkDuplicates.jar"; \
		    elif [ $(which picard 2> /dev/null | wc -l || echo) == "1" ]; then \
		      MARKDUP="$(ls $(dirname $(which picard))/../share/picard*/MarkDuplicates.jar 2> /dev/null || echo)"; \
		    fi

# SYS command. line 368

 if [ -f "${MARKDUP}" ]; then \
			java -Xmx4G -jar ${MARKDUP} \
				INPUT="/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370912_1/out/align/ctl1/ss_50M_2000_GRCh38.filt.bam" OUTPUT="/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370912_1/out/align/ctl1/ss_50M_2000_GRCh38.dupmark.bam" \
				METRICS_FILE="/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370912_1/out/qc/ctl1/ss_50M_2000_GRCh38.dup.qc" VALIDATION_STRINGENCY=LENIENT \
				ASSUME_SORTED=true REMOVE_DUPLICATES=false; \
		    else \
			picard MarkDuplicates \
				INPUT="/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370912_1/out/align/ctl1/ss_50M_2000_GRCh38.filt.bam" OUTPUT="/srv/scratch/shared/