From 5d4fabb8f2a31ad87fccfb0b8ee426b464867f54 Mon Sep 17 00:00:00 2001
From: Calvin Morrison <mutantturkey@gmail.com>
Date: Tue, 23 Jul 2013 16:37:19 -0400
Subject: added a nucleobase counter

---
 count_nucleobases.c | 41 +++++++++++++++++++++++++++++++++++++++++
 1 file changed, 41 insertions(+)
 create mode 100755 count_nucleobases.c

(limited to 'count_nucleobases.c')

diff --git a/count_nucleobases.c b/count_nucleobases.c
new file mode 100755
index 0000000..6edf5e8
--- /dev/null
+++ b/count_nucleobases.c
@@ -0,0 +1,41 @@
+#include <stdio.h>
+#include <stdlib.h>
+#include <unistd.h>
+#include <stdint.h>
+main() {
+
+long unsigned long a = 0;
+long unsigned long c = 0;
+long unsigned long g = 0;
+long unsigned long t = 0;
+
+FILE *fh = fopen("dna.txt", "r" );
+  if(fh == NULL) {
+    fprintf(stderr, "could not open dna.txt\n");
+    exit(EXIT_FAILURE);
+  }
+
+  unsigned int ch = 0;
+  while ((ch = getc(fh)) != EOF) {
+    switch(ch) {
+      case 'A':
+      case 'a':
+        a++;
+        break;
+      case 'C':
+      case 'c':
+        c++;
+        break;
+      case 'G':
+      case 'g':
+        g++;
+        break;
+      case 'T':
+      case 't':
+        t++;
+        break;
+    }
+  }
+
+  printf("A:%lu\nC:%lu\nG:%lu\nT:%lu\n", a, c, g, t);
+}
-- 
cgit v1.2.3