X-Git-Url: https://git.donarmstrong.com/?p=samtools.git;a=blobdiff_plain;f=bcftools%2Fprob1.c;h=e3d6b5eb42fee41bb3b2e8f392b5ae38d1934bc5;hp=12894371c2eee9a6c5e05a5def471da413735f9b;hb=da45f88bac09817fcd87916d3e104f0e4bc53874;hpb=faf76ae18e5ddf2d03253e2ed77f444bcdf4914c diff --git a/bcftools/prob1.c b/bcftools/prob1.c index 1289437..e3d6b5e 100644 --- a/bcftools/prob1.c +++ b/bcftools/prob1.c @@ -401,7 +401,7 @@ int call_multiallelic_gt(bcf1_t *b, bcf_p1aux_t *ma, double threshold) lk = -log(1-lk/lk_sum)/0.2302585; if ( idp>=0 && ((uint16_t*)b->gi[idp].data)[isample]==0 ) { - ((uint8_t*)b->gi[old_n_gi].data)[isample] = als | 1<<7; + ((uint8_t*)b->gi[old_n_gi].data)[isample] = 1<<7; ((uint8_t*)b->gi[old_n_gi+1].data)[isample] = 0; continue; } @@ -447,9 +447,9 @@ int call_multiallelic_gt(bcf1_t *b, bcf_p1aux_t *ma, double threshold) if ( a.is_tested) ksprintf(&s, ";PV4=%.2g,%.2g,%.2g,%.2g", a.p[0], a.p[1], a.p[2], a.p[3]); ksprintf(&s, ";DP4=%d,%d,%d,%d;MQ=%d", a.d[0], a.d[1], a.d[2], a.d[3], a.mq); } + kputc('\0', &s); rm_info(&s, "I16="); rm_info(&s, "QS="); - kputc('\0', &s); } kputs(b->fmt, &s); kputc('\0', &s); free(b->str);