March 29 2023 4:41:55.534 PM ss_gg_align_test(): FORTRAN90 version Test ss_gg_align(). TEST01: SS_GG_FSQ - Forward score quadratic; SS_GG_FSL - Forward score linear; GAP_OPEN penalty = -2.00000 GAP_EXTEND penalty = -0.500000 Matching scores by PAM120_SCORE. Verify that the FSQ and FSL tables agree. Sequences A and B: 1 G C 2 C C 3 T G 4 G T 5 A G 6 T A 7 A T 8 T T 9 A A 10 G G 11 C C 12 T T Matching Scores: 0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1 0.0 -5.0 -5.0 5.0 -1.0 5.0 1.0 -1.0 -1.0 1.0 5.0 -5.0 -1.0 2 0.0 9.0 9.0 -5.0 -3.0 -5.0 -3.0 -3.0 -3.0 -3.0 -5.0 9.0 -3.0 3 0.0 -3.0 -3.0 -1.0 4.0 -1.0 1.0 4.0 4.0 1.0 -1.0 -3.0 4.0 4 0.0 -5.0 -5.0 5.0 -1.0 5.0 1.0 -1.0 -1.0 1.0 5.0 -5.0 -1.0 5 0.0 -3.0 -3.0 1.0 1.0 1.0 3.0 1.0 1.0 3.0 1.0 -3.0 1.0 6 0.0 -3.0 -3.0 -1.0 4.0 -1.0 1.0 4.0 4.0 1.0 -1.0 -3.0 4.0 7 0.0 -3.0 -3.0 1.0 1.0 1.0 3.0 1.0 1.0 3.0 1.0 -3.0 1.0 8 0.0 -3.0 -3.0 -1.0 4.0 -1.0 1.0 4.0 4.0 1.0 -1.0 -3.0 4.0 9 0.0 -3.0 -3.0 1.0 1.0 1.0 3.0 1.0 1.0 3.0 1.0 -3.0 1.0 10 0.0 -5.0 -5.0 5.0 -1.0 5.0 1.0 -1.0 -1.0 1.0 5.0 -5.0 -1.0 11 0.0 9.0 9.0 -5.0 -3.0 -5.0 -3.0 -3.0 -3.0 -3.0 -5.0 9.0 -3.0 12 0.0 -3.0 -3.0 -1.0 4.0 -1.0 1.0 4.0 4.0 1.0 -1.0 -3.0 4.0 SS_GG_FSQ: 0 1 2 3 4 5 6 7 8 9 10 11 12 SF: 0 0.0 -2.5 -3.0 -3.5 -4.0 -4.5 -5.0 -5.5 -6.0 -6.5 -7.0 -7.5 -8.0 1 -2.5 -5.0 -5.5 2.0 -0.5 1.0 -1.5 -2.0 -2.5 -3.0 -1.5 -4.0 -4.5 2 -3.0 6.5 4.0 3.5 3.0 2.5 2.0 1.5 1.0 0.5 0.0 7.5 5.0 3 -3.5 4.0 3.5 3.0 7.5 5.0 4.5 6.0 5.5 3.0 2.5 5.0 11.5 4 -4.0 3.5 1.0 8.5 6.0 12.5 10.0 9.5 9.0 8.5 8.0 7.5 9.0 5 -4.5 3.0 0.5 6.0 9.5 10.0 15.5 13.0 12.5 12.0 11.5 11.0 10.5 6 -5.0 2.5 0.0 5.5 10.0 9.5 13.0 19.5 17.0 16.5 16.0 15.5 15.0 7 -5.5 2.0 -0.5 5.0 7.5 11.0 12.5 17.0 20.5 20.0 17.5 17.0 16.5 8 -6.0 1.5 -1.0 4.5 9.0 8.5 12.0 16.5 21.0 21.5 19.0 18.5 21.0 9 -6.5 1.0 -1.5 4.0 6.5 10.0 11.5 16.0 18.5 24.0 22.5 21.0 20.5 10 -7.0 0.5 -2.0 3.5 6.0 11.5 11.0 15.5 18.0 21.5 29.0 26.5 26.0 11 -7.5 2.0 9.5 7.0 6.5 9.0 10.5 15.0 17.5 21.0 26.5 38.0 35.5 12 -8.0 -0.5 7.0 8.5 11.0 8.5 10.0 14.5 19.0 20.5 26.0 35.5 42.0 EF: 0 0.0 -2.5 -3.0 -3.5 -4.0 -4.5 -5.0 -5.5 -6.0 -6.5 -7.0 -7.5 -8.0 1 -4.5 -5.0 -5.5 -6.0 -0.5 -1.0 -1.5 -2.0 -2.5 -3.0 -3.5 -4.0 -4.5 2 -5.0 -5.5 4.0 3.5 3.0 2.5 2.0 1.5 1.0 0.5 0.0 -0.5 5.0 3 -5.5 -6.0 1.5 1.0 0.5 5.0 4.5 4.0 3.5 3.0 2.5 2.0 2.5 4 -6.0 -6.5 1.0 0.5 6.0 5.5 10.0 9.5 9.0 8.5 8.0 7.5 7.0 5 -6.5 -7.0 0.5 0.0 3.5 7.0 7.5 13.0 12.5 12.0 11.5 11.0 10.5 6 -7.0 -7.5 0.0 -0.5 3.0 7.5 7.0 10.5 17.0 16.5 16.0 15.5 15.0 7 -7.5 -8.0 -0.5 -1.0 2.5 5.0 8.5 10.0 14.5 18.0 17.5 17.0 16.5 8 -8.0 -8.5 -1.0 -1.5 2.0 6.5 6.0 9.5 14.0 18.5 19.0 18.5 18.0 9 -8.5 -9.0 -1.5 -2.0 1.5 4.0 7.5 9.0 13.5 16.0 21.5 21.0 20.5 10 -9.0 -9.5 -2.0 -2.5 1.0 3.5 9.0 8.5 13.0 15.5 19.0 26.5 26.0 11 -9.5-10.0 -0.5 7.0 6.5 6.0 6.5 8.0 12.5 15.0 18.5 24.0 35.5 12-10.0-10.5 -3.0 4.5 6.0 8.5 8.0 7.5 12.0 16.5 18.0 23.5 33.0 FF: 0 0.0 -4.5 -5.0 -5.5 -6.0 -6.5 -7.0 -7.5 -8.0 -8.5 -9.0 -9.5-10.0 1 -2.5 -5.0 -5.5 -6.0 -6.5 -7.0 -7.5 -8.0 -8.5 -9.0 -9.5-10.0-10.5 2 -3.0 -5.5 -6.0 -0.5 -3.0 -1.5 -4.0 -4.5 -5.0 -5.5 -4.0 -6.5 -7.0 3 -3.5 4.0 1.5 1.0 0.5 0.0 -0.5 -1.0 -1.5 -2.0 -2.5 5.0 2.5 4 -4.0 3.5 1.0 0.5 5.0 2.5 2.0 3.5 3.0 0.5 0.0 4.5 9.0 5 -4.5 3.0 0.5 6.0 4.5 10.0 7.5 7.0 6.5 6.0 5.5 5.0 8.5 6 -5.0 2.5 0.0 5.5 7.0 9.5 13.0 10.5 10.0 9.5 9.0 8.5 8.0 7 -5.5 2.0 -0.5 5.0 7.5 9.0 12.5 17.0 14.5 14.0 13.5 13.0 12.5 8 -6.0 1.5 -1.0 4.5 7.0 8.5 12.0 16.5 18.0 17.5 15.0 14.5 14.0 9 -6.5 1.0 -1.5 4.0 6.5 8.0 11.5 16.0 18.5 19.0 16.5 16.0 18.5 10 -7.0 0.5 -2.0 3.5 6.0 7.5 11.0 15.5 18.0 21.5 20.0 18.5 18.0 11 -7.5 0.0 -2.5 3.0 5.5 9.0 10.5 15.0 17.5 21.0 26.5 24.0 23.5 12 -8.0 -0.5 7.0 4.5 5.0 8.5 10.0 14.5 17.0 20.5 26.0 35.5 33.0 TF: 0 0 2 2 2 2 2 2 2 2 2 2 2 2 1 4 7 6 1 2 1 2 2 2 2 1 2 2 2 4 1 3 2 2 2 2 2 2 2 2 1 2 3 4 4 1 1 1 2 2 1 1 2 2 4 1 4 4 4 6 1 2 1 2 2 2 2 3 2 4 5 4 4 7 4 1 4 1 2 2 3 2 2 2 6 4 4 7 4 1 4 4 1 3 2 2 2 3 7 4 4 7 4 4 1 5 4 1 1 3 2 3 8 4 4 7 4 1 4 5 5 1 1 3 2 1 9 4 4 7 4 4 1 5 4 4 1 1 2 2 10 4 4 6 5 4 1 5 4 4 4 1 2 2 11 4 1 1 2 2 4 4 4 4 4 4 1 2 12 4 4 4 1 1 6 5 5 1 4 4 4 1 SS_GG_FSL: 0 1 2 3 4 5 6 7 8 9 10 11 12 SF: 0 0.0 -2.5 -3.0 -3.5 -4.0 -4.5 -5.0 -5.5 -6.0 -6.5 -7.0 -7.5 -8.0 1 -2.5 -5.0 -5.5 2.0 -0.5 1.0 -1.5 -2.0 -2.5 -3.0 -1.5 -4.0 -4.5 2 -3.0 6.5 4.0 3.5 3.0 2.5 2.0 1.5 1.0 0.5 0.0 7.5 5.0 3 -3.5 4.0 3.5 3.0 7.5 5.0 4.5 6.0 5.5 3.0 2.5 5.0 11.5 4 -4.0 3.5 1.0 8.5 6.0 12.5 10.0 9.5 9.0 8.5 8.0 7.5 9.0 5 -4.5 3.0 0.5 6.0 9.5 10.0 15.5 13.0 12.5 12.0 11.5 11.0 10.5 6 -5.0 2.5 0.0 5.5 10.0 9.5 13.0 19.5 17.0 16.5 16.0 15.5 15.0 7 -5.5 2.0 -0.5 5.0 7.5 11.0 12.5 17.0 20.5 20.0 17.5 17.0 16.5 8 -6.0 1.5 -1.0 4.5 9.0 8.5 12.0 16.5 21.0 21.5 19.0 18.5 21.0 9 -6.5 1.0 -1.5 4.0 6.5 10.0 11.5 16.0 18.5 24.0 22.5 21.0 20.5 10 -7.0 0.5 -2.0 3.5 6.0 11.5 11.0 15.5 18.0 21.5 29.0 26.5 26.0 11 -7.5 2.0 9.5 7.0 6.5 9.0 10.5 15.0 17.5 21.0 26.5 38.0 35.5 12 -8.0 -0.5 7.0 8.5 11.0 8.5 10.0 14.5 19.0 20.5 26.0 35.5 42.0 EF: 0 0.0 -2.5 -3.0 -3.5 -4.0 -4.5 -5.0 -5.5 -6.0 -6.5 -7.0 -7.5 -8.0 1 -4.5 -5.0 -5.5 -6.0 -0.5 -1.0 -1.5 -2.0 -2.5 -3.0 -3.5 -4.0 -4.5 2 -5.0 -5.5 4.0 3.5 3.0 2.5 2.0 1.5 1.0 0.5 0.0 -0.5 5.0 3 -5.5 -6.0 1.5 1.0 0.5 5.0 4.5 4.0 3.5 3.0 2.5 2.0 2.5 4 -6.0 -6.5 1.0 0.5 6.0 5.5 10.0 9.5 9.0 8.5 8.0 7.5 7.0 5 -6.5 -7.0 0.5 0.0 3.5 7.0 7.5 13.0 12.5 12.0 11.5 11.0 10.5 6 -7.0 -7.5 0.0 -0.5 3.0 7.5 7.0 10.5 17.0 16.5 16.0 15.5 15.0 7 -7.5 -8.0 -0.5 -1.0 2.5 5.0 8.5 10.0 14.5 18.0 17.5 17.0 16.5 8 -8.0 -8.5 -1.0 -1.5 2.0 6.5 6.0 9.5 14.0 18.5 19.0 18.5 18.0 9 -8.5 -9.0 -1.5 -2.0 1.5 4.0 7.5 9.0 13.5 16.0 21.5 21.0 20.5 10 -9.0 -9.5 -2.0 -2.5 1.0 3.5 9.0 8.5 13.0 15.5 19.0 26.5 26.0 11 -9.5-10.0 -0.5 7.0 6.5 6.0 6.5 8.0 12.5 15.0 18.5 24.0 35.5 12-10.0-10.5 -3.0 4.5 6.0 8.5 8.0 7.5 12.0 16.5 18.0 23.5 33.0 FF: 0 0.0 -4.5 -5.0 -5.5 -6.0 -6.5 -7.0 -7.5 -8.0 -8.5 -9.0 -9.5-10.0 1 -2.5 -5.0 -5.5 -6.0 -6.5 -7.0 -7.5 -8.0 -8.5 -9.0 -9.5-10.0-10.5 2 -3.0 -5.5 -6.0 -0.5 -3.0 -1.5 -4.0 -4.5 -5.0 -5.5 -4.0 -6.5 -7.0 3 -3.5 4.0 1.5 1.0 0.5 0.0 -0.5 -1.0 -1.5 -2.0 -2.5 5.0 2.5 4 -4.0 3.5 1.0 0.5 5.0 2.5 2.0 3.5 3.0 0.5 0.0 4.5 9.0 5 -4.5 3.0 0.5 6.0 4.5 10.0 7.5 7.0 6.5 6.0 5.5 5.0 8.5 6 -5.0 2.5 0.0 5.5 7.0 9.5 13.0 10.5 10.0 9.5 9.0 8.5 8.0 7 -5.5 2.0 -0.5 5.0 7.5 9.0 12.5 17.0 14.5 14.0 13.5 13.0 12.5 8 -6.0 1.5 -1.0 4.5 7.0 8.5 12.0 16.5 18.0 17.5 15.0 14.5 14.0 9 -6.5 1.0 -1.5 4.0 6.5 8.0 11.5 16.0 18.5 19.0 16.5 16.0 18.5 10 -7.0 0.5 -2.0 3.5 6.0 7.5 11.0 15.5 18.0 21.5 20.0 18.5 18.0 11 -7.5 0.0 -2.5 3.0 5.5 9.0 10.5 15.0 17.5 21.0 26.5 24.0 23.5 12 -8.0 -0.5 7.0 4.5 5.0 8.5 10.0 14.5 17.0 20.5 26.0 35.5 33.0 TF: 0 0 2 2 2 2 2 2 2 2 2 2 2 2 1 4 7 6 1 2 1 2 2 2 2 1 2 2 2 4 1 3 2 2 2 2 2 2 2 2 1 2 3 4 4 1 1 1 2 2 1 1 2 2 4 1 4 4 4 6 1 2 1 2 2 2 2 3 2 4 5 4 4 7 4 1 4 1 2 2 3 2 2 2 6 4 4 7 4 1 4 4 1 3 2 2 2 3 7 4 4 7 4 4 1 5 4 1 1 3 2 3 8 4 4 7 4 1 4 5 5 1 1 3 2 1 9 4 4 7 4 4 1 5 4 4 1 1 2 2 10 4 4 6 5 4 1 5 4 4 4 1 2 2 11 4 1 1 2 2 4 4 4 4 4 4 1 2 12 4 4 4 1 1 6 5 5 1 4 4 4 1 TEST02: SS_GG_BSQ - Backward score quadratic; SS_GG_BSL - Backward score linear. GAP_OPEN penalty = -2.00000 GAP_EXTEND penalty = -0.500000 Matching scores by PAM120_SCORE. Verify that the BSQ and BSL tables agree. Sequences A and B: 1 G C 2 C C 3 T G 4 G T 5 A G 6 T A 7 A T 8 T T 9 A A 10 G G 11 C C 12 T T SS_GG_BSQ: 0 1 2 3 4 5 6 7 8 9 10 11 12 SB: 0 42.0 42.5 41.0 35.5 31.0 25.5 22.0 19.5 15.0 11.5 6.0 -3.5 -8.0 1 44.5 45.0 33.5 36.0 31.5 26.0 22.5 20.0 15.5 12.0 6.5 -3.0 -7.5 2 35.0 35.5 36.0 38.5 32.0 26.5 25.0 20.5 16.0 12.5 7.0 -2.5 -7.0 3 30.5 31.0 31.5 32.0 34.5 27.0 23.5 21.0 16.5 13.0 7.5 -2.0 -6.5 4 25.0 25.5 28.0 26.5 28.0 29.5 24.0 21.5 17.0 13.5 8.0 -1.5 -6.0 5 23.5 24.0 24.5 27.0 24.5 27.0 26.5 22.0 17.5 14.0 8.5 -1.0 -5.5 6 21.0 21.5 22.0 22.5 23.0 25.5 26.0 22.5 18.0 14.5 9.0 -0.5 -5.0 7 19.5 20.0 20.5 21.0 21.5 22.0 22.5 25.0 18.5 15.0 9.5 0.0 -4.5 8 15.0 15.5 16.0 16.5 17.0 17.5 18.0 18.5 21.0 15.5 10.0 0.5 -4.0 9 11.5 12.0 12.5 13.0 13.5 14.0 14.5 15.0 15.5 18.0 10.5 1.0 -3.5 10 6.0 6.5 7.0 7.5 8.0 8.5 9.0 9.5 10.0 10.5 13.0 1.5 -3.0 11 -3.5 -3.0 -2.5 -2.0 -1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 4.0 -2.5 12 -8.0 -7.5 -7.0 -6.5 -6.0 -5.5 -5.0 -4.5 -4.0 -3.5 -3.0 -2.5 0.0 EB: 0 40.0 38.5 33.0 28.5 23.0 19.5 17.0 12.5 9.0 3.5 -6.0-10.5-10.0 1 42.5 33.0 33.5 29.0 23.5 20.0 17.5 13.0 9.5 4.0 -5.5-10.0 -9.5 2 35.0 35.5 36.0 29.5 24.0 22.5 18.0 13.5 10.0 4.5 -5.0 -9.5 -9.0 3 30.5 31.0 31.5 32.0 24.5 21.0 18.5 14.0 10.5 5.0 -4.5 -9.0 -8.5 4 25.0 25.5 26.0 26.5 27.0 21.5 19.0 14.5 11.0 5.5 -4.0 -8.5 -8.0 5 23.5 24.0 24.5 24.0 24.5 24.0 19.5 15.0 11.5 6.0 -3.5 -8.0 -7.5 6 21.0 21.5 22.0 22.5 23.0 23.5 20.0 15.5 12.0 6.5 -3.0 -7.5 -7.0 7 19.5 20.0 20.5 21.0 21.5 22.0 22.5 16.0 12.5 7.0 -2.5 -7.0 -6.5 8 15.0 15.5 16.0 16.5 17.0 17.5 18.0 18.5 13.0 7.5 -2.0 -6.5 -6.0 9 11.5 12.0 12.5 13.0 13.5 14.0 14.5 15.0 15.5 8.0 -1.5 -6.0 -5.5 10 6.0 6.5 7.0 7.5 8.0 8.5 9.0 9.5 10.0 10.5 -1.0 -5.5 -5.0 11 -3.5 -3.0 -2.5 -2.0 -1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 -5.0 -4.5 12 -8.0 -7.5 -7.0 -6.5 -6.0 -5.5 -5.0 -4.5 -4.0 -3.5 -3.0 -2.5 0.0 FB: 0 42.0 42.5 33.0 35.5 31.0 25.5 22.0 19.5 15.0 11.5 6.0 -3.5 -8.0 1 32.5 33.0 33.5 36.0 31.5 26.0 22.5 20.0 15.5 12.0 6.5 -3.0 -7.5 2 28.0 28.5 29.0 29.5 32.0 26.5 23.0 20.5 16.0 12.5 7.0 -2.5 -7.0 3 22.5 23.0 25.5 24.0 25.5 27.0 23.5 21.0 16.5 13.0 7.5 -2.0 -6.5 4 21.0 21.5 22.0 24.5 22.0 24.5 24.0 21.5 17.0 13.5 8.0 -1.5 -6.0 5 18.5 19.0 19.5 20.0 20.5 23.0 23.5 22.0 17.5 14.0 8.5 -1.0 -5.5 6 17.0 17.5 18.0 18.5 19.0 19.5 20.0 22.5 18.0 14.5 9.0 -0.5 -5.0 7 12.5 13.0 13.5 14.0 14.5 15.0 15.5 16.0 18.5 15.0 9.5 0.0 -4.5 8 9.0 9.5 10.0 10.5 11.0 11.5 12.0 12.5 13.0 15.5 10.0 0.5 -4.0 9 3.5 4.0 4.5 5.0 5.5 6.0 6.5 7.0 7.5 8.0 10.5 1.0 -3.5 10 -6.0 -5.5 -5.0 -4.5 -4.0 -3.5 -3.0 -2.5 -2.0 -1.5 -1.0 1.5 -3.0 11-10.5-10.0 -9.5 -9.0 -8.5 -8.0 -7.5 -7.0 -6.5 -6.0 -5.5 -5.0 -2.5 12-10.0 -9.5 -9.0 -8.5 -8.0 -7.5 -7.0 -6.5 -6.0 -5.5 -5.0 -4.5 0.0 TB: 0 4 4 1 4 5 4 4 4 4 5 4 4 4 1 1 1 7 4 4 4 4 4 4 4 5 4 4 2 2 2 2 1 4 4 1 5 4 4 4 5 4 3 2 2 3 2 1 4 4 4 4 5 4 4 4 4 2 2 1 2 1 1 4 4 5 4 4 4 4 5 2 2 2 1 3 1 1 5 4 4 4 5 4 6 2 2 3 3 3 1 1 4 5 4 4 4 4 7 2 2 2 3 2 2 3 1 4 4 4 5 4 8 2 2 2 2 2 3 2 2 1 4 4 4 4 9 2 2 3 2 3 2 2 2 2 1 4 4 4 10 3 3 2 2 2 2 2 2 2 2 1 4 4 11 2 2 2 3 2 2 3 3 2 2 2 1 4 12 2 2 2 2 2 2 2 2 2 2 2 2 0 SS_GG_BSL: 0 1 2 3 4 5 6 7 8 9 10 11 12 SB: 0 42.0 42.5 41.0 35.5 31.0 25.5 22.0 19.5 15.0 11.5 6.0 -3.5 -8.0 1 44.5 45.0 33.5 36.0 31.5 26.0 22.5 20.0 15.5 12.0 6.5 -3.0 -7.5 2 35.0 35.5 36.0 38.5 32.0 26.5 25.0 20.5 16.0 12.5 7.0 -2.5 -7.0 3 30.5 31.0 31.5 32.0 34.5 27.0 23.5 21.0 16.5 13.0 7.5 -2.0 -6.5 4 25.0 25.5 28.0 26.5 28.0 29.5 24.0 21.5 17.0 13.5 8.0 -1.5 -6.0 5 23.5 24.0 24.5 27.0 24.5 27.0 26.5 22.0 17.5 14.0 8.5 -1.0 -5.5 6 21.0 21.5 22.0 22.5 23.0 25.5 26.0 22.5 18.0 14.5 9.0 -0.5 -5.0 7 19.5 20.0 20.5 21.0 21.5 22.0 22.5 25.0 18.5 15.0 9.5 0.0 -4.5 8 15.0 15.5 16.0 16.5 17.0 17.5 18.0 18.5 21.0 15.5 10.0 0.5 -4.0 9 11.5 12.0 12.5 13.0 13.5 14.0 14.5 15.0 15.5 18.0 10.5 1.0 -3.5 10 6.0 6.5 7.0 7.5 8.0 8.5 9.0 9.5 10.0 10.5 13.0 1.5 -3.0 11 -3.5 -3.0 -2.5 -2.0 -1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 4.0 -2.5 12 -8.0 -7.5 -7.0 -6.5 -6.0 -5.5 -5.0 -4.5 -4.0 -3.5 -3.0 -2.5 0.0 EB: 0 40.0 38.5 33.0 28.5 23.0 19.5 17.0 12.5 9.0 3.5 -6.0-10.5-10.0 1 42.5 33.0 33.5 29.0 23.5 20.0 17.5 13.0 9.5 4.0 -5.5-10.0 -9.5 2 35.0 35.5 36.0 29.5 24.0 22.5 18.0 13.5 10.0 4.5 -5.0 -9.5 -9.0 3 30.5 31.0 31.5 32.0 24.5 21.0 18.5 14.0 10.5 5.0 -4.5 -9.0 -8.5 4 25.0 25.5 26.0 26.5 27.0 21.5 19.0 14.5 11.0 5.5 -4.0 -8.5 -8.0 5 23.5 24.0 24.5 24.0 24.5 24.0 19.5 15.0 11.5 6.0 -3.5 -8.0 -7.5 6 21.0 21.5 22.0 22.5 23.0 23.5 20.0 15.5 12.0 6.5 -3.0 -7.5 -7.0 7 19.5 20.0 20.5 21.0 21.5 22.0 22.5 16.0 12.5 7.0 -2.5 -7.0 -6.5 8 15.0 15.5 16.0 16.5 17.0 17.5 18.0 18.5 13.0 7.5 -2.0 -6.5 -6.0 9 11.5 12.0 12.5 13.0 13.5 14.0 14.5 15.0 15.5 8.0 -1.5 -6.0 -5.5 10 6.0 6.5 7.0 7.5 8.0 8.5 9.0 9.5 10.0 10.5 -1.0 -5.5 -5.0 11 -3.5 -3.0 -2.5 -2.0 -1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 -5.0 -4.5 12 -8.0 -7.5 -7.0 -6.5 -6.0 -5.5 -5.0 -4.5 -4.0 -3.5 -3.0 -2.5 0.0 FB: 0 42.0 42.5 33.0 35.5 31.0 25.5 22.0 19.5 15.0 11.5 6.0 -3.5 -8.0 1 32.5 33.0 33.5 36.0 31.5 26.0 22.5 20.0 15.5 12.0 6.5 -3.0 -7.5 2 28.0 28.5 29.0 29.5 32.0 26.5 23.0 20.5 16.0 12.5 7.0 -2.5 -7.0 3 22.5 23.0 25.5 24.0 25.5 27.0 23.5 21.0 16.5 13.0 7.5 -2.0 -6.5 4 21.0 21.5 22.0 24.5 22.0 24.5 24.0 21.5 17.0 13.5 8.0 -1.5 -6.0 5 18.5 19.0 19.5 20.0 20.5 23.0 23.5 22.0 17.5 14.0 8.5 -1.0 -5.5 6 17.0 17.5 18.0 18.5 19.0 19.5 20.0 22.5 18.0 14.5 9.0 -0.5 -5.0 7 12.5 13.0 13.5 14.0 14.5 15.0 15.5 16.0 18.5 15.0 9.5 0.0 -4.5 8 9.0 9.5 10.0 10.5 11.0 11.5 12.0 12.5 13.0 15.5 10.0 0.5 -4.0 9 3.5 4.0 4.5 5.0 5.5 6.0 6.5 7.0 7.5 8.0 10.5 1.0 -3.5 10 -6.0 -5.5 -5.0 -4.5 -4.0 -3.5 -3.0 -2.5 -2.0 -1.5 -1.0 1.5 -3.0 11-10.5-10.0 -9.5 -9.0 -8.5 -8.0 -7.5 -7.0 -6.5 -6.0 -5.5 -5.0 -2.5 12-10.0 -9.5 -9.0 -8.5 -8.0 -7.5 -7.0 -6.5 -6.0 -5.5 -5.0 -4.5 0.0 TB: 0 4 4 1 4 5 4 4 4 4 5 4 4 4 1 1 1 7 4 4 4 4 4 4 4 5 4 4 2 2 2 2 1 4 4 1 5 4 4 4 5 4 3 2 2 3 2 1 4 4 4 4 5 4 4 4 4 2 2 1 2 1 1 4 4 5 4 4 4 4 5 2 2 2 1 3 1 1 5 4 4 4 5 4 6 2 2 3 3 3 1 1 4 5 4 4 4 4 7 2 2 2 3 2 2 3 1 4 4 4 5 4 8 2 2 2 2 2 3 2 2 1 4 4 4 4 9 2 2 3 2 3 2 2 2 2 1 4 4 4 10 3 3 2 2 2 2 2 2 2 2 1 4 4 11 2 2 2 3 2 2 3 3 2 2 2 1 4 12 2 2 2 2 2 2 2 2 2 2 2 2 0 TEST03: SS_GG_FSQ - Forward score quadratic; SS_GG_FSL - Forward score linear; GAP_OPEN penalty = -2.00000 GAP_EXTEND penalty = -0.500000 Matching scores by PAM120_SCORE. Verify that the FSQ and FSL tables agree. Sequences A and B: 1 G G 2 C G 3 T G 4 T Matching scores: 0 0.0 0.0 0.0 0.0 0.0 1 0.0 5.0 5.0 5.0 -1.0 2 0.0 -5.0 -5.0 -5.0 -3.0 3 0.0 -1.0 -1.0 -1.0 4.0 SS_GG_FSQ: 0 1 2 3 4 SF: 0 0.0 -2.5 -3.0 -3.5 -4.0 1 -2.5 5.0 2.5 2.0 1.5 2 -3.0 2.5 0.0 -0.5 -1.0 3 -3.5 2.0 1.5 -1.0 3.5 EF: 0 0.0 -2.5 -3.0 -3.5 -4.0 1 -4.5 -5.0 2.5 2.0 1.5 2 -5.0 -5.5 0.0 -0.5 -1.0 3 -5.5 -6.0 -0.5 -1.0 -1.5 FF: 0 0.0 -4.5 -5.0 -5.5 -6.0 1 -2.5 -5.0 -5.5 -6.0 -6.5 2 -3.0 2.5 0.0 -0.5 -1.0 3 -3.5 2.0 -0.5 -1.0 -1.5 TF: 0 0 2 2 2 2 1 4 1 3 3 2 2 4 4 7 6 7 3 4 4 1 7 1 SS_GG_FSL: 0 1 2 3 4 SF: 0 0.0 -2.5 -3.0 -3.5 -4.0 1 -2.5 5.0 2.5 2.0 1.5 2 -3.0 2.5 0.0 -0.5 -1.0 3 -3.5 2.0 1.5 -1.0 3.5 EF: 0 0.0 -2.5 -3.0 -3.5 -4.0 1 -4.5 -5.0 2.5 2.0 1.5 2 -5.0 -5.5 0.0 -0.5 -1.0 3 -5.5 -6.0 -0.5 -1.0 -1.5 FF: 0 0.0 -4.5 -5.0 -5.5 -6.0 1 -2.5 -5.0 -5.5 -6.0 -6.5 2 -3.0 2.5 0.0 -0.5 -1.0 3 -3.5 2.0 -0.5 -1.0 -1.5 TF: 0 0 2 2 2 2 1 4 1 3 3 2 2 4 4 7 6 7 3 4 4 1 7 1 TEST04: SS_GG_BSQ - Backward score quadratic; SS_GG_BSL - Backward score linear. GAP_OPEN penalty = -2.00000 GAP_EXTEND penalty = -0.500000 Matching scores by PAM120_SCORE. Verify that the BSQ and BSL tables agree. Sequences A and B: 1 G G 2 C G 3 T G 4 T SS_GG_BSQ: 0 1 2 3 4 SB: 0 3.5 4.0 6.5 1.0 -3.5 1 -2.0 -1.5 -1.0 1.5 -3.0 2 0.5 1.0 1.5 4.0 -2.5 3 -4.0 -3.5 -3.0 -2.5 0.0 EB: 0 3.5 4.0 -1.5 -6.0 -5.5 1 -2.0 -1.5 -1.0 -5.5 -5.0 2 0.5 1.0 1.5 -5.0 -4.5 3 -4.0 -3.5 -3.0 -2.5 0.0 FB: 0 -2.5 -2.0 -1.5 1.0 -3.5 1 -2.0 -1.5 -1.0 1.5 -3.0 2 -6.5 -6.0 -5.5 -5.0 -2.5 3 -6.0 -5.5 -5.0 -4.5 0.0 TB: 0 3 3 1 4 4 1 6 6 7 4 4 2 2 2 2 1 4 3 2 2 2 2 0 SS_GG_BSL: 0 1 2 3 4 SB: 0 3.5 4.0 6.5 1.0 -3.5 1 -2.0 -1.5 -1.0 1.5 -3.0 2 0.5 1.0 1.5 4.0 -2.5 3 -4.0 -3.5 -3.0 -2.5 0.0 EB: 0 3.5 4.0 -1.5 -6.0 -5.5 1 -2.0 -1.5 -1.0 -5.5 -5.0 2 0.5 1.0 1.5 -5.0 -4.5 3 -4.0 -3.5 -3.0 -2.5 0.0 FB: 0 -2.5 -2.0 -1.5 1.0 -3.5 1 -2.0 -1.5 -1.0 1.5 -3.0 2 -6.5 -6.0 -5.5 -5.0 -2.5 3 -6.0 -5.5 -5.0 -4.5 0.0 TB: 0 3 3 1 4 4 1 6 6 7 4 4 2 2 2 2 1 4 3 2 2 2 2 0 TEST05: SS_GG_BSQ - Backward score quadratic; SS_GG_FSQ - Forward score quadratic. GAP_OPEN penalty = -2.00000 GAP_EXTEND penalty = -0.500000 Matching scores by PAM120_SCORE. Verify that the FSQ and (reversed) BSQ tables agree. Sequences A and B: 1 G G 2 C G 3 T G 4 T SS_GG_FSQ: 0 1 2 3 4 SF: 0 0.0 -2.5 -3.0 -3.5 -4.0 1 -2.5 5.0 2.5 2.0 1.5 2 -3.0 2.5 0.0 -0.5 -1.0 3 -3.5 2.0 1.5 -1.0 3.5 EF: 0 0.0 -2.5 -3.0 -3.5 -4.0 1 -4.5 -5.0 2.5 2.0 1.5 2 -5.0 -5.5 0.0 -0.5 -1.0 3 -5.5 -6.0 -0.5 -1.0 -1.5 FF: 0 0.0 -4.5 -5.0 -5.5 -6.0 1 -2.5 -5.0 -5.5 -6.0 -6.5 2 -3.0 2.5 0.0 -0.5 -1.0 3 -3.5 2.0 -0.5 -1.0 -1.5 TF: 0 0 2 2 2 2 1 4 1 3 3 2 2 4 4 7 6 7 3 4 4 1 7 1 Reversed sequences A and B: 1 T T 2 C G 3 G G 4 G SS_GG_BSQ: (The tables are printed in reverse order) 4 3 2 1 0 SB: 3 0.0 -2.5 -3.0 -3.5 -4.0 2 -2.5 5.0 2.5 2.0 1.5 1 -3.0 2.5 0.0 -0.5 -1.0 0 -3.5 2.0 1.5 -1.0 3.5 EB: 3 0.0 -2.5 -3.0 -3.5 -4.0 2 -4.5 -5.0 2.5 2.0 1.5 1 -5.0 -5.5 0.0 -0.5 -1.0 0 -5.5 -6.0 -0.5 -1.0 -1.5 FB: 3 0.0 -4.5 -5.0 -5.5 -6.0 2 -2.5 -5.0 -5.5 -6.0 -6.5 1 -3.0 2.5 0.0 -0.5 -1.0 0 -3.5 2.0 -0.5 -1.0 -1.5 TB: 3 0 2 2 2 2 2 4 1 3 3 2 1 4 4 7 6 7 0 4 4 1 7 1 TEST06: SS_GG_FSL - Forward linear; SS_GG_BSL - Backward linear. GAP_OPEN penalty = -2.00000 GAP_EXTEND penalty = -0.500000 Matching scores by PAM120_SCORE. Do a score computation on a sub-block. Verify that the (M2,N2) entry for FSL agrees with the (M1,N1) entry for BSL. Sequences A and B: 1 G C 2 C C 3 T G 4 G T 5 A G 6 T A 7 A T 8 T T 9 A A 10 G G 11 C C 12 T T SS_GG_FSL: 3 4 5 6 7 8 9 4 0.0 -2.5 -3.0 -3.5 -4.0 -4.5 -5.0 5 -2.5 1.0 -1.5 0.0 -2.5 -3.0 -1.5 6 -3.0 1.5 0.0 -0.5 4.0 1.5 1.0 7 -3.5 -1.0 2.5 3.0 1.5 5.0 4.5 8 -4.0 0.5 0.0 3.5 7.0 5.5 6.0 9 -4.5 -2.0 1.5 3.0 4.5 8.0 8.5 10 -5.0 -2.5 3.0 2.5 4.0 5.5 9.0 SS_GG_BSL: 3 4 5 6 7 8 9 4 9.0 7.5 9.0 4.5 1.0 -1.5 -5.0 5 10.5 8.0 6.5 6.0 3.5 -3.0 -4.5 6 6.0 6.5 9.0 5.5 2.0 -0.5 -4.0 7 2.5 3.0 3.5 6.0 4.5 -2.0 -3.5 8 2.0 -1.5 1.0 -0.5 2.0 0.5 -3.0 9 -1.5 1.0 -2.5 -2.0 -1.5 1.0 -2.5 10 -5.0 -4.5 -4.0 -3.5 -3.0 -2.5 0.0 TEST07: SS_GG_FSQ - forward score quadratic SS_GG_FPQ - forward path quadratic GAP_OPEN penalty = -2.00000 GAP_EXTEND penalty = -0.500000 Matching scores by PAM120_SCORE. Sequences A and B: 1 G C 2 C C 3 T G 4 G T 5 A G 6 T A 7 A T 8 T T 9 A A 10 G G 11 C C 12 T T SS_GG_FSQ reports optimal matching score is 42.0000 Matching path: 1 0 0 2 1 -1 3 2 1 4 -1 2 5 -1 3 6 3 4 7 4 5 8 5 6 9 6 7 10 7 -1 11 8 8 12 9 9 13 10 10 14 11 11 15 12 12 Sequence/sequence matching, Affine gap penalty: # A B # Increm Score 0 0 0.00 1 G | -2.50 -2.50 2 C = C 1 9.00 6.50 | C 2 -2.50 4.00 | G 3 -0.50 3.50 3 T = T 4 4.00 7.50 4 G = G 5 5.00 12.50 5 A = A 6 3.00 15.50 6 T = T 7 4.00 19.50 7 A | -2.50 17.00 8 T = T 8 4.00 21.00 9 A = A 9 3.00 24.00 10 G = G 10 5.00 29.00 11 C = C 11 9.00 38.00 12 T = T 12 4.00 42.00 TEST08: SS_GG_BSQ - backward score quadratic SS_GG_BPQ - backward path quadratic GAP_OPEN penalty = -2.00000 GAP_EXTEND penalty = -0.500000 Matching scores by PAM120_SCORE. Compare path and score with FSQ/FPQ computation. Sequences A and B: 1 G C 2 C C 3 T G 4 G T 5 A G 6 T A 7 A T 8 T T 9 A A 10 G G 11 C C 12 T T SS_GG_BSQ reports optimal matching score is 42.0000 Matching path: 1 0 0 2 1 -1 3 2 1 4 -1 2 5 -1 3 6 3 4 7 4 5 8 5 6 9 6 7 10 7 -1 11 8 8 12 9 9 13 10 10 14 11 11 15 12 12 Sequence/sequence matching, Affine gap penalty: # A B # Increm Score 0 0 0.00 1 G | -2.50 -2.50 2 C = C 1 9.00 6.50 | C 2 -2.50 4.00 | G 3 -0.50 3.50 3 T = T 4 4.00 7.50 4 G = G 5 5.00 12.50 5 A = A 6 3.00 15.50 6 T = T 7 4.00 19.50 7 A | -2.50 17.00 8 T = T 8 4.00 21.00 9 A = A 9 3.00 24.00 10 G = G 10 5.00 29.00 11 C = C 11 9.00 38.00 12 T = T 12 4.00 42.00 TEST09: SS_GG_RPL - path routine; GAP_OPEN penalty = -2.00000 GAP_EXTEND penalty = -0.500000 Matching scores by PAM120_SCORE. Verify that RPL agrees with FSQ/FPQ and BSQ/BPQ. Sequences A and B: 1 G C 2 C C 3 T G 4 G T 5 A G 6 T A 7 A T 8 T T 9 A A 10 G G 11 C C 12 T T Matching path: 1 0 0 2 1 -1 3 2 1 4 -1 2 5 -1 3 6 3 4 7 4 5 8 5 6 9 6 7 10 7 -1 11 8 8 12 9 9 13 10 10 14 11 11 15 12 12 Sequence/sequence matching, Affine gap penalty: # A B # Increm Score 0 0 0.00 1 G | -2.50 -2.50 2 C = C 1 9.00 6.50 | C 2 -2.50 4.00 | G 3 -0.50 3.50 3 T = T 4 4.00 7.50 4 G = G 5 5.00 12.50 5 A = A 6 3.00 15.50 6 T = T 7 4.00 19.50 7 A | -2.50 17.00 8 T = T 8 4.00 21.00 9 A = A 9 3.00 24.00 10 G = G 10 5.00 29.00 11 C = C 11 9.00 38.00 12 T = T 12 4.00 42.00 TEST12: SS_GG_FSQ - Forward score quadratic; SS_GG_FPQ - Forward path quadratic. GAP_OPEN penalty = -2.00000 GAP_EXTEND penalty = -0.500000 Matching scores by PAM120_SCORE. Sequence/sequence matching, Affine gap penalty: # A B # Increm Score 0 0 0.00 1 M = M 1 8.00 8.00 2 M = M 2 8.00 16.00 3 A = A 3 3.00 19.00 4 A = A 4 3.00 22.00 5 E = E 5 5.00 27.00 6 A = A 6 3.00 30.00 7 G = G 7 5.00 35.00 8 G = G 8 5.00 40.00 9 E | -2.50 37.50 10 E | -0.50 37.00 11 G | -0.50 36.50 12 G | -0.50 36.00 13 P = P 9 6.00 42.00 14 V = V 10 5.00 47.00 15 T = T 11 4.00 51.00 16 A = A 12 3.00 54.00 17 G = G 13 5.00 59.00 18 A = A 14 3.00 62.00 19 A = A 15 3.00 65.00 20 G = G 16 5.00 70.00 21 G = G 17 5.00 75.00 22 G = G 18 5.00 80.00 23 A = A 19 3.00 83.00 24 A = A 20 3.00 86.00 | C 21 -2.50 83.50 | C 22 -0.50 83.00 | C 23 -0.50 82.50 25 A = A 24 3.00 85.50 26 A = A 25 3.00 88.50 27 S = S 26 3.00 91.50 28 G = G 27 5.00 96.50 29 A = A 28 3.00 99.50 30 Y = Y 29 8.00 107.50 31 P = P 30 6.00 113.50 32 A = A 31 3.00 116.50 33 V = V 32 5.00 121.50 34 C = C 33 9.00 130.50 35 R = R 34 6.00 136.50 36 V = V 35 5.00 141.50 37 K = K 36 5.00 146.50 38 I = I 37 6.00 152.50 39 P = P 38 6.00 158.50 40 A = A 39 3.00 161.50 41 A = A 40 3.00 164.50 42 L = L 41 5.00 169.50 43 P = P 42 6.00 175.50 44 V = V 43 5.00 180.50 45 A = A 44 3.00 183.50 46 A = A 45 3.00 186.50 47 A = A 46 3.00 189.50 48 A = A 47 3.00 192.50 49 P = P 48 6.00 198.50 50 F = F 49 8.00 206.50 51 P = P 50 6.00 212.50 52 G = G 51 5.00 217.50 53 L = L 52 5.00 222.50 54 A = A 53 3.00 225.50 55 E = E 54 5.00 230.50 56 A = A 55 3.00 233.50 57 G = G 56 5.00 238.50 58 V = V 57 5.00 243.50 59 A = A 58 3.00 246.50 60 A = A 59 3.00 249.50 TEST13: SS_GG_RPL - path routine; GAP_OPEN penalty = -2.00000 GAP_EXTEND penalty = -0.500000 Matching scores by PAM120_SCORE. Compare with FSQ/FPQ calculation. Sequence/sequence matching, Affine gap penalty: # A B # Increm Score 0 0 0.00 1 M = M 1 8.00 8.00 2 M = M 2 8.00 16.00 3 A = A 3 3.00 19.00 4 A = A 4 3.00 22.00 5 E = E 5 5.00 27.00 6 A = A 6 3.00 30.00 7 G | -2.50 27.50 8 G | -0.50 27.00 9 E | -0.50 26.50 10 E | -0.50 26.00 11 G = G 7 5.00 31.00 12 G = G 8 5.00 36.00 13 P = P 9 6.00 42.00 14 V = V 10 5.00 47.00 15 T = T 11 4.00 51.00 16 A = A 12 3.00 54.00 17 G = G 13 5.00 59.00 18 A = A 14 3.00 62.00 19 A = A 15 3.00 65.00 20 G = G 16 5.00 70.00 21 G = G 17 5.00 75.00 22 G = G 18 5.00 80.00 23 A = A 19 3.00 83.00 24 A = A 20 3.00 86.00 | C 21 -2.50 83.50 | C 22 -0.50 83.00 | C 23 -0.50 82.50 25 A = A 24 3.00 85.50 26 A = A 25 3.00 88.50 27 S = S 26 3.00 91.50 28 G = G 27 5.00 96.50 29 A = A 28 3.00 99.50 30 Y = Y 29 8.00 107.50 31 P = P 30 6.00 113.50 32 A = A 31 3.00 116.50 33 V = V 32 5.00 121.50 34 C = C 33 9.00 130.50 35 R = R 34 6.00 136.50 36 V = V 35 5.00 141.50 37 K = K 36 5.00 146.50 38 I = I 37 6.00 152.50 39 P = P 38 6.00 158.50 40 A = A 39 3.00 161.50 41 A = A 40 3.00 164.50 42 L = L 41 5.00 169.50 43 P = P 42 6.00 175.50 44 V = V 43 5.00 180.50 45 A = A 44 3.00 183.50 46 A = A 45 3.00 186.50 47 A = A 46 3.00 189.50 48 A = A 47 3.00 192.50 49 P = P 48 6.00 198.50 50 F = F 49 8.00 206.50 51 P = P 50 6.00 212.50 52 G = G 51 5.00 217.50 53 L = L 52 5.00 222.50 54 A = A 53 3.00 225.50 55 E = E 54 5.00 230.50 56 A = A 55 3.00 233.50 57 G = G 56 5.00 238.50 58 V = V 57 5.00 243.50 59 A = A 58 3.00 246.50 60 A = A 59 3.00 249.50 TEST14: SS_GG_FSQ - Forward score quadratic; SS_GG_FSL - Forward score linear; GAP_OPEN penalty = -2.00000 GAP_EXTEND penalty = -0.500000 Matching scores by SIMPLE_SCORE. Verify that FSQ agrees with Myers/Miller calculation. Verify that FSQ and FSL agree. Sequences A and B: 1 A A 2 G A 3 T G 4 A 5 C Matching scores: 0 0.0 0.0 0.0 0.0 1 0.0 0.0 0.0 -1.0 2 0.0 -1.0 -1.0 0.0 3 0.0 -1.0 -1.0 -1.0 4 0.0 0.0 0.0 -1.0 5 0.0 -1.0 -1.0 -1.0 SS_GG_FSQ: 0 1 2 3 SF: 0 0.0 -2.5 -3.0 -3.5 1 -2.5 0.0 -2.5 -3.0 2 -3.0 -2.5 -1.0 -2.5 3 -3.5 -3.0 -3.5 -2.0 4 -4.0 -3.5 -3.0 -4.5 5 -4.5 -4.0 -4.5 -4.0 EF: 0 0.0 -2.5 -3.0 -3.5 1 -4.5 -5.0 -2.5 -3.0 2 -5.0 -5.5 -5.0 -3.5 3 -5.5 -6.0 -5.5 -6.0 4 -6.0 -6.5 -6.0 -5.5 5 -6.5 -7.0 -6.5 -7.0 FF: 0 0.0 -4.5 -5.0 -5.5 1 -2.5 -5.0 -5.5 -6.0 2 -3.0 -2.5 -5.0 -5.5 3 -3.5 -3.0 -3.5 -5.0 4 -4.0 -3.5 -4.0 -4.5 5 -4.5 -4.0 -4.5 -5.0 TF: 0 0 2 2 2 1 4 1 3 2 2 4 4 1 1 3 4 4 5 1 4 4 5 1 5 5 4 4 5 1 SS_GG_FSL: 0 1 2 3 SF: 0 0.0 -2.5 -3.0 -3.5 1 -2.5 0.0 -2.5 -3.0 2 -3.0 -2.5 -1.0 -2.5 3 -3.5 -3.0 -3.5 -2.0 4 -4.0 -3.5 -3.0 -4.5 5 -4.5 -4.0 -4.5 -4.0 EF: 0 0.0 -2.5 -3.0 -3.5 1 -4.5 -5.0 -2.5 -3.0 2 -5.0 -5.5 -5.0 -3.5 3 -5.5 -6.0 -5.5 -6.0 4 -6.0 -6.5 -6.0 -5.5 5 -6.5 -7.0 -6.5 -7.0 FF: 0 0.0 -4.5 -5.0 -5.5 1 -2.5 -5.0 -5.5 -6.0 2 -3.0 -2.5 -5.0 -5.5 3 -3.5 -3.0 -3.5 -5.0 4 -4.0 -3.5 -4.0 -4.5 5 -4.5 -4.0 -4.5 -5.0 TF: 0 0 2 2 2 1 4 1 3 2 2 4 4 1 1 3 4 4 5 1 4 4 5 1 5 5 4 4 5 1 TEST15: SS_GG_BSQ - Backward score quadratic; SS_GG_BSL - Backward score linear; GAP_OPEN penalty = -2.00000 GAP_EXTEND penalty = -0.500000 Matching scores by SIMPLE_SCORE. Verify that BSQ agrees with Myers/Miller calculation. Compare BSQ and BSL results. Sequences A and B: 1 A A 2 G A 3 T G 4 A 5 C SS_GG_BSQ: 0 1 2 3 SF: 0 -4.0 -3.5 -5.0 -4.5 1 -4.5 -4.0 -3.5 -4.0 2 -2.0 -3.5 -4.0 -3.5 3 -3.5 -1.0 -3.5 -3.0 4 -4.0 -3.5 -1.0 -2.5 5 -3.5 -3.0 -2.5 0.0 EF: 0 -6.0 -7.5 -7.0 -6.5 1 -6.5 -6.0 -6.5 -6.0 2 -6.0 -6.5 -6.0 -5.5 3 -3.5 -6.0 -5.5 -5.0 4 -4.0 -3.5 -5.0 -4.5 5 -3.5 -3.0 -2.5 0.0 FF: 0 -5.0 -4.5 -5.0 -4.5 1 -4.5 -4.0 -4.5 -4.0 2 -6.0 -3.5 -4.0 -3.5 3 -6.5 -6.0 -3.5 -3.0 4 -6.0 -5.5 -5.0 -2.5 5 -5.5 -5.0 -4.5 0.0 TF: 0 1 1 5 4 1 5 4 1 4 2 1 4 5 4 3 3 1 5 4 4 3 3 1 4 5 2 2 2 0 SS_GG_BSL: 0 1 2 3 SF: 0 -4.0 -3.5 -5.0 -4.5 1 -4.5 -4.0 -3.5 -4.0 2 -2.0 -3.5 -4.0 -3.5 3 -3.5 -1.0 -3.5 -3.0 4 -4.0 -3.5 -1.0 -2.5 5 -3.5 -3.0 -2.5 0.0 EF: 0 -6.0 -7.5 -7.0 -6.5 1 -6.5 -6.0 -6.5 -6.0 2 -6.0 -6.5 -6.0 -5.5 3 -3.5 -6.0 -5.5 -5.0 4 -4.0 -3.5 -5.0 -4.5 5 -3.5 -3.0 -2.5 0.0 FF: 0 -5.0 -4.5 -5.0 -4.5 1 -4.5 -4.0 -4.5 -4.0 2 -6.0 -3.5 -4.0 -3.5 3 -6.5 -6.0 -3.5 -3.0 4 -6.0 -5.5 -5.0 -2.5 5 -5.5 -5.0 -4.5 0.0 TF: 0 1 1 5 4 1 5 4 1 4 2 1 4 5 4 3 3 1 5 4 4 3 3 1 4 5 2 2 2 0 TEST16: SS_GG_FSQ - forward score quadratic SS_GG_FPQ - forward path quadratic GAP_OPEN penalty = -2.00000 GAP_EXTEND penalty = -0.500000 Matching scores by SIMPLE_SCORE. Compare results with Myers/Miller calculation Sequences A and B: 1 A A 2 G A 3 T G 4 A 5 C Matching path: 1 0 0 2 1 1 3 2 -1 4 3 -1 5 4 2 6 5 3 Sequence/sequence matching, Affine gap penalty: # A B # Increm Score 0 0 0.00 1 A = A 1 0.00 0.00 2 G | -2.50 -2.50 3 T | -0.50 -3.00 4 A = A 2 0.00 -3.00 5 C - G 3 -1.00 -4.00 TEST17: SS_GG_RPL - forward path linear GAP_OPEN penalty = -2.00000 GAP_EXTEND penalty = -0.500000 Matching scores by SIMPLE_SCORE. Compare with FPQ results. Compare results with Myers/Miller calculation Sequences A and B: 1 A A 2 G A 3 T G 4 A 5 C Matching path: 1 0 0 2 1 1 3 2 -1 4 3 -1 5 4 2 6 5 3 Sequence/sequence matching, Affine gap penalty: # A B # Increm Score 0 0 0.00 1 A = A 1 0.00 0.00 2 G | -2.50 -2.50 3 T | -0.50 -3.00 4 A = A 2 0.00 -3.00 5 C - G 3 -1.00 -4.00 TEST18: SS_GG_RPL - forward path linear GAP_OPEN penalty = -2.00000 GAP_EXTEND penalty = -0.500000 Matching scores by SIMPLE_SCORE. We simply switched the two sequences. Compare with unswitched results. Sequences A and B: 1 A A 2 A G 3 G T 4 A 5 C Matching path: 1 0 0 2 1 1 3 -1 2 4 -1 3 5 2 4 6 3 5 Sequence/sequence matching, Affine gap penalty: # A B # Increm Score 0 0 0.00 1 A = A 1 0.00 0.00 | G 2 -2.50 -2.50 | T 3 -0.50 -3.00 2 A = A 4 0.00 -3.00 3 G - C 5 -1.00 -4.00 TEST19: SS_GG_BSQ - backward score quadratic SS_GG_BPQ - backward path quadratic GAP_OPEN penalty = -2.00000 GAP_EXTEND penalty = -0.500000 Matching scores by SIMPLE_SCORE. Compare path and score with FSQ/FPQ computation. Sequences A and B: 1 A A 2 G A 3 T G 4 A 5 C SS_GG_BSQ reports optimal matching score is -4.00000 Matching path: 1 0 0 2 1 1 3 2 -1 4 3 -1 5 4 2 6 5 3 Sequence/sequence matching, Affine gap penalty: # A B # Increm Score 0 0 0.00 1 A = A 1 0.00 0.00 2 G | -2.50 -2.50 3 T | -0.50 -3.00 4 A = A 2 0.00 -3.00 5 C - G 3 -1.00 -4.00 TEST20: SS_GG_FSQ - Forward score quadratic; SS_GG_FPQ - Forward path quadratic; SS_GG_RPL - Recursive path linear; GAP_OPEN penalty = -2.00000 GAP_EXTEND penalty = -0.500000 Matching scores by PAM120_SCORE. Compare the score computed by FSQ/FPQ with the score associated with the path determined by RPL. If the scores don't match, the paths differ, and presumably, the RPL algorithm has failed. The test is carried out by comparing a fixed sequence with many "mutated" variations. Sample test number 10 Sequence/sequence matching, Affine gap penalty: # A B # Increm Score 0 0 0.00 | W 1 -2.50 -2.50 1 G = G 2 5.00 2.50 2 C = C 3 9.00 11.50 3 T = T 4 4.00 15.50 4 G - P 5 -2.00 13.50 5 A = A 6 3.00 16.50 6 T = T 7 4.00 20.50 7 A = A 8 3.00 23.50 8 T = T 9 4.00 27.50 9 A = A 10 3.00 30.50 | D 11 -2.50 28.00 10 G = G 12 5.00 33.00 11 C = C 13 9.00 42.00 12 T = T 14 4.00 46.00 Sample test number 20 Sequence/sequence matching, Affine gap penalty: # A B # Increm Score 0 0 0.00 1 G | -2.50 -2.50 2 C | -0.50 -3.00 3 T = T 1 4.00 1.00 | X 2 -2.50 -1.50 4 G = G 3 5.00 3.50 5 A = A 4 3.00 6.50 6 T = T 5 4.00 10.50 7 A = A 6 3.00 13.50 8 T = T 7 4.00 17.50 9 A | -2.50 15.00 10 G | -0.50 14.50 11 C = C 8 9.00 23.50 12 T = T 9 4.00 27.50 Sample test number 30 Sequence/sequence matching, Affine gap penalty: # A B # Increm Score 0 0 0.00 1 G = G 1 5.00 5.00 2 C = C 2 9.00 14.00 3 T = T 3 4.00 18.00 4 G = G 4 5.00 23.00 5 A = A 5 3.00 26.00 6 T = T 6 4.00 30.00 7 A - M 7 -2.00 28.00 8 T = T 8 4.00 32.00 9 A = A 9 3.00 35.00 10 G = G 10 5.00 40.00 11 C = C 11 9.00 49.00 12 T = T 12 4.00 53.00 Sample test number 40 Sequence/sequence matching, Affine gap penalty: # A B # Increm Score 0 0 0.00 1 G = G 1 5.00 5.00 2 C = C 2 9.00 14.00 3 T = T 3 4.00 18.00 4 G = G 4 5.00 23.00 5 A | -2.50 20.50 6 T | -0.50 20.00 7 A | -0.50 19.50 8 T | -0.50 19.00 9 A | -0.50 18.50 10 G | -0.50 18.00 11 C | -0.50 17.50 12 T | -0.50 17.00 Sample test number 50 Sequence/sequence matching, Affine gap penalty: # A B # Increm Score 0 0 0.00 1 G = G 1 5.00 5.00 2 C = C 2 9.00 14.00 3 T = T 3 4.00 18.00 | Z 4 -2.50 15.50 | S 5 -0.50 15.00 | V 6 -0.50 14.50 4 G = G 7 5.00 19.50 5 A = A 8 3.00 22.50 6 T = T 9 4.00 26.50 7 A = A 10 3.00 29.50 | R 11 -2.50 27.00 | W 12 -0.50 26.50 8 T = T 13 4.00 30.50 9 A = A 14 3.00 33.50 10 G = G 15 5.00 38.50 11 C | -2.50 36.00 12 T | -0.50 35.50 Sample test number 60 Sequence/sequence matching, Affine gap penalty: # A B # Increm Score 0 0 0.00 1 G | -2.50 -2.50 2 C | -0.50 -3.00 3 T | -0.50 -3.50 4 G | -0.50 -4.00 5 A = A 1 3.00 -1.00 6 T = T 2 4.00 3.00 7 A = A 3 3.00 6.00 8 T = T 4 4.00 10.00 9 A = A 5 3.00 13.00 10 G = G 6 5.00 18.00 11 C = C 7 9.00 27.00 12 T = T 8 4.00 31.00 Sample test number 70 Sequence/sequence matching, Affine gap penalty: # A B # Increm Score 0 0 0.00 1 G = G 1 5.00 5.00 2 C = C 2 9.00 14.00 3 T = T 3 4.00 18.00 | P 4 -2.50 15.50 4 G = G 5 5.00 20.50 5 A - G 6 1.00 21.50 6 T - A 7 1.00 22.50 7 A = A 8 3.00 25.50 8 T = T 9 4.00 29.50 9 A = A 10 3.00 32.50 10 G = G 11 5.00 37.50 11 C = C 12 9.00 46.50 12 T = T 13 4.00 50.50 | Y 14 -2.50 48.00 Sample test number 80 Sequence/sequence matching, Affine gap penalty: # A B # Increm Score 0 0 0.00 1 G | -2.50 -2.50 2 C | -0.50 -3.00 3 T | -0.50 -3.50 4 G = G 1 5.00 1.50 5 A = A 2 3.00 4.50 6 T = T 3 4.00 8.50 7 A = A 4 3.00 11.50 8 T = T 5 4.00 15.50 9 A = A 6 3.00 18.50 10 G = G 7 5.00 23.50 11 C | -2.50 21.00 | B 8 -2.50 18.50 | W 9 -0.50 18.00 12 T = T 10 4.00 22.00 Sample test number 90 Sequence/sequence matching, Affine gap penalty: # A B # Increm Score 0 0 0.00 1 G = G 1 5.00 5.00 2 C = C 2 9.00 14.00 | A 3 -2.50 11.50 | F 4 -0.50 11.00 | N 5 -0.50 10.50 3 T = T 6 4.00 14.50 4 G = G 7 5.00 19.50 5 A = A 8 3.00 22.50 6 T = T 9 4.00 26.50 7 A = A 10 3.00 29.50 8 T = T 11 4.00 33.50 9 A = A 12 3.00 36.50 10 G = G 13 5.00 41.50 11 C = C 14 9.00 50.50 12 T = T 15 4.00 54.50 Sample test number 100 Sequence/sequence matching, Affine gap penalty: # A B # Increm Score 0 0 0.00 1 G = G 1 5.00 5.00 2 C = C 2 9.00 14.00 3 T = T 3 4.00 18.00 4 G = G 4 5.00 23.00 5 A = A 5 3.00 26.00 6 T = T 6 4.00 30.00 7 A = A 7 3.00 33.00 8 T = T 8 4.00 37.00 9 A = A 9 3.00 40.00 10 G | -2.50 37.50 11 C | -0.50 37.00 12 T | -0.50 36.50 Number of score agreements = 100 Number of score disagreements = 0 SS_GG_ALIGN_TEST Normal end of execution. March 29 2023 4:41:55.540 PM