-----------------------
Start of variable list:
-----------------------
Dataset_DIR:
/restricted/projectnb/waxmanlab/kkarri/G160_Samples
Dataset_Label:
G160
bin_mode:
Off
SCRIPT_DIR:
/restricted/projectnb/waxmanlab/kkarri/G160_Samples/Scripts/10_bamCorrelate
OUTPUT_DIR:
/restricted/projectnb/waxmanlab/kkarri/G160_Samples/Scripts/10_bamCorrelate/Output
JOB_COUNTER:
03
zMin:
0.6
zMax:
1.0
BED_file_name:
RefSeq_GeneBody
-----------------------
End of variable list
-----------------------
==========================================================
Starting on : Thu Aug 24 07:49:22 EDT 2017
Running on node : scc-pf2
Current directory : /restricted/projectnb/waxmanlab/kkarri/G160_Samples/Scripts/10_bamCorrelate
Current job ID : 383585
Current job name : Step_10_bamCor_03
Task index number : undefined
Parameter for multiple cores : 16
==========================================================

Change dir to scratch directory


Print scratch directory location:

/scratch/383585.1.p8

Loading required modules...

Getting BED file line count...
Before ENCODE_Blacklist filter:
665744
Overlap BED file with mm9-blacklist.bed and filter out overlap
RefSeq_GeneBody
mm9-blacklist
----------------------------------------------------------------------------
cleaning input
sorting input files
classify common or unique peaks
getting the union of the peak sets
getting the merged common peaks
getting distance of nearest peak
counting peaks
Running overlap.R script to get peak overlap, peak width, and peak proximity statistics ...
Creating summary file
Done!
----------------------------------------------------------------------------

Getting BED file line count...
After ENCODE_Blacklist filter:
660786

Copy input sample BAM files to Input folder


Save Input file lists to a variable

Input_List:

Labels_List:


List files in scratch directory:

total 37M
drwxr-xr-x   4 kkarri pulmseq 4.0K Aug 24 07:49 .
drwxrwxrwt. 72 root   root    4.0K Aug 24 07:49 ..
drwxr-xr-x   2 kkarri pulmseq 4.0K Aug 24 07:49 Input
-rwxr-xr-x   1 kkarri pulmseq  205 Aug 24 07:49 Input_Samples.txt
-rwxr-xr-x   1 kkarri pulmseq  36M Aug 24 07:49 RefSeq_GeneBody.bed
drwxr-xr-x   2 kkarri pulmseq 4.0K Aug 24 07:49 RefSeq_GeneBody_mm9-blacklist_Output
-rwxr-xr-x   1 kkarri pulmseq  71K Aug 24 07:49 mm9-blacklist.bed
-rwxr-xr-x   1 kkarri pulmseq 7.6K Aug 24 07:49 overlap.R
-rwxr-xr-x   1 kkarri pulmseq 6.0K Aug 24 07:49 overlap.sh
-rw-r--r--   1 kkarri pulmseq  80K Aug 24 07:49 peak1_dump.bed

Starting to run my commands


Starting bamCorrelate command


Ending bamCorrelate command


Starting plotCorrelation command


Ending plotCorrelation command


List files:

total 37M
drwxr-xr-x   4 kkarri pulmseq 4.0K Aug 24 07:49 .
drwxrwxrwt. 72 root   root    4.0K Aug 24 07:49 ..
drwxr-xr-x   2 kkarri pulmseq 4.0K Aug 24 07:49 Input
-rwxr-xr-x   1 kkarri pulmseq  205 Aug 24 07:49 Input_Samples.txt
-rwxr-xr-x   1 kkarri pulmseq  36M Aug 24 07:49 RefSeq_GeneBody.bed
drwxr-xr-x   2 kkarri pulmseq 4.0K Aug 24 07:49 RefSeq_GeneBody_mm9-blacklist_Output
-rwxr-xr-x   1 kkarri pulmseq  71K Aug 24 07:49 mm9-blacklist.bed
-rwxr-xr-x   1 kkarri pulmseq 7.6K Aug 24 07:49 overlap.R
-rwxr-xr-x   1 kkarri pulmseq 6.0K Aug 24 07:49 overlap.sh
-rw-r--r--   1 kkarri pulmseq  80K Aug 24 07:49 peak1_dump.bed
==========================================================
Finished on : Thu Aug 24 07:49:51 EDT 2017
0 hours, 0 minutes and 29 seconds elapsed.
==========================================================