new algorithm based on PairHMM/impa 13939992d