94d6890bdac9d6c53a6070eef9bcb82611c7d25f
chmalee
  Fri Feb 14 11:32:08 2020 -0800
Adding mouseOver to missense constraint track and finishing up track description, refs #20394

diff --git src/hg/makeDb/gnomad/aaToGenomic.py src/hg/makeDb/gnomad/aaToGenomic.py
index 734efc4..17bb40e 100755
--- src/hg/makeDb/gnomad/aaToGenomic.py
+++ src/hg/makeDb/gnomad/aaToGenomic.py
@@ -259,33 +259,33 @@
     ret = bedInfo(transcript.chrom, chromStart, chromEnd, transcript.name, "%d" % int(pscore),
         transcript.strand, chromStart, chromEnd, "255,0,0", cdnaEndExon-cdnaStartExon+1,
         outSizes, outStarts)
     checkAaLenToBlockLen(ret, transcript, cdnaEndPos, cdnaStartPos, aaStart, aaEnd)
     sanityCheckBed12(ret, transcript)
     return ret
 
 def aaToBedByStrand(aaStart, aaEnd, score, observed, expected, chi2, gene, transcript, lineNum):
     """Dispatch to write function depending on strand."""
     ret = []
     if transcript.strand == "+":
         ret = aaToBedPosStrand(aaStart, aaEnd, score, transcript, lineNum)
     else:
         ret = aaToBedNegStrand(aaStart, aaEnd, score, transcript, lineNum)
     # now add the extra scoring information
-    color = findColor(score, chi2)
+    color = "160,160,160" # findColor(score, chi2)
     ret = ret._replace(itemRgb=color)
-    print("%s\t%s\t%d\t%.3f\t%.3f\t%.3f" % (bedToStr(ret), gene, observed, expected, score, chi2))
+    print("%s\t%s\t%d\t%.3f\t%.3f\t%.3f\tO/E: %0.3f" % (bedToStr(ret), gene, observed, expected, score, chi2, score))
 
 def aaToBed(aaFh, verbose):
     """
     Using the transcript dict, get the coordinates for an amino acid range in bed12 format.
     Example input line:
     ENST00000337907.3   RERE    1   1-507   8716356 8424825 97  197.9807    0.489947    51.505535   RERE_1
     Example output line:
     chr1    8424824 8716356 ENST00000337907.3   0.49    -   8424824 8716356 224,165,8   13  74,163,81,99,100,125,49,105,97,106,126,71,325   0,1047,57962,101160,130298,132640,143861,176448,191709,192652,249795,259544,291207 RERE 97 197.98 0.48 51.50
     Example output line:
     """
     lineCount = 1
     for line in aaFh:
         if line.startswith('transcript'):
             continue
         splitLine = line.split('\t')