# HG changeset patch # User davidvanzessen # Date 1429191379 14400 # Node ID ac9a4307861a53b688fc2149fe06384aa2146ae8 # Parent 2a7343e4be5a91376e607d9c9862461aebaa2fb7 Uploaded diff -r 2a7343e4be5a -r ac9a4307861a mutation_analysis.py --- a/mutation_analysis.py Thu Apr 09 08:58:59 2015 -0400 +++ b/mutation_analysis.py Thu Apr 16 09:36:19 2015 -0400 @@ -168,4 +168,4 @@ o.write("ID\tRGYWC\tWRCY\tWA\tTW\n") first = False continue - o.write(ID + "\t" + str(int(RGYWCount[ID])) + "\t" + str(int(WRCYCount[ID])) + "\t" + str(int(WACount[ID])) + "\t" + str(int(TWCount[ID])) + "\n") + 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") diff -r 2a7343e4be5a -r ac9a4307861a mutation_analysis.r --- a/mutation_analysis.r Thu Apr 09 08:58:59 2015 -0400 +++ b/mutation_analysis.r Thu Apr 16 09:36:19 2015 -0400 @@ -284,7 +284,7 @@ result = data.frame(matrx) -row.names(result) = c("Number of Mutations (%)", "Transition (%)", "Transversions (%)", "Transitions at G C (%)", "Targeting of C.G (%)", "FR S/R (ratio)", "CDR S/R (ratio)") +row.names(result) = c("Number of Mutations (%)", "Transition (%)", "Transversions (%)", "Transitions at G C (%)", "Targeting of C.G (%)", "FR R/S (ratio)", "CDR R/S (ratio)") write.table(x=result, file="mutations.txt", sep=",",quote=F,row.names=T,col.names=F)