# HG changeset patch # User davidvanzessen # Date 1458303428 14400 # Node ID 64e6a7803e07ec0c515fb6213acd474d85d84db3 # Parent 1558fbbc19efbb04dff442ff0706052d0a2415a3 Uploaded diff -r 1558fbbc19ef -r 64e6a7803e07 gene_identification.py --- a/gene_identification.py Fri Mar 18 08:02:50 2016 -0400 +++ b/gene_identification.py Fri Mar 18 08:17:08 2016 -0400 @@ -199,7 +199,7 @@ else: #cg4 gene o.write(ID + "\tcg4\t" + str(int(cg4hits / varsInCG * 100)) + "\t" + str(int(cghits / possiblecg * 100)) + "\t" + start_location[ID + "_cg"] + "\n") else: #its a cm gene - o.write(ID + "\tcm\t0\t" + str(int(cmhits / possiblecm * 100)) + "\t" + start_location[ID + "_cg"] + "\n") + o.write(ID + "\tcm\t100\t" + str(int(cmhits / possiblecm * 100)) + "\t" + start_location[ID + "_cg"] + "\n") seq_write_count += 1 print "Time: %i" % (int(time.time() * 1000) - starttime)