X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=bam2bcf.c;h=bd5503d7e55cef2479c1a8d12d3dc0855eab06f2;hb=27eec65490fb3461f0af054868d9a0949e1843ed;hp=98acb5dd2ff051eacb86f8efd72c00c370798617;hpb=29e08d7455e9b283ef4650af7fde6fe0310d8158;p=samtools.git diff --git a/bam2bcf.c b/bam2bcf.c index 98acb5d..bd5503d 100644 --- a/bam2bcf.c +++ b/bam2bcf.c @@ -144,8 +144,8 @@ int bcf_call_combine(int n, const bcf_callret1_t *calls, int ref_base /*4-bit*/, x = call->n_alleles * (call->n_alleles + 1) / 2; // get the possible genotypes for (i = z = 0; i < call->n_alleles; ++i) - for (j = i; j < call->n_alleles; ++j) - g[z++] = call->a[i] * 5 + call->a[j]; + for (j = 0; j <= i; ++j) + g[z++] = call->a[j] * 5 + call->a[i]; for (i = 0; i < n; ++i) { uint8_t *PL = call->PL + x * i; const bcf_callret1_t *r = calls + i;