a64bc5a6cfe9756104a86f633b0f0a01772c4b83 lrnassar Wed Oct 6 16:55:59 2021 -0700 Updating a few things on otto. For the most part these are changes to the wrapper emails, which had not been updated in git, and the big changes are were done by engineers directly to the hive copy without updating the tree version. The big updates for clinvar were Max, and decipher was Braney I believe. To clarify, these now updated versions are the ones that have been running daily/weekly, so all of it has been in action already. diff --git src/hg/utils/otto/omim/omimGene2ToBigBed.py src/hg/utils/otto/omim/omimGene2ToBigBed.py index 852a6dd..c68ae4b 100644 --- src/hg/utils/otto/omim/omimGene2ToBigBed.py +++ src/hg/utils/otto/omim/omimGene2ToBigBed.py @@ -64,31 +64,32 @@ for line in open(inFname): row = line.rstrip("\n").split("\t") # 1 chr12 # 2 4477392 # 3 4488878 # 4 605380 # 5 0 # 6 . # 7 4477392 # 8 4488878 # 9 color #10 FGF23, ADHR, HPDR2, PHPTC, HFTC2 #11 #12 Tumoral calcinosis, hyperphosphatemic, familial, 2|3|$Hypophosphatemic rickets, autosomal dominant|3|Autosomal dominant - chrom, start, end, mimId, score, strand, thickStart, thickEnd, dummyColor, syms, oldDisorderStr, phenoStr = row + #13 PhenoMapKey + chrom, start, end, mimId, score, strand, thickStart, thickEnd, dummyColor, syms, oldDisorderStr, phenoStr,phenoMapKey = row mapKeys = [] if phenoStr=="": newPhenoStr = "" newPhenoPlEnding = "" else: phenoList = phenoStr.split("$") newPhenos = [] for phenoPart in phenoList: name, mapKey, inhMode = phenoPart.split("|") inhCode = inhModeToCode(inhMode) phenoLabels = [] phenoLabels.append(name) if inhMode!="": @@ -123,29 +124,29 @@ color = "190,190,190" else: mapKeys.sort() maxKey = mapKeys[-1] if maxKey=="4": color = "105,50,155" elif maxKey=="3": color = "0,85,0" elif maxKey=="2": color = "102,150,102" elif maxKey=="1": color = "170,196,170" else: assert(False) - newRow = (chrom, start, end, mimId, score, strand, thickStart, thickEnd, color, mainSym, mouseOver) + newRow = (chrom, start, end, mimId, score, strand, thickStart, thickEnd, color, mainSym, phenoMapKey, mouseOver) ofh.write("\t".join(newRow)) ofh.write("\n") #ofh.flush() ofh.close() cmd = "bedSort %s %s" % (tmpFname, tmpFname) os.system(cmd) cmd = "bedToBigBed %s %s %s -tab -extraIndex=name -as=omimGene2.as -type=bed9+" % (tmpFname, chromSizesFname, outFname) os.system(cmd) main()