-----------------------
Start of variable list:
-----------------------
Sample_ID:
G168_G169_M5
Dataset_DIR:
/restricted/projectnb/waxmanlab/kkarri/G168_G169/G168/TCPOBOP_frStrand_Env/TCPOBOP
Sample_Labels_DIR:
/restricted/projectnb/waxmanlab/kkarri/G168_G169/G168/TCPOBOP_frStrand_Env/TCPOBOP/Scripts/00_Setup_Pipeline
STRAND_RULE:
1+-,1-+,2++,2--
BU_User:
kkarri
VM_DIR_UCSC:
/net/waxman-server/mnt/data/waxmanlabvm_home/kkarri/G168_TCPO_Env
SCRIPT_DIR:
/restricted/projectnb/waxmanlab/kkarri/G168_G169/G168/TCPOBOP_frStrand_Env/TCPOBOP/Scripts/11_UCSC_BigWig
-----------------------
End of variable list
-----------------------
==========================================================
Starting on : Wed Mar 13 10:52:59 EDT 2019
Running on node : scc-kb4
Current directory : /restricted/projectnb/waxmanlab/kkarri/G168_G169/G168/TCPOBOP_frStrand_Env/TCPOBOP/Scripts/11_UCSC_BigWig
Current job ID : 3566202
Current job name : Step_11_G168_G169_M5
Task index number : undefined
Parameter for multiple cores : 1
==========================================================

Change dir to scratch directory


Print scratch directory location:

/scratch/3566202.1.linga

Loading required modules...


INPUT_BAM:
G168_G169_M5_primary_unique.bam


List files in scratch directory:

total 1.9G
drwx------   4 kkarri waxmanlab 4.0K Mar 13 10:53 .
drwxrwxrwt. 31 root   root      164K Mar 13 10:52 ..
-rw-r--r--   1 kkarri waxmanlab 1.9G Mar 13 10:53 G168_G169_M5_primary_unique.bam
drwxr-xr-x   2 kkarri waxmanlab 4.0K Mar 13 10:53 UCSC_BigWig
drwxr-xr-x   2 kkarri waxmanlab 4.0K Mar 13 10:53 UCSC_Tools
-rwxr-xr-x   1 kkarri waxmanlab  584 Mar 13 10:53 mm9.chrom.sizes

Starting to run my commands


Starting samtools sort and index


Finished samtools sort and index


Starting bam2wig.py command

Printing command:
bam2wig.py --input-file=G168_G169_M5_sorted.bam --chromSize=mm9.chrom.sizes --out-prefix=G168_G169_M5 --wigsum=1000000000 --strand=1+-,1-+,2++,2--

Skip multi-hits:False

Finished bam2wig.py command


Sample from wiggle files:

==> G168_G169_M5.Forward.wig <==
variableStep chrom=chrY
90049	0.17
90050	0.17
90051	0.17
90052	0.17
90053	0.17
90054	0.17
90055	0.17
90056	0.17
90057	0.17

==> G168_G169_M5.Reverse.wig <==
variableStep chrom=chrY
434726	-0.34
434727	-0.34
434728	-0.34
434729	-0.34
434730	-0.34
434731	-0.34
434732	-0.34
434733	-0.34
434734	-0.34

#--------------------------------
G168_G169_M5.Forward.wig
#--------------------------------
Line count where 2nd column does not equal "-":
79676635
Line count where 2nd column is only a negative sign:
0
Copy G168_G169_M5.Forward.wig to UCSC_Tools folder
Copy mm9.chrom.sizes to UCSC_Tools folder
Start wigToBigWig
End wigToBigWig

#--------------------------------
G168_G169_M5.Reverse.wig
#--------------------------------
Line count where 2nd column does not equal "-":
77655142
Line count where 2nd column is only a negative sign:
0
Copy G168_G169_M5.Reverse.wig to UCSC_Tools folder
Copy mm9.chrom.sizes to UCSC_Tools folder
Start wigToBigWig
End wigToBigWig

List of files in UCSC_Tools folder:

-rw-r--r-- 1 kkarri waxmanlab 1141511305 Mar 13 12:42 ./UCSC_Tools/G168_G169_M5.Forward.wig
-rw-r--r-- 1 kkarri waxmanlab 1190718275 Mar 13 12:46 ./UCSC_Tools/G168_G169_M5.Reverse.wig
-rwxr-xr-x 1 kkarri waxmanlab        584 Mar 13 12:46 ./UCSC_Tools/mm9.chrom.sizes
-rwx------ 1 kkarri waxmanlab    1012225 Mar 13 10:53 ./UCSC_Tools/wigToBigWig

Copy output to storage dir


List files in scratch

total 5.9G
drwx------   4 kkarri waxmanlab 4.0K Mar 13 11:13 .
drwxrwxrwt. 42 root   root      164K Mar 13 12:45 ..
-rw-r--r--   1 kkarri waxmanlab 1.1G Mar 13 12:40 G168_G169_M5.Forward.wig
-rw-r--r--   1 kkarri waxmanlab 1.2G Mar 13 12:40 G168_G169_M5.Reverse.wig
-rw-r--r--   1 kkarri waxmanlab 1.9G Mar 13 10:53 G168_G169_M5_primary_unique.bam
-rw-r--r--   1 kkarri waxmanlab 1.9G Mar 13 11:07 G168_G169_M5_sorted.bam
-rw-r--r--   1 kkarri waxmanlab 3.2M Mar 13 11:08 G168_G169_M5_sorted.bam.bai
drwxr-xr-x   2 kkarri waxmanlab 4.0K Mar 13 12:47 UCSC_BigWig
drwxr-xr-x   2 kkarri waxmanlab 4.0K Mar 13 12:46 UCSC_Tools
-rwxr-xr-x   1 kkarri waxmanlab  584 Mar 13 10:53 mm9.chrom.sizes
==========================================================
Finished on : Wed Mar 13 12:48:38 EDT 2019
1 hours, 55 minutes and 39 seconds elapsed.
==========================================================