a823a26a0ed014c44482f50eb1c4072b3860f9de max Mon Oct 31 14:27:59 2016 -0700 small fix to doLocusName to handle genes with with spaces in names diff --git src/utils/doLocusName src/utils/doLocusName index f2c0fa3..8de63da 100755 --- src/utils/doLocusName +++ src/utils/doLocusName @@ -63,28 +63,28 @@ # only known Gene has no bin field if geneTableName!="knownGene": cmdParts.append("cut -f2- |") # break genes into exons cmdParts.append("genePredToBed stdin stdout | bedToExons stdin stdout | ") # replace transcript with symbol, skip dupes cmdParts.append("tabRepl %(tempFname)s 3 /dev/stdin | sort -u | ") # merge adjacent exons cmdParts.append("bedSort stdin stdout | bedOverlapMerge /dev/stdin /dev/stdout |") # annotate regions cmdParts.append("bedBetween stdin /dev/stdout -a -s %(chromSizesFname)s -m 100 |") # and sort cmdParts.append("bedSort stdin %(outTmp)s") runCmd(("".join(cmdParts)) % locals()) - cmd = "hgLoadBed %s locusName %s" % (db, outTmp) + cmd = "hgLoadBed -tab %s locusName %s" % (db, outTmp) runCmd(cmd) os.remove(outTmp) # ----------- main -------------- def main(): args, options = parseArgs() db, geneTable = args doLocusName(db, geneTable) main()