Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 942,274 A→G 7.0% L1116L (CTT→CTC PP_0806 ← surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947942,2740AG 7.0% 163.8 / 2.0 57L1116L (CTT→CTCPP_0806surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base A (29/24);  new base G (2/2);  total (31/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

GTGCTGCCGTCGTTCACGGCAGGATTCAGGGTGACGGTGAAGGTGCCATCCGCGGCGACGGTACCGGTCGCGAGGATGTTGCCGGCCGCGTCACGCACCTGCACGGTAGTGCCGGCTTCGCCTCGACCACTGAGCTGGCTGCCGGCAAGGCCGATCACCAGATTGGTCGGTGCTTCAGGCGGGGTGGTATCGCCGCCATCGATATCCGGCGCGGTGACCGAGCTGGTAACGGAAGTGTTGCCAGCGGCATCGGTCGCGGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGCGGGT  >  NC_002947/942128‑942422
                                                                                                                                                  |                                                                                                                                                    
gTGCTGCCGTCGTTCACGGCAGGATTCAGGGTGACGGTGAAGGTGCCATCCGCGGCGACGGTACCGGTCGCGAGGATGTTGCCGGCCGCGTCACGCACCTGCACGGTAGTGCCGGCTTCGCCTCGACCACTGAGCTGGCTGCCGGCAAg                                                                                                                                                    <  1:171867/149‑1 (MQ=255)
    tgcCGTCGTTCACGGCAGGATTCAGGGTGACGGTGAAGGTGCCATCCGCGGCGACGGTACCGGTCGCGAGGATGTTGCCGGCCGCGTCACGCACCTGCACGGTAGTGCCGGCTTCGCCTCGACCACTGAGCTGGCTGCCGGCAAGGcc                                                                                                                                                 >  4:296755/1‑148 (MQ=255)
      ccGGCGTTCTCGGCAGGATTCAGGGTGACGGTGAAGGTGCCATCCGCGGCGACGGTACCGGTCGCGAGGATGTTGCCGGCCGCGTCACGCACCTGCACGGTAGTGCCGGCTTCGCCTCGACCTCTGAGCTGGCTGGAGGCAAGGCCg                                                                                                                                                <  4:63324/147‑1 (MQ=17)
       cgtcgtTCACGGCAGGATTCAGGGTGACGGTGAAGGTGCCATCCGCGGCGACGGTACCGGTCGCGAGGATGTTGCCGGCCGCGTCACGCACCTGCACGGTAGTGCCGGCTTCGCCTCGACCACTGAGCTGGCTGCCGGCAAGGCCGAt                                                                                                                                              >  2:260317/1‑148 (MQ=255)
           gtTCACGGCAGGATTCAGGGTGACGGTGAAGGTGCCATCCGCGGCGACGGTACCGGTCGCGAGGATGTTGCCGGCCGCGTCACGCACCTGCACGGTAGTGCCGGCTTCGCCTCGACCACTGAGCTGGCTGCCGGCAAGGCCGATCACCa                                                                                                                                         <  5:228964/149‑1 (MQ=255)
              cACGGCAGGATTCAGGGTGACGGTGAAGGTGCCATCCGCGGCGACGGTACCGGTCGCGAGGATGTTGCCGGCCGCGTCACGCACCTGCACGGTAGTGCCGGCTTCGCCTCGACCACTGAGCTGGCTGCCGGCAAGGCCGATCACCAGa                                                                                                                                       <  7:224676/148‑1 (MQ=255)
                    aGGATTCAGGGTGACGGTGAAGGTGCCATCCGCGGCGACGGTACCGGTCGCGAGGATGTTGCCGGCCGCGTCACGCACCTGCACGGTAGTGCCGGCTTCGCCTCGACCACTGAGCTGGCTGCCGGCAAGGCCGATCACCAGATTGGTCg                                                                                                                                <  7:165187/149‑1 (MQ=255)
                         tCAGGGTGACGGTGAAGGTGCCATCCGCGGCGACGGTACCGGTCGCGAGGATGTTGCCGGCCGCGTCACGCACCTGCACGGTAGTGCCGGCTTCGCCTCGACCACTGAGCTGGCTGCCGGCAAGGCCGATCACCAGATTGGTCGGTGCt                                                                                                                           >  1:349677/1‑149 (MQ=255)
                                    gTGAAGGTGCCATCCGCGGCGACGGTACCGGTCGCGAGGATGTTGCCGGCCGCGTCACGCACCTGCACGGTAGTGCCGGCTTCGCCTCGACCACTGAGCTGGCTGCCGGCAAGGCCGATCACCAGATTGGTCGGTGCTTCAGGCGGggt                                                                                                                >  8:275873/1‑149 (MQ=255)
                                     tGAAGGTGCCATCCGCGGCGACGGTACCGGTCGCGAGGATGTTGCCGGCCGCGTCACGCACCTGCACGGTAGTGCCGGCTTCGCCTCGACCACTGAGCTGGCTGCCGGCAAGGCCGATCACCAGATTGGTCGGt                                                                                                                              >  6:22903/1‑134 (MQ=255)
                                     tGAAGGTGCCATCCGCGGCGACGGTACCGGTCGCGAGGATGTTGCCGGCCGCGTCACGCACCTGCACGGTAGTGCCGGCTTCGCCTCGACCACTGAGCTGGCTGCCGGCAAGGCCGATCACCAGATTGGTCGGt                                                                                                                              <  5:22903/134‑1 (MQ=255)
                                      gAAGGTGCCATCCGCGGCGACGGTACCGGTCGCGAGGATGTTGCCGGCCGCGTCACGCACCTGCACGGTAGTGCCGGCTTCGCCTCGACCACTGAGCTGGCTGCCGGCAAGGCCGATCACCAGATTGGTCGGTGCTTCAGGCGGggtgg                                                                                                              >  7:142628/1‑149 (MQ=255)
                                         ggTGCCATCCGCGGCGACGGTACCGGTCGCGAGGATGTTGCCGGCCGCGTCACGCACCTGCACGGTAGTGCCGGCTTCGCCTCGCCCACTGAGCTGGCTGCCGGCAAGGCCGATCACCAGATTGGTCGGTGCTTCAGGCGGGGTGGTAt                                                                                                           >  2:130085/1‑149 (MQ=255)
                                         ggTGCCATCCGCGGCGACGGTACCGGTCGCGAGGATGTTGCCGGCCGCGTCACGCACCTGCACGGTAGTGCCGGCTTCGCCTCGACCACTGAGCTGGCTGCCGGCAAGGCCGATCACCAGATTGGTCGGTGCTTCAGGCGGGGTGGTAt                                                                                                           >  8:63779/1‑149 (MQ=255)
                                          gTGCCATCCGCGGCGACGGTACCGGTCGCGAGGATGTTGCCGGCCGCGTCACGCACCTGCACGGTAGTGCCGGCTTCGCCTCGACCACTGAGCTGGCTGCCGGCAAGGCCGATCACCAGATTGGTCGGTGCTTCAGGCgggg                                                                                                                 <  3:113840/142‑1 (MQ=255)
                                          gTGCCATCCGCGGCGACGGTACCGGTCGCGAGGATGTTGCCGGCCGCGTCACGCACCTGCACGGTAGTGCCGGCTTCGCCTCGACCACTGAGCTGGCTGCCGGCAAGGCCGATCACCAGATTGGTCGGTGCTTCAGGCgggg                                                                                                                 >  4:113840/1‑142 (MQ=255)
                                                tCCGCGGCGACGGTACCGGTCGCGAGGATGTTGCCGGCCGCGTCACGCACCTGCACGGTAGTGCCGGCTTCGCCTCGACCACTGAGCTGGCTGCCGGCAAGGCCGATCACCAGATTGGTCg                                                                                                                                >  8:165038/1‑121 (MQ=255)
                                                tCCGCGGCGACGGTACCGGTCGCGAGGATGTTGCCGGCCGCGTCACGCACCTGCACGGTAGTGCCGGCTTCGCCTCGACCACTGAGCTGGCTGCCGGCAAGGCCGATCACCAGATTGGTCg                                                                                                                                <  7:165038/121‑1 (MQ=255)
                                                   gcggcgACGGTACCGGTCGCGAGGATGTTGCCGGCCGCGTCACGCACCTGCACGGTAGTGCCGGCTTCGCCTCGACCACTGAGCTGGCTGCCGGCAAGGCCGATCACCAGATTGGTCGGTGCTTCAGGCGGGGTGGTATCGCCGCCATc                                                                                                 <  1:260317/149‑1 (MQ=255)
                                                                cGGTCGCGAGGATGTTGCCGGCCGCGTCACGCACCTGCACGGTAGTGCCGGCTTCGCCTAGACCACAGAGCTGGCTGCCGGCAATGCCGATCACCAGATTGGTCGGTGc                                                                                                                            >  4:216766/1‑109 (MQ=255)
                                                                cGGTCGAGAGGATGTTGCCGGCCGCGTCACGCACCTGCACGGTAGTGCCGGCTTCGCCTCGACCACTGAGCTGGCTGCCGGCAAGGCCGATCACCAGATTGGTCGGTGc                                                                                                                            <  3:216766/109‑1 (MQ=255)
                                                                  gTCGCGAGGATGTTGCCGGCCGCGTCACGCACCTGCACGGTAGTGCCGGCTTCGCCTCGACCACTGAGCTGGCTGCCGGCAAGGCCGATCACCAGATTGGTCGGTGCTTCAGGCGGGGTGGTATCGCCGCCATCGATATCCGGCGCGGt                                                                                  >  5:119620/1‑149 (MQ=255)
                                                                   tCGCGAGGATGTTGCCGGCCGCGTCACGCACCTGCACGGTAGTGCCGGCTTCGCCTCGACCACTGAGCTGGCTGCCGGCAAGGCCGATCACCAGATTGGTCGGTGCTTCAGGCGGGGTGGTATCGCCGCCATCGATATCCGGCGCGGt                                                                                  >  7:222741/1‑148 (MQ=255)
                                                                     gcgAGGATGTTGCCGGCCGCGTCACGCACCTGCACGGTAGTGCCGGCTTCGCCTCGACCACTGAGCTGGCTGCCGGCAAGGCCGATCACCAGATTGGTCGGTGCTTCAGGCGGggtggt                                                                                                             <  1:167968/119‑1 (MQ=255)
                                                                     gcgAGGATGTTGCCGGCCGCGTCACGCACCTGCACGGTAGTGCCGGCTTCGCCTCGACCACTGAGCTGGCTGCCGGCAAGGCCGATCACCAGATTGGTCGGTGCTTCAGGCGGggtggt                                                                                                             >  2:167968/1‑119 (MQ=255)
                                                                      cgAGGATGTTGCCGGCCGCGTCACGCACCTGCACGGTAGTGCCGGCTTCGCCTCGACCACTGAGCTGGCTGCCGGCAAGGCCGATCACCAGATTGGTCGGTGCTTCAGGCGGGGTGGTATCGCCGCCATCGATATCCGGCGCGGTGAcc                                                                              <  7:63779/149‑1 (MQ=255)
                                                                                   ggccgCGTCACGCACCTGCACGGTAGTGCCGGCTTCGCCTCGACCCCTGAGCTGGCTGCCGGCAAGGCCGCTCCCCAGAGTGGTCGGTGCTTCCGGCGGGGGGGTATCGCCGCCATCGATATCCGGCGCGGTGCCCGAGCTGGTAACg                                                                  >  3:216133/1‑148 (MQ=18)
                                                                                       gcgTCACGCACCTGCACGGTAGTGCCGGCTTCGCCTCGACCACTGAGCTGGCTGCCGGCAAGGCCGATCACCAGATTGGTCGGTGCTTCAGGCGGGGTGGTATCGCCGCCATCGATATCCGGCGCGGTGAcc                                                                              <  5:261475/132‑1 (MQ=255)
                                                                                       gcgTCACGCACCTGCACGGTAGTGCCGGCTTCGCCTCGACCACTGAGCTGGCTGCCGGCAAGGCCGATCACCAGATTGGTCGGTGCTTCAGGCGGGGTGGTATCGCCGCCATCGATATCCGGCGCGGTGAcc                                                                              >  6:261475/1‑132 (MQ=255)
                                                                                       gcgTCACGCACCTGCACGGTAGTGCCGGCTTCGCCTCGACCACTGAGCTGGCTGCCGGCAAGGCCGATCACCAGATTGGTCGGTGCTTCAGGCGGGGTGGTATCGCCGCCATCGATATCCGGCGCGGTGACCGAGCTGGTAACGGAAgt                                                             <  6:119620/149‑1 (MQ=255)
                                                                                       gcgTCACGCACCTGCACGGTAGTGCCGGCTTCGCCTCGACCACTGAGCTGGCTGCCGGCAAGGCCGATCACCAGATTGGTCGGTGCTTCAGGCGGGGTGGTATCGCCGCCATCGATATCCGGCGCGGTGACCGAGCTGGCAAcacttcc                                                             >  8:194217/1‑143 (MQ=255)
                                                                                         gTCACGCACCTGCACGGTAGTGCCGGCTTCGCCTCGACCACTGAGCTGGCTGCCGGCAAGGCCGATCACCAGATTGGTCgg                                                                                                                               <  5:196846/81‑1 (MQ=255)
                                                                                         gTCACGCACCTGCACGGTAGTGCCGGCTTCGCCTCGACCACTGAGCTGGCTGCCGGCAAGGCCGATCACCAGATTGGTCgg                                                                                                                               >  6:196846/1‑81 (MQ=255)
                                                                                            aCGCACCTGCACGGTAGTGCCGGCTTCGCCTCGACCACTGAGCTGGCTGCCGGCAAGGCCGATCACCAGATTGGTCGGTGCTTCAGGCGGGGTGGTATCGCCGCCATCGATATCCGGCGCGGTGACCGAGCTGGTAACGGAAGTGTTGc                                                        <  7:275873/149‑1 (MQ=255)
                                                                                                aCCTGCACGGTAGTGCCGGCTTCGCCTCGACCACTGAGCTGGCTGCCGGCAAGGCCGATCACCAGATTGGTCGGTGCTTCAGGCGGGGTGGTATCGCCGCCATCGATATCCGGCGCGGTGACCGAGCTGGTAACGGAAGTGTTGCCAGc                                                    >  7:338243/1‑149 (MQ=255)
                                                                                                  cTGCACGGTAGTGCCGGCTTCGCCTCGACCACTGAGCTGGCTGCCGGCGAGGCCGATCACCAGATTGGTCGGTGCTTCAGGCGGGGTGGTATCGCCGCCATCGATATCCGGCGCGGTGACCGAGCTGGTAACGGAAGTGTTGCCAGCgg                                                  >  2:335779/1‑149 (MQ=255)
                                                                                                  cTGCACGGTAGTGCCGGCTTCGCCTCGACCACTGAGCTGGCTGCCGGCGAGGCCGATCACCAGATTGGTCGGTGCTTCAGGCGGGGTGGTATCGCCGCCATCGATATCCGGCGCGGTGACCGAGCTGGTAACGGAAGTGTTGCCAGCgg                                                  >  6:24449/1‑149 (MQ=255)
                                                                                                    gCACGGTAGTGCCGGCTTCGCCTCGACCACTGAGCTGGCTGCCGGCAAGGCCGATCACCAGATTGGTCGGTGCTTCAGGCGGGGTGGTATCGCCGCCATCGATATCCGGCGCGGTGACCGAGCTGGTAACGGAAGTGTTGCCAGCGGCa                                                >  8:217247/1‑149 (MQ=255)
                                                                                                      aCGGTAGTGCCGGCTTCGCCTCGACCACTGAGCTGGCTGCCGGCAAGGCCGATCACCAGATTGGTCGGTGCTTCAGGCGGGGTGGTATCGCCGCCATCGATATCCGGCGCGGTGACCGAGCTGGTAACGGAAGTGTTGCCAGCGGCATc                                              <  7:217247/149‑1 (MQ=255)
                                                                                                                      cGCCTCGACCACTGAGCTGGCTGCCGGCAAGGCCGATCACCAGATTGGTCGGTGCTTCAGGCGGGGTGGTATCGCCGCCATCGATATCCGGCGCGGTGACCGAGCTGGTAACGGAAGTGTTGCCAGCGGCATCGGTCGCGGTGGCTTg                               <  2:349677/148‑1 (MQ=255)
                                                                                                                         cTCGACCACTGAGCTGGCTGCCGGCAAGGCCGATCACCAGATTGGTCGGTGCTTCAGGCGGGGTGGTATCGCCGCCATCGATATCCGGCGCGGTGACCGAGCTGGTAACGGAAGTGTTGCCAGCGGCATCGGTCGCGGTGGCTTgc                              <  4:78565/146‑1 (MQ=255)
                                                                                                                         cTCGACCACTGAGCTGGCTGCCGGCAAGGCCGATCACCAGATTGGTCGGTGCTTCAGGCGGGGTGGTATCGCCGCCATCGATATCCGGCGCGGTGACCGAGCTGGTAACGGAAGTGTTGCCAGCGGCATCGGTCGCGGTGGCTTgc                              >  3:78565/1‑146 (MQ=255)
                                                                                                                          tCGACCACTGAGCTGGCTGCCGGCAAGGCCGATCACCAGATTGGTCGGTGCTTCAGGCGGGGTGGTATCGCCGCCATCGATATCCGGCGCGGTGACCGAGCTGGTAACGGAAGTGTTGCCAGCGGCATCGGTCGCGGTGGCTTgcagca                          >  3:147019/1‑149 (MQ=255)
                                                                                                                               cccTGCGCGGGCGGCCGGCAAGGCCGATCACCCGATTGGTCGGTGCGTCAGGCGGGGTGGTATCGCCGCCATCGATATCCGGCGCGGTGACCGAGCTGGTAACGGAAGTGTTGCCAGCGGCATCGGTCGCGGTGGCTTGCAGCACTTCg                     <  4:216133/147‑1 (MQ=37)
                                                                                                                               cACTGAGCTGGCTGCCGGCAAGGCCGATCACCAGATTGGTCGGTGCTTCAGGCGGGGTGGTATCGCCGCCATCGATATCCGGCGCGGTGACCGAGCTGGt                                                                      >  6:185034/1‑100 (MQ=255)
                                                                                                                               cACTGAGCTGGCTGCCGGCAAGGCCGATCACCAGATTGGTCGGTGCTTCAGGCGGGGTGGTATCGCCGCCATCGATATCCGGCGCGGTGACCGAGCTGGt                                                                      <  5:185034/100‑1 (MQ=255)
                                                                                                                                    aGCTGGCTGCCGGCAAGGCCGATCACCAGATTGGTCGGTGCTTCAGGCGGGGTGGTATCGCCGCCATCGATATCCGGCGCGGTGACCGAGCTGGTAACGGAAGTGTTGCCAGCGGCATCGGTCGCGGTGGCTTGCAGCACTTCGCCGt                 <  3:296755/148‑1 (MQ=255)
                                                                                                                                      ctggctgCCGGCGAGGCCGATCACCAGATTGGTCGGTGCTTCAGGCGGGGTGGTATCGCCGCCATCGATATCCGGCGCGGTGACCGAGCTGGTAACGGAAGTGTTGCCAGCGGCATCGGTCGCGGTGGCTTGCAGCACTTCGCCGTCaa              <  1:335779/149‑1 (MQ=255)
                                                                                                                                      ctggctgCCGGCGAGGCCGATCACCAGATTGGTCGGTGCTTCAGGCGGGGTGGTATCGCCGCCATCGATATCCGGCGCGGTGACCGAGCTGGTAACGGAAGTGTTGCCAGCGGCATCGGTCGCGGTGGCTTGCAGCACTTCGCCGTCaa              <  5:24449/149‑1 (MQ=255)
                                                                                                                                      ctggctgCCGGCAAGGCCGATCACCAGATTGGTCGGTGCTTCAGGCGGGGTGGTATCGCCGCCATCGATATCCGGCGCGGTGACCGAGCTGGTAACGGAAGTGTTGCCAGCGGCATCGGTCGCGGTGGCTTGCAGCACTTCGCCGTCaa              >  3:33486/1‑149 (MQ=255)
                                                                                                                                        ggctgCCGGCAAGGCCGATCACCAGATTGGACGGTGCTTCAGGCGGGGTGGTATCGCCGCCATCGATATCCGGCGCGGTGACCGAGCTGGTAACGGAAGTGTTGCCAGCGGCATCGGTCGCGGTGGCTTGCAGCACTTCGCCGTCAACt            >  1:31584/1‑149 (MQ=255)
                                                                                                                                          ctgCCGGCAAGGCCGATCACCAGATTGGTCGGTGCTTCAGGCGGGGTGGTATCGCCGCCATCGATATCCGGCGCGGTGACCGAGCTGGTAACGGAAGTGTTGCCAGCGGCATCGGTCGCGGTGGCTTGCAGCACTTCGCCGTCAACTtg          >  5:82529/1‑149 (MQ=255)
                                                                                                                                              cGGCAAGGCCGATCACCAGATTGGTCGGTGCTTCAGGCGGGGTGGTATCGCCGCCATCGATATCCGGCGCGGTGACCGAGCTGGTAACGGAAGTGTTGCCAGCGGCATCGGTCGCGGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGc      >  2:104080/1‑149 (MQ=255)
                                                                                                                                               ggCAAGGCCGATCACCAGATTGGTCGGTGCTTCAGGCGGGGTGGTATCGCCGCCATCGATATCCGGCGCGGTGACCGAGCTGGTAACGGAAGTGTTGCCAGCGGCATCGGTCGCGGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGCg     >  1:44730/1‑149 (MQ=255)
                                                                                                                                                gCAAGGCCGATCACCAGATTGGTCGGTGCTTCAGGCGGGGTGGTATCGCCGCCATCGATATCCGGCGCGGTGACCGAGCTGGTAACGGAAGTGTTGCCAGCGGCATCGGTCGCGGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGCgg    <  7:128238/149‑1 (MQ=255)
                                                                                                                                                 cAAGGCCGATCACCAGATTGGTCGGTGCTTCAGGCGGGGTTGTATCGCCGCCATCGATATCCGGCGCGGTGACCGAGCTGGTAACGGAAGTGTTGCCAGCGGCATCGGTCGCGGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGCggg   <  8:222741/149‑1 (MQ=255)
                                                                                                                                                 cAAGGCCGATCACCAGATTGGTCGGTGCTTCAGGCGGGGTGGTATCGCCGCCATCGATATCCGGCGCGGTGACCGAGCTGGTAACGGAAGTGTTGCCAGCGGCATCGGTCGCGGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGCggg   >  1:241366/1‑149 (MQ=255)
                                                                                                                                                  aaGGCCGATCACCAGATTGGTCGGTGCTTCAGGCGGGGTGGTATCGCCGCCATCGATATCCGGCGCGGTGACCGAGCTGGTAACGGAAGTGTTGCCAGCGGCATCGGTCGCGGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGCGggt  >  2:318342/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GTGCTGCCGTCGTTCACGGCAGGATTCAGGGTGACGGTGAAGGTGCCATCCGCGGCGACGGTACCGGTCGCGAGGATGTTGCCGGCCGCGTCACGCACCTGCACGGTAGTGCCGGCTTCGCCTCGACCACTGAGCTGGCTGCCGGCAAGGCCGATCACCAGATTGGTCGGTGCTTCAGGCGGGGTGGTATCGCCGCCATCGATATCCGGCGCGGTGACCGAGCTGGTAACGGAAGTGTTGCCAGCGGCATCGGTCGCGGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGCGGGT  >  NC_002947/942128‑942422

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: