Mercurial > repos > davidvanzessen > shm_csr

comparison shm_csr.py @ 16:4e596473c25c draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author davidvanzessen
date Wed, 23 Nov 2016 09:05:40 -0500
parents ad9be244b104
children b95fa7e426c3
comparison
equal deleted inserted replaced
15:61d0a6318711 16:4e596473c25c
282 # for testing 282 # for testing
283 seq_motif_file = outfile[:outfile.rindex("/")] + "/motif_per_seq.txt" 283 seq_motif_file = outfile[:outfile.rindex("/")] + "/motif_per_seq.txt"
284 with open(seq_motif_file, 'w') as o: 284 with open(seq_motif_file, 'w') as o:
285 o.write("ID\tRGYWC\tWRCY\tWA\tTW\n") 285 o.write("ID\tRGYWC\tWRCY\tWA\tTW\n")
286 for ID in IDlist: 286 for ID in IDlist:
287 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") 287 #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")
288 o.write(ID + "\t" + str(RGYWCount[ID]) + "\t" + str(WRCYCount[ID]) + "\t" + str(WACount[ID]) + "\t" + str(TWCount[ID]) + "\n")