Script file	/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/src/TF_chipseq_pipeline/chipseq.bds
Program ID	chipseq.bds.20160929_153859_619
Start time	2016-09-29 15:38:59
Run time	14:21:17.504
Tasks executed	8
Tasks failed	0
Tasks failed names
Arguments*	[-title, ZF11, -nth, 5, -no_pseudo_rep, -macs2_for_pooled_rep_only, -species, mm10, -url_base, http://mitra.stanford.edu/kundaje/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out, -fastq, /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/RawReads/ZF11_allReps.fastq.gz, -ctl_fastq, /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/InputData/r120621A_l1i7_aB180__wt__bioBirAonly.fastq.gz, $]
System*	sge
Cpus*	-1
Exit value	0

Thread ID	Parent	Tasks
thread_Root	None	chipseq.bds.20160929_153859_619/task.callpeak_spp.spp_rep1.line_54.id_16 chipseq.bds.20160929_153859_619/task.graphviz.report.line_97.id_17
thread_42	thread_Root	chipseq.bds.20160929_153859_619_parallel_42/task.align_bwa.bwa_aln_ctl1.line_108.id_10 chipseq.bds.20160929_153859_619_parallel_42/task.align_bwa.bwa_sam_ctl1.line_140.id_11 chipseq.bds.20160929_153859_619_parallel_42/task.postalign_bam.dedup_bam_1_ctl1.line_95.id_12 chipseq.bds.20160929_153859_619_parallel_42/task.postalign_bam.markdup_bam_picard_ctl1.line_347.id_13 chipseq.bds.20160929_153859_619_parallel_42/task.postalign_bam.dedup_bam_2_ctl1.line_144.id_14 chipseq.bds.20160929_153859_619_parallel_42/task.postalign_bam.bam_to_tag_ctl1.line_504.id_15
thread_41	thread_Root

Task details

Task

Execution

Time

Dependencies

Task program, Errors, StdOut / StdErr

Num	1
ID	task.align_bwa.bwa_aln_ctl1.line_108.id_10
Name	bwa_aln ctl1
Thread	thread_42
PID	9817

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

Start	2016-09-29 15:39:01
End	2016-09-29 16:47:43
Elapsed	01:08:42
Timeout	00:00:-1
Wall Timeout	100 days

Input files	/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/InputData/r120621A_l1i7_aB180__wt__bioBirAonly.fastq.gz
Output files	/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/align/ctl1/r120621A_l1i7_aB180__wt__bioBirAonly.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 2 /mnt/data/bds_pipeline_genome_data/mm10/bwa_index/GRCm38.p4.genome.fa /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/InputData/r120621A_l1i7_aB180__wt__bioBirAonly.fastq.gz > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/align/ctl1/r120621A_l1i7_aB180__wt__bioBirAonly.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--------------------
28557 (process ID) old priority 0, new priority 10

 
--------------------Stderr--------------------
discarding /software/miniconda3/bin from PATH
prepending /software/miniconda3/envs/aquas_chipseq/bin to PATH
[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] 1.2% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 76.59 sec
[bwa_aln_core] write to the disk... 0.10 sec
[bwa_aln_core] 262144 sequences have been processed.
[bwa_read_seq] 1.3% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 65.45 sec
[bwa_aln_core] write to the disk... 0.06 sec
[bwa_aln_core] 524288 sequences have been processed.
[bwa_read_seq] 1.2% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 70.92 sec
[bwa_aln_core] write to the disk... 0.06 sec
[bwa_aln_core] 786432 sequences have been processed.
[bwa_read_seq] 1.2% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 68.30 sec
[bwa_aln_core] write to the disk... 0.06 sec
[bwa_aln_core] 1048576 sequences have been processed.
[bwa_read_seq] 1.1% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 63.21 sec
[bwa_aln_core] write to the disk... 0.06 sec
[bwa_aln_core] 1310720 sequences have been processed.
[bwa_read_seq] 1.1% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 63.19 sec
[bwa_aln_core] write to the disk... 0.06 sec
[bwa_aln_core] 1572864 sequences have been processed.
[bwa_read_seq] 1.2% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 62.58 sec
[bwa_aln_core] write to the disk... 0.05 sec
[bwa_aln_core] 1835008 sequences have been processed.
[bwa_read_seq] 1.5% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 61.00 sec
[bwa_aln_core] write to the disk... 0.05 sec
[bwa_aln_core] 2097152 sequences have been processed.
[bwa_read_seq] 2.7% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 61.99 sec
[bwa_aln_core] write to the disk... 0.04 sec
[bwa_aln_core] 2359296 sequences have been processed.
[bwa_read_seq] 6.1% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 61.98 sec
[bwa_aln_core] write to the disk... 0.04 sec
[bwa_aln_core] 2621440 sequences have been processed.
[bwa_read_seq] 10.0% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 65.17 sec
[bwa_aln_core] write to the disk... 0.05 sec
[bwa_aln_core] 2883584 sequences have been processed.
[bwa_read_seq] 12.0% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 93.93 sec
[bwa_aln_core] write to the disk... 0.06 sec
[bwa_aln_core] 3145728 sequences have been processed.
[bwa_read_seq] 12.3% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 70.46 sec
[bwa_aln_core] write to the disk... 0.08 sec
[bwa_aln_core] 3407872 sequences have been processed.
[bwa_read_seq] 11.2% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 67.28 sec
[bwa_aln_core] write to the disk... 0.05 sec
[bwa_aln_core] 3670016 sequences have been processed.
[bwa_read_seq] 9.3% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 84.68 sec
[bwa_aln_core] write to the disk... 0.05 sec
[bwa_aln_core] 3932160 sequences have been processed.
[bwa_read_seq] 6.8% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 61.89 sec
[bwa_aln_core] write to the disk... 0.05 sec
[bwa_aln_core] 4194304 sequences have been processed.
[bwa_read_seq] 4.9% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 55.48 sec
[bwa_aln_core] write to the disk... 0.05 sec
[bwa_aln_core] 4456448 sequences have been processed.
[bwa_read_seq] 3.4% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 55.72 sec
[bwa_aln_core] write to the disk... 0.06 sec
[bwa_aln_core] 4718592 sequences have been processed.
[bwa_read_seq] 2.5% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 60.76 sec
[bwa_aln_core] write to the disk... 0.05 sec
[bwa_aln_core] 4980736 sequences have been processed.
[bwa_read_seq] 1.9% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 60.95 sec
[bwa_aln_core] write to the disk... 0.03 sec
[bwa_aln_core] 5242880 sequences have been processed.
[bwa_read_seq] 1.5% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 57.18 sec
[bwa_aln_core] write to the disk... 0.05 sec
[bwa_aln_core] 5505024 sequences have been processed.
[bwa_read_seq] 1.5% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 53.60 sec
[bwa_aln_core] write to the disk... 0.05 sec
[bwa_aln_core] 5767168 sequences have been processed.
[bwa_read_seq] 1.3% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 55.66 sec
[bwa_aln_core] write to the disk... 0.06 sec
[bwa_aln_core] 6029312 sequences have been processed.
[bwa_read_seq] 1.3% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 59.33 sec
[bwa_aln_core] write to the disk... 0.03 sec
[bwa_aln_core] 6291456 sequences have been processed.
[bwa_read_seq] 1.2% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 58.52 sec
[bwa_aln_core] write to the disk... 0.05 sec
[bwa_aln_core] 6553600 sequences have been processed.
[bwa_read_seq] 1.2% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 57.29 sec
[bwa_aln_core] write to the disk... 0.03 sec
[bwa_aln_core] 6815744 sequences have been processed.
[bwa_read_seq] 1.2% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 54.01 sec
[bwa_aln_core] write to the disk... 0.03 sec
[bwa_aln_core] 7077888 sequences have been processed.
[bwa_read_seq] 1.5% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 50.15 sec
[bwa_aln_core] write to the disk... 0.05 sec
[bwa_aln_core] 7340032 sequences have been processed.
[bwa_read_seq] 2.5% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 52.92 sec
[bwa_aln_core] write to the disk... 0.05 sec
[bwa_aln_core] 7602176 sequences have been processed.
[bwa_read_seq] 6.2% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 51.15 sec
[bwa_aln_core] write to the disk... 0.05 sec
[bwa_aln_core] 7864320 sequences have been processed.
[bwa_read_seq] 11.1% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 58.39 sec
[bwa_aln_core] write to the disk... 0.05 sec
[bwa_aln_core] 8126464 sequences have been processed.
[bwa_read_seq] 13.3% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 68.18 sec
[bwa_aln_core] write to the disk... 0.04 sec
[bwa_aln_core] 8388608 sequences have been processed.
[bwa_read_seq] 13.6% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 65.24 sec
[bwa_aln_core] write to the disk... 0.04 sec
[bwa_aln_core] 8650752 sequences have been processed.
[bwa_read_seq] 12.7% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 59.92 sec
[bwa_aln_core] write to the disk... 0.04 sec
[bwa_aln_core] 8912896 sequences have been processed.
[bwa_read_seq] 11.1% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 63.15 sec
[bwa_aln_core] write to the disk... 0.04 sec
[bwa_aln_core] 9175040 sequences have been processed.
[bwa_read_seq] 9.0% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 58.07 sec
[bwa_aln_core] write to the disk... 0.04 sec
[bwa_aln_core] 9437184 sequences have been processed.
[bwa_read_seq] 6.9% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 51.81 sec
[bwa_aln_core] write to the disk... 0.06 sec
[bwa_aln_core] 9699328 sequences have been processed.
[bwa_read_seq] 5.1% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 50.19 sec
[bwa_aln_core] write to the disk... 0.06 sec
[bwa_aln_core] 9961472 sequences have been processed.
[bwa_read_seq] 3.8% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 57.17 sec
[bwa_aln_core] write to the disk... 0.05 sec
[bwa_aln_core] 10223616 sequences have been processed.
[bwa_read_seq] 2.9% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 56.37 sec
[bwa_aln_core] write to the disk... 0.06 sec
[bwa_aln_core] 10485760 sequences have been processed.
[bwa_read_seq] 2.2% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 58.70 sec
[bwa_aln_core] write to the disk... 0.05 sec
[bwa_aln_core] 10747904 sequences have been processed.
[bwa_read_seq] 1.2% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 54.92 sec
[bwa_aln_core] write to the disk... 0.04 sec
[bwa_aln_core] 11010048 sequences have been processed.
[bwa_read_seq] 1.2% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 55.35 sec
[bwa_aln_core] write to the disk... 0.06 sec
[bwa_aln_core] 11272192 sequences have been processed.
[bwa_read_seq] 1.2% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 56.26 sec
[bwa_aln_core] write to the disk... 0.06 sec
[bwa_aln_core] 11534336 sequences have been processed.
[bwa_read_seq] 1.2% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 59.71 sec
[bwa_aln_core] write to the disk... 0.05 sec
[bwa_aln_core] 11796480 sequences have been processed.
[bwa_read_seq] 1.1% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 56.03 sec
[bwa_aln_core] write to the disk... 0.05 sec
[bwa_aln_core] 12058624 sequences have been processed.
[bwa_read_seq] 1.1% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 57.33 sec
[bwa_aln_core] write to the disk... 0.04 sec
[bwa_aln_core] 12320768 sequences have been processed.
[bwa_read_seq] 1.2% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 57.97 sec
[bwa_aln_core] write to the disk... 0.04 sec
[bwa_aln_core] 12582912 sequences have been processed.
[bwa_read_seq] 1.8% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 57.26 sec
[bwa_aln_core] write to the disk... 0.05 sec
[bwa_aln_core] 12845056 sequences have been processed.
[bwa_read_seq] 4.1% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 56.46 sec
[bwa_aln_core] write to the disk... 0.05 sec
[bwa_aln_core] 13107200 sequences have been processed.
[bwa_read_seq] 8.7% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 60.54 sec
[bwa_aln_core] write to the disk... 0.06 sec
[bwa_aln_core] 13369344 sequences have been processed.
[bwa_read_seq] 11.9% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 68.77 sec
[bwa_aln_core] write to the disk... 0.04 sec
[bwa_aln_core] 13631488 sequences have been processed.
[bwa_read_seq] 12.5% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 70.58 sec
[bwa_aln_core] write to the disk... 0.03 sec
[bwa_aln_core] 13893632 sequences have been processed.
[bwa_read_seq] 12.1% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 63.55 sec
[bwa_aln_core] write to the disk... 0.05 sec
[bwa_aln_core] 14155776 sequences have been processed.
[bwa_read_seq] 10.5% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 65.06 sec
[bwa_aln_core] write to the disk... 0.05 sec
[bwa_aln_core] 14417920 sequences have been processed.
[bwa_read_seq] 8.4% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 61.28 sec
[bwa_aln_core] write to the disk... 0.06 sec
[bwa_aln_core] 14680064 sequences have been processed.
[bwa_read_seq] 6.1% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 74.02 sec
[bwa_aln_core] write to the disk... 0.07 sec
[bwa_aln_core] 14942208 sequences have been processed.
[bwa_read_seq] 4.3% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 61.62 sec
[bwa_aln_core] write to the disk... 0.05 sec
[bwa_aln_core] 15204352 sequences have been processed.
[bwa_read_seq] 3.0% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 58.57 sec
[bwa_aln_core] write to the disk... 0.05 sec
[bwa_aln_core] 15466496 sequences have been processed.
[bwa_read_seq] 2.4% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 61.07 sec
[bwa_aln_core] write to the disk... 0.03 sec
[bwa_aln_core] 15728640 sequences have been processed.
[bwa_read_seq] 2.0% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 61.10 sec
[bwa_aln_core] write to the disk... 0.06 sec
[bwa_aln_core] 15990784 sequences have been processed.
[bwa_read_seq] 1.6% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 62.41 sec
[bwa_aln_core] write to the disk... 0.06 sec
[bwa_aln_core] 16252928 sequences have been processed.
[bwa_read_seq] 1.3% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 63.00 sec
[bwa_aln_core] write to the disk... 0.06 sec
[bwa_aln_core] 16515072 sequences have been processed.
[bwa_read_seq] 1.3% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 61.49 sec
[bwa_aln_core] write to the disk... 0.06 sec
[bwa_aln_core] 16777216 sequences have been processed.
[bwa_read_seq] 1.8% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 62.60 sec
[bwa_aln_core] write to the disk... 0.05 sec
[bwa_aln_core] 17039360 sequences have been processed.
[bwa_read_seq] 1.7% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 63.62 sec
[bwa_aln_core] write to the disk... 0.04 sec
[bwa_aln_core] 17301504 sequences have been processed.
[bwa_read_seq] 1.3% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 70.05 sec
[bwa_aln_core] write to the disk... 0.04 sec
[bwa_aln_core] 17563648 sequences have been processed.
[bwa_read_seq] 1.3% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 63.41 sec
[bwa_aln_core] write to the disk... 0.06 sec
[bwa_aln_core] 17825792 sequences have been processed.
[bwa_read_seq] 1.4% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 63.65 sec
[bwa_aln_core] write to the disk... 0.04 sec
[bwa_aln_core] 18087936 sequences have been processed.
[bwa_read_seq] 2.2% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 63.25 sec
[bwa_aln_core] write to the disk... 0.05 sec
[bwa_aln_core] 18350080 sequences have been processed.
[bwa_read_seq] 4.8% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 66.50 sec
[bwa_aln_core] write to the disk... 0.06 sec
[bwa_aln_core] 18612224 sequences have been processed.
[bwa_read_seq] 9.5% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 73.95 sec
[bwa_aln_core] write to the disk... 0.04 sec
[bwa_aln_core] 18874368 sequences have been processed.
[bwa_read_seq] 12.1% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 71.87 sec
[bwa_aln_core] write to the disk... 0.05 sec
[bwa_aln_core] 19136512 sequences have been processed.
[bwa_read_seq] 12.5% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 73.71 sec
[bwa_aln_core] write to the disk... 0.06 sec
[bwa_aln_core] 19398656 sequences have been processed.
[bwa_read_seq] 11.4% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 70.67 sec
[bwa_aln_core] write to the disk... 0.04 sec
[bwa_aln_core] 19660800 sequences have been processed.
[bwa_read_seq] 9.7% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 82.67 sec
[bwa_aln_core] write to the disk... 0.05 sec
[bwa_aln_core] 19922944 sequences have been processed.
[bwa_read_seq] 7.5% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 72.51 sec
[bwa_aln_core] write to the disk... 0.05 sec
[bwa_aln_core] 20185088 sequences have been processed.
[bwa_read_seq] 5.6% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 67.41 sec
[bwa_aln_core] write to the disk... 0.05 sec
[bwa_aln_core] 20447232 sequences have been processed.
[bwa_read_seq] 3.9% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 62.39 sec
[bwa_aln_core] write to the disk... 0.03 sec
[bwa_aln_core] 20709376 sequences have been processed.
[bwa_read_seq] 2.4% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 61.10 sec
[bwa_aln_core] write to the disk... 0.06 sec
[bwa_aln_core] 20971520 sequences have been processed.
[bwa_read_seq] 1.9% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 61.40 sec
[bwa_aln_core] write to the disk... 0.05 sec
[bwa_aln_core] 21233664 sequences have been processed.
[bwa_read_seq] 1.6% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 67.76 sec
[bwa_aln_core] write to the disk... 0.04 sec
[bwa_aln_core] 21495808 sequences have been processed.
[bwa_read_seq] 1.5% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 61.32 sec
[bwa_aln_core] write to the disk... 0.03 sec
[bwa_aln_core] 21757952 sequences have been processed.
[bwa_read_seq] 1.5% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 62.44 sec
[bwa_aln_core] write to the disk... 0.04 sec
[bwa_aln_core] 22020096 sequences have been processed.
[bwa_read_seq] 1.4% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 62.03 sec
[bwa_aln_core] write to the disk... 0.05 sec
[bwa_aln_core] 22282240 sequences have been processed.
[bwa_read_seq] 1.4% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 64.51 sec
[bwa_aln_core] write to the disk... 0.07 sec
[bwa_aln_core] 22544384 sequences have been processed.
[bwa_read_seq] 1.3% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 60.30 sec
[bwa_aln_core] write to the disk... 0.07 sec
[bwa_aln_core] 22806528 sequences have been processed.
[bwa_read_seq] 1.5% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 61.04 sec
[bwa_aln_core] write to the disk... 0.05 sec
[bwa_aln_core] 23068672 sequences have been processed.
[bwa_read_seq] 1.7% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 59.71 sec
[bwa_aln_core] write to the disk... 0.05 sec
[bwa_aln_core] 23330816 sequences have been processed.
[bwa_read_seq] 2.9% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 58.02 sec
[bwa_aln_core] write to the disk... 0.06 sec
[bwa_aln_core] 23592960 sequences have been processed.
[bwa_read_seq] 6.8% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 61.03 sec
[bwa_aln_core] write to the disk... 0.06 sec
[bwa_aln_core] 23855104 sequences have been processed.
[bwa_read_seq] 11.4% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 65.42 sec
[bwa_aln_core] write to the disk... 0.06 sec
[bwa_aln_core] 24117248 sequences have been processed.
[bwa_read_seq] 13.5% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 72.24 sec
[bwa_aln_core] write to the disk... 0.06 sec
[bwa_aln_core] 24379392 sequences have been processed.
[bwa_read_seq] 14.4% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 73.29 sec
[bwa_aln_core] write to the disk... 0.05 sec
[bwa_aln_core] 24641536 sequences have been processed.
[bwa_read_seq] 13.5% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 74.67 sec
[bwa_aln_core] write to the disk... 0.07 sec
[bwa_aln_core] 24903680 sequences have been processed.
[bwa_read_seq] 11.4% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 68.65 sec
[bwa_aln_core] write to the disk... 0.05 sec
[bwa_aln_core] 25165824 sequences have been processed.
[bwa_read_seq] 9.6% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 73.74 sec
[bwa_aln_core] write to the disk... 0.04 sec
[bwa_aln_core] 25427968 sequences have been processed.
[bwa_read_seq] 7.5% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 65.81 sec
[bwa_aln_core] write to the disk... 0.04 sec
[bwa_aln_core] 25690112 sequences have been processed.
[bwa_read_seq] 5.5% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 69.95 sec
[bwa_aln_core] write to the disk... 0.04 sec
[bwa_aln_core] 25952256 sequences have been processed.
[bwa_read_seq] 3.7% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 71.11 sec
[bwa_aln_core] write to the disk... 0.04 sec
[bwa_aln_core] 26214400 sequences have been processed.
[bwa_read_seq] 2.9% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 74.97 sec
[bwa_aln_core] write to the disk... 0.05 sec
[bwa_aln_core] 26476544 sequences have been processed.
[bwa_read_seq] 2.1% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 65.22 sec
[bwa_aln_core] write to the disk... 0.05 sec
[bwa_aln_core] 26738688 sequences have been processed.
[bwa_read_seq] 1.3% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 72.74 sec
[bwa_aln_core] write to the disk... 0.05 sec
[bwa_aln_core] 27000832 sequences have been processed.
[bwa_read_seq] 1.4% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 75.11 sec
[bwa_aln_core] write to the disk... 0.06 sec
[bwa_aln_core] 27262976 sequences have been processed.
[bwa_read_seq] 1.3% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 67.00 sec
[bwa_aln_core] write to the disk... 0.04 sec
[bwa_aln_core] 27525120 sequences have been processed.
[bwa_read_seq] 1.2% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 72.13 sec
[bwa_aln_core] write to the disk... 0.05 sec
[bwa_aln_core] 27787264 sequences have been processed.
[bwa_read_seq] 1.3% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 72.32 sec
[bwa_aln_core] write to the disk... 0.03 sec
[bwa_aln_core] 28049408 sequences have been processed.
[bwa_read_seq] 1.2% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 69.20 sec
[bwa_aln_core] write to the disk... 0.05 sec
[bwa_aln_core] 28311552 sequences have been processed.
[bwa_read_seq] 1.4% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 65.79 sec
[bwa_aln_core] write to the disk... 0.04 sec
[bwa_aln_core] 28573696 sequences have been processed.
[bwa_read_seq] 2.2% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 72.60 sec
[bwa_aln_core] write to the disk... 0.04 sec
[bwa_aln_core] 28835840 sequences have been processed.
[bwa_read_seq] 5.0% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 69.83 sec
[bwa_aln_core] write to the disk... 0.06 sec
[bwa_aln_core] 29097984 sequences have been processed.
[bwa_read_seq] 9.8% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 66.96 sec
[bwa_aln_core] write to the disk... 0.08 sec
[bwa_aln_core] 29360128 sequences have been processed.
[bwa_read_seq] 12.2% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 65.54 sec
[bwa_aln_core] write to the disk... 0.04 sec
[bwa_aln_core] 29622272 sequences have been processed.
[bwa_read_seq] 12.7% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 67.69 sec
[bwa_aln_core] write to the disk... 0.05 sec
[bwa_aln_core] 29884416 sequences have been processed.
[bwa_read_seq] 11.6% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 59.61 sec
[bwa_aln_core] write to the disk... 0.05 sec
[bwa_aln_core] 30146560 sequences have been processed.
[bwa_read_seq] 9.8% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 60.70 sec
[bwa_aln_core] write to the disk... 0.05 sec
[bwa_aln_core] 30408704 sequences have been processed.
[bwa_read_seq] 7.7% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 54.54 sec
[bwa_aln_core] write to the disk... 0.06 sec
[bwa_aln_core] 30670848 sequences have been processed.
[bwa_read_seq] 5.8% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 58.94 sec
[bwa_aln_core] write to the disk... 0.05 sec
[bwa_aln_core] 30932992 sequences have been processed.
[bwa_read_seq] 4.4% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 59.18 sec
[bwa_aln_core] write to the disk... 0.04 sec
[bwa_aln_core] 31195136 sequences have been processed.
[bwa_read_seq] 3.0% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 63.43 sec
[bwa_aln_core] write to the disk... 0.06 sec
[bwa_aln_core] 31457280 sequences have been processed.
[bwa_read_seq] 2.4% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 62.95 sec
[bwa_aln_core] write to the disk... 0.06 sec
[bwa_aln_core] 31719424 sequences have been processed.
[bwa_read_seq] 1.9% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 57.21 sec
[bwa_aln_core] write to the disk... 0.05 sec
[bwa_aln_core] 31981568 sequences have been processed.
[bwa_read_seq] 2.0% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 10.02 sec
[bwa_aln_core] write to the disk... 0.01 sec
[bwa_aln_core] 32022635 sequences have been processed.
[main] Version: 0.7.13-r1126
[main] CMD: bwa aln -q 5 -l 32 -k 2 -t 2 /mnt/data/bds_pipeline_genome_data/mm10/bwa_index/GRCm38.p4.genome.fa /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/InputData/r120621A_l1i7_aB180__wt__bioBirAonly.fastq.gz
[main] Real time: 4109.083 sec; CPU: 7934.439 sec

Num	2
ID	task.align_bwa.bwa_sam_ctl1.line_140.id_11
Name	bwa_sam ctl1
Thread	thread_42
PID	9830

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

Start	2016-09-29 16:47:44
End	2016-09-29 17:25:21
Elapsed	00:37:36
Timeout	00:00:-1
Wall Timeout	100 days

Input files	/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/InputData/r120621A_l1i7_aB180__wt__bioBirAonly.fastq.gz /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/align/ctl1/r120621A_l1i7_aB180__wt__bioBirAonly.sai
Output files	/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/align/ctl1/r120621A_l1i7_aB180__wt__bioBirAonly.bam /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/qc/ctl1/r120621A_l1i7_aB180__wt__bioBirAonly.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/mm10/bwa_index/GRCm38.p4.genome.fa /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/align/ctl1/r120621A_l1i7_aB180__wt__bioBirAonly.sai /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/InputData/r120621A_l1i7_aB180__wt__bioBirAonly.fastq.gz | samtools view -Su /dev/stdin \
			| sambamba sort -t 2 /dev/stdin -o /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/align/ctl1/r120621A_l1i7_aB180__wt__bioBirAonly.bam

# SYS command. line 149

 sambamba flagstat -t 2 /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/align/ctl1/r120621A_l1i7_aB180__wt__bioBirAonly.bam > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/qc/ctl1/r120621A_l1i7_aB180__wt__bioBirAonly.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--------------------
45113 (process ID) old priority 0, new priority 10

 
--------------------Stderr--------------------
discarding /software/miniconda3/bin from PATH
prepending /software/miniconda3/envs/aquas_chipseq/bin to PATH
[bwa_read_seq] 1.2% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 101.10 sec
[bwa_aln_core] refine gapped alignments... 9.97 sec
[bwa_aln_core] print alignments... [samopen] SAM header is present: 96 sequences.
0.57 sec
[bwa_aln_core] 262144 sequences have been processed.
[bwa_read_seq] 1.3% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 46.54 sec
[bwa_aln_core] refine gapped alignments... 5.73 sec
[bwa_aln_core] print alignments... 0.35 sec
[bwa_aln_core] 524288 sequences have been processed.
[bwa_read_seq] 1.2% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 46.94 sec
[bwa_aln_core] refine gapped alignments... 6.13 sec
[bwa_aln_core] print alignments... 0.41 sec
[bwa_aln_core] 786432 sequences have been processed.
[bwa_read_seq] 1.2% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 45.14 sec
[bwa_aln_core] refine gapped alignments... 5.24 sec
[bwa_aln_core] print alignments... 0.32 sec
[bwa_aln_core] 1048576 sequences have been processed.
[bwa_read_seq] 1.1% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 45.94 sec
[bwa_aln_core] refine gapped alignments... 5.32 sec
[bwa_aln_core] print alignments... 0.37 sec
[bwa_aln_core] 1310720 sequences have been processed.
[bwa_read_seq] 1.1% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 48.32 sec
[bwa_aln_core] refine gapped alignments... 6.04 sec
[bwa_aln_core] print alignments... 0.44 sec
[bwa_aln_core] 1572864 sequences have been processed.
[bwa_read_seq] 1.2% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 50.19 sec
[bwa_aln_core] refine gapped alignments... 8.01 sec
[bwa_aln_core] print alignments... 0.29 sec
[bwa_aln_core] 1835008 sequences have been processed.
[bwa_read_seq] 1.5% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 40.27 sec
[bwa_aln_core] refine gapped alignments... 5.90 sec
[bwa_aln_core] print alignments... 0.64 sec
[bwa_aln_core] 2097152 sequences have been processed.
[bwa_read_seq] 2.7% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 38.56 sec
[bwa_aln_core] refine gapped alignments... 6.61 sec
[bwa_aln_core] print alignments... 0.62 sec
[bwa_aln_core] 2359296 sequences have been processed.
[bwa_read_seq] 6.1% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 37.61 sec
[bwa_aln_core] refine gapped alignments... 6.24 sec
[bwa_aln_core] print alignments... 0.65 sec
[bwa_aln_core] 2621440 sequences have been processed.
[bwa_read_seq] 10.0% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 38.98 sec
[bwa_aln_core] refine gapped alignments... 6.50 sec
[bwa_aln_core] print alignments... 0.65 sec
[bwa_aln_core] 2883584 sequences have been processed.
[bwa_read_seq] 12.0% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 49.73 sec
[bwa_aln_core] refine gapped alignments... 6.32 sec
[bwa_aln_core] print alignments... 0.74 sec
[bwa_aln_core] 3145728 sequences have been processed.
[bwa_read_seq] 12.3% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 39.27 sec
[bwa_aln_core] refine gapped alignments... 8.02 sec
[bwa_aln_core] print alignments... 0.69 sec
[bwa_aln_core] 3407872 sequences have been processed.
[bwa_read_seq] 11.2% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 35.91 sec
[bwa_aln_core] refine gapped alignments... 5.92 sec
[bwa_aln_core] print alignments... 0.67 sec
[bwa_aln_core] 3670016 sequences have been processed.
[bwa_read_seq] 9.3% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 36.62 sec
[bwa_aln_core] refine gapped alignments... 7.67 sec
[bwa_aln_core] print alignments... 0.69 sec
[bwa_aln_core] 3932160 sequences have been processed.
[bwa_read_seq] 6.8% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 51.41 sec
[bwa_aln_core] refine gapped alignments... 9.54 sec
[bwa_aln_core] print alignments... 0.49 sec
[bwa_aln_core] 4194304 sequences have been processed.
[bwa_read_seq] 4.9% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 55.46 sec
[bwa_aln_core] refine gapped alignments... 9.01 sec
[bwa_aln_core] print alignments... 0.71 sec
[bwa_aln_core] 4456448 sequences have been processed.
[bwa_read_seq] 3.4% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 55.66 sec
[bwa_aln_core] refine gapped alignments... 9.24 sec
[bwa_aln_core] print alignments... 0.68 sec
[bwa_aln_core] 4718592 sequences have been processed.
[bwa_read_seq] 2.5% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 28.19 sec
[bwa_aln_core] refine gapped alignments... 0.41 sec
[bwa_aln_core] print alignments... 0.31 sec
[bwa_aln_core] 4980736 sequences have been processed.
[bwa_read_seq] 1.9% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 14.17 sec
[bwa_aln_core] refine gapped alignments... 0.50 sec
[bwa_aln_core] print alignments... 0.47 sec
[bwa_aln_core] 5242880 sequences have been processed.
[bwa_read_seq] 1.5% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 9.52 sec
[bwa_aln_core] refine gapped alignments... 0.44 sec
[bwa_aln_core] print alignments... 0.44 sec
[bwa_aln_core] 5505024 sequences have been processed.
[bwa_read_seq] 1.5% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 2.86 sec
[bwa_aln_core] refine gapped alignments... 0.39 sec
[bwa_aln_core] print alignments... 0.36 sec
[bwa_aln_core] 5767168 sequences have been processed.
[bwa_read_seq] 1.3% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 12.82 sec
[bwa_aln_core] refine gapped alignments... 0.46 sec
[bwa_aln_core] print alignments... 0.32 sec
[bwa_aln_core] 6029312 sequences have been processed.
[bwa_read_seq] 1.3% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 3.28 sec
[bwa_aln_core] refine gapped alignments... 0.43 sec
[bwa_aln_core] print alignments... 0.32 sec
[bwa_aln_core] 6291456 sequences have been processed.
[bwa_read_seq] 1.2% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 3.32 sec
[bwa_aln_core] refine gapped alignments... 0.48 sec
[bwa_aln_core] print alignments... 0.42 sec
[bwa_aln_core] 6553600 sequences have been processed.
[bwa_read_seq] 1.2% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 2.70 sec
[bwa_aln_core] refine gapped alignments... 0.40 sec
[bwa_aln_core] print alignments... 0.32 sec
[bwa_aln_core] 6815744 sequences have been processed.
[bwa_read_seq] 1.2% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 2.69 sec
[bwa_aln_core] refine gapped alignments... 0.39 sec
[bwa_aln_core] print alignments... 0.37 sec
[bwa_aln_core] 7077888 sequences have been processed.
[bwa_read_seq] 1.5% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 3.38 sec
[bwa_aln_core] refine gapped alignments... 0.47 sec
[bwa_aln_core] print alignments... 0.50 sec
[bwa_aln_core] 7340032 sequences have been processed.
[bwa_read_seq] 2.5% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 2.67 sec
[bwa_aln_core] refine gapped alignments... 0.39 sec
[bwa_aln_core] print alignments... 0.27 sec
[bwa_aln_core] 7602176 sequences have been processed.
[bwa_read_seq] 6.2% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 2.66 sec
[bwa_aln_core] refine gapped alignments... 0.40 sec
[bwa_aln_core] print alignments... 0.31 sec
[bwa_aln_core] 7864320 sequences have been processed.
[bwa_read_seq] 11.1% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 3.19 sec
[bwa_aln_core] refine gapped alignments... 0.40 sec
[bwa_aln_core] print alignments... 0.31 sec
[bwa_aln_core] 8126464 sequences have been processed.
[bwa_read_seq] 13.3% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 3.81 sec
[bwa_aln_core] refine gapped alignments... 0.40 sec
[bwa_aln_core] print alignments... 0.30 sec
[bwa_aln_core] 8388608 sequences have been processed.
[bwa_read_seq] 13.6% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 4.42 sec
[bwa_aln_core] refine gapped alignments... 0.40 sec
[bwa_aln_core] print alignments... 0.39 sec
[bwa_aln_core] 8650752 sequences have been processed.
[bwa_read_seq] 12.7% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 5.02 sec
[bwa_aln_core] refine gapped alignments... 0.39 sec
[bwa_aln_core] print alignments... 0.40 sec
[bwa_aln_core] 8912896 sequences have been processed.
[bwa_read_seq] 11.1% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 5.74 sec
[bwa_aln_core] refine gapped alignments... 0.39 sec
[bwa_aln_core] print alignments... 0.49 sec
[bwa_aln_core] 9175040 sequences have been processed.
[bwa_read_seq] 9.0% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 6.54 sec
[bwa_aln_core] refine gapped alignments... 0.39 sec
[bwa_aln_core] print alignments... 0.44 sec
[bwa_aln_core] 9437184 sequences have been processed.
[bwa_read_seq] 6.9% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 7.27 sec
[bwa_aln_core] refine gapped alignments... 0.39 sec
[bwa_aln_core] print alignments... 0.46 sec
[bwa_aln_core] 9699328 sequences have been processed.
[bwa_read_seq] 5.1% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 7.17 sec
[bwa_aln_core] refine gapped alignments... 0.39 sec
[bwa_aln_core] print alignments... 0.40 sec
[bwa_aln_core] 9961472 sequences have been processed.
[bwa_read_seq] 3.8% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 7.95 sec
[bwa_aln_core] refine gapped alignments... 0.39 sec
[bwa_aln_core] print alignments... 0.40 sec
[bwa_aln_core] 10223616 sequences have been processed.
[bwa_read_seq] 2.9% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 12.95 sec
[bwa_aln_core] refine gapped alignments... 0.39 sec
[bwa_aln_core] print alignments... 0.30 sec
[bwa_aln_core] 10485760 sequences have been processed.
[bwa_read_seq] 2.2% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 13.00 sec
[bwa_aln_core] refine gapped alignments... 0.38 sec
[bwa_aln_core] print alignments... 0.37 sec
[bwa_aln_core] 10747904 sequences have been processed.
[bwa_read_seq] 1.2% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 12.53 sec
[bwa_aln_core] refine gapped alignments... 0.39 sec
[bwa_aln_core] print alignments... 0.35 sec
[bwa_aln_core] 11010048 sequences have been processed.
[bwa_read_seq] 1.2% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 12.55 sec
[bwa_aln_core] refine gapped alignments... 0.38 sec
[bwa_aln_core] print alignments... 0.33 sec
[bwa_aln_core] 11272192 sequences have been processed.
[bwa_read_seq] 1.2% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 12.55 sec
[bwa_aln_core] refine gapped alignments... 0.38 sec
[bwa_aln_core] print alignments... 0.29 sec
[bwa_aln_core] 11534336 sequences have been processed.
[bwa_read_seq] 1.2% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 12.64 sec
[bwa_aln_core] refine gapped alignments... 0.38 sec
[bwa_aln_core] print alignments... 0.34 sec
[bwa_aln_core] 11796480 sequences have been processed.
[bwa_read_seq] 1.1% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 13.15 sec
[bwa_aln_core] refine gapped alignments... 0.38 sec
[bwa_aln_core] print alignments... 0.26 sec
[bwa_aln_core] 12058624 sequences have been processed.
[bwa_read_seq] 1.1% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 13.44 sec
[bwa_aln_core] refine gapped alignments... 0.39 sec
[bwa_aln_core] print alignments... 0.34 sec
[bwa_aln_core] 12320768 sequences have been processed.
[bwa_read_seq] 1.2% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 13.73 sec
[bwa_aln_core] refine gapped alignments... 0.38 sec
[bwa_aln_core] print alignments... 0.34 sec
[bwa_aln_core] 12582912 sequences have been processed.
[bwa_read_seq] 1.8% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 13.27 sec
[bwa_aln_core] refine gapped alignments... 0.38 sec
[bwa_aln_core] print alignments... 0.27 sec
[bwa_aln_core] 12845056 sequences have been processed.
[bwa_read_seq] 4.1% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 13.27 sec
[bwa_aln_core] refine gapped alignments... 0.38 sec
[bwa_aln_core] print alignments... 0.38 sec
[bwa_aln_core] 13107200 sequences have been processed.
[bwa_read_seq] 8.7% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 10.86 sec
[bwa_aln_core] refine gapped alignments... 0.38 sec
[bwa_aln_core] print alignments... 0.28 sec
[bwa_aln_core] 13369344 sequences have been processed.
[bwa_read_seq] 11.9% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 9.31 sec
[bwa_aln_core] refine gapped alignments... 0.39 sec
[bwa_aln_core] print alignments... 0.40 sec
[bwa_aln_core] 13631488 sequences have been processed.
[bwa_read_seq] 12.5% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 9.52 sec
[bwa_aln_core] refine gapped alignments... 0.39 sec
[bwa_aln_core] print alignments... 0.36 sec
[bwa_aln_core] 13893632 sequences have been processed.
[bwa_read_seq] 12.1% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 9.30 sec
[bwa_aln_core] refine gapped alignments... 0.39 sec
[bwa_aln_core] print alignments... 0.40 sec
[bwa_aln_core] 14155776 sequences have been processed.
[bwa_read_seq] 10.5% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 9.34 sec
[bwa_aln_core] refine gapped alignments... 0.39 sec
[bwa_aln_core] print alignments... 0.29 sec
[bwa_aln_core] 14417920 sequences have been processed.
[bwa_read_seq] 8.4% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 9.23 sec
[bwa_aln_core] refine gapped alignments... 0.41 sec
[bwa_aln_core] print alignments... 0.37 sec
[bwa_aln_core] 14680064 sequences have been processed.
[bwa_read_seq] 6.1% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 9.26 sec
[bwa_aln_core] refine gapped alignments... 0.38 sec
[bwa_aln_core] print alignments... 0.38 sec
[bwa_aln_core] 14942208 sequences have been processed.
[bwa_read_seq] 4.3% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 9.49 sec
[bwa_aln_core] refine gapped alignments... 0.38 sec
[bwa_aln_core] print alignments... 0.37 sec
[bwa_aln_core] 15204352 sequences have been processed.
[bwa_read_seq] 3.0% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 10.69 sec
[bwa_aln_core] refine gapped alignments... 0.38 sec
[bwa_aln_core] print alignments... 0.38 sec
[bwa_aln_core] 15466496 sequences have been processed.
[bwa_read_seq] 2.4% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 10.76 sec
[bwa_aln_core] refine gapped alignments... 0.39 sec
[bwa_aln_core] print alignments... 0.30 sec
[bwa_aln_core] 15728640 sequences have been processed.
[bwa_read_seq] 2.0% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 11.46 sec
[bwa_aln_core] refine gapped alignments... 0.38 sec
[bwa_aln_core] print alignments... 0.36 sec
[bwa_aln_core] 15990784 sequences have been processed.
[bwa_read_seq] 1.6% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 11.78 sec
[bwa_aln_core] refine gapped alignments... 0.38 sec
[bwa_aln_core] print alignments... 0.28 sec
[bwa_aln_core] 16252928 sequences have been processed.
[bwa_read_seq] 1.3% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 12.52 sec
[bwa_aln_core] refine gapped alignments... 0.38 sec
[bwa_aln_core] print alignments... 0.36 sec
[bwa_aln_core] 16515072 sequences have been processed.
[bwa_read_seq] 1.3% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 12.98 sec
[bwa_aln_core] refine gapped alignments... 0.38 sec
[bwa_aln_core] print alignments... 0.36 sec
[bwa_aln_core] 16777216 sequences have been processed.
[bwa_read_seq] 1.8% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 13.16 sec
[bwa_aln_core] refine gapped alignments... 0.38 sec
[bwa_aln_core] print alignments... 0.30 sec
[bwa_aln_core] 17039360 sequences have been processed.
[bwa_read_seq] 1.7% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 12.67 sec
[bwa_aln_core] refine gapped alignments... 0.38 sec
[bwa_aln_core] print alignments... 0.42 sec
[bwa_aln_core] 17301504 sequences have been processed.
[bwa_read_seq] 1.3% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 12.60 sec
[bwa_aln_core] refine gapped alignments... 0.38 sec
[bwa_aln_core] print alignments... 0.34 sec
[bwa_aln_core] 17563648 sequences have been processed.
[bwa_read_seq] 1.3% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 12.64 sec
[bwa_aln_core] refine gapped alignments... 0.38 sec
[bwa_aln_core] print alignments... 0.39 sec
[bwa_aln_core] 17825792 sequences have been processed.
[bwa_read_seq] 1.4% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 12.61 sec
[bwa_aln_core] refine gapped alignments... 0.38 sec
[bwa_aln_core] print alignments... 0.37 sec
[bwa_aln_core] 18087936 sequences have been processed.
[bwa_read_seq] 2.2% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 12.75 sec
[bwa_aln_core] refine gapped alignments... 0.39 sec
[bwa_aln_core] print alignments... 0.29 sec
[bwa_aln_core] 18350080 sequences have been processed.
[bwa_read_seq] 4.8% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 12.78 sec
[bwa_aln_core] refine gapped alignments... 0.39 sec
[bwa_aln_core] print alignments... 0.28 sec
[bwa_aln_core] 18612224 sequences have been processed.
[bwa_read_seq] 9.5% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 12.72 sec
[bwa_aln_core] refine gapped alignments... 0.39 sec
[bwa_aln_core] print alignments... 0.33 sec
[bwa_aln_core] 18874368 sequences have been processed.
[bwa_read_seq] 12.1% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 12.69 sec
[bwa_aln_core] refine gapped alignments... 0.39 sec
[bwa_aln_core] print alignments... 0.42 sec
[bwa_aln_core] 19136512 sequences have been processed.
[bwa_read_seq] 12.5% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 12.92 sec
[bwa_aln_core] refine gapped alignments... 0.39 sec
[bwa_aln_core] print alignments... 0.40 sec
[bwa_aln_core] 19398656 sequences have been processed.
[bwa_read_seq] 11.4% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 12.96 sec
[bwa_aln_core] refine gapped alignments... 0.39 sec
[bwa_aln_core] print alignments... 0.40 sec
[bwa_aln_core] 19660800 sequences have been processed.
[bwa_read_seq] 9.7% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 13.09 sec
[bwa_aln_core] refine gapped alignments... 0.39 sec
[bwa_aln_core] print alignments... 0.40 sec
[bwa_aln_core] 19922944 sequences have been processed.
[bwa_read_seq] 7.5% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 13.18 sec
[bwa_aln_core] refine gapped alignments... 0.39 sec
[bwa_aln_core] print alignments... 0.48 sec
[bwa_aln_core] 20185088 sequences have been processed.
[bwa_read_seq] 5.6% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 13.21 sec
[bwa_aln_core] refine gapped alignments... 0.39 sec
[bwa_aln_core] print alignments... 0.30 sec
[bwa_aln_core] 20447232 sequences have been processed.
[bwa_read_seq] 3.9% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 10.94 sec
[bwa_aln_core] refine gapped alignments... 0.39 sec
[bwa_aln_core] print alignments... 0.29 sec
[bwa_aln_core] 20709376 sequences have been processed.
[bwa_read_seq] 2.4% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 3.72 sec
[bwa_aln_core] refine gapped alignments... 0.38 sec
[bwa_aln_core] print alignments... 0.27 sec
[bwa_aln_core] 20971520 sequences have been processed.
[bwa_read_seq] 1.9% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 3.76 sec
[bwa_aln_core] refine gapped alignments... 0.39 sec
[bwa_aln_core] print alignments... 0.27 sec
[bwa_aln_core] 21233664 sequences have been processed.
[bwa_read_seq] 1.6% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 3.76 sec
[bwa_aln_core] refine gapped alignments... 0.38 sec
[bwa_aln_core] print alignments... 0.29 sec
[bwa_aln_core] 21495808 sequences have been processed.
[bwa_read_seq] 1.5% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 3.67 sec
[bwa_aln_core] refine gapped alignments... 0.38 sec
[bwa_aln_core] print alignments... 0.27 sec
[bwa_aln_core] 21757952 sequences have been processed.
[bwa_read_seq] 1.5% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 3.73 sec
[bwa_aln_core] refine gapped alignments... 0.38 sec
[bwa_aln_core] print alignments... 0.26 sec
[bwa_aln_core] 22020096 sequences have been processed.
[bwa_read_seq] 1.4% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 3.73 sec
[bwa_aln_core] refine gapped alignments... 0.38 sec
[bwa_aln_core] print alignments... 0.35 sec
[bwa_aln_core] 22282240 sequences have been processed.
[bwa_read_seq] 1.4% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 3.92 sec
[bwa_aln_core] refine gapped alignments... 0.38 sec
[bwa_aln_core] print alignments... 0.41 sec
[bwa_aln_core] 22544384 sequences have been processed.
[bwa_read_seq] 1.3% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 3.94 sec
[bwa_aln_core] refine gapped alignments... 0.38 sec
[bwa_aln_core] print alignments... 0.40 sec
[bwa_aln_core] 22806528 sequences have been processed.
[bwa_read_seq] 1.5% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 3.84 sec
[bwa_aln_core] refine gapped alignments... 0.38 sec
[bwa_aln_core] print alignments... 0.26 sec
[bwa_aln_core] 23068672 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.38 sec
[bwa_aln_core] print alignments... 0.28 sec
[bwa_aln_core] 23330816 sequences have been processed.
[bwa_read_seq] 2.9% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 3.88 sec
[bwa_aln_core] refine gapped alignments... 0.38 sec
[bwa_aln_core] print alignments... 0.28 sec
[bwa_aln_core] 23592960 sequences have been processed.
[bwa_read_seq] 6.8% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 3.84 sec
[bwa_aln_core] refine gapped alignments... 0.38 sec
[bwa_aln_core] print alignments... 0.30 sec
[bwa_aln_core] 23855104 sequences have been processed.
[bwa_read_seq] 11.4% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 3.92 sec
[bwa_aln_core] refine gapped alignments... 0.44 sec
[bwa_aln_core] print alignments... 0.53 sec
[bwa_aln_core] 24117248 sequences have been processed.
[bwa_read_seq] 13.5% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 3.88 sec
[bwa_aln_core] refine gapped alignments... 0.39 sec
[bwa_aln_core] print alignments... 0.37 sec
[bwa_aln_core] 24379392 sequences have been processed.
[bwa_read_seq] 14.4% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 3.91 sec
[bwa_aln_core] refine gapped alignments... 0.39 sec
[bwa_aln_core] print alignments... 0.47 sec
[bwa_aln_core] 24641536 sequences have been processed.
[bwa_read_seq] 13.5% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 3.85 sec
[bwa_aln_core] refine gapped alignments... 0.39 sec
[bwa_aln_core] print alignments... 0.45 sec
[bwa_aln_core] 24903680 sequences have been processed.
[bwa_read_seq] 11.4% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 3.99 sec
[bwa_aln_core] refine gapped alignments... 0.38 sec
[bwa_aln_core] print alignments... 0.50 sec
[bwa_aln_core] 25165824 sequences have been processed.
[bwa_read_seq] 9.6% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 4.25 sec
[bwa_aln_core] refine gapped alignments... 0.38 sec
[bwa_aln_core] print alignments... 0.44 sec
[bwa_aln_core] 25427968 sequences have been processed.
[bwa_read_seq] 7.5% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 4.21 sec
[bwa_aln_core] refine gapped alignments... 0.38 sec
[bwa_aln_core] print alignments... 0.44 sec
[bwa_aln_core] 25690112 sequences have been processed.
[bwa_read_seq] 5.5% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 4.24 sec
[bwa_aln_core] refine gapped alignments... 0.38 sec
[bwa_aln_core] print alignments... 0.40 sec
[bwa_aln_core] 25952256 sequences have been processed.
[bwa_read_seq] 3.7% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 4.35 sec
[bwa_aln_core] refine gapped alignments... 0.38 sec
[bwa_aln_core] print alignments... 0.44 sec
[bwa_aln_core] 26214400 sequences have been processed.
[bwa_read_seq] 2.9% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 4.28 sec
[bwa_aln_core] refine gapped alignments... 0.38 sec
[bwa_aln_core] print alignments... 0.41 sec
[bwa_aln_core] 26476544 sequences have been processed.
[bwa_read_seq] 2.1% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 4.38 sec
[bwa_aln_core] refine gapped alignments... 0.38 sec
[bwa_aln_core] print alignments... 0.42 sec
[bwa_aln_core] 26738688 sequences have been processed.
[bwa_read_seq] 1.3% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 4.33 sec
[bwa_aln_core] refine gapped alignments... 0.37 sec
[bwa_aln_core] print alignments... 0.35 sec
[bwa_aln_core] 27000832 sequences have been processed.
[bwa_read_seq] 1.4% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 4.29 sec
[bwa_aln_core] refine gapped alignments... 0.38 sec
[bwa_aln_core] print alignments... 0.35 sec
[bwa_aln_core] 27262976 sequences have been processed.
[bwa_read_seq] 1.3% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 4.32 sec
[bwa_aln_core] refine gapped alignments... 0.38 sec
[bwa_aln_core] print alignments... 0.29 sec
[bwa_aln_core] 27525120 sequences have been processed.
[bwa_read_seq] 1.2% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 4.27 sec
[bwa_aln_core] refine gapped alignments... 0.38 sec
[bwa_aln_core] print alignments... 0.28 sec
[bwa_aln_core] 27787264 sequences have been processed.
[bwa_read_seq] 1.3% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 4.21 sec
[bwa_aln_core] refine gapped alignments... 0.37 sec
[bwa_aln_core] print alignments... 0.33 sec
[bwa_aln_core] 28049408 sequences have been processed.
[bwa_read_seq] 1.2% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 4.32 sec
[bwa_aln_core] refine gapped alignments... 0.37 sec
[bwa_aln_core] print alignments... 0.37 sec
[bwa_aln_core] 28311552 sequences have been processed.
[bwa_read_seq] 1.4% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 4.36 sec
[bwa_aln_core] refine gapped alignments... 0.37 sec
[bwa_aln_core] print alignments... 0.36 sec
[bwa_aln_core] 28573696 sequences have been processed.
[bwa_read_seq] 2.2% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 4.37 sec
[bwa_aln_core] refine gapped alignments... 0.42 sec
[bwa_aln_core] print alignments... 0.49 sec
[bwa_aln_core] 28835840 sequences have been processed.
[bwa_read_seq] 5.0% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 4.47 sec
[bwa_aln_core] refine gapped alignments... 0.39 sec
[bwa_aln_core] print alignments... 0.38 sec
[bwa_aln_core] 29097984 sequences have been processed.
[bwa_read_seq] 9.8% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 4.30 sec
[bwa_aln_core] refine gapped alignments... 0.38 sec
[bwa_aln_core] print alignments... 0.29 sec
[bwa_aln_core] 29360128 sequences have been processed.
[bwa_read_seq] 12.2% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 4.30 sec
[bwa_aln_core] refine gapped alignments... 0.39 sec
[bwa_aln_core] print alignments... 0.32 sec
[bwa_aln_core] 29622272 sequences have been processed.
[bwa_read_seq] 12.7% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 4.32 sec
[bwa_aln_core] refine gapped alignments... 0.38 sec
[bwa_aln_core] print alignments... 0.30 sec
[bwa_aln_core] 29884416 sequences have been processed.
[bwa_read_seq] 11.6% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 4.31 sec
[bwa_aln_core] refine gapped alignments... 0.43 sec
[bwa_aln_core] print alignments... 0.53 sec
[bwa_aln_core] 30146560 sequences have been processed.
[bwa_read_seq] 9.8% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 4.40 sec
[bwa_aln_core] refine gapped alignments... 0.38 sec
[bwa_aln_core] print alignments... 0.29 sec
[bwa_aln_core] 30408704 sequences have been processed.
[bwa_read_seq] 7.7% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 4.37 sec
[bwa_aln_core] refine gapped alignments... 0.38 sec
[bwa_aln_core] print alignments... 0.38 sec
[bwa_aln_core] 30670848 sequences have been processed.
[bwa_read_seq] 5.8% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 5.42 sec
[bwa_aln_core] refine gapped alignments... 0.43 sec
[bwa_aln_core] print alignments... 0.52 sec
[bwa_aln_core] 30932992 sequences have been processed.
[bwa_read_seq] 4.4% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 4.83 sec
[bwa_aln_core] refine gapped alignments... 0.38 sec
[bwa_aln_core] print alignments... 0.35 sec
[bwa_aln_core] 31195136 sequences have been processed.
[bwa_read_seq] 3.0% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 4.86 sec
[bwa_aln_core] refine gapped alignments... 0.38 sec
[bwa_aln_core] print alignments... 0.39 sec
[bwa_aln_core] 31457280 sequences have been processed.
[bwa_read_seq] 2.4% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 4.97 sec
[bwa_aln_core] refine gapped alignments... 0.38 sec
[bwa_aln_core] print alignments... 0.28 sec
[bwa_aln_core] 31719424 sequences have been processed.
[bwa_read_seq] 1.9% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 4.96 sec
[bwa_aln_core] refine gapped alignments... 0.38 sec
[bwa_aln_core] print alignments... 0.41 sec
[bwa_aln_core] 31981568 sequences have been processed.
[bwa_read_seq] 2.0% bases are trimmed.
[bwa_aln_core] convert to sequence coordinate... 4.24 sec
[bwa_aln_core] refine gapped alignments... 0.30 sec
[bwa_aln_core] print alignments... 0.05 sec
[bwa_aln_core] 32022635 sequences have been processed.
[main] Version: 0.7.13-r1126
[main] CMD: bwa samse /mnt/data/bds_pipeline_genome_data/mm10/bwa_index/GRCm38.p4.genome.fa /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/align/ctl1/r120621A_l1i7_aB180__wt__bioBirAonly.sai /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/InputData/r120621A_l1i7_aB180__wt__bioBirAonly.fastq.gz
[main] Real time: 2022.299 sec; CPU: 1937.196 sec

Num	3
ID	task.postalign_bam.dedup_bam_1_ctl1.line_95.id_12
Name	dedup_bam_1 ctl1
Thread	thread_42
PID	9846

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

Start	2016-09-29 17:25:22
End	2016-09-29 17:29:32
Elapsed	00:04:10
Timeout	00:00:-1
Wall Timeout	100 days

Input files	/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/align/ctl1/r120621A_l1i7_aB180__wt__bioBirAonly.bam
Output files	/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/align/ctl1/r120621A_l1i7_aB180__wt__bioBirAonly.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 2 /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/align/ctl1/r120621A_l1i7_aB180__wt__bioBirAonly.bam -n -o /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/align/ctl1/r120621A_l1i7_aB180__wt__bioBirAonly.qnmsrt.bam; \
				samtools view -h /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/align/ctl1/r120621A_l1i7_aB180__wt__bioBirAonly.qnmsrt.bam | $(which assign_multimappers.py) -k 0 | \
				samtools view -F 1804 -Su /dev/stdin | \
				sambamba sort -t 2 /dev/stdin -o /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/align/ctl1/r120621A_l1i7_aB180__wt__bioBirAonly.filt.bam; \
				rm -f /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/align/ctl1/r120621A_l1i7_aB180__wt__bioBirAonly.qnmsrt.bam; \
			    else \
			    	samtools view -F 1804 -q 30 -u /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/align/ctl1/r120621A_l1i7_aB180__wt__bioBirAonly.bam | \
			    	sambamba sort -t 2 /dev/stdin -o /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/align/ctl1/r120621A_l1i7_aB180__wt__bioBirAonly.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--------------------
46668 (process ID) old priority 0, new priority 10

 
--------------------Stderr--------------------
discarding /software/miniconda3/bin from PATH
prepending /software/miniconda3/envs/aquas_chipseq/bin to PATH

Num	4
ID	task.postalign_bam.markdup_bam_picard_ctl1.line_347.id_13
Name	markdup_bam_picard ctl1
Thread	thread_42
PID	9851

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

Start	2016-09-29 17:29:33
End	2016-09-29 17:33:15
Elapsed	00:03:41
Timeout	00:00:-1
Wall Timeout	100 days

Input files	/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/align/ctl1/r120621A_l1i7_aB180__wt__bioBirAonly.filt.bam
Output files	/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/align/ctl1/r120621A_l1i7_aB180__wt__bioBirAonly.dupmark.bam /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/qc/ctl1/r120621A_l1i7_aB180__wt__bioBirAonly.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/MouseMutantData/ZF11/out/align/ctl1/r120621A_l1i7_aB180__wt__bioBirAonly.filt.bam" OUTPUT="/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/align/ctl1/r120621A_l1i7_aB180__wt__bioBirAonly.dupmark.bam" \
				METRICS_FILE="/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/qc/ctl1/r120621A_l1i7_aB180__wt__bioBirAonly.dup.qc" VALIDATION_STRINGENCY=LENIENT \
				ASSUME_SORTED=true REMOVE_DUPLICATES=false; \
		    else \
			picard MarkDuplicates \
				INPUT="/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/align/ctl1/r120621A_l1i7_aB180__wt__bioBirAonly.filt.bam" OUTPUT="/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/align/ctl1/r120621A_l1i7_aB180__wt__bioBirAonly.dupmark.bam" \
				METRICS_FILE="/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/qc/ctl1/r120621A_l1i7_aB180__wt__bioBirAonly.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--------------------
25052 (process ID) old priority 0, new priority 10

 
--------------------Stderr--------------------
discarding /software/miniconda3/bin from PATH
prepending /software/miniconda3/envs/aquas_chipseq/bin to PATH
Picked up _JAVA_OPTIONS: -Xms256M -Xmx4G -XX:ParallelGCThreads=1
[Thu Sep 29 17:33:54 PDT 2016] net.sf.picard.sam.MarkDuplicates INPUT=[/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/align/ctl1/r120621A_l1i7_aB180__wt__bioBirAonly.filt.bam] OUTPUT=/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/align/ctl1/r120621A_l1i7_aB180__wt__bioBirAonly.dupmark.bam METRICS_FILE=/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/qc/ctl1/r120621A_l1i7_aB180__wt__bioBirAonly.dup.qc REMOVE_DUPLICATES=false ASSUME_SORTED=true VALIDATION_STRINGENCY=LENIENT    PROGRAM_RECORD_ID=MarkDuplicates PROGRAM_GROUP_NAME=MarkDuplicates MAX_SEQUENCES_FOR_DISK_READ_ENDS_MAP=50000 MAX_FILE_HANDLES_FOR_READ_ENDS_MAP=8000 SORTING_COLLECTION_SIZE_RATIO=0.25 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
[Thu Sep 29 17:33:54 PDT 2016] Executing as imk1@kadru on Linux 4.4.0-22-generic amd64; OpenJDK 64-Bit Server VM 1.6.0-77-b77; Picard version: 1.97(Unversioned directory)
INFO	2016-09-29 17:33:54	MarkDuplicates	Start of doWork freeMemory: 255389744; totalMemory: 257294336; maxMemory: 3817799680
INFO	2016-09-29 17:33:54	MarkDuplicates	Reading input file and constructing read end information.
INFO	2016-09-29 17:33:54	MarkDuplicates	Will retain up to 15149998 data points before spilling to disk.
INFO	2016-09-29 17:33:59	MarkDuplicates	Read     1,000,000 records.  Elapsed time: 00:00:05s.  Time for last 1,000,000:    5s.  Last read position: chr1:169,651,309
INFO	2016-09-29 17:33:59	MarkDuplicates	Tracking 0 as yet unmatched pairs. 0 records in RAM.
INFO	2016-09-29 17:34:02	MarkDuplicates	Read     2,000,000 records.  Elapsed time: 00:00:08s.  Time for last 1,000,000:    3s.  Last read position: chr2:120,600,340
INFO	2016-09-29 17:34:02	MarkDuplicates	Tracking 0 as yet unmatched pairs. 0 records in RAM.
INFO	2016-09-29 17:34:06	MarkDuplicates	Read     3,000,000 records.  Elapsed time: 00:00:11s.  Time for last 1,000,000:    3s.  Last read position: chr3:96,659,413
INFO	2016-09-29 17:34:06	MarkDuplicates	Tracking 0 as yet unmatched pairs. 0 records in RAM.
INFO	2016-09-29 17:34:11	MarkDuplicates	Read     4,000,000 records.  Elapsed time: 00:00:16s.  Time for last 1,000,000:    4s.  Last read position: chr4:108,195,091
INFO	2016-09-29 17:34:11	MarkDuplicates	Tracking 0 as yet unmatched pairs. 0 records in RAM.
INFO	2016-09-29 17:34:16	MarkDuplicates	Read     5,000,000 records.  Elapsed time: 00:00:22s.  Time for last 1,000,000:    5s.  Last read position: chr5:103,594,339
INFO	2016-09-29 17:34:16	MarkDuplicates	Tracking 0 as yet unmatched pairs. 0 records in RAM.
INFO	2016-09-29 17:34:20	MarkDuplicates	Read     6,000,000 records.  Elapsed time: 00:00:26s.  Time for last 1,000,000:    4s.  Last read position: chr6:98,592,333
INFO	2016-09-29 17:34:20	MarkDuplicates	Tracking 0 as yet unmatched pairs. 0 records in RAM.
INFO	2016-09-29 17:34:23	MarkDuplicates	Read     7,000,000 records.  Elapsed time: 00:00:29s.  Time for last 1,000,000:    2s.  Last read position: chr7:114,109,505
INFO	2016-09-29 17:34:23	MarkDuplicates	Tracking 0 as yet unmatched pairs. 0 records in RAM.
INFO	2016-09-29 17:34:26	MarkDuplicates	Read     8,000,000 records.  Elapsed time: 00:00:32s.  Time for last 1,000,000:    3s.  Last read position: chr8:121,515,352
INFO	2016-09-29 17:34:26	MarkDuplicates	Tracking 0 as yet unmatched pairs. 0 records in RAM.
INFO	2016-09-29 17:34:29	MarkDuplicates	Read     9,000,000 records.  Elapsed time: 00:00:34s.  Time for last 1,000,000:    2s.  Last read position: chr10:14,536,964
INFO	2016-09-29 17:34:29	MarkDuplicates	Tracking 0 as yet unmatched pairs. 0 records in RAM.
INFO	2016-09-29 17:34:33	MarkDuplicates	Read    10,000,000 records.  Elapsed time: 00:00:38s.  Time for last 1,000,000:    3s.  Last read position: chr11:32,924,585
INFO	2016-09-29 17:34:33	MarkDuplicates	Tracking 0 as yet unmatched pairs. 0 records in RAM.
INFO	2016-09-29 17:34:36	MarkDuplicates	Read    11,000,000 records.  Elapsed time: 00:00:41s.  Time for last 1,000,000:    2s.  Last read position: chr12:41,050,039
INFO	2016-09-29 17:34:36	MarkDuplicates	Tracking 0 as yet unmatched pairs. 0 records in RAM.
INFO	2016-09-29 17:34:38	MarkDuplicates	Read    12,000,000 records.  Elapsed time: 00:00:43s.  Time for last 1,000,000:    2s.  Last read position: chr13:80,878,598
INFO	2016-09-29 17:34:38	MarkDuplicates	Tracking 0 as yet unmatched pairs. 0 records in RAM.
INFO	2016-09-29 17:34:44	MarkDuplicates	Read    13,000,000 records.  Elapsed time: 00:00:49s.  Time for last 1,000,000:    5s.  Last read position: chr15:12,027,335
INFO	2016-09-29 17:34:44	MarkDuplicates	Tracking 0 as yet unmatched pairs. 0 records in RAM.
INFO	2016-09-29 17:34:46	MarkDuplicates	Read    14,000,000 records.  Elapsed time: 00:00:51s.  Time for last 1,000,000:    2s.  Last read position: chr16:56,024,589
INFO	2016-09-29 17:34:46	MarkDuplicates	Tracking 0 as yet unmatched pairs. 0 records in RAM.
INFO	2016-09-29 17:34:49	MarkDuplicates	Read    15,000,000 records.  Elapsed time: 00:00:54s.  Time for last 1,000,000:    3s.  Last read position: chr18:10,030,498
INFO	2016-09-29 17:34:49	MarkDuplicates	Tracking 0 as yet unmatched pairs. 0 records in RAM.
INFO	2016-09-29 17:34:58	MarkDuplicates	Read    16,000,000 records.  Elapsed time: 00:01:03s.  Time for last 1,000,000:    9s.  Last read position: chrX:36,565,109
INFO	2016-09-29 17:34:58	MarkDuplicates	Tracking 0 as yet unmatched pairs. 0 records in RAM.
INFO	2016-09-29 17:34:59	MarkDuplicates	Read 16572675 records. 0 pairs never matched.
INFO	2016-09-29 17:34:59	MarkDuplicates	After buildSortedReadEndLists freeMemory: 2535064936; totalMemory: 2622160896; maxMemory: 3817799680
INFO	2016-09-29 17:34:59	MarkDuplicates	Will retain up to 119306240 duplicate indices before spilling to disk.
INFO	2016-09-29 17:35:03	MarkDuplicates	Traversing read pair information and detecting duplicates.
INFO	2016-09-29 17:35:03	MarkDuplicates	Traversing fragment information and detecting duplicates.
INFO	2016-09-29 17:35:07	MarkDuplicates	Sorting list of duplicate records.
INFO	2016-09-29 17:35:07	MarkDuplicates	After generateDuplicateIndexes freeMemory: 2019147976; totalMemory: 2995716096; maxMemory: 3817799680
INFO	2016-09-29 17:35:07	MarkDuplicates	Marking 13188493 records as duplicates.
INFO	2016-09-29 17:35:07	MarkDuplicates	Found 0 optical duplicate clusters.
INFO	2016-09-29 17:36:30	MarkDuplicates	Written    10,000,000 records.  Elapsed time: 00:01:22s.  Time for last 10,000,000:   82s.  Last read position: chr11:32,924,585
INFO	2016-09-29 17:37:25	MarkDuplicates	Before output close freeMemory: 2984774448; totalMemory: 3012886528; maxMemory: 3817799680
INFO	2016-09-29 17:37:26	MarkDuplicates	After output close freeMemory: 2960918504; totalMemory: 2989031424; maxMemory: 3817799680
[Thu Sep 29 17:37:26 PDT 2016] net.sf.picard.sam.MarkDuplicates done. Elapsed time: 3.53 minutes.
Runtime.totalMemory()=2989031424

Num	5
ID	task.postalign_bam.dedup_bam_2_ctl1.line_144.id_14
Name	dedup_bam_2 ctl1
Thread	thread_42
PID	9855

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

Start	2016-09-29 17:33:17
End	2016-09-29 17:36:32
Elapsed	00:03:14
Timeout	00:00:-1
Wall Timeout	100 days

Input files	/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/align/ctl1/r120621A_l1i7_aB180__wt__bioBirAonly.dupmark.bam
Output files	/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/align/ctl1/r120621A_l1i7_aB180__wt__bioBirAonly.nodup.bam /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/qc/ctl1/r120621A_l1i7_aB180__wt__bioBirAonly.nodup.flagstat.qc /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/qc/ctl1/r120621A_l1i7_aB180__wt__bioBirAonly.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; 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 150

 samtools view -F 1804 -b /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/align/ctl1/r120621A_l1i7_aB180__wt__bioBirAonly.dupmark.bam > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/align/ctl1/r120621A_l1i7_aB180__wt__bioBirAonly.nodup.bam

# SYS command. line 153

 sambamba index -t 2 /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/align/ctl1/r120621A_l1i7_aB180__wt__bioBirAonly.nodup.bam

# SYS command. line 155

 sambamba flagstat -t 2 /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/align/ctl1/r120621A_l1i7_aB180__wt__bioBirAonly.nodup.bam > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/qc/ctl1/r120621A_l1i7_aB180__wt__bioBirAonly.nodup.flagstat.qc

# SYS command. line 165

 bedtools bamtobed -i /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/align/ctl1/r120621A_l1i7_aB180__wt__bioBirAonly.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/MouseMutantData/ZF11/out/qc/ctl1/r120621A_l1i7_aB180__wt__bioBirAonly.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--------------------
47025 (process ID) old priority 0, new priority 10

 
--------------------Stderr--------------------
discarding /software/miniconda3/bin from PATH
prepending /software/miniconda3/envs/aquas_chipseq/bin to PATH

Num	6
ID	task.postalign_bam.bam_to_tag_ctl1.line_504.id_15
Name	bam_to_tag ctl1
Thread	thread_42
PID	9859

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

Start	2016-09-29 17:36:33
End	2016-09-29 17:37:43
Elapsed	00:01:09
Timeout	00:00:-1
Wall Timeout	100 days

Input files	/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/align/ctl1/r120621A_l1i7_aB180__wt__bioBirAonly.nodup.bam
Output files	/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/align/ctl1/r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign.gz
Dependencies

 
# SYS command. line 506

 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 509

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

# SYS command. line 511

 echo

# SYS command. line 513

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

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

Waiting for 50 seconds.

 
--------------------Stderr--------------------
discarding /software/miniconda3/bin from PATH
prepending /software/miniconda3/envs/aquas_chipseq/bin to PATH

Num	7
ID	task.callpeak_spp.spp_rep1.line_54.id_16
Name	spp rep1
Thread	thread_Root
PID	9860

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

Start	2016-09-29 17:37:44
End	2016-09-30 06:00:14
Elapsed	12:22:29
Timeout	00:00:-1
Wall Timeout	100 days

Input files	/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/align/rep1/ZF11_allReps.nodup.tagAlign.gz /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/align/ctl1/r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign.gz
Output files	/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/peak/spp/rep1/ZF11_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign.regionPeak.gz /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/peak/spp/rep1/ZF11_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign.ccscore /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/peak/spp/rep1/ZF11_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign.pdf
Dependencies

 
# SYS command. line 56

 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 59

 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 63

 Rscript ${RUN_SPP} -c=/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/align/rep1/ZF11_allReps.nodup.tagAlign.gz -p=1 -i=/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/align/ctl1/r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign.gz \
			-npeak=300000 -odir=/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/peak/spp/rep1 -speak=125 -savr -savp -rf -out=/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/peak/spp/rep1/ZF11_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign.ccscore

# SYS command. line 67

 zcat /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/peak/spp/rep1/ZF11_allReps.nodup.tagAlign_VS_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign.regionPeak.gz | awk 'BEGIN{OFS="\t"}{ if ($2<0) $2=0; print $1,int($2),int($3),$4,$5,$6,$7,$8,$9,$10;}' | gzip -f -c > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/peak/spp/rep1/ZF11_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign.regionPeak.gz

# SYS command. line 69

 rm -f /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/peak/spp/rep1/ZF11_allReps.nodup.tagAlign_VS_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign.regionPeak.gz

# SYS command. line 71

 mv /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/peak/spp/rep1/ZF11_allReps.nodup.tagAlign.pdf /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/peak/spp/rep1/ZF11_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign.pdf

# SYS command. line 74

 if [ $(zcat /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/peak/spp/rep1/ZF11_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign.regionPeak.gz | wc -l ) == "0" ]; then rm -f /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/peak/spp/rep1/ZF11_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign.regionPeak.gz; fi

# SYS command. line 77

 if [ ! -f /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/peak/spp/rep1/ZF11_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign.regionPeak.gz ]; then error_in_spp_output_peak_does_not_exist; fi

# SYS command. line 79

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

   
--------------------Stdout--------------------
25325 (process ID) old priority 0, new priority 10
################
ChIP data: /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/align/rep1/ZF11_allReps.nodup.tagAlign.gz 
Control data: /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/align/ctl1/r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign.gz 
strandshift(min): -500 
strandshift(step): 5 
strandshift(max) 1500 
user-defined peak shift 125 
exclusion(min): 10 
exclusion(max): NaN 
num parallel nodes: 1 
FDR threshold: 0.01 
NumPeaks Threshold: 3e+05 
Output Directory: /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/peak/spp/rep1 
narrowPeak output file name: NA 
regionPeak output file name: /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/peak/spp/rep1/ZF11_allReps.nodup.tagAlign_VS_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign.regionPeak 
Rdata filename: NA 
plot pdf filename: /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/peak/spp/rep1/ZF11_allReps.nodup.tagAlign.pdf 
result filename: /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/peak/spp/rep1/ZF11_allReps.nodup.tagAlign_x_r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign.ccscore 
Overwrite files?: TRUE

Decompressing ChIP file
Decompressing control file
Reading ChIP tagAlign/BAM file /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/align/rep1/ZF11_allReps.nodup.tagAlign.gz 
opened /tmp/9860.1.q/RtmpaqF0c9/ZF11_allReps.nodup.tagAlign632f5d1e4570
done. read 30207286 fragments
ChIP data read length 50 
[1] TRUE
Reading Control tagAlign/BAM file /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/align/ctl1/r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign.gz 
opened /tmp/9860.1.q/RtmpaqF0c9/r120621A_l1i7_aB180__wt__bioBirAonly.nodup.tagAlign632f22f19de
done. read 3384182 fragments
Control data read length 50 
Calculating peak characteristics
Minimum cross-correlation value 0.3042838 
Minimum cross-correlation shift 1500 
Top 3 cross-correlation values 0.450586479744112 
Top 3 estimates for fragment length 125 
Window half size 375 
Phantom peak location 55 
Phantom peak Correlation 0.4259703 
Normalized Strand cross-correlation coefficient (NSC) 1.48081 
Relative Strand Cross correlation Coefficient (RSC) 1.202292 
Phantom Peak Quality Tag 1 
null device 
          1 
Removing read stacks
Finding peaks
finding background exclusion regions ... done
determining peaks on provided 1 control datasets:
using reversed signal for FDR calculations
bg.weight= 0.1456462  excluding systematic background anomalies ... done
determining peaks on real data:
bg.weight= 6.865953  excluding systematic background anomalies ... done
calculating statistical thresholds
FDR 0.99 threshold= 2.000083 
Detected 1433504 peaks 

 
--------------------Stderr--------------------
discarding /software/miniconda3/bin from PATH
prepending /software/miniconda3/envs/aquas_chipseq/bin to PATH
Loading required package: caTools
There were 30 warnings (use warnings() to see them)

Num	8
ID	task.graphviz.report.line_97.id_17
Name	report
Thread	thread_Root
PID	43042

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

Start	2016-09-30 06:00:15
End	2016-09-30 06:00:17
Elapsed	00:00:01
Timeout	00:00:-1
Wall Timeout	100 days

Input files
Output files
Dependencies

 
# SYS command. line 98

 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 99

 dot -Tsvg /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/report/ZF11_workflow.dot > /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/report/ZF11_workflow.svg 2> /dev/null || echo "svg: /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out/report/ZF11_workflow.svg" # to suppress dot font error (exit code=1)

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

 
--------------------Stderr--------------------
discarding /software/miniconda3/bin from PATH
prepending /software/miniconda3/envs/aquas_chipseq/bin to PATH

Type	Name	Value
string	_	/usr/bin/bds
bool	allowEmpty	false
string[]	args	[-title, ZF11, -nth, 5, -no_pseudo_rep, -macs2_for_pooled_rep_only, -species, mm10, -url_base, http://mitra.stanford.edu/kundaje/imk1/TFBindingPredictionProject/MouseMutantData/ZF11/out, -fastq, /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/RawReads/ZF11_allReps.fastq.gz, -ctl_fastq, /srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/InputData/r120621A_l1i7_aB180__wt__bioBirAonly.fastq.gz, $]
bool	canFail	false
int	cpus	-1
int	cpusLocal	56
int	day	86400
string	DEEPLIFT_DIR	/srv/scratch/imk1/TFBindingPredictionProject/src/deeplift_private/deeplift
string	DISPLAY	localhost:19.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
int	K	1024
string	KERAS_DIR	/srv/scratch/imk1/TFBindingPredictionProject/src/keras
string	KRB5CCNAME	FILE:/tmp/krb5cc_1048_eHVWt0
string	LANG	en_US.UTF-8
string	LOADEDMODULES
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	MODULESHOME	/software/env_module/3.2.10
string	NLSPATH	/usr/dt/lib/nls/msg/%L/%N.cat
string	node
int	P	1125899906842624
string	PATH	/users/imk1/perl5/bin:/srv/scratch/imk1/TFBindingPredictionProject/src/bedtools2/bin:/usr/local/cuda/bin:/srv/persistent/imk1/TFBindingPredictionProject/src/rcade/:/usr/local/sbin:/usr/local/bin:/usr/sbin:/usr/bin:/sbin:/bin:/usr/games:/usr/local/games:/software/miniconda3/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/MouseMutantData/ZF11
string	programName	chipseq.bds
string	programPath	/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/src/TF_chipseq_pipeline/chipseq.bds
string	PWD	/srv/scratch/shared/surya/imk1/TFBindingPredictionProject/MouseMutantData/ZF11
string	PYTHONPATH	/software/python_anaconda/2.2.0/bin/python:/srv/scratch/imk1/TFBindingPredictionProject/src/pybedtools/
string	queue
int	retry	0
string	RULEFITBASE	/srv/scratch/imk1/TFBindingPredictionProject/src/RuleFit/
string	SGE_CELL	default
string	SGE_ROOT	/var/lib/gridengine
string	SHELL	/bin/bash
string	SHLVL	3
string	SSH_CLIENT	171.65.77.8 58475 22
string	SSH_CONNECTION	171.65.77.8 58475 171.65.76.63 22
string	SSH_TTY	/dev/pts/38
string	STY	38311.ZF11.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	418
string	XFILESEARCHPATH	/usr/dt/app-defaults/%L/Dt

BigDataScript report: chipseq.bds

Timeline

Task Graph

Parallel (threads) details

Task details

Global scope