TimD1/vcfdist

Fix incredibly rare edge case when assigning credit to variants at same position

TimD1 opened this issue · 1 comments

TimD1 commented

When comparing PAV to T2T, there are 3 cases (out of 7 million variants) where adjacent INS+DEL variants cause issues with credit assignment. This happens at

  • sc 124112, chr5:80,538,408
  • sc 77422, chr11:41,966,547
  • sc 1741, chr17:1,131,664

To view these errors, uncomment the "Non-zero edit distance..." warning on line 1235 of dist.cpp.

TimD1 commented

Closing, since v2.5.0 updates the behavior of sync points and fixes these issues.