BigDataScript report: chipseq.bds

Script file /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/src/chipseq_pipeline/chipseq.bds
Program ID chipseq.bds.20161113_130805_528
Start time 2016-11-13 13:08:05
Run time 07:31:54.874
Tasks executed 10
Tasks failed 1
Tasks failed names 
xcor rep1
Arguments* [-title, SRR1370897_1, -nth, 5, -macs2_for_pooled_rep_only, -species, hg38, -url_base, http://mitra.stanford.edu/kundaje/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out, -fastq, /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/RawReadsFiles/SRR1370897_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_42 thread_Root 
 chipseq.bds.20161113_130805_528_parallel_42/task.postalign_bam.dedup_bam_2_ctl1.line_144.id_11
chipseq.bds.20161113_130805_528_parallel_42/task.postalign_bam.bam_to_tag_ctl1.line_505.id_13
thread_41 thread_Root 
 chipseq.bds.20161113_130805_528_parallel_41/task.align_bwa.bwa_sam_rep1.line_140.id_10
chipseq.bds.20161113_130805_528_parallel_41/task.postalign_bam.dedup_bam_1_rep1.line_95.id_12
chipseq.bds.20161113_130805_528_parallel_41/task.postalign_bam.markdup_bam_picard_rep1.line_348.id_14
chipseq.bds.20161113_130805_528_parallel_41/task.postalign_bam.dedup_bam_2_rep1.line_144.id_15
chipseq.bds.20161113_130805_528_parallel_41/task.postalign_bam.bam_to_tag_rep1.line_505.id_16
chipseq.bds.20161113_130805_528_parallel_41/task.postalign_bed.subsample_tag_rep1.line_45.id_17
chipseq.bds.20161113_130805_528_parallel_41/task.postalign_bed.spr_rep1.line_255.id_18
chipseq.bds.20161113_130805_528_parallel_41/task.postalign_bed.xcor_rep1.line_212.id_19

Task details

Task Execution Time Dependencies Task program, Errors, StdOut / StdErr
Num 1
ID task.align_bwa.bwa_sam_rep1.line_140.id_10
Name bwa_sam rep1
Thread thread_41
PID 9030
OK true
Exit Code 0
Retries
State FINISHED
Dep. OK
Cpus 1
Mem
Start 2016-11-13 13:08:08
End 2016-11-13 13:35:25
Elapsed 00:27:17
Timeout 00:00:-1
Wall Timeout 100 days
Input files /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/RawReadsFiles/SRR1370897_1.fastq.gz /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/align/rep1/SRR1370897_1.sai
Output files /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/align/rep1/SRR1370897_1.bam /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/qc/rep1/SRR1370897_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; 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 147

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

# SYS command. line 149

 sambamba flagstat -t 1 /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/align/rep1/SRR1370897_1.bam > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/qc/rep1/SRR1370897_1.flagstat.qc

# SYS command. line 151

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

 
--------------------Stderr--------------------
[bwa_read_seq] 0.6% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 4.25 sec
[bwa_aln_core] refine gapped alignments... 0.55 sec
[bwa_aln_core] print alignments... [samopen] SAM header is present: 493 sequences.
0.65 sec
[bwa_aln_core] 262144 sequences have been processed.
[bwa_read_seq] 0.5% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 4.01 sec
[bwa_aln_core] refine gapped alignments... 0.51 sec
[bwa_aln_core] print alignments... 0.39 sec
[bwa_aln_core] 524288 sequences have been processed.
[bwa_read_seq] 0.5% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 3.99 sec
[bwa_aln_core] refine gapped alignments... 0.51 sec
[bwa_aln_core] print alignments... 0.38 sec
[bwa_aln_core] 786432 sequences have been processed.
[bwa_read_seq] 0.5% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 3.99 sec
[bwa_aln_core] refine gapped alignments... 0.50 sec
[bwa_aln_core] print alignments... 0.65 sec
[bwa_aln_core] 1048576 sequences have been processed.
[bwa_read_seq] 0.5% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 3.76 sec
[bwa_aln_core] refine gapped alignments... 0.51 sec
[bwa_aln_core] print alignments... 0.57 sec
[bwa_aln_core] 1310720 sequences have been processed.
[bwa_read_seq] 0.8% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 3.77 sec
[bwa_aln_core] refine gapped alignments... 0.52 sec
[bwa_aln_core] print alignments... 0.55 sec
[bwa_aln_core] 1572864 sequences have been processed.
[bwa_read_seq] 0.9% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 3.61 sec
[bwa_aln_core] refine gapped alignments... 0.50 sec
[bwa_aln_core] print alignments... 0.58 sec
[bwa_aln_core] 1835008 sequences have been processed.
[bwa_read_seq] 0.8% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 3.68 sec
[bwa_aln_core] refine gapped alignments... 0.52 sec
[bwa_aln_core] print alignments... 0.54 sec
[bwa_aln_core] 2097152 sequences have been processed.
[bwa_read_seq] 0.7% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 4.03 sec
[bwa_aln_core] refine gapped alignments... 0.52 sec
[bwa_aln_core] print alignments... 0.55 sec
[bwa_aln_core] 2359296 sequences have been processed.
[bwa_read_seq] 0.6% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 3.79 sec
[bwa_aln_core] refine gapped alignments... 0.50 sec
[bwa_aln_core] print alignments... 0.42 sec
[bwa_aln_core] 2621440 sequences have been processed.
[bwa_read_seq] 0.5% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 3.77 sec
[bwa_aln_core] refine gapped alignments... 0.51 sec
[bwa_aln_core] print alignments... 0.54 sec
[bwa_aln_core] 2883584 sequences have been processed.
[bwa_read_seq] 0.5% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 3.63 sec
[bwa_aln_core] refine gapped alignments... 0.50 sec
[bwa_aln_core] print alignments... 0.47 sec
[bwa_aln_core] 3145728 sequences have been processed.
[bwa_read_seq] 0.6% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 3.93 sec
[bwa_aln_core] refine gapped alignments... 0.52 sec
[bwa_aln_core] print alignments... 0.51 sec
[bwa_aln_core] 3407872 sequences have been processed.
[bwa_read_seq] 0.6% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 3.75 sec
[bwa_aln_core] refine gapped alignments... 0.50 sec
[bwa_aln_core] print alignments... 0.40 sec
[bwa_aln_core] 3670016 sequences have been processed.
[bwa_read_seq] 0.6% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 3.82 sec
[bwa_aln_core] refine gapped alignments... 0.50 sec
[bwa_aln_core] print alignments... 0.54 sec
[bwa_aln_core] 3932160 sequences have been processed.
[bwa_read_seq] 0.9% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 3.54 sec
[bwa_aln_core] refine gapped alignments... 0.56 sec
[bwa_aln_core] print alignments... 0.52 sec
[bwa_aln_core] 4194304 sequences have been processed.
[bwa_read_seq] 0.6% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 5.24 sec
[bwa_aln_core] refine gapped alignments... 0.50 sec
[bwa_aln_core] print alignments... 0.39 sec
[bwa_aln_core] 4456448 sequences have been processed.
[bwa_read_seq] 0.5% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 4.16 sec
[bwa_aln_core] refine gapped alignments... 0.51 sec
[bwa_aln_core] print alignments... 0.49 sec
[bwa_aln_core] 4718592 sequences have been processed.
[bwa_read_seq] 0.4% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 3.65 sec
[bwa_aln_core] refine gapped alignments... 0.50 sec
[bwa_aln_core] print alignments... 0.51 sec
[bwa_aln_core] 4980736 sequences have been processed.
[bwa_read_seq] 0.4% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 4.13 sec
[bwa_aln_core] refine gapped alignments... 0.50 sec
[bwa_aln_core] print alignments... 0.56 sec
[bwa_aln_core] 5242880 sequences have been processed.
[bwa_read_seq] 0.5% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 3.72 sec
[bwa_aln_core] refine gapped alignments... 0.50 sec
[bwa_aln_core] print alignments... 0.47 sec
[bwa_aln_core] 5505024 sequences have been processed.
[bwa_read_seq] 0.8% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 4.18 sec
[bwa_aln_core] refine gapped alignments... 0.53 sec
[bwa_aln_core] print alignments... 0.49 sec
[bwa_aln_core] 5767168 sequences have been processed.
[bwa_read_seq] 1.0% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 3.58 sec
[bwa_aln_core] refine gapped alignments... 0.51 sec
[bwa_aln_core] print alignments... 0.54 sec
[bwa_aln_core] 6029312 sequences have been processed.
[bwa_read_seq] 1.0% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 4.04 sec
[bwa_aln_core] refine gapped alignments... 0.52 sec
[bwa_aln_core] print alignments... 0.62 sec
[bwa_aln_core] 6291456 sequences have been processed.
[bwa_read_seq] 0.9% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 3.57 sec
[bwa_aln_core] refine gapped alignments... 0.51 sec
[bwa_aln_core] print alignments... 0.48 sec
[bwa_aln_core] 6553600 sequences have been processed.
[bwa_read_seq] 0.8% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 4.09 sec
[bwa_aln_core] refine gapped alignments... 0.53 sec
[bwa_aln_core] print alignments... 0.43 sec
[bwa_aln_core] 6815744 sequences have been processed.
[bwa_read_seq] 0.7% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 4.07 sec
[bwa_aln_core] refine gapped alignments... 0.52 sec
[bwa_aln_core] print alignments... 0.49 sec
[bwa_aln_core] 7077888 sequences have been processed.
[bwa_read_seq] 0.7% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 3.65 sec
[bwa_aln_core] refine gapped alignments... 0.51 sec
[bwa_aln_core] print alignments... 0.54 sec
[bwa_aln_core] 7340032 sequences have been processed.
[bwa_read_seq] 0.8% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 4.86 sec
[bwa_aln_core] refine gapped alignments... 0.50 sec
[bwa_aln_core] print alignments... 0.48 sec
[bwa_aln_core] 7602176 sequences have been processed.
[bwa_read_seq] 1.2% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 5.68 sec
[bwa_aln_core] refine gapped alignments... 0.51 sec
[bwa_aln_core] print alignments... 0.56 sec
[bwa_aln_core] 7864320 sequences have been processed.
[bwa_read_seq] 1.3% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 4.84 sec
[bwa_aln_core] refine gapped alignments... 0.51 sec
[bwa_aln_core] print alignments... 0.52 sec
[bwa_aln_core] 8126464 sequences have been processed.
[bwa_read_seq] 0.9% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 4.23 sec
[bwa_aln_core] refine gapped alignments... 0.50 sec
[bwa_aln_core] print alignments... 0.59 sec
[bwa_aln_core] 8388608 sequences have been processed.
[bwa_read_seq] 0.8% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 4.75 sec
[bwa_aln_core] refine gapped alignments... 0.50 sec
[bwa_aln_core] print alignments... 0.56 sec
[bwa_aln_core] 8650752 sequences have been processed.
[bwa_read_seq] 0.6% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 4.22 sec
[bwa_aln_core] refine gapped alignments... 0.51 sec
[bwa_aln_core] print alignments... 0.49 sec
[bwa_aln_core] 8912896 sequences have been processed.
[bwa_read_seq] 0.6% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 3.60 sec
[bwa_aln_core] refine gapped alignments... 0.50 sec
[bwa_aln_core] print alignments... 0.51 sec
[bwa_aln_core] 9175040 sequences have been processed.
[bwa_read_seq] 0.8% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 4.27 sec
[bwa_aln_core] refine gapped alignments... 0.52 sec
[bwa_aln_core] print alignments... 0.39 sec
[bwa_aln_core] 9437184 sequences have been processed.
[bwa_read_seq] 0.7% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 4.11 sec
[bwa_aln_core] refine gapped alignments... 0.52 sec
[bwa_aln_core] print alignments... 0.49 sec
[bwa_aln_core] 9699328 sequences have been processed.
[bwa_read_seq] 0.8% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 4.97 sec
[bwa_aln_core] refine gapped alignments... 0.51 sec
[bwa_aln_core] print alignments... 0.61 sec
[bwa_aln_core] 9961472 sequences have been processed.
[bwa_read_seq] 1.0% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 7.05 sec
[bwa_aln_core] refine gapped alignments... 0.54 sec
[bwa_aln_core] print alignments... 0.53 sec
[bwa_aln_core] 10223616 sequences have been processed.
[bwa_read_seq] 0.8% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 4.41 sec
[bwa_aln_core] refine gapped alignments... 0.50 sec
[bwa_aln_core] print alignments... 0.61 sec
[bwa_aln_core] 10485760 sequences have been processed.
[bwa_read_seq] 0.8% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 3.67 sec
[bwa_aln_core] refine gapped alignments... 0.51 sec
[bwa_aln_core] print alignments... 0.34 sec
[bwa_aln_core] 10747904 sequences have been processed.
[bwa_read_seq] 0.6% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 4.08 sec
[bwa_aln_core] refine gapped alignments... 0.50 sec
[bwa_aln_core] print alignments... 0.60 sec
[bwa_aln_core] 11010048 sequences have been processed.
[bwa_read_seq] 0.7% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 12.86 sec
[bwa_aln_core] refine gapped alignments... 0.53 sec
[bwa_aln_core] print alignments... 0.53 sec
[bwa_aln_core] 11272192 sequences have been processed.
[bwa_read_seq] 0.7% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 15.49 sec
[bwa_aln_core] refine gapped alignments... 0.52 sec
[bwa_aln_core] print alignments... 0.58 sec
[bwa_aln_core] 11534336 sequences have been processed.
[bwa_read_seq] 1.0% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 18.39 sec
[bwa_aln_core] refine gapped alignments... 0.78 sec
[bwa_aln_core] print alignments... 0.49 sec
[bwa_aln_core] 11796480 sequences have been processed.
[bwa_read_seq] 1.7% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 20.66 sec
[bwa_aln_core] refine gapped alignments... 0.52 sec
[bwa_aln_core] print alignments... 0.52 sec
[bwa_aln_core] 12058624 sequences have been processed.
[bwa_read_seq] 1.1% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 20.01 sec
[bwa_aln_core] refine gapped alignments... 0.52 sec
[bwa_aln_core] print alignments... 0.44 sec
[bwa_aln_core] 12320768 sequences have been processed.
[bwa_read_seq] 0.6% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 21.69 sec
[bwa_aln_core] refine gapped alignments... 0.50 sec
[bwa_aln_core] print alignments... 0.62 sec
[bwa_aln_core] 12582912 sequences have been processed.
[bwa_read_seq] 0.5% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 22.96 sec
[bwa_aln_core] refine gapped alignments... 0.50 sec
[bwa_aln_core] print alignments... 0.69 sec
[bwa_aln_core] 12845056 sequences have been processed.
[bwa_read_seq] 0.5% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 29.73 sec
[bwa_aln_core] refine gapped alignments... 1.84 sec
[bwa_aln_core] print alignments... 0.58 sec
[bwa_aln_core] 13107200 sequences have been processed.
[bwa_read_seq] 0.5% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 28.24 sec
[bwa_aln_core] refine gapped alignments... 0.50 sec
[bwa_aln_core] print alignments... 0.40 sec
[bwa_aln_core] 13369344 sequences have been processed.
[bwa_read_seq] 0.5% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 12.43 sec
[bwa_aln_core] refine gapped alignments... 0.49 sec
[bwa_aln_core] print alignments... 0.58 sec
[bwa_aln_core] 13631488 sequences have been processed.
[bwa_read_seq] 0.7% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 21.76 sec
[bwa_aln_core] refine gapped alignments... 0.54 sec
[bwa_aln_core] print alignments... 0.44 sec
[bwa_aln_core] 13893632 sequences have been processed.
[bwa_read_seq] 0.5% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 14.24 sec
[bwa_aln_core] refine gapped alignments... 0.50 sec
[bwa_aln_core] print alignments... 0.53 sec
[bwa_aln_core] 14155776 sequences have been processed.
[bwa_read_seq] 0.8% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 10.32 sec
[bwa_aln_core] refine gapped alignments... 0.49 sec
[bwa_aln_core] print alignments... 0.54 sec
[bwa_aln_core] 14417920 sequences have been processed.
[bwa_read_seq] 0.7% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 23.55 sec
[bwa_aln_core] refine gapped alignments... 0.50 sec
[bwa_aln_core] print alignments... 0.54 sec
[bwa_aln_core] 14680064 sequences have been processed.
[bwa_read_seq] 0.6% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 28.32 sec
[bwa_aln_core] refine gapped alignments... 0.50 sec
[bwa_aln_core] print alignments... 0.61 sec
[bwa_aln_core] 14942208 sequences have been processed.
[bwa_read_seq] 0.5% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 26.01 sec
[bwa_aln_core] refine gapped alignments... 0.50 sec
[bwa_aln_core] print alignments... 0.41 sec
[bwa_aln_core] 15204352 sequences have been processed.
[bwa_read_seq] 0.5% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 18.75 sec
[bwa_aln_core] refine gapped alignments... 0.50 sec
[bwa_aln_core] print alignments... 0.49 sec
[bwa_aln_core] 15466496 sequences have been processed.
[bwa_read_seq] 0.5% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 18.43 sec
[bwa_aln_core] refine gapped alignments... 0.50 sec
[bwa_aln_core] print alignments... 0.49 sec
[bwa_aln_core] 15728640 sequences have been processed.
[bwa_read_seq] 0.6% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 23.50 sec
[bwa_aln_core] refine gapped alignments... 0.50 sec
[bwa_aln_core] print alignments... 0.48 sec
[bwa_aln_core] 15990784 sequences have been processed.
[bwa_read_seq] 0.6% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 16.44 sec
[bwa_aln_core] refine gapped alignments... 0.50 sec
[bwa_aln_core] print alignments... 0.45 sec
[bwa_aln_core] 16252928 sequences have been processed.
[bwa_read_seq] 0.7% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 19.80 sec
[bwa_aln_core] refine gapped alignments... 2.08 sec
[bwa_aln_core] print alignments... 0.44 sec
[bwa_aln_core] 16515072 sequences have been processed.
[bwa_read_seq] 0.6% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 24.37 sec
[bwa_aln_core] refine gapped alignments... 0.55 sec
[bwa_aln_core] print alignments... 0.56 sec
[bwa_aln_core] 16777216 sequences have been processed.
[bwa_read_seq] 0.6% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 25.05 sec
[bwa_aln_core] refine gapped alignments... 0.50 sec
[bwa_aln_core] print alignments... 0.47 sec
[bwa_aln_core] 17039360 sequences have been processed.
[bwa_read_seq] 0.4% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 19.32 sec
[bwa_aln_core] refine gapped alignments... 0.50 sec
[bwa_aln_core] print alignments... 0.34 sec
[bwa_aln_core] 17301504 sequences have been processed.
[bwa_read_seq] 0.5% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 12.96 sec
[bwa_aln_core] refine gapped alignments... 0.50 sec
[bwa_aln_core] print alignments... 0.34 sec
[bwa_aln_core] 17563648 sequences have been processed.
[bwa_read_seq] 0.6% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 9.12 sec
[bwa_aln_core] refine gapped alignments... 0.50 sec
[bwa_aln_core] print alignments... 0.47 sec
[bwa_aln_core] 17825792 sequences have been processed.
[bwa_read_seq] 0.8% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 13.74 sec
[bwa_aln_core] refine gapped alignments... 0.49 sec
[bwa_aln_core] print alignments... 0.51 sec
[bwa_aln_core] 18087936 sequences have been processed.
[bwa_read_seq] 0.7% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 14.35 sec
[bwa_aln_core] refine gapped alignments... 0.51 sec
[bwa_aln_core] print alignments... 0.45 sec
[bwa_aln_core] 18350080 sequences have been processed.
[bwa_read_seq] 1.0% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 8.64 sec
[bwa_aln_core] refine gapped alignments... 0.49 sec
[bwa_aln_core] print alignments... 0.51 sec
[bwa_aln_core] 18612224 sequences have been processed.
[bwa_read_seq] 0.9% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 4.25 sec
[bwa_aln_core] refine gapped alignments... 0.50 sec
[bwa_aln_core] print alignments... 0.55 sec
[bwa_aln_core] 18874368 sequences have been processed.
[bwa_read_seq] 0.8% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 13.19 sec
[bwa_aln_core] refine gapped alignments... 0.49 sec
[bwa_aln_core] print alignments... 0.58 sec
[bwa_aln_core] 19136512 sequences have been processed.
[bwa_read_seq] 0.7% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 12.86 sec
[bwa_aln_core] refine gapped alignments... 0.50 sec
[bwa_aln_core] print alignments... 0.48 sec
[bwa_aln_core] 19398656 sequences have been processed.
[bwa_read_seq] 0.7% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 8.77 sec
[bwa_aln_core] refine gapped alignments... 0.51 sec
[bwa_aln_core] print alignments... 0.47 sec
[bwa_aln_core] 19660800 sequences have been processed.
[bwa_read_seq] 0.8% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 8.25 sec
[bwa_aln_core] refine gapped alignments... 0.51 sec
[bwa_aln_core] print alignments... 0.38 sec
[bwa_aln_core] 19922944 sequences have been processed.
[bwa_read_seq] 1.2% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 21.59 sec
[bwa_aln_core] refine gapped alignments... 0.50 sec
[bwa_aln_core] print alignments... 0.48 sec
[bwa_aln_core] 20185088 sequences have been processed.
[bwa_read_seq] 1.0% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 14.69 sec
[bwa_aln_core] refine gapped alignments... 0.56 sec
[bwa_aln_core] print alignments... 0.59 sec
[bwa_aln_core] 20447232 sequences have been processed.
[bwa_read_seq] 1.0% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 12.72 sec
[bwa_aln_core] refine gapped alignments... 0.50 sec
[bwa_aln_core] print alignments... 0.57 sec
[bwa_aln_core] 20709376 sequences have been processed.
[bwa_read_seq] 0.9% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 6.27 sec
[bwa_aln_core] refine gapped alignments... 0.54 sec
[bwa_aln_core] print alignments... 0.62 sec
[bwa_aln_core] 20971520 sequences have been processed.
[bwa_read_seq] 0.7% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 4.13 sec
[bwa_aln_core] refine gapped alignments... 0.50 sec
[bwa_aln_core] print alignments... 0.45 sec
[bwa_aln_core] 21233664 sequences have been processed.
[bwa_read_seq] 0.7% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 5.40 sec
[bwa_aln_core] refine gapped alignments... 0.50 sec
[bwa_aln_core] print alignments... 0.53 sec
[bwa_aln_core] 21495808 sequences have been processed.
[bwa_read_seq] 0.7% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 4.61 sec
[bwa_aln_core] refine gapped alignments... 0.50 sec
[bwa_aln_core] print alignments... 0.67 sec
[bwa_aln_core] 21757952 sequences have been processed.
[bwa_read_seq] 0.7% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 5.00 sec
[bwa_aln_core] refine gapped alignments... 0.50 sec
[bwa_aln_core] print alignments... 0.51 sec
[bwa_aln_core] 22020096 sequences have been processed.
[bwa_read_seq] 0.7% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 15.41 sec
[bwa_aln_core] refine gapped alignments... 0.51 sec
[bwa_aln_core] print alignments... 0.55 sec
[bwa_aln_core] 22282240 sequences have been processed.
[bwa_read_seq] 0.8% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 5.97 sec
[bwa_aln_core] refine gapped alignments... 0.50 sec
[bwa_aln_core] print alignments... 0.48 sec
[bwa_aln_core] 22544384 sequences have been processed.
[bwa_read_seq] 0.9% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 5.82 sec
[bwa_aln_core] refine gapped alignments... 0.55 sec
[bwa_aln_core] print alignments... 0.61 sec
[bwa_aln_core] 22806528 sequences have been processed.
[bwa_read_seq] 0.8% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 5.14 sec
[bwa_aln_core] refine gapped alignments... 0.50 sec
[bwa_aln_core] print alignments... 0.58 sec
[bwa_aln_core] 23068672 sequences have been processed.
[bwa_read_seq] 0.8% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 4.61 sec
[bwa_aln_core] refine gapped alignments... 0.50 sec
[bwa_aln_core] print alignments... 0.56 sec
[bwa_aln_core] 23330816 sequences have been processed.
[bwa_read_seq] 0.7% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 3.96 sec
[bwa_aln_core] refine gapped alignments... 0.51 sec
[bwa_aln_core] print alignments... 0.45 sec
[bwa_aln_core] 23592960 sequences have been processed.
[bwa_read_seq] 0.7% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 6.27 sec
[bwa_aln_core] refine gapped alignments... 0.50 sec
[bwa_aln_core] print alignments... 0.56 sec
[bwa_aln_core] 23855104 sequences have been processed.
[bwa_read_seq] 1.1% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 12.17 sec
[bwa_aln_core] refine gapped alignments... 0.51 sec
[bwa_aln_core] print alignments... 0.56 sec
[bwa_aln_core] 24117248 sequences have been processed.
[bwa_read_seq] 1.6% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 5.95 sec
[bwa_aln_core] refine gapped alignments... 0.52 sec
[bwa_aln_core] print alignments... 0.54 sec
[bwa_aln_core] 24379392 sequences have been processed.
[bwa_read_seq] 1.7% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 3.83 sec
[bwa_aln_core] refine gapped alignments... 0.47 sec
[bwa_aln_core] print alignments... 0.31 sec
[bwa_aln_core] 24555168 sequences have been processed.
[main] Version: 0.7.13-r1126
[main] CMD: bwa samse /mnt/data/bds_pipeline_genome_data/hg38/bwa_index/GRCh38.p3.genome.fa /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/align/rep1/SRR1370897_1.sai /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/RawReadsFiles/SRR1370897_1.fastq.gz
[main] Real time: 1198.265 sec; CPU: 1059.756 sec
Num 2
ID task.postalign_bam.dedup_bam_2_ctl1.line_144.id_11
Name dedup_bam_2 ctl1
Thread thread_42
PID 9031
OK true
Exit Code 0
Retries
State FINISHED
Dep. OK
Cpus 4
Mem
Start 2016-11-13 13:08:08
End 2016-11-13 13:51:33
Elapsed 00:43:24
Timeout 00:00:-1
Wall Timeout 100 days
Input files /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/align/ctl1/ss_50M_2000_GRCh38.dupmark.bam
Output files /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/align/ctl1/ss_50M_2000_GRCh38.nodup.bam /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/qc/ctl1/ss_50M_2000_GRCh38.nodup.flagstat.qc /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/qc/ctl1/ss_50M_2000_GRCh38.nodup.pbc.qc
Dependencies 
 
# SYS command. line 146

 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 150

 samtools view -F 1804 -b /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/align/ctl1/ss_50M_2000_GRCh38.dupmark.bam > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/align/ctl1/ss_50M_2000_GRCh38.nodup.bam

# SYS command. line 153

 sambamba index -t 4 /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/align/ctl1/ss_50M_2000_GRCh38.nodup.bam

# SYS command. line 155

 sambamba flagstat -t 4 /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/align/ctl1/ss_50M_2000_GRCh38.nodup.bam > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/qc/ctl1/ss_50M_2000_GRCh38.nodup.flagstat.qc

# SYS command. line 165

 bedtools bamtobed -i /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/align/ctl1/ss_50M_2000_GRCh38.dupmark.bam | \
				awk 'BEGIN{OFS="\t"}{print $1,$2,$3,$6}' | \
				grep -v 'chrM' | sort | uniq -c | \
				awk 'BEGIN{mt=0;m0=0;m1=0;m2=0} ($1==1){m1=m1+1} ($1==2){m2=m2+1} {m0=m0+1} {mt=mt+$1} END{m1_m2=-1.0; if(m2>0) m1_m2=m1/m2; printf "%d\t%d\t%d\t%d\t%f\t%f\t%f\n",mt,m0,m1,m2,m0/mt,m1/m0,m1_m2}' > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/qc/ctl1/ss_50M_2000_GRCh38.nodup.pbc.qc

# SYS command. line 169

 TASKTIME=$[$(date +%s)-${STARTTIME}]; if [ ${TASKTIME} -lt 60 ]; then echo "Waiting for $[60-${TASKTIME}] seconds."; sleep $[60-${TASKTIME}]; fi
 
   
--------------------Stdout--------------------
10247 (process ID) old priority 0, new priority 10
Num 3
ID task.postalign_bam.dedup_bam_1_rep1.line_95.id_12
Name dedup_bam_1 rep1
Thread thread_41
PID 9316
OK true
Exit Code 0
Retries
State FINISHED
Dep. OK
Cpus
Mem
Start 2016-11-13 13:35:27
End 2016-11-13 15:40:31
Elapsed 02:05:03
Timeout 00:00:-1
Wall Timeout 100 days
Input files /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/align/rep1/SRR1370897_1.bam
Output files /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/align/rep1/SRR1370897_1.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; 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 105

 if [[ 0 > 0 ]]; then \
				sambamba sort -t 1 /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/align/rep1/SRR1370897_1.bam -n -o /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/align/rep1/SRR1370897_1.qnmsrt.bam; \
				samtools view -h /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/align/rep1/SRR1370897_1.qnmsrt.bam | $(which assign_multimappers.py) -k 0 | \
				samtools view -F 1804 -Su /dev/stdin | \
				sambamba sort -t 1 /dev/stdin -o /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/align/rep1/SRR1370897_1.filt.bam; \
				rm -f /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/align/rep1/SRR1370897_1.qnmsrt.bam; \
			    else \
			    	samtools view -F 1804 -q 30 -u /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/align/rep1/SRR1370897_1.bam | \
			    	sambamba sort -t 1 /dev/stdin -o /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/align/rep1/SRR1370897_1.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--------------------
44090 (process ID) old priority 0, new priority 10
Num 4
ID task.postalign_bam.bam_to_tag_ctl1.line_505.id_13
Name bam_to_tag ctl1
Thread thread_42
PID 9335
OK true
Exit Code 0
Retries
State FINISHED
Dep. OK
Cpus
Mem
Start 2016-11-13 13:51:52
End 2016-11-13 16:00:03
Elapsed 02:08:11
Timeout 00:00:-1
Wall Timeout 100 days
Input files /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/align/ctl1/ss_50M_2000_GRCh38.nodup.bam
Output files /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/align/ctl1/ss_50M_2000_GRCh38.nodup.tagAlign.gz
Dependencies 
 
# SYS command. line 507

 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 510

 bedtools bamtobed -i /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/align/ctl1/ss_50M_2000_GRCh38.nodup.bam | awk 'BEGIN{OFS="\t"}{$4="N";$5="1000";print $0}' | gzip -nc > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/align/ctl1/ss_50M_2000_GRCh38.nodup.tagAlign.gz

# SYS command. line 512

 echo

# SYS command. line 514

 TASKTIME=$[$(date +%s)-${STARTTIME}]; if [ ${TASKTIME} -lt 60 ]; then echo "Waiting for $[60-${TASKTIME}] seconds."; sleep $[60-${TASKTIME}]; fi
 
   
--------------------Stdout--------------------
2237 (process ID) old priority 0, new priority 10
Num 5
ID task.postalign_bam.markdup_bam_picard_rep1.line_348.id_14
Name markdup_bam_picard rep1
Thread thread_41
PID 9519
OK true
Exit Code 0
Retries
State FINISHED
Dep. OK
Cpus 1
Mem
Start 2016-11-13 15:40:32
End 2016-11-13 18:04:19
Elapsed 02:23:47
Timeout 00:00:-1
Wall Timeout 100 days
Input files /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/align/rep1/SRR1370897_1.filt.bam
Output files /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/align/rep1/SRR1370897_1.dupmark.bam /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/qc/rep1/SRR1370897_1.dup.qc
Dependencies 
 
# SYS command. line 350

 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 354

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

# SYS command. line 360

 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 369

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

# SYS command. line 381

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

 
--------------------Stderr--------------------
Picked up _JAVA_OPTIONS: -Xms256M -Xmx12G -XX:ParallelGCThreads=1
[Sun Nov 13 17:58:08 PST 2016] picard.sam.markduplicates.MarkDuplicates INPUT=[/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/align/rep1/SRR1370897_1.filt.bam] OUTPUT=/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/align/rep1/SRR1370897_1.dupmark.bam METRICS_FILE=/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/qc/rep1/SRR1370897_1.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
[Sun Nov 13 17:58:08 PST 2016] Executing as imk1@surya 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-13 17:58:09	MarkDuplicates	Start of doWork freeMemory: 254879344; totalMemory: 257425408; maxMemory: 11453595648
INFO	2016-11-13 17:58:09	MarkDuplicates	Reading input file and constructing read end information.
INFO	2016-11-13 17:58:09	MarkDuplicates	Will retain up to 44052290 data points before spilling to disk.
INFO	2016-11-13 17:58:21	MarkDuplicates	Read     1,000,000 records.  Elapsed time: 00:00:09s.  Time for last 1,000,000:    9s.  Last read position: chr1:117,209,917
INFO	2016-11-13 17:58:21	MarkDuplicates	Tracking 0 as yet unmatched pairs. 0 records in RAM.
INFO	2016-11-13 17:58:33	MarkDuplicates	Read     2,000,000 records.  Elapsed time: 00:00:20s.  Time for last 1,000,000:   11s.  Last read position: chr2:9,424,442
INFO	2016-11-13 17:58:33	MarkDuplicates	Tracking 0 as yet unmatched pairs. 0 records in RAM.
INFO	2016-11-13 17:58:43	MarkDuplicates	Read     3,000,000 records.  Elapsed time: 00:00:31s.  Time for last 1,000,000:   10s.  Last read position: chr2:199,996,972
INFO	2016-11-13 17:58:43	MarkDuplicates	Tracking 0 as yet unmatched pairs. 0 records in RAM.
INFO	2016-11-13 17:58:49	MarkDuplicates	Read     4,000,000 records.  Elapsed time: 00:00:36s.  Time for last 1,000,000:    5s.  Last read position: chr3:129,701,387
INFO	2016-11-13 17:58:49	MarkDuplicates	Tracking 0 as yet unmatched pairs. 0 records in RAM.
INFO	2016-11-13 17:59:02	MarkDuplicates	Read     5,000,000 records.  Elapsed time: 00:00:49s.  Time for last 1,000,000:   13s.  Last read position: chr4:121,546,638
INFO	2016-11-13 17:59:02	MarkDuplicates	Tracking 0 as yet unmatched pairs. 0 records in RAM.
INFO	2016-11-13 17:59:07	MarkDuplicates	Read     6,000,000 records.  Elapsed time: 00:00:54s.  Time for last 1,000,000:    5s.  Last read position: chr5:82,791,477
INFO	2016-11-13 17:59:07	MarkDuplicates	Tracking 0 as yet unmatched pairs. 0 records in RAM.
INFO	2016-11-13 17:59:17	MarkDuplicates	Read     7,000,000 records.  Elapsed time: 00:01:05s.  Time for last 1,000,000:   10s.  Last read position: chr6:19,061,998
INFO	2016-11-13 17:59:17	MarkDuplicates	Tracking 0 as yet unmatched pairs. 0 records in RAM.
INFO	2016-11-13 17:59:21	MarkDuplicates	Read     8,000,000 records.  Elapsed time: 00:01:09s.  Time for last 1,000,000:    4s.  Last read position: chr6:144,216,450
INFO	2016-11-13 17:59:21	MarkDuplicates	Tracking 0 as yet unmatched pairs. 0 records in RAM.
INFO	2016-11-13 17:59:31	MarkDuplicates	Read     9,000,000 records.  Elapsed time: 00:01:18s.  Time for last 1,000,000:    9s.  Last read position: chr7:132,409,342
INFO	2016-11-13 17:59:31	MarkDuplicates	Tracking 0 as yet unmatched pairs. 0 records in RAM.
INFO	2016-11-13 17:59:39	MarkDuplicates	Read    10,000,000 records.  Elapsed time: 00:01:26s.  Time for last 1,000,000:    7s.  Last read position: chr8:142,402,266
INFO	2016-11-13 17:59:39	MarkDuplicates	Tracking 0 as yet unmatched pairs. 0 records in RAM.
INFO	2016-11-13 17:59:42	MarkDuplicates	Read    11,000,000 records.  Elapsed time: 00:01:29s.  Time for last 1,000,000:    3s.  Last read position: chr10:49,252,224
INFO	2016-11-13 17:59:42	MarkDuplicates	Tracking 0 as yet unmatched pairs. 0 records in RAM.
INFO	2016-11-13 17:59:47	MarkDuplicates	Read    12,000,000 records.  Elapsed time: 00:01:35s.  Time for last 1,000,000:    5s.  Last read position: chr11:64,356,620
INFO	2016-11-13 17:59:47	MarkDuplicates	Tracking 0 as yet unmatched pairs. 0 records in RAM.
INFO	2016-11-13 18:00:13	MarkDuplicates	Read    13,000,000 records.  Elapsed time: 00:02:00s.  Time for last 1,000,000:   25s.  Last read position: chr12:76,377,956
INFO	2016-11-13 18:00:13	MarkDuplicates	Tracking 0 as yet unmatched pairs. 0 records in RAM.
INFO	2016-11-13 18:00:28	MarkDuplicates	Read    14,000,000 records.  Elapsed time: 00:02:16s.  Time for last 1,000,000:   15s.  Last read position: chr13:96,116,659
INFO	2016-11-13 18:00:28	MarkDuplicates	Tracking 0 as yet unmatched pairs. 0 records in RAM.
INFO	2016-11-13 18:00:32	MarkDuplicates	Read    15,000,000 records.  Elapsed time: 00:02:19s.  Time for last 1,000,000:    3s.  Last read position: chr15:73,329,627
INFO	2016-11-13 18:00:32	MarkDuplicates	Tracking 0 as yet unmatched pairs. 0 records in RAM.
INFO	2016-11-13 18:00:37	MarkDuplicates	Read    16,000,000 records.  Elapsed time: 00:02:24s.  Time for last 1,000,000:    5s.  Last read position: chr17:39,679,448
INFO	2016-11-13 18:00:37	MarkDuplicates	Tracking 0 as yet unmatched pairs. 0 records in RAM.
INFO	2016-11-13 18:00:40	MarkDuplicates	Read    17,000,000 records.  Elapsed time: 00:02:27s.  Time for last 1,000,000:    2s.  Last read position: chr19:14,800,929
INFO	2016-11-13 18:00:40	MarkDuplicates	Tracking 0 as yet unmatched pairs. 0 records in RAM.
INFO	2016-11-13 18:00:46	MarkDuplicates	Read    18,000,000 records.  Elapsed time: 00:02:34s.  Time for last 1,000,000:    6s.  Last read position: chr21:35,633,624
INFO	2016-11-13 18:00:46	MarkDuplicates	Tracking 0 as yet unmatched pairs. 0 records in RAM.
INFO	2016-11-13 18:00:49	MarkDuplicates	Read    19,000,000 records.  Elapsed time: 00:02:37s.  Time for last 1,000,000:    2s.  Last read position: chrX:88,376,830
INFO	2016-11-13 18:00:49	MarkDuplicates	Tracking 0 as yet unmatched pairs. 0 records in RAM.
INFO	2016-11-13 18:00:52	MarkDuplicates	Read 19622310 records. 0 pairs never matched.
INFO	2016-11-13 18:00:56	MarkDuplicates	After buildSortedReadEndLists freeMemory: 4271229816; totalMemory: 5958533120; maxMemory: 11453595648
INFO	2016-11-13 18:00:56	MarkDuplicates	Will retain up to 357924864 duplicate indices before spilling to disk.
INFO	2016-11-13 18:00:57	MarkDuplicates	Traversing read pair information and detecting duplicates.
INFO	2016-11-13 18:00:57	MarkDuplicates	Traversing fragment information and detecting duplicates.
INFO	2016-11-13 18:00:58	MarkDuplicates	Sorting list of duplicate records.
INFO	2016-11-13 18:01:25	MarkDuplicates	After generateDuplicateIndexes freeMemory: 5690098496; totalMemory: 8630829056; maxMemory: 11453595648
INFO	2016-11-13 18:01:25	MarkDuplicates	Marking 2797721 records as duplicates.
INFO	2016-11-13 18:01:25	MarkDuplicates	Found 0 optical duplicate clusters.
INFO	2016-11-13 18:02:51	MarkDuplicates	Written    10,000,000 records.  Elapsed time: 00:01:26s.  Time for last 10,000,000:   86s.  Last read position: chr8:142,402,266
INFO	2016-11-13 18:04:18	MarkDuplicates	Before output close freeMemory: 8613964736; totalMemory: 8695316480; maxMemory: 11453595648
INFO	2016-11-13 18:04:18	MarkDuplicates	After output close freeMemory: 8746073672; totalMemory: 8827437056; maxMemory: 11453595648
[Sun Nov 13 18:04:18 PST 2016] picard.sam.markduplicates.MarkDuplicates done. Elapsed time: 6.17 minutes.
Runtime.totalMemory()=8827437056
Num 6
ID task.postalign_bam.dedup_bam_2_rep1.line_144.id_15
Name dedup_bam_2 rep1
Thread thread_41
PID 9677
OK true
Exit Code 0
Retries
State FINISHED
Dep. OK
Cpus
Mem
Start 2016-11-13 18:04:21
End 2016-11-13 20:12:24
Elapsed 02:08:02
Timeout 00:00:-1
Wall Timeout 100 days
Input files /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/align/rep1/SRR1370897_1.dupmark.bam
Output files /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/align/rep1/SRR1370897_1.nodup.bam /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/qc/rep1/SRR1370897_1.nodup.flagstat.qc /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/qc/rep1/SRR1370897_1.nodup.pbc.qc
Dependencies 
 
# SYS command. line 146

 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 150

 samtools view -F 1804 -b /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/align/rep1/SRR1370897_1.dupmark.bam > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/align/rep1/SRR1370897_1.nodup.bam

# SYS command. line 153

 sambamba index -t 1 /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/align/rep1/SRR1370897_1.nodup.bam

# SYS command. line 155

 sambamba flagstat -t 1 /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/align/rep1/SRR1370897_1.nodup.bam > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/qc/rep1/SRR1370897_1.nodup.flagstat.qc

# SYS command. line 165

 bedtools bamtobed -i /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/align/rep1/SRR1370897_1.dupmark.bam | \
				awk 'BEGIN{OFS="\t"}{print $1,$2,$3,$6}' | \
				grep -v 'chrM' | sort | uniq -c | \
				awk 'BEGIN{mt=0;m0=0;m1=0;m2=0} ($1==1){m1=m1+1} ($1==2){m2=m2+1} {m0=m0+1} {mt=mt+$1} END{m1_m2=-1.0; if(m2>0) m1_m2=m1/m2; printf "%d\t%d\t%d\t%d\t%f\t%f\t%f\n",mt,m0,m1,m2,m0/mt,m1/m0,m1_m2}' > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/qc/rep1/SRR1370897_1.nodup.pbc.qc

# SYS command. line 169

 TASKTIME=$[$(date +%s)-${STARTTIME}]; if [ ${TASKTIME} -lt 60 ]; then echo "Waiting for $[60-${TASKTIME}] seconds."; sleep $[60-${TASKTIME}]; fi
 
   
--------------------Stdout--------------------
9481 (process ID) old priority 0, new priority 10
Num 7
ID task.postalign_bam.bam_to_tag_rep1.line_505.id_16
Name bam_to_tag rep1
Thread thread_41
PID 9855
OK true
Exit Code 0
Retries
State FINISHED
Dep. OK
Cpus
Mem
Start 2016-11-13 20:12:27
End 2016-11-13 20:27:56
Elapsed 00:15:28
Timeout 00:00:-1
Wall Timeout 100 days
Input files /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/align/rep1/SRR1370897_1.nodup.bam
Output files /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/align/rep1/SRR1370897_1.nodup.tagAlign.gz
Dependencies 
 
# SYS command. line 507

 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 510

 bedtools bamtobed -i /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/align/rep1/SRR1370897_1.nodup.bam | awk 'BEGIN{OFS="\t"}{$4="N";$5="1000";print $0}' | gzip -nc > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/align/rep1/SRR1370897_1.nodup.tagAlign.gz

# SYS command. line 512

 echo

# SYS command. line 514

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

Waiting for 10 seconds.
Num 8
ID task.postalign_bed.subsample_tag_rep1.line_45.id_17
Name subsample_tag rep1
Thread thread_41
PID 9928
OK true
Exit Code 0
Retries
State FINISHED
Dep. OK
Cpus
Mem
Start 2016-11-13 20:27:57
End 2016-11-13 20:37:07
Elapsed 00:09:10
Timeout 00:00:-1
Wall Timeout 100 days
Input files /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/align/rep1/SRR1370897_1.nodup.tagAlign.gz
Output files /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/align/rep1/SRR1370897_1.nodup.15M.tagAlign.gz
Dependencies 
 
# SYS command. line 47

 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 50

 zcat /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/align/rep1/SRR1370897_1.nodup.tagAlign.gz | \
			grep -v "chrM" | shuf -n 15000000 --random-source=/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/align/rep1/SRR1370897_1.nodup.tagAlign.gz | gzip -nc > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/align/rep1/SRR1370897_1.nodup.15M.tagAlign.gz

# SYS command. line 53

 TASKTIME=$[$(date +%s)-${STARTTIME}]; if [ ${TASKTIME} -lt 60 ]; then echo "Waiting for $[60-${TASKTIME}] seconds."; sleep $[60-${TASKTIME}]; fi
 
   
--------------------Stdout--------------------
1253 (process ID) old priority 0, new priority 10
Waiting for 24 seconds.
Num 9
ID task.postalign_bed.spr_rep1.line_255.id_18
Name spr rep1
Thread thread_41
PID 9929
OK true
Exit Code 0
Retries
State FINISHED
Dep. OK
Cpus
Mem
Start 2016-11-13 20:27:58
End 2016-11-13 20:39:36
Elapsed 00:11:38
Timeout 00:00:-1
Wall Timeout 100 days
Input files /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/align/rep1/SRR1370897_1.nodup.tagAlign.gz
Output files /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/align/pseudo_reps/rep1/pr1/SRR1370897_1.nodup.pr1.tagAlign.gz /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/align/pseudo_reps/rep1/pr2/SRR1370897_1.nodup.pr2.tagAlign.gz
Dependencies 
 
# SYS command. line 257

 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 260

 nlines=$( zcat /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/align/rep1/SRR1370897_1.nodup.tagAlign.gz | wc -l )

# SYS command. line 261

 nlines=$(( (nlines + 1) / 2 ))

# SYS command. line 265

 zcat /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/align/rep1/SRR1370897_1.nodup.tagAlign.gz | shuf --random-source=/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/align/rep1/SRR1370897_1.nodup.tagAlign.gz | split -d -l $((nlines)) - /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/align/pseudo_reps/rep1/pr1/SRR1370897_1.nodup.

# SYS command. line 268

 gzip -nc /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/align/pseudo_reps/rep1/pr1/SRR1370897_1.nodup.00 > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/align/pseudo_reps/rep1/pr1/SRR1370897_1.nodup.pr1.tagAlign.gz

# SYS command. line 269

 rm -f /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/align/pseudo_reps/rep1/pr1/SRR1370897_1.nodup.00

# SYS command. line 270

 gzip -nc /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/align/pseudo_reps/rep1/pr1/SRR1370897_1.nodup.01 > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/align/pseudo_reps/rep1/pr2/SRR1370897_1.nodup.pr2.tagAlign.gz

# SYS command. line 271

 rm -f /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/align/pseudo_reps/rep1/pr1/SRR1370897_1.nodup.01

# SYS command. line 273

 TASKTIME=$[$(date +%s)-${STARTTIME}]; if [ ${TASKTIME} -lt 60 ]; then echo "Waiting for $[60-${TASKTIME}] seconds."; sleep $[60-${TASKTIME}]; fi
 
   
--------------------Stdout--------------------
12379 (process ID) old priority 0, new priority 10
Num 10
ID task.postalign_bed.xcor_rep1.line_212.id_19
Name xcor rep1
Thread thread_41
PID 9955
OK false
Exit Code 1
Retries
State ERROR
Dep. ERROR
Cpus
Mem
Start 2016-11-13 20:39:37
End 2016-11-13 20:39:37
Elapsed 00:00:00
Timeout 00:00:-1
Wall Timeout 100 days
Input files /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/align/rep1/SRR1370897_1.nodup.15M.tagAlign.gz
Output files /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/qc/rep1/SRR1370897_1.nodup.15M.cc.qc /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/qc/rep1/SRR1370897_1.nodup.15M.cc.plot.pdf
Dependencies 
 
# SYS command. line 214

 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 217

 if [[ $(which run_spp_nodups.R 2> /dev/null | wc -l || echo) == "1" ]]; then RUN_SPP=$(which run_spp_nodups.R); \
		    else RUN_SPP=$(which run_spp.R); \
		    fi

# SYS command. line 223

 Rscript ${RUN_SPP} -rf \
			-c=/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/align/rep1/SRR1370897_1.nodup.15M.tagAlign.gz -p=1 \
			-filtchr=chrM -savp=/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/qc/rep1/SRR1370897_1.nodup.15M.cc.plot.pdf -out=/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/qc/rep1/SRR1370897_1.nodup.15M.cc.qc 

# SYS command. line 226

 sed -r 's/,[^\t]+//g' /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/qc/rep1/SRR1370897_1.nodup.15M.cc.qc > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/qc/rep1/SRR1370897_1.nodup.15M.cc.qc.tmp

# SYS command. line 227

 mv /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/qc/rep1/SRR1370897_1.nodup.15M.cc.qc.tmp /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out/qc/rep1/SRR1370897_1.nodup.15M.cc.qc

# SYS command. line 229

 TASKTIME=$[$(date +%s)-${STARTTIME}]; if [ ${TASKTIME} -lt 60 ]; then echo "Waiting for $[60-${TASKTIME}] seconds."; sleep $[60-${TASKTIME}]; fi
 
     
--------------------Post mortem info--------------------
==============================================================
job_number:                 9955
exec_file:                  job_scripts/9955
submission_time:            Sun Nov 13 20:39:37 2016
owner:                      imk1
uid:                        1048
group:                      users
gid:                        100
sge_o_home:                 /users/imk1/
sge_o_log_name:             imk1
sge_o_path:                 /software/ascp/:/software/sratoolkit/2.5.2/bin/:/users/imk1/anaconda2/bin:/users/imk1/perl5/bin:/srv/scratch/imk1/TFBindingPredictionProject/src/bedtools2/bin:/usr/local/cuda/bin:/srv/persistent/imk1/TFBindingPredictionProject/src/rcade/:/srv/scratch/imk1/TFBindingPredictionProject/src/R-3.2.3/bin:/usr/local/sbin:/usr/local/bin:/usr/sbin:/usr/bin:/sbin:/bin:/usr/games:/usr/local/games:/snap/bin:/software/miniconda3/bin:/usr/lib/jvm/java-8-oracle/bin:/usr/lib/jvm/java-8-oracle/db/bin:/usr/lib/jvm/java-8-oracle/jre/bin:/users/imk1/edirect:/users/imk1/edirect
sge_o_shell:                /bin/bash
sge_o_workdir:              /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1
sge_o_host:                 surya
account:                    sge
stderr_path_list:           NONE:NONE:/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/chipseq.bds.20161113_130805_528_parallel_41/task.postalign_bed.xcor_rep1.line_212.id_19.stderr.cluster
mail_list:                  imk1@surya
notify:                     FALSE
job_name:                   STDIN
stdout_path_list:           NONE:NONE:/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/chipseq.bds.20161113_130805_528_parallel_41/task.postalign_bed.xcor_rep1.line_212.id_19.stdout.cluster
jobshare:                   0
env_list:                   MAIL=/var/mail/imk1,SSH_CLIENT=171.65.77.8 55144 22,USER=imk1,J2SDKDIR=/usr/lib/jvm/java-8-oracle,SHLVL=3,PERL_LOCAL_LIB_ROOT=/users/imk1/perl5,J2REDIR=/usr/lib/jvm/java-8-oracle/jre,HOME=/users/imk1/,DEEPLIFT_DIR=/srv/scratch/imk1/TFBindingPredictionProject/src/deeplift_private/deeplift,SSH_TTY=/dev/pts/4,ASCPKEY=/software/ascp/etc/asperaweb_id_dsa.openssh,_LMFILES_=/software/modulefiles/modsappsdir/modsappsdir:/software/modulefiles/sratoolkit/2.5.2:/software/modulefiles/ascp/3.5.6,LOGNAME=imk1,_=/usr/bin/bds,XDG_SESSION_ID=20,TERM=screen,PERL_MB_OPT=--install_base "/users/imk1/perl5",KERAS_DIR=/users/imk1/.local/lib/python2.7/site-packages/keras/,SGE_ROOT=/var/lib/gridengine,PATH=/software/ascp/:/software/sratoolkit/2.5.2/bin/:/users/imk1/anaconda2/bin:/users/imk1/perl5/bin:/srv/scratch/imk1/TFBindingPredictionProject/src/bedtools2/bin:/usr/local/cuda/bin:/srv/persistent/imk1/TFBindingPredictionProject/src/rcade/:/srv/scratch/imk1/TFBindingPredictionProject/src/R-3.2.3/bin:/usr/local/sbin:/usr/local/bin:/usr/sbin:/usr/bin:/sbin:/bin:/usr/games:/usr/local/games:/snap/bin:/software/miniconda3/bin:/usr/lib/jvm/java-8-oracle/bin:/usr/lib/jvm/java-8-oracle/db/bin:/usr/lib/jvm/java-8-oracle/jre/bin:/users/imk1/edirect:/users/imk1/edirect,DERBY_HOME=/usr/lib/jvm/java-8-oracle/db,PERL5LIB=/users/imk1/perl5/lib/perl5,XDG_RUNTIME_DIR=/run/user/1048,DISPLAY=localhost:12.0,SGE_CELL=default,STY=9560.SRR1370897_1.BDS,LANG=en_US.UTF-8,SHELL=/bin/bash,KRB5CCNAME=FILE:/tmp/krb5cc_1048_N5SeQf,XFILESEARCHPATH=/usr/dt/app-defaults/%L/Dt,RULEFITBASE=/srv/scratch/imk1/TFBindingPredictionProject/src/RuleFit/,PERL_MM_OPT=INSTALL_BASE=/users/imk1/perl5,UTIL_SCRIPTS_DIR=/srv/scratch/imk1/TFBindingPredictionProject/src/av_scripts/,MODULE_VERSION=3.2.10,MODULE_VERSION_STACK=3.2.10,WINDOW=0,NLSPATH=/usr/dt/lib/nls/msg/%L/%N.cat,PWD=/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1,JAVA_HOME=/usr/lib/jvm/java-8-oracle,LOADEDMODULES=modsappsdir/modsappsdir:sratoolkit/2.5.2:ascp/3.5.6,SSH_CONNECTION=171.65.77.8 55144 171.65.76.63 22,TERMCAP=SC|screen|VT 100/ANSI X3.64 virtual terminal:\,PYTHONPATH=/users/imk1/.local/lib/python2.7/site-packages/keras/:/srv/scratch/imk1/TFBindingPredictionProject/src/pybedtools/,ENHANCER_SCRIPTS_DIR=/srv/scratch/imk1/TFBindingPredictionProject/src/enhancer_prediction_code/,MODULEPATH=/usr/local/Modules/versions				:/usr/local/Modules/$MODULE_VERSION/modulefiles	:/modules/				:/software/modulefiles,MODULESAPPSDIR=/software,MODULESHOME=/software/env_module/3.2.10
script_file:                STDIN
usage    1:                 cpu=00:00:01, mem=0.09915 GBs, io=0.06815, vmem=51.738M, maxvmem=303.250M
scheduling info:            queue instance "q@kali" dropped because it is full
                            queue instance "q@nandi" dropped because it is full
                            queue instance "q@surya" dropped because it is full
                            queue instance "q@kadru" dropped because it is full

Global scope

Note: Global scope when program finished execution.
Type Name Value
string _ /usr/bin/bds
string _LMFILES_ /software/modulefiles/modsappsdir/modsappsdir:/software/modulefiles/sratoolkit/2.5.2:/software/modulefiles/ascp/3.5.6
bool allowEmpty false
string[] args [-title, SRR1370897_1, -nth, 5, -macs2_for_pooled_rep_only, -species, hg38, -url_base, http://mitra.stanford.edu/kundaje/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1/out, -fastq, /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/RawReadsFiles/SRR1370897_1.fastq.gz, -ctl_bam, /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/EncodeOtherZnfData/Encode3Datahg38/InputData/ss_50M_2000_GRCh38.bam, -final_stage, xcor, $]
string ASCPKEY /software/ascp/etc/asperaweb_id_dsa.openssh
bool canFail false
int cpus -1
int cpusLocal 56
int day 86400
string DEEPLIFT_DIR /srv/scratch/imk1/TFBindingPredictionProject/src/deeplift_private/deeplift
string DERBY_HOME /usr/lib/jvm/java-8-oracle/db
string DISPLAY localhost:12.0
real E 2.718281828459045
string ENHANCER_SCRIPTS_DIR /srv/scratch/imk1/TFBindingPredictionProject/src/enhancer_prediction_code/
int G 1073741824
string HOME /users/imk1/
int hour 3600
string J2REDIR /usr/lib/jvm/java-8-oracle/jre
string J2SDKDIR /usr/lib/jvm/java-8-oracle
string JAVA_HOME /usr/lib/jvm/java-8-oracle
int K 1024
string KERAS_DIR /users/imk1/.local/lib/python2.7/site-packages/keras/
string KRB5CCNAME FILE:/tmp/krb5cc_1048_N5SeQf
string LANG en_US.UTF-8
string LOADEDMODULES modsappsdir/modsappsdir:sratoolkit/2.5.2:ascp/3.5.6
string LOGNAME imk1
int M 1048576
string MAIL /var/mail/imk1
int mem -1
int minute 60
string MODULE_VERSION 3.2.10
string MODULE_VERSION_STACK 3.2.10
string MODULEPATH /usr/local/Modules/versions\t\t\t\t:/usr/local/Modules/$MODULE_VERSION/modulefiles\t:/modules/\t\t\t\t:/software/modulefiles
string MODULESAPPSDIR /software
string MODULESHOME /software/env_module/3.2.10
string NLSPATH /usr/dt/lib/nls/msg/%L/%N.cat
string node
int P 1125899906842624
string PATH /software/ascp/:/software/sratoolkit/2.5.2/bin/:/users/imk1/anaconda2/bin:/users/imk1/perl5/bin:/srv/scratch/imk1/TFBindingPredictionProject/src/bedtools2/bin:/usr/local/cuda/bin:/srv/persistent/imk1/TFBindingPredictionProject/src/rcade/:/srv/scratch/imk1/TFBindingPredictionProject/src/R-3.2.3/bin:/usr/local/sbin:/usr/local/bin:/usr/sbin:/usr/bin:/sbin:/bin:/usr/games:/usr/local/games:/snap/bin:/software/miniconda3/bin:/usr/lib/jvm/java-8-oracle/bin:/usr/lib/jvm/java-8-oracle/db/bin:/usr/lib/jvm/java-8-oracle/jre/bin:/users/imk1/edirect:/users/imk1/edirect
string PERL5LIB /users/imk1/perl5/lib/perl5
string PERL_LOCAL_LIB_ROOT /users/imk1/perl5
string PERL_MB_OPT --install_base \"/users/imk1/perl5\"
string PERL_MM_OPT INSTALL_BASE=/users/imk1/perl5
real PI 3.141592653589793
string ppwd /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1
string programName chipseq.bds
string programPath /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/src/chipseq_pipeline/chipseq.bds
string PWD /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/HughesZnfData/SRR1370897_1
string PYTHONPATH /users/imk1/.local/lib/python2.7/site-packages/keras/:/srv/scratch/imk1/TFBindingPredictionProject/src/pybedtools/
string queue
int retry 0
string RULEFITBASE /srv/scratch/imk1/TFBindingPredictionProject/src/RuleFit/
string SHELL /bin/bash
string SHLVL 3
string SSH_CLIENT 171.65.77.8 55144 22
string SSH_CONNECTION 171.65.77.8 55144 171.65.76.63 22
string SSH_TTY /dev/pts/4
string STY 9560.SRR1370897_1.BDS
string system sge
int T 1099511627776
string TERM screen
string TERMCAP SC|screen|VT 100/ANSI X3.64 virtual terminal:\\\n\t:DO=\\E[%dB:LE=\\E[%dD:RI=\\E[%dC:UP=\\E[%dA:bs:bt=\\E[Z:\\\n\t:cd=\\E[J:ce=\\E[K:cl=\\E[H\\E[J:cm=\\E[%i%d;%dH:ct=\\E[3g:\\\n\t:do=^J:nd=\\E[C:pt:rc=\\E8:rs=\\Ec:sc=\\E7:st=\\EH:up=\\EM:\\\n\t:le=^H:bl=^G:cr=^M:it#8:ho=\\E[H:nw=\\EE:ta=^I:is=\\E)0:\\\n\t:li#24:co#80:am:xn:xv:LP:sr=\\EM:al=\\E[L:AL=\\E[%dL:\\\n\t:cs=\\E[%i%d;%dr:dl=\\E[M:DL=\\E[%dM:dc=\\E[P:DC=\\E[%dP:\\\n\t:im=\\E[4h:ei=\\E[4l:mi:IC=\\E[%d@:ks=\\E[?1h\\E=:\\\n\t:ke=\\E[?1l\\E>:vi=\\E[?25l:ve=\\E[34h\\E[?25h:vs=\\E[34l:\\\n\t:ti=\\E[?1049h:te=\\E[?1049l:Km=\\E[M:k0=\\E[10~:k1=\\EOP:\\\n\t:k2=\\EOQ:k3=\\EOR:k4=\\EOS:k5=\\E[15~:k6=\\E[17~:k7=\\E[18~:\\\n\t:k8=\\E[19~:k9=\\E[20~:k;=\\E[21~:F1=\\E[23~:F2=\\E[24~:\\\n\t:kh=\\E[1~:@1=\\E[1~:kH=\\E[4~:@7=\\E[4~:kN=\\E[6~:kP=\\E[5~:\\\n\t:kI=\\E[2~:kD=\\E[3~:ku=\\EOA:kd=\\EOB:kr=\\EOC:kl=\\EOD:
int timeout -1
string USER imk1
string UTIL_SCRIPTS_DIR /srv/scratch/imk1/TFBindingPredictionProject/src/av_scripts/
int walltimeout 8640000
int week 604800
string WINDOW 0
string XDG_RUNTIME_DIR /run/user/1048
string XDG_SESSION_ID 20
string XFILESEARCHPATH /usr/dt/app-defaults/%L/Dt