aboutsummaryrefslogtreecommitdiff
path: root/src/melting_range.c
diff options
context:
space:
mode:
authorCalvin Morrison <mutantturkey@gmail.com>2014-01-29 11:53:30 -0500
committerCalvin Morrison <mutantturkey@gmail.com>2014-01-29 11:53:30 -0500
commit94d04a1e503121a98b403f882c18a4f0799267d7 (patch)
tree0d2cf5586b31bddc9bca99b4b07ebb4b993f1130 /src/melting_range.c
parent73531da5cdf33f9bde7d4db0e4ce96f1e41f581b (diff)
add filtering based on consecutive mer lengths
Diffstat (limited to 'src/melting_range.c')
-rw-r--r--src/melting_range.c58
1 files changed, 0 insertions, 58 deletions
diff --git a/src/melting_range.c b/src/melting_range.c
deleted file mode 100644
index 2c89195..0000000
--- a/src/melting_range.c
+++ /dev/null
@@ -1,58 +0,0 @@
-#include <stdlib.h>
-#include <stdio.h>
-#include <string.h>
-#include <unistd.h>
-
-float melting_temperature(char *mer) {
-
- float a = 0;
- float c = 0;
- float g = 0;
- float t = 0;
- int i = 0;
-
- for(i = 0; i < strlen(mer); i++) {
- switch(mer[i]) {
- case 'A':
- a++;
- break;
- case 'C':
- c++;
- break;
- case 'G':
- g++;
- break;
- case 'T':
- t++;
- break;
- default:
- break;
- }
- }
-
- if(strlen(mer) < 13)
- return ((a+t) * 2) + ((c+g) * 4);
- else
- return 64.9 + 41.0*(g+c-16.4)/(a+t+g+c);
-}
-
-int main(int argc, char **argv){
-
- if(argc < 3) {
- printf("please supply the min and max as stdargs");
- exit(EXIT_FAILURE);
- }
- float min = atof(argv[1]);
- float max = atof(argv[2]);
-
- char mer[24] = { 0 };
- int count = 0;
-
- while(fscanf(stdin, "%s\t%d\n", &mer, &count) == 2) {
- float temp = melting_temperature(mer);
- if( (temp > min) && (temp < max) )
- printf("%s\t%d\n", mer, count);
- }
-
- exit(EXIT_SUCCESS);
-}