bdc06272b74428b4b206745d07b0ac4c527444f8 max Mon May 18 07:01:01 2020 -0700 fixing the Orf1ab problem, refs #25114 as per email from Jason diff --git src/utils/genbankToBigGenePred src/utils/genbankToBigGenePred index fb0af70..100fee7 100755 --- src/utils/genbankToBigGenePred +++ src/utils/genbankToBigGenePred @@ -40,30 +40,34 @@ global doneProducts strand = "+" if ft.strand==-1: strand = "-" elif ft.strand==0: strand = "." start = ft.location.start end = ft.location.end size = ft.location.end-ft.location.start geneList = ft.qualifiers.get("gene") geneId = "" if geneList and len(geneList)==1: geneId = geneList[0] + # requested by Jason Fernandes, email Apr 29, 2020, 9:34 PM + if chromName=="NC_045512v2" and start==265 and end==13483: + geneId = "Orf1a" + prodList = ft.qualifiers.get("product") prod = "" if prodList and len(prodList)==1: prod = prodList[0] if nameField=="product" and prod in doneProducts: return doneProducts.add(prod) # set the BED name field if nameField=="gene": name = geneId elif nameField=="product": name = prod