363d60920206cd70ab8ec35a67a3e25c03eb8a46
braney
  Thu Jun 16 17:36:21 2022 -0700
add 30way sequence logo on hg38

diff --git src/hg/makeDb/doc/hg38/multiz30way.txt src/hg/makeDb/doc/hg38/multiz30way.txt
index e8658d6..a04cb3e 100644
--- src/hg/makeDb/doc/hg38/multiz30way.txt
+++ src/hg/makeDb/doc/hg38/multiz30way.txt
@@ -2394,15 +2394,50 @@
 set ylabel " Relative Frequency"
 set y2label " Cumulative Relative Frequency (CRF)"
 set y2range [0:1]
 set y2tics
 set xrange [-5:1.5]
 set yrange [0:0.04]
 
 plot "histogram.data" using 2:5 title " RelFreq" with impulses, \
         "histogram.data" using 2:7 axes x1y2 title " CRF" with lines
 ' | gnuplot > histo.png
 
     # verify it looks sane
     display histo.png &
 
 #############################################################################
+# sequence logo
+    mkdir /hive/data/genomes/hg38/bed/multiz30way/logo
+    cd /hive/data/genomes/hg38/bed/multiz30way/logo
+    mkdir wigs bw 
+
+    ls ../maf/* | awk '{printf "./makeWig.sh %s ../consPhyloP/all/phyloP30way.bw ../../../chrom.sizes\n", $1}' > jobs
+
+    cat << '_EOF_' > makeWig.sh
+maf=$1
+scale=$2
+chromSizes=$3
+outRoot=`basename $1 .maf`
+mafCounts -scale=$scale $maf wigs/$outRoot
+for i in A C G T
+do
+    wigToBigWig -keepAllChromosomes -fixedSummaries wigs/$outRoot.$i.wig $chromSizes bw/$outRoot.$i.bw
+done
+'_EOF_'
+
+    chmod +x makeWig.sh
+
+    ssh ku "cd /hive/data/genomes/hg38/bed/multiz30way/logo; para make jobs"
+
+    for i in A C G T
+    do
+        bigWigCat multiz30Logo.$i.bw bw/*.$i.bw &
+    done
+
+    wait
+
+    for i in A C G T
+    do
+        ln -s `pwd`/multiz30Logo.$i.bw /gbdb/hg38/multiz30way/
+    done
+