blob: b6f66b583b7b64a87354f6febcf5299471e1f922 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
|
#include <stdint.h>
#include <stdio.h>
#include <string.h>
#include <unistd.h>
// This function takes a char array containing sequences and converts it
// into a kmer index (essentially a base 4 radix converstion to base 10, with
// some mumbo-jumbo of taking each of the characters we want (ACGT) and turning
// them into a radix-4 number we can use.
//
// kind of convoluted but it works.
//
// Arguemnts:
// char *str - a NULL terminated character array
// long kmer - how long of a index value you want to return
// long error_pos - what index to return for a non ACGT character
//
long convert_kmer_to_index(char *str, long kmer, long error_pos) {
char **ptr = NULL;
char vals[kmer];
long i = 0;
long ret = 0;
for(i = 0; i < kmer; i++) {
int val = 0;
switch(str[i]) {
case 'a':
case 'A':
val = 48;
break;
case 'c':
case 'C':
val = 49;
break;
case 'g':
case 'G':
val = 50;
break;
case 't':
case 'T':
val = 51;
break;
default:
return error_pos;
}
vals[i] = val;
}
ret = strtol(vals, ptr, 4);
return ret;
}
|