7 return c == '-' || c == '.';
\r
10 const char *WordToStrAmino(unsigned Word, unsigned WordLength)
\r
12 static char Str[32];
\r
13 for (unsigned i = 0; i < WordLength; ++i)
\r
15 unsigned Letter = Word%20;
\r
16 Str[WordLength-i-1] = g_LetterToCharAmino[Letter];
\r
19 Str[WordLength] = 0;
\r
23 const char *WordToStrAmino2(unsigned Word, unsigned WordLength, char *Str)
\r
25 for (unsigned i = 0; i < WordLength; ++i)
\r
27 unsigned Letter = Word%20;
\r
28 Str[WordLength-i-1] = g_LetterToCharAmino[Letter];
\r
31 Str[WordLength] = 0;
\r
35 const char *WordToStrNucleo(unsigned Word, unsigned WordLength)
\r
37 static char Str[32];
\r
38 for (unsigned i = 0; i < WordLength; ++i)
\r
40 unsigned Letter = Word%4;
\r
41 Str[WordLength-i-1] = g_LetterToCharNucleo[Letter];
\r
44 Str[WordLength] = 0;
\r
48 const char *WordToStr(unsigned Word, unsigned WordLength, bool Nucleo)
\r
50 return (Nucleo ? WordToStrNucleo : WordToStrAmino)(Word, WordLength);
\r
53 byte *RevCompAlloc(const byte *Seq, unsigned L)
\r
55 byte *RCSeq = MYALLOC(byte, L, Alpha);
\r
57 for (unsigned i = 0; i < L; ++i)
\r
58 RCSeq[L-i-1] = g_CharToCompChar[Seq[i]];
\r
63 void RevCompInPlace(byte *Seq, unsigned L)
\r
65 unsigned L1 = L - 1;
\r
67 for (unsigned i = 0; i < L2; ++i)
\r
69 unsigned j = L1 - i;
\r
70 unsigned ci = Seq[i];
\r
71 unsigned cj = Seq[j];
\r
73 unsigned ri = g_CharToCompChar[ci];
\r
74 unsigned rj = g_CharToCompChar[cj];
\r
81 Seq[L2] = g_CharToCompChar[Seq[L2]];
\r
84 void RevComp(const byte *Seq, unsigned L, byte *RCSeq)
\r
86 for (unsigned i = 0; i < L; ++i)
\r
87 RCSeq[L-i-1] = g_CharToCompChar[Seq[i]];
\r
90 unsigned char GetAminoCharFrom3NucChars(unsigned char c1, unsigned char c2,
\r
93 unsigned Letter1 = g_CharToLetterNucleo[c1];
\r
94 unsigned Letter2 = g_CharToLetterNucleo[c2];
\r
95 unsigned Letter3 = g_CharToLetterNucleo[c3];
\r
96 unsigned Word = Letter1*(4*4) + Letter2*4 + Letter3;
\r
98 unsigned Letter = g_CodonWordToAminoLetter[Word];
\r
99 return g_LetterToCharAmino[Letter];
\r