# HG changeset patch
# User davidvanzessen
# Date 1480346842 18000
# Node ID 0bea8c187a90be2cc5eb472d137a13e46526b304
# Parent  c281993e515b1b932bf1bbe348b17441310bab5a
Uploaded

diff -r c281993e515b -r 0bea8c187a90 aa_histogram.r
--- a/aa_histogram.r	Mon Nov 28 10:24:25 2016 -0500
+++ b/aa_histogram.r	Mon Nov 28 10:27:22 2016 -0500
@@ -50,9 +50,6 @@
 
         print("---------------- write/print ----------------")
 
-        png(filename=paste(outdir, "/aa_histogram_", gene, ".png", sep=""), width=1280, height=720)
-        print(m)
-        dev.off()
 
         dat.sums = data.frame(index=1:length(mutations.at.position), mutations.at.position=mutations.at.position, aa.at.position=aa.at.position)
 
@@ -60,4 +57,8 @@
         write.table(mutations.by.id.gene, paste(outdir, "/aa_histogram_count_", gene, ".txt", sep=""), sep="\t",quote=F,row.names=F,col.names=T)
         write.table(absent.aa.by.id.gene, paste(outdir, "/aa_histogram_absent_", gene, ".txt", sep=""), sep="\t",quote=F,row.names=F,col.names=T)
         write.table(dat_dt, paste(outdir, "/aa_histogram_", gene, ".txt", sep=""), sep="\t",quote=F,row.names=F,col.names=T)
+        
+        png(filename=paste(outdir, "/aa_histogram_", gene, ".png", sep=""), width=1280, height=720)
+        print(m)
+        dev.off()
 }