-
Notifications
You must be signed in to change notification settings - Fork 2
/
Copy pathbase.cpp
80 lines (74 loc) · 1.71 KB
/
base.cpp
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
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
#include "base.h"
#include "utils.h"
namespace BASE {
base char2base(char c)
{
switch (c) {
case 'A': return A;
case 'C': return C;
case 'G': return G;
case 'T': return T;
//case 'N': return N;
case 'N': return A;
default: panic("Invalid %s: %u\n", __func__, c);
}
}
char base2char(base b)
{
switch(b) {
case A: return 'A';
case C: return 'C';
case G: return 'G';
case T: return 'T';
//case N: return 'N';
default: panic("Invalid %s: %u\n", __func__, b);
}
}
base inv_base(base b)
{
switch (b) {
case A: return T;
case C: return G;
case G: return C;
case T: return A;
//case N: return N;
default: panic("Invalid %s: %u\n", __func__, b);
}
}
char inv_base(char c)
{
switch (c) {
case 'A': return 'T';
case 'C': return 'G';
case 'G': return 'C';
case 'T': return 'A';
//case 'N': return 'N';
default: panic("Invalid %s: %u\n", __func__, c);
}
}
bool valid_base(base b)
{
switch (b) {
case A:
case C:
case G:
case T:
//case N:
return true;
default:
return false;
}
}
bool valid_base(char c)
{
switch (c) {
case 'A':
case 'C':
case 'G':
case 'T':
return true;
default:
return false;
}
}
}