# HG changeset patch
# User davidvanzessen
# Date 1480001059 18000
# Node ID b95fa7e426c3a634e7f7e54f0ce7d502424bb185
# Parent  4e596473c25c16a27999f79090bbf862940be7ef
Uploaded

diff -r 4e596473c25c -r b95fa7e426c3 shm_csr.py
--- a/shm_csr.py	Wed Nov 23 09:05:40 2016 -0500
+++ b/shm_csr.py	Thu Nov 24 10:24:19 2016 -0500
@@ -247,12 +247,9 @@
 	
 
 def get_xyz(lst, gene, f, fname):
-	
-	x = int(round(f(lst)))
+	x = round(round(f(lst), 1))
 	y = valuedic[gene + "_" + fname]
 	z = str(round(x / float(y) * 100, 1)) if y != 0 else "0"
-	if gene == "unmatched":
-		print x, y, z
 	return (str(x), str(y), z)
 
 dic = {"RGYW": RGYWCount, "WRCY": WRCYCount, "WA": WACount, "TW": TWCount}
@@ -282,7 +279,7 @@
 # for testing
 seq_motif_file = outfile[:outfile.rindex("/")] + "/motif_per_seq.txt"
 with open(seq_motif_file, 'w') as o:
-	o.write("ID\tRGYWC\tWRCY\tWA\tTW\n")
+	o.write("ID\tRGYW\tWRCY\tWA\tTW\n")
 	for ID in IDlist:
 		#o.write(ID + "\t" + str(round(RGYWCount[ID], 2)) + "\t" + str(round(WRCYCount[ID], 2)) + "\t" + str(round(WACount[ID], 2)) + "\t" + str(round(TWCount[ID], 2)) + "\n")
 		o.write(ID + "\t" + str(RGYWCount[ID]) + "\t" + str(WRCYCount[ID]) + "\t" + str(WACount[ID]) + "\t" + str(TWCount[ID]) + "\n")