f5a1b930828e2737ee09d3dc71b3a0c9c2068452 max Wed Mar 27 08:21:03 2024 -0700 otto job for decipher failed, otto email diff --git src/hg/utils/otto/decipher/processDecipher.py src/hg/utils/otto/decipher/processDecipher.py index 9bf61ca..7479478 100755 --- src/hg/utils/otto/decipher/processDecipher.py +++ src/hg/utils/otto/decipher/processDecipher.py @@ -70,31 +70,31 @@ if len(sys.argv) != 6: print ("Usage: makeUcscBed.py ") sys.exit(0) cnvOutput = open(sys.argv[3], "w") snvOutput = open(sys.argv[4], "w") snvRawOutput = open(sys.argv[5], "w") geneFields = namedtuple("geneFields", ["chrom", "chromStart", "chromEnd", "name"]) geneList = [] with open(sys.argv[2]) as geneFile: for line in geneFile: fields = line.strip('\n').split('\t') geneList.append(geneFields(chrom=fields[0],chromStart=int(fields[1]),chromEnd=int(fields[2]),name=fields[3])) - with open(sys.argv[1], "r") as inputFile: + with open(sys.argv[1], "r", encoding="utf8") as inputFile: header = inputFile.readline() for line in inputFile: fields = line.strip('\n').split('\t') parsedJson = json.loads(fields[13]) varType = parsedJson["variant_type"]; if (varType == "SNV") : # Make the basic table file result = "\t".join(fields[:4]) print (result, file=snvOutput) # Make the "raw" data file with more fields # format: id, chrom, start, end, ref, alt, transcript, gene, genotype, inheritance, pathgenicity, # contribution, phenotypes fields[0] = fields[0].replace("chr", "") # strip off initial "chr" for these names