--- a/sequence_overview.r	Mon Nov 07 03:04:07 2016 -0500
+++ b/sequence_overview.r	Mon Nov 07 06:52:53 2016 -0500
@@ -279,7 +279,9 @@
 NTresult = data.frame(nt=c("A", "C", "T", "G"))
 
 for(clazz in gene.classes){
+	print(paste("class:", clazz))
 	NToverview.sub = NToverview[grepl(paste("^", clazz, sep=""), NToverview$best_match),]
+	print(paste("nrow:", nrow(NToverview.sub)))
 	new.col.x = c(sum(NToverview.sub$A), sum(NToverview.sub$C), sum(NToverview.sub$T), sum(NToverview.sub$G))
 	new.col.y = sum(new.col.x)
 	new.col.z = round(new.col.x / new.col.y * 100, 2)