- if (ref4 < 4) {
- call->a[0] = ref4;
- for (i = 3, j = 1; i >= 0; --i) {
- if ((sum[i]&3) != ref4) {
- if (sum[i]>>2 != 0) call->a[j++] = sum[i]&3;
- else break;
- }
- }
- if (j < 4 && i >= 0) call->unseen = j, call->a[j++] = sum[i]&3;
- call->n_alleles = j;
- } else {
- for (i = 3, j = 0; i >= 0; --i)
+ call->a[0] = ref4;
+ for (i = 3, j = 1; i >= 0; --i) {
+ if ((sum[i]&3) != ref4) {