Script file	/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/src/TF_chipseq_pipeline/chipseq.bds
Program ID	chipseq.bds.20161109_175125_591
Start time	2016-11-09 17:51:25
Run time	03:32:40.956
Tasks executed	4
Tasks failed	0
Tasks failed names
Arguments*	[-title, SRR1370890_1, -nth, 5, -macs2_for_pooled_rep_only, -species, hg38, -url_base, http://mitra.stanford.edu/kundaje/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370890_1/out, -fastq, /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/RawReadsFiles/SRR1370890_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

Thread ID	Parent	Tasks
thread_Root	None
thread_42	thread_Root	chipseq.bds.20161109_175125_591_parallel_42/task.postalign_bam.dedup_bam_1_ctl1.line_95.id_10 chipseq.bds.20161109_175125_591_parallel_42/task.postalign_bam.markdup_bam_picard_ctl1.line_347.id_12
thread_41	thread_Root	chipseq.bds.20161109_175125_591_parallel_41/task.align_bwa.bwa_aln_rep1.line_108.id_11 chipseq.bds.20161109_175125_591_parallel_41/task.align_bwa.bwa_sam_rep1.line_140.id_13

Task details

Task

Execution

Time

Dependencies

Task program, Errors, StdOut / StdErr

Num	1
ID	task.postalign_bam.dedup_bam_1_ctl1.line_95.id_10
Name	dedup_bam_1 ctl1
Thread	thread_42
PID	5970

OK	true
Exit Code	0
Retries
State	FINISHED
Dep.	OK
Cpus	4
Mem

Start	2016-11-09 17:51:29
End	2016-11-09 18:33:47
Elapsed	00:42: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/SRR1370890_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/SRR1370890_1/out/align/ctl1/ss_50M_2000_GRCh38.qnmsrt.bam; \
				samtools view -h /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370890_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/SRR1370890_1/out/align/ctl1/ss_50M_2000_GRCh38.filt.bam; \
				rm -f /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370890_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/SRR1370890_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--------------------
52287 (process ID) old priority 0, new priority 10

Num	2
ID	task.align_bwa.bwa_aln_rep1.line_108.id_11
Name	bwa_aln rep1
Thread	thread_41
PID	5979

OK	true
Exit Code	0
Retries
State	FINISHED
Dep.	OK
Cpus
Mem

Start	2016-11-09 17:51:30
End	2016-11-09 20:52:02
Elapsed	03:00:32
Timeout	00:00:-1
Wall Timeout	100 days

Input files	/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/RawReadsFiles/SRR1370890_1.fastq.gz
Output files	/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370890_1/out/align/rep1/SRR1370890_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/SRR1370890_1.fastq.gz > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370890_1/out/align/rep1/SRR1370890_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--------------------
23753 (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.8% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 40.25 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 262144 sequences have been processed.
[bwa_read_seq] 0.8% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 38.76 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 524288 sequences have been processed.
[bwa_read_seq] 0.8% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 39.71 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 786432 sequences have been processed.
[bwa_read_seq] 0.9% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 39.16 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 1048576 sequences have been processed.
[bwa_read_seq] 0.9% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 39.29 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 1310720 sequences have been processed.
[bwa_read_seq] 1.0% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 39.09 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 1572864 sequences have been processed.
[bwa_read_seq] 1.1% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 40.81 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 1835008 sequences have been processed.
[bwa_read_seq] 1.0% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 41.29 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 2097152 sequences have been processed.
[bwa_read_seq] 0.7% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 39.18 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 2359296 sequences have been processed.
[bwa_read_seq] 0.7% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 39.09 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 2621440 sequences have been processed.
[bwa_read_seq] 0.7% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 39.18 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 2883584 sequences have been processed.
[bwa_read_seq] 0.8% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 39.37 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 3145728 sequences have been processed.
[bwa_read_seq] 0.9% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 38.55 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 3407872 sequences have been processed.
[bwa_read_seq] 1.0% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 39.53 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 3670016 sequences have been processed.
[bwa_read_seq] 1.1% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 40.21 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 3932160 sequences have been processed.
[bwa_read_seq] 0.7% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 39.31 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 4194304 sequences have been processed.
[bwa_read_seq] 0.6% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 39.07 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 4456448 sequences have been processed.
[bwa_read_seq] 0.6% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 38.54 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 4718592 sequences have been processed.
[bwa_read_seq] 0.6% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 36.58 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 4980736 sequences have been processed.
[bwa_read_seq] 0.6% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 39.08 sec
[bwa_aln_core] write to the disk... 0.03 sec
[bwa_aln_core] 5242880 sequences have been processed.
[bwa_read_seq] 0.6% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 40.14 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 5505024 sequences have been processed.
[bwa_read_seq] 0.7% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 41.02 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 5767168 sequences have been processed.
[bwa_read_seq] 0.7% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 40.59 sec
[bwa_aln_core] write to the disk... 0.03 sec
[bwa_aln_core] 6029312 sequences have been processed.
[bwa_read_seq] 0.9% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 40.96 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 6291456 sequences have been processed.
[bwa_read_seq] 0.9% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 41.40 sec
[bwa_aln_core] write to the disk... 0.03 sec
[bwa_aln_core] 6553600 sequences have been processed.
[bwa_read_seq] 0.5% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 40.58 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 6815744 sequences have been processed.
[bwa_read_seq] 0.5% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 40.35 sec
[bwa_aln_core] write to the disk... 0.03 sec
[bwa_aln_core] 7077888 sequences have been processed.
[bwa_read_seq] 0.6% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 41.33 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 7340032 sequences have been processed.
[bwa_read_seq] 0.6% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 41.37 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 7602176 sequences have been processed.
[bwa_read_seq] 0.6% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 40.53 sec
[bwa_aln_core] write to the disk... 0.03 sec
[bwa_aln_core] 7864320 sequences have been processed.
[bwa_read_seq] 0.6% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 40.54 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 8126464 sequences have been processed.
[bwa_read_seq] 0.7% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 40.48 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 8388608 sequences have been processed.
[bwa_read_seq] 0.8% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 38.83 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 8650752 sequences have been processed.
[bwa_read_seq] 0.9% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 39.39 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 8912896 sequences have been processed.
[bwa_read_seq] 0.7% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 39.70 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 9175040 sequences have been processed.
[bwa_read_seq] 0.7% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 39.61 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 9437184 sequences have been processed.
[bwa_read_seq] 0.7% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 38.93 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 9699328 sequences have been processed.
[bwa_read_seq] 0.8% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 39.85 sec
[bwa_aln_core] write to the disk... 0.03 sec
[bwa_aln_core] 9961472 sequences have been processed.
[bwa_read_seq] 0.8% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 39.47 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 10223616 sequences have been processed.
[bwa_read_seq] 0.9% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 40.15 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 10485760 sequences have been processed.
[bwa_read_seq] 1.0% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 39.82 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 10747904 sequences have been processed.
[bwa_read_seq] 1.1% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 40.66 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 11010048 sequences have been processed.
[bwa_read_seq] 0.6% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 39.02 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 11272192 sequences have been processed.
[bwa_read_seq] 0.7% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 39.08 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 11534336 sequences have been processed.
[bwa_read_seq] 0.7% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 38.65 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 11796480 sequences have been processed.
[bwa_read_seq] 0.8% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 39.21 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 12058624 sequences have been processed.
[bwa_read_seq] 0.9% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 39.49 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 12320768 sequences have been processed.
[bwa_read_seq] 1.0% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 40.00 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 12582912 sequences have been processed.
[bwa_read_seq] 1.2% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 40.89 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 12845056 sequences have been processed.
[bwa_read_seq] 1.2% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 40.18 sec
[bwa_aln_core] write to the disk... 0.03 sec
[bwa_aln_core] 13107200 sequences have been processed.
[bwa_read_seq] 0.9% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 42.02 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 13369344 sequences have been processed.
[bwa_read_seq] 0.5% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 40.14 sec
[bwa_aln_core] write to the disk... 0.03 sec
[bwa_aln_core] 13631488 sequences have been processed.
[bwa_read_seq] 0.5% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 40.18 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 13893632 sequences have been processed.
[bwa_read_seq] 0.6% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 40.39 sec
[bwa_aln_core] write to the disk... 0.03 sec
[bwa_aln_core] 14155776 sequences have been processed.
[bwa_read_seq] 0.6% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 40.04 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 14417920 sequences have been processed.
[bwa_read_seq] 0.6% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 40.55 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 14680064 sequences have been processed.
[bwa_read_seq] 0.7% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 40.68 sec
[bwa_aln_core] write to the disk... 0.03 sec
[bwa_aln_core] 14942208 sequences have been processed.
[bwa_read_seq] 0.7% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 41.16 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 15204352 sequences have been processed.
[bwa_read_seq] 0.8% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 41.25 sec
[bwa_aln_core] write to the disk... 0.03 sec
[bwa_aln_core] 15466496 sequences have been processed.
[bwa_read_seq] 0.9% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 41.31 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 15728640 sequences have been processed.
[bwa_read_seq] 0.6% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 40.39 sec
[bwa_aln_core] write to the disk... 0.03 sec
[bwa_aln_core] 15990784 sequences have been processed.
[bwa_read_seq] 0.5% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 39.82 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 16252928 sequences have been processed.
[bwa_read_seq] 0.5% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 39.02 sec
[bwa_aln_core] write to the disk... 0.03 sec
[bwa_aln_core] 16515072 sequences have been processed.
[bwa_read_seq] 0.5% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 39.09 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 16777216 sequences have been processed.
[bwa_read_seq] 0.6% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 39.44 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 17039360 sequences have been processed.
[bwa_read_seq] 0.6% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 39.46 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 17301504 sequences have been processed.
[bwa_read_seq] 0.7% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 40.00 sec
[bwa_aln_core] write to the disk... 0.03 sec
[bwa_aln_core] 17563648 sequences have been processed.
[bwa_read_seq] 0.7% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 38.90 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 17825792 sequences have been processed.
[bwa_read_seq] 0.8% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 40.22 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 18087936 sequences have been processed.
[bwa_read_seq] 0.9% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 39.33 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 18350080 sequences have been processed.
[bwa_read_seq] 1.0% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 34.77 sec
[bwa_aln_core] write to the disk... 0.02 sec
[bwa_aln_core] 18572862 sequences have been processed.
[main] Version: 0.7.13-r1126
[main] CMD: 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/SRR1370890_1.fastq.gz
[main] Real time: 2866.992 sec; CPU: 2864.888 sec

Num	3
ID	task.postalign_bam.markdup_bam_picard_ctl1.line_347.id_12
Name	markdup_bam_picard ctl1
Thread	thread_42
PID	6267

OK	false
Exit Code	0
Retries
State	RUNNING
Dep.	WAIT
Cpus	1
Mem

Start	2016-11-09 18:33:48
End	2016-11-09 18:33:48
Elapsed	00:00:00
Timeout	00:00:-1
Wall Timeout	100 days

Input files	/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370890_1/out/align/ctl1/ss_50M_2000_GRCh38.filt.bam
Output files	/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370890_1/out/align/ctl1/ss_50M_2000_GRCh38.dupmark.bam /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370890_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/SRR1370890_1/out/align/ctl1/ss_50M_2000_GRCh38.filt.bam" OUTPUT="/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370890_1/out/align/ctl1/ss_50M_2000_GRCh38.dupmark.bam" \
				METRICS_FILE="/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370890_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/SRR1370890_1/out/align/ctl1/ss_50M_2000_GRCh38.filt.bam" OUTPUT="/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370890_1/out/align/ctl1/ss_50M_2000_GRCh38.dupmark.bam" \
				METRICS_FILE="/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370890_1/out/qc/ctl1/ss_50M_2000_GRCh38.dup.qc" VALIDATION_STRINGENCY=LENIENT \
				ASSUME_SORTED=true REMOVE_DUPLICATES=false; \
			fi

# SYS command. line 380

 TASKTIME=$[$(date +%s)-${STARTTIME}]; if [ ${TASKTIME} -lt 60 ]; then echo "Waiting for $[60-${TASKTIME}] seconds."; sleep $[60-${TASKTIME}]; fi

   
--------------------Stdout--------------------
62101 (process ID) old priority 0, new priority 10

 
--------------------Stderr--------------------
Picked up _JAVA_OPTIONS: -Xms256M -Xmx4G -XX:ParallelGCThreads=1
[Wed Nov 09 21:21:34 PST 2016] picard.sam.markduplicates.MarkDuplicates INPUT=[/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370890_1/out/align/ctl1/ss_50M_2000_GRCh38.filt.bam] OUTPUT=/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370890_1/out/align/ctl1/ss_50M_2000_GRCh38.dupmark.bam METRICS_FILE=/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370890_1/out/qc/ctl1/ss_50M_2000_GRCh38.dup.qc REMOVE_DUPLICATES=false ASSUME_SORTED=true VALIDATION_STRINGENCY=LENIENT    MAX_SEQUENCES_FOR_DISK_READ_ENDS_MAP=50000 MAX_FILE_HANDLES_FOR_READ_ENDS_MAP=8000 SORTING_COLLECTION_SIZE_RATIO=0.25 PROGRAM_RECORD_ID=MarkDuplicates PROGRAM_GROUP_NAME=MarkDuplicates DUPLICATE_SCORING_STRATEGY=SUM_OF_BASE_QUALITIES READ_NAME_REGEX=[a-zA-Z0-9]+:[0-9]:([0-9]+):([0-9]+):([0-9]+).* OPTICAL_DUPLICATE_PIXEL_DISTANCE=100 VERBOSITY=INFO QUIET=false COMPRESSION_LEVEL=5 MAX_RECORDS_IN_RAM=500000 CREATE_INDEX=false CREATE_MD5_FILE=false
[Wed Nov 09 21:21:34 PST 2016] Executing as imk1@wotan on Linux 4.4.0-45-generic amd64; OpenJDK 64-Bit Server VM 1.8.0_92-b15; Picard version: 1.126(4691ee611ac205d4afe2a1b7a2ea975a6f997426_1417447214) IntelDeflater
INFO	2016-11-09 21:21:34	MarkDuplicates	Start of doWork freeMemory: 254877736; totalMemory: 257425408; maxMemory: 3817865216
INFO	2016-11-09 21:21:34	MarkDuplicates	Reading input file and constructing read end information.
INFO	2016-11-09 21:21:34	MarkDuplicates	Will retain up to 14684096 data points before spilling to disk.
INFO	2016-11-09 21:21:49	MarkDuplicates	Read     1,000,000 records.  Elapsed time: 00:00:14s.  Time for last 1,000,000:   14s.  Last read position: chr1:22,298,551
INFO	2016-11-09 21:21:49	MarkDuplicates	Tracking 2 as yet unmatched pairs. 2 records in RAM.
INFO	2016-11-09 21:22:01	MarkDuplicates	Read     2,000,000 records.  Elapsed time: 00:00:25s.  Time for last 1,000,000:   11s.  Last read position: chr1:41,895,424
INFO	2016-11-09 21:22:01	MarkDuplicates	Tracking 8 as yet unmatched pairs. 8 records in RAM.
INFO	2016-11-09 21:22:12	MarkDuplicates	Read     3,000,000 records.  Elapsed time: 00:00:37s.  Time for last 1,000,000:   11s.  Last read position: chr1:63,396,527
INFO	2016-11-09 21:22:12	MarkDuplicates	Tracking 0 as yet unmatched pairs. 0 records in RAM.
INFO	2016-11-09 21:22:24	MarkDuplicates	Read     4,000,000 records.  Elapsed time: 00:00:48s.  Time for last 1,000,000:   11s.  Last read position: chr1:89,147,707
INFO	2016-11-09 21:22:24	MarkDuplicates	Tracking 4 as yet unmatched pairs. 4 records in RAM.
INFO	2016-11-09 21:22:41	MarkDuplicates	Read     5,000,000 records.  Elapsed time: 00:01:06s.  Time for last 1,000,000:   17s.  Last read position: chr1:115,615,553
INFO	2016-11-09 21:22:41	MarkDuplicates	Tracking 4 as yet unmatched pairs. 4 records in RAM.
INFO	2016-11-09 21:23:01	MarkDuplicates	Read     6,000,000 records.  Elapsed time: 00:01:26s.  Time for last 1,000,000:   19s.  Last read position: chr1:164,468,682
INFO	2016-11-09 21:23:01	MarkDuplicates	Tracking 6 as yet unmatched pairs. 6 records in RAM.
INFO	2016-11-09 21:23:15	MarkDuplicates	Read     7,000,000 records.  Elapsed time: 00:01:40s.  Time for last 1,000,000:   14s.  Last read position: chr1:189,102,252
INFO	2016-11-09 21:23:15	MarkDuplicates	Tracking 0 as yet unmatched pairs. 0 records in RAM.
INFO	2016-11-09 21:23:26	MarkDuplicates	Read     8,000,000 records.  Elapsed time: 00:01:51s.  Time for last 1,000,000:   11s.  Last read position: chr1:213,979,760
INFO	2016-11-09 21:23:26	MarkDuplicates	Tracking 4 as yet unmatched pairs. 4 records in RAM.
INFO	2016-11-09 21:23:38	MarkDuplicates	Read     9,000,000 records.  Elapsed time: 00:02:03s.  Time for last 1,000,000:   11s.  Last read position: chr1:231,814,268
INFO	2016-11-09 21:23:38	MarkDuplicates	Tracking 2 as yet unmatched pairs. 2 records in RAM.

Num	4
ID	task.align_bwa.bwa_sam_rep1.line_140.id_13
Name	bwa_sam rep1
Thread	thread_41
PID	6406

OK	false
Exit Code	0
Retries
State	RUNNING
Dep.	WAIT
Cpus	1
Mem

Start	2016-11-09 20:52:03
End	2016-11-09 20:52:03
Elapsed	00:00:00
Timeout	00:00:-1
Wall Timeout	100 days

Input files	/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/RawReadsFiles/SRR1370890_1.fastq.gz /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370890_1/out/align/rep1/SRR1370890_1.sai
Output files	/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370890_1/out/align/rep1/SRR1370890_1.bam /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370890_1/out/qc/rep1/SRR1370890_1.flagstat.qc
Dependencies

 
# SYS command. line 142

 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 147

 bwa samse /mnt/data/bds_pipeline_genome_data/hg38/bwa_index/GRCh38.p3.genome.fa /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370890_1/out/align/rep1/SRR1370890_1.sai /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/RawReadsFiles/SRR1370890_1.fastq.gz | samtools view -Su /dev/stdin \
			| sambamba sort -t 1 /dev/stdin -o /srv/scratch/shared/surya/imk1/TFBindingPrediction

BigDataScript report: chipseq.bds

Timeline

Task Graph

Parallel (threads) details

Task details