50c0049649c7268bd12ab0689cda7f8261b46a63 chmalee Mon May 6 16:11:42 2019 -0700 DGV Gold track mostly ready for QA, refs #23371 diff --git src/hg/utils/automation/translateDgvGold.py src/hg/utils/automation/translateDgvGold.py index 9e72afc..c412ff3 100755 --- src/hg/utils/automation/translateDgvGold.py +++ src/hg/utils/automation/translateDgvGold.py @@ -41,31 +41,31 @@ extraString = "\t".join(info) if doHeader: global header extraHeader.remove("Name") extraHeader.remove("outer_start") extraHeader.remove("inner_start") extraHeader.remove("inner_end") extraHeader.remove("outer_end") header += "\t" + "\t".join(extraHeader) return name, subType, outerStart, innerStart, innerEnd, outerEnd, extraString ### main ### -if len(sys.argv) != 1: +if len(sys.argv) != 2: print("Usage: translateDgvGold.py input.gff | sort -k1,1 -k2,2n | uniq > out.bed12\n"\ "\nWrites a bed12 to stdout.\n") sys.exit(1) doHeader = True for line in open(sys.argv[1]): gffFields = line.strip().split("\t") extraInfo = gffFields[8].split(";") chrom = gffFields[0] name, subType, outerStart, innerStart, innerEnd, chromEnd, extraString = parseExtraInfo(extraInfo, doHeader) if subType not in colors: sys.exit("CNV subtype: %s not recognized.\n" % subType) chromStart = str(int(outerStart) - 1)