38 for (
int i = 0; i < size; i++) {
39 int current[6] = { 0, 0, 0, 0, 0, 0 };
42 switch (sequence[i]) {
139 int *cs = consensus[i].i;
140 double sum = (double) (cs[0] + cs[1] + cs[2] + cs[3] + cs[4] + cs[5]);
142 for (
int j = 0; j < 6; j++) {
143 int currentj = current[j];
145 if (cs[j] > currentj) {
146 cr.
conformity += (double) (cs[j] - currentj) / sum;
161 for (
int i = 0; i < size; i++) {
162 int *cs = consensus[i].i;
163 switch (sequence[i]) {
254 fprintf(stderr,
"Illegal character '%c'", sequence[i]);
void SQ_add_sequence(const char *sequence) OVERRIDE
consensus_result SQ_calc_consensus(const char *sequence) const OVERRIDE
void SQ_add_sequence(const char *sequence) OVERRIDE
consensus_result SQ_calc_consensus(const char *sequence) const OVERRIDE