diff options
| author | Calvin Morrison <mutantturkey@gmail.com> | 2014-01-29 11:53:30 -0500 | 
|---|---|---|
| committer | Calvin Morrison <mutantturkey@gmail.com> | 2014-01-29 11:53:30 -0500 | 
| commit | 94d04a1e503121a98b403f882c18a4f0799267d7 (patch) | |
| tree | 0d2cf5586b31bddc9bca99b4b07ebb4b993f1130 /src/melting_range.c | |
| parent | 73531da5cdf33f9bde7d4db0e4ce96f1e41f581b (diff) | |
add filtering based on consecutive mer lengths
Diffstat (limited to 'src/melting_range.c')
| -rw-r--r-- | src/melting_range.c | 58 | 
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); -} | 
