3c1d42b5846d46612b495f1aeece9a20857480ec max Fri May 25 16:22:04 2018 -0700 crispr track for ce11, refs #21517 diff --git src/utils/doLocusName src/utils/doLocusName index 077e516..c5db8aa 100755 --- src/utils/doLocusName +++ src/utils/doLocusName @@ -35,30 +35,32 @@ print ("Running: "+ cmd) ret = os.system(cmd) if ret!=0 and mustRun: print "Could not run command %s" % cmd sys.exit(1) return ret def doLocusName(db, geneTableName, options): " create and fill the locusName table for db " if geneTableName=="ensGene": query = "SELECT * from ensemblToGeneName" elif geneTableName=="xenoRefGene": query = "SELECT name,name2 from xenoRefGene" elif geneTableName=="knownGene": query = "SELECT kgID, geneSymbol from kgXref where geneSymbol<>'' and geneSymbol not like '% %'" + elif geneTableName.startswith("ws"): + query = "SELECT name,name from %s" % geneTableName else: assert(False) # no support for this gene table yet tempFh = tempfile.NamedTemporaryFile() #tempFh = open("temp.tmp", "w") for row in runSql(db, query): tempFh.write("\t".join(row)+"\n") tempFh.flush() tempFname = tempFh.name chromSizesFname = "/cluster/data/%s/chrom.sizes" % db outTmp = "locusName.tmp.bed" cmdParts = [] # get genePred