- f0 = acnt[ma->ref] + acnt[ma->alt] == 0? -1.
- : (double)acnt[ref] / (acnt[ma->ref] + acnt[ma->alt]);
+ if (ma->rcnt[ma->ref] + ma->rcnt[ma->alt]) // never used...
+ fr = (double)ma->rcnt[ref] / (ma->rcnt[ma->ref] + ma->rcnt[ma->alt]);
+ if (acnt[ma->ref] + acnt[ma->alt])
+ f0 = (double)acnt[ref] / (acnt[ma->ref] + acnt[ma->alt]);