# HG changeset patch
# User davidvanzessen
# Date 1426599865 14400
# Node ID 3f4b4ef46c7fc62691d481542b0d398f555f9a12
# Parent  cb7c65e3e43fb154ff66c8d426e4eda9fae4d594
Uploaded

diff -r cb7c65e3e43f -r 3f4b4ef46c7f gene_identification.py
--- a/gene_identification.py	Wed Mar 11 08:56:37 2015 -0400
+++ b/gene_identification.py	Tue Mar 17 09:44:25 2015 -0400
@@ -200,7 +200,7 @@
 				elif cg3hits > cg1hits and cg3hits > cg2hits and cg3hits > cg4hits: #cg3 gene
 					o.write(line.rstrip() + "\tcg3\t" + str(int(cg3hits / varsInCG * 100)) + "\t" + str(int(cghits / possiblecg * 100)) + "\t" + start_location[ID + "_cg"] + "\n")
 				else: #cg4 gene
-					o.write(line.rstrip() + "\tcg3\t" + str(int(cg4hits / varsInCG * 100)) + "\t" + str(int(cghits / possiblecg * 100)) + "\t" + start_location[ID + "_cg"] + "\n")
+					o.write(line.rstrip() + "\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(line.rstrip() + "\tcm\t0\t" + str(int(cmhits / possiblecm * 100)) + "\t" + start_location[ID + "_cg"] + "\n")