> produceTSSmaps > 0. Running on Wed Jun 9 14:08:32 2021 > 1. Request Memory to Operating System Loading the R color schema > 2. Reading Chromosome Sizes Dictionary (hash): chr8 | 7 chr9 | 8 chr10 | 9 chr11 | 10 chr12 | 11 chrX | 20 chr13 | 12 chrY | 21 chr14 | 13 chr15 | 14 chr16 | 15 chr17 | 16 chr18 | 17 chr19 | 18 chr1 | 0 chr2 | 1 chr3 | 2 chr4 | 3 chrM | 19 chr5 | 4 chr6 | 5 chr7 | 6 Chromosome sizes: 0 | 197195432 1 | 181748087 2 | 159599783 3 | 155630120 4 | 152537259 5 | 149517037 6 | 152524553 7 | 131738871 8 | 124076172 9 | 129993255 10 | 121843856 11 | 121257530 12 | 120284312 13 | 125194864 14 | 103494974 15 | 98319150 16 | 95272651 17 | 90772031 18 | 61342430 19 | 16299 20 | 166650296 21 | 15902555 Size was successfully acquired for 22 chromosomes > 3. Prepare the Output Folder The prefix folder ./ exists The prefix folder ./catalina/ exists The prefix folder ./catalina/results/ exists Creating the output folder ./catalina/results/H3K27ac_mESC_TSSmap_5000 Filename for the HEAT matrix ./catalina/results/H3K27ac_mESC_TSSmap_5000/HEATmatrix_H3K27ac_mESC_5000.txt Filename for the processed HEAT matrix ./catalina/results/H3K27ac_mESC_TSSmap_5000/HEATmatrix_H3K27ac_mESC_5000_processed.txt Filename for the HEAT gene rank ./catalina/results/H3K27ac_mESC_TSSmap_5000/HEATmatrix_H3K27ac_mESC_5000_genes.txt PDF file for the HEAT plot file ./catalina/results/H3K27ac_mESC_TSSmap_5000/PlotHEATmap_H3K27ac_mESC_5000.pdf > 4. Reading Refgene file data/mm9/refGene.txt: 12.61 Mb 46881 transcripts of RefSeq were successfully acquired > 5. Processing the BAM/SAM file samples/mouse/H3K27ac_mESC.bam: 1988.53 Mb Reading:  [ 0 lines] [ 100000 lines] [ 200000 lines] [ 300000 lines] [ 400000 lines] [ 500000 lines] [ 600000 lines] [ 700000 lines] [ 800000 lines] [ 900000 lines] [ 1000000 lines] [ 1100000 lines] [ 1200000 lines] [ 1300000 lines] [ 1400000 lines] [ 1500000 lines] [ 1600000 lines] [ 1700000 lines] [ 1800000 lines] [ 1900000 lines] [ 2000000 lines] [ 2100000 lines] [ 2200000 lines] [ 2300000 lines] [ 2400000 lines] [ 2500000 lines] [ 2600000 lines] [ 2700000 lines] [ 2800000 lines] [ 2900000 lines] [ 3000000 lines] [ 3100000 lines] [ 3200000 lines] [ 3300000 lines] [ 3400000 lines] [ 3500000 lines] [ 3600000 lines] [ 3700000 lines] [ 3800000 lines] [ 3900000 lines] [ 4000000 lines] [ 4100000 lines] [ 4200000 lines] [ 4300000 lines] [ 4400000 lines] [ 4500000 lines] [ 4600000 lines] [ 4700000 lines] [ 4800000 lines] [ 4900000 lines] [ 5000000 lines] [ 5100000 lines] [ 5200000 lines] [ 5300000 lines] [ 5400000 lines] [ 5500000 lines] [ 5600000 lines] [ 5700000 lines] [ 5800000 lines] [ 5900000 lines] [ 6000000 lines] [ 6100000 lines] [ 6200000 lines] [ 6300000 lines] [ 6400000 lines] [ 6500000 lines] [ 6600000 lines] [ 6700000 lines] [ 6800000 lines] [ 6900000 lines] [ 7000000 lines] [ 7100000 lines] [ 7200000 lines] [ 7300000 lines] [ 7400000 lines] [ 7500000 lines] [ 7600000 lines] [ 7700000 lines] [ 7800000 lines] [ 7900000 lines] [ 8000000 lines] [ 8100000 lines] [ 8200000 lines] [ 8300000 lines] [ 8400000 lines] [ 8500000 lines] [ 8600000 lines] [ 8700000 lines] [ 8800000 lines] [ 8900000 lines] [ 9000000 lines] [ 9100000 lines] [ 9200000 lines] [ 9300000 lines] [ 9400000 lines] [ 9500000 lines] [ 9600000 lines] [ 9700000 lines] [ 9800000 lines] [ 9900000 lines] [ 10000000 lines] [ 10100000 lines] [ 10200000 lines] [ 10300000 lines] [ 10400000 lines] [ 10500000 lines] [ 10600000 lines] [ 10700000 lines] [ 10800000 lines] [ 10900000 lines] [ 11000000 lines] [ 11100000 lines] [ 11200000 lines] [ 11300000 lines] [ 11400000 lines] [ 11500000 lines] [ 11600000 lines] [ 11700000 lines] [ 11800000 lines] [ 11900000 lines] [ 12000000 lines] [ 12100000 lines] [ 12200000 lines] [ 12300000 lines] [ 12400000 lines] [ 12500000 lines] [ 12600000 lines] [ 12700000 lines] [ 12800000 lines] [ 12900000 lines] [ 13000000 lines] [ 13100000 lines] [ 13200000 lines] [ 13300000 lines] [ 13400000 lines] [ 13500000 lines] [ 13600000 lines] [ 13700000 lines] [ 13800000 lines] [ 13900000 lines] [ 14000000 lines] [ 14100000 lines] [ 14200000 lines] [ 14300000 lines] [ 14400000 lines] [ 14500000 lines] [ 14600000 lines] [ 14700000 lines] [ 14800000 lines] [ 14900000 lines] [ 15000000 lines] [ 15100000 lines] [ 15200000 lines] [ 15300000 lines] [ 15400000 lines] [ 15500000 lines] [ 15600000 lines] [ 15700000 lines] [ 15800000 lines] [ 15900000 lines] [ 16000000 lines] [ 16100000 lines] [ 16200000 lines] [ 16300000 lines] [ 16400000 lines] [ 16500000 lines] [ 16600000 lines] [ 16700000 lines] [ 16800000 lines] [ 16900000 lines] [ 17000000 lines] [ 17100000 lines] [ 17200000 lines] [ 17300000 lines] [ 17400000 lines] [ 17500000 lines] [ 17600000 lines] [ 17700000 lines] [ 17800000 lines] [ 17900000 lines] [ 18000000 lines] [ 18100000 lines] [ 18200000 lines] [ 18300000 lines] [ 18400000 lines] [ 18500000 lines] [ 18600000 lines] [ 18700000 lines] [ 18800000 lines] [ 18900000 lines] [ 19000000 lines] [ 19100000 lines] [ 19200000 lines] [ 19300000 lines] [ 19400000 lines] [ 19500000 lines] [ 19600000 lines] [ 19700000 lines] [ 19800000 lines] [ 19900000 lines] [ 20000000 lines] [ 20100000 lines] [ 20200000 lines] [ 20300000 lines] [ 20400000 lines] [ 20500000 lines] [ 20600000 lines] [ 20700000 lines] [ 20800000 lines] [ 20900000 lines] [ 21000000 lines] [ 21100000 lines] [ 21200000 lines] [ 21300000 lines] [ 21400000 lines] [ 21500000 lines] [ 21600000 lines] [ 21700000 lines] [ 21800000 lines] [ 21900000 lines] [ 22000000 lines] [ 22100000 lines] [ 22200000 lines] [ 22300000 lines] [ 22400000 lines] [ 22500000 lines] [ 22600000 lines] [ 22700000 lines] [ 22800000 lines] [ 22900000 lines] [ 23000000 lines] [ 23100000 lines] [ 23200000 lines] [ 23300000 lines] [ 23400000 lines] [ 23500000 lines] [ 23600000 lines] [ 23700000 lines] [ 23800000 lines] [ 23900000 lines] [ 24000000 lines] [ 24100000 lines] [ 24200000 lines] [ 24300000 lines] [ 24400000 lines] [ 24500000 lines] [ 24600000 lines] [ 24700000 lines] [ 24800000 lines] [ 24900000 lines] [ 25000000 lines] [ 25100000 lines] [ 25200000 lines] [ 25300000 lines] [ 25400000 lines] [ 25500000 lines] [ 25600000 lines] [ 25700000 lines] [ 25800000 lines] [ 25900000 lines] [ 26000000 lines] [ 26100000 lines] [ 26200000 lines] [ 26300000 lines] [ 26400000 lines] [ 26500000 lines] [ 26600000 lines] [ 26700000 lines] [ 26800000 lines] [ 26900000 lines] [ 27000000 lines] [ 27100000 lines] [ 27200000 lines] [ 27300000 lines] [ 27400000 lines] [ 27500000 lines] [ 27600000 lines] [ 27700000 lines] [ 27800000 lines] [ 27900000 lines] [ 28000000 lines] [ 28100000 lines] [ 28200000 lines] [ 28300000 lines] [ 28400000 lines] [ 28500000 lines] [ 28600000 lines] [ 28700000 lines] [ 28800000 lines] [ 28900000 lines] [ 29000000 lines] [ 29100000 lines] [ 29200000 lines] [ 29300000 lines] [ 29400000 lines] [ 29500000 lines] [ 29600000 lines] [ 29700000 lines] [ 29800000 lines] [ 29900000 lines] [ 30000000 lines] [ 30100000 lines] [ 30200000 lines] [ 30300000 lines] [ 30400000 lines] [ 30500000 lines] [ 30600000 lines] [ 30700000 lines] [ 30800000 lines] [ 30900000 lines] [ 31000000 lines] [ 31100000 lines] [ 31200000 lines] [ 31300000 lines] [ 31400000 lines] [ 31500000 lines] [ 31600000 lines] [ 31700000 lines] [ 31800000 lines] [ 31900000 lines] [ 32000000 lines] [ 32100000 lines] [ 32200000 lines] [ 32300000 lines] [ 32400000 lines] [ 32500000 lines] [ 32600000 lines] [ 32700000 lines] [ 32800000 lines] [ 32900000 lines] [ 33000000 lines] [ 33100000 lines] [ 33200000 lines] [ 33300000 lines] [ 33400000 lines] [ 33500000 lines] [ 33600000 lines] [ 33700000 lines] [ 33800000 lines] [ 33900000 lines] [ 34000000 lines] 34078546 reads of the SAM file processed 17039279 forward reads acquired 17039401 reverse reads acquired > 6. Processing the List of target genes List of non-available genes: 1700074H08Rik 2010315B03Rik 4930526I15Rik AF067063 Aida B020031M17Rik BC147527 Brox Ccdc82 Ccl28 Csprs Defb33 Dhrsx Disp1 Gm20736 Gm20738 Gm20747 Gm20806 Gm20809 Gm20815 Gm20816 Gm20819 Gm20822 Gm20823 Gm20824 Gm20826 Gm20831 Gm20852 Gm20854 Gm20857 Gm20858 Gm20865 Gm20867 Gm20868 Gm20871 Gm20877 Gm20879 Gm20901 Gm20909 Gm20917 Gm20924 Gm20931 Gm21064 Gm21242 Gm21258 Gm21292 Gm21302 Gm21308 Gm21316 Gm21330 Gm21379 Gm21394 Gm21440 Gm21454 Gm21469 Gm21679 Gm21724 Gm21745 Gm21751 Gm21753 Gm21763 Gm21783 Gm21805 Gm21810 Gm21818 Gm21828 Gm21838 Gm21843 Gm21908 Gm21918 Gm21943 Gm21944 Gm28093 Gm28371 Gm28459 Gm28840 Gm29317 Gm31186 Gm41077 Hhipl2 Hmgcs1 Jrkl LOC100041034 LOC100041877 LOC100861691 LOC102638047 LOC664787 Mia3 Mucl1 Mucl2 Nim1k Nlrp4g Orly Ppp2r3d Rbm31y Sly Spry3 Srsy Ssty1 Ssty2 Taf1a Tcstv1 Tcstv3 Tmlhe Ulbp1 Vamp7 Vmn2r123 Zf12 Zfp131 26074 target genes were successfully processed Extracting gene avg. values in ./catalina/results/H3K27ac_mESC_TSSmap_5000/HEATmatrix_H3K27ac_mESC_5000_genes.txt cut -f 1,102 ./catalina/results/H3K27ac_mESC_TSSmap_5000/HEATmatrix_H3K27ac_mESC_5000.txt | sort -nrk 2 > ./catalina/results/H3K27ac_mESC_TSSmap_5000/HEATmatrix_H3K27ac_mESC_5000_genes.txt Ranking the heat file: ./catalina/results/H3K27ac_mESC_TSSmap_5000/HEATmatrix_H3K27ac_mESC_5000_processed.txt sort -nk 102 ./catalina/results/H3K27ac_mESC_TSSmap_5000/HEATmatrix_H3K27ac_mESC_5000.txt | cut -f 2-101 > ./catalina/results/H3K27ac_mESC_TSSmap_5000/HEATmatrix_H3K27ac_mESC_5000_processed.txt > 7. Running R to produce the plots > Preparing the Rscript file ./catalina/results/H3K27ac_mESC_TSSmap_5000/RscriptTSSmap_H3K27ac_mESC_5000.txt R CMD BATCH ./catalina/results/H3K27ac_mESC_TSSmap_5000/RscriptTSSmap_H3K27ac_mESC_5000.txt Removing the out R file rm -f *.Rout __________________________________________________________________________ Memory usage: 111.98 Megabytes CPU time: 27.11 secs Total time: 52 secs (0.87 mins)