Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,326,724 A→G S15P (TCT→CCT)  mreD ← cell wall structural complex MreBCD transmembrane component MreD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,326,7240AG98.8% 1013.3 / 0.4 251S15P (TCT→CCT) mreDcell wall structural complex MreBCD transmembrane component MreD
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base G (137/111);  total (138/113)
Fisher's exact test for biased strand distribution p-value = 5.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.56e-01

Only 100 of 272 total aligned reads displayed.

ACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCAGGCGGTT  >  REL606/3326651‑3326797
                                                                         |                                                                         
acCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGACaggga                                                                         >  1:2125917/1‑76 (MQ=255)
acCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAaggga                                                                         <  2:848508/76‑1 (MQ=255)
acCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAaggga                                                                         <  2:2781974/76‑1 (MQ=255)
acCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAaggga                                                                         <  2:6956901/76‑1 (MQ=255)
acCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAaggga                                                                         <  1:3454891/76‑1 (MQ=255)
 cccAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAagggag                                                                        <  2:3662796/76‑1 (MQ=255)
   cAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGGGAGcc                                                                      >  1:6848734/1‑76 (MQ=255)
   cAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGGGAGcc                                                                      >  1:8032336/1‑76 (MQ=255)
   cAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGGGAGcc                                                                      <  1:8874926/76‑1 (MQ=255)
     gTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGGGAGCCAg                                                                    <  1:3205770/76‑1 (MQ=255)
     gTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGGGAGCCAg                                                                    >  1:8403443/1‑76 (MQ=255)
      tttGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGGGAGCCAGa                                                                   >  2:5509054/1‑76 (MQ=255)
      tttGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGGGAGCCAGa                                                                   <  2:5955816/76‑1 (MQ=255)
      tttGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGGGAGCCAGa                                                                   >  2:1927607/1‑76 (MQ=255)
       ttGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGGGAGCCAGAt                                                                  >  1:463357/1‑76 (MQ=255)
       ttGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGGGAGCCAGAt                                                                  <  1:5885943/76‑1 (MQ=255)
       ttGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGGGAGCCAGAt                                                                  <  2:8222072/76‑1 (MQ=255)
       ttGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGGGAGCCAGAt                                                                  >  2:7098853/1‑76 (MQ=255)
       ttGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGGGAGCCAGAt                                                                  <  2:3184905/76‑1 (MQ=255)
        tGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGGGAGCCAGAtt                                                                 >  1:2155015/1‑76 (MQ=255)
         ggCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGGGAGCCAGATTa                                                                >  2:7845775/1‑76 (MQ=255)
         ggCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGGGAGCCAGATTa                                                                >  2:5092718/1‑76 (MQ=255)
         ggCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGGGAGCCAGATTa                                                                <  1:7034313/76‑1 (MQ=255)
          gCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGGGAGCCAGATTAc                                                               <  1:8531951/76‑1 (MQ=255)
             ggAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGGGAGCCAGATTACCCa                                                            <  1:5586544/76‑1 (MQ=255)
              gAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGGGAGCCAGATTACCCAg                                                           <  2:2998952/76‑1 (MQ=255)
              gAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGGGAGCCAGATTACCCAg                                                           >  1:1920687/1‑76 (MQ=255)
               aaaaCAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGGGAGCCAGATTACCCAGc                                                          >  2:1081428/1‑76 (MQ=255)
                aaaCAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCCATGAGGAAAGGGAGCCAGATTACCCAGCg                                                         >  2:2684030/1‑76 (MQ=255)
                   cAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGGGAGCCAGATTACCCAGCGTcc                                                      <  2:1226132/76‑1 (MQ=255)
                   cAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGGGAGCCAGATTACCCAGCGTcc                                                      >  1:1091085/1‑76 (MQ=255)
                     aTCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGGGAGCCAGATTACCCAGCGTCCCt                                                    <  1:1496730/76‑1 (MQ=255)
                     aTCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGGGAGCCAGATTACCCAGCGTCCCt                                                    >  1:7379414/1‑76 (MQ=255)
                     aTCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGGGAGCCAGATTACCCAGCGTCCCt                                                    <  1:6965582/76‑1 (MQ=255)
                      tCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGGGAGCCAGATTACCCAGCGTCCCTg                                                   <  2:662395/76‑1 (MQ=255)
                       cAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGGGAGCCAGATTACCCAGCGTCCCTgg                                                  >  1:2867018/1‑76 (MQ=255)
                       cAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGGGAGCCAGATTACCCAGCGTCCCTgg                                                  >  1:6671108/1‑76 (MQ=255)
                         ggTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGGGAGCCAGATTACCCAGCGTCCCTGGCt                                                <  2:5016740/76‑1 (MQ=255)
                          gTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGGGAGCCAGATTACCCAGCGTCCCTGGCTa                                               <  2:8895418/76‑1 (MQ=255)
                          gTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGGGAGCCAGATTACCCAGCGTCCCTGGCTa                                               <  2:6765198/76‑1 (MQ=255)
                           ttATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGGGAGCCAGATTACCCAGCGTCCTTGGCTAc                                              >  1:5729497/1‑76 (MQ=255)
                           ttATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGGGAGCCAGATTACCCAGCGTCCCTGGCTAc                                              <  2:7767748/76‑1 (MQ=255)
                           ttATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGGGAGCCAGATTACCCAGCGTCCCTGGCTAc                                              <  1:3047978/76‑1 (MQ=255)
                           ttATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGGGAGCCAGATTACCCAGCGTCCCTGGCTAc                                              >  2:8739795/1‑76 (MQ=255)
                            tATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGGGAGCCAGATTACCCAGCGTCCCTGGCTACg                                             <  2:6229096/76‑1 (MQ=255)
                             aTCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGGGAGCCAGATTACCCAGCGTCCCTGGCTACGa                                            >  2:6475195/1‑76 (MQ=255)
                              tCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGGGAGCCAGATTACCCAGCGTCCCGGGCTACGAt                                           <  2:2118480/76‑1 (MQ=255)
                               ccGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGGGAGCCAGATTACCCAGCGTCCCTGGCTACGATa                                          >  1:4523728/1‑76 (MQ=255)
                                cGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGGGAGCCAGATTACCCAGCGTCCCTGGCTACGATAg                                         >  1:6747228/1‑76 (MQ=255)
                                cGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGGGAGCCAGATTACCCAGCGTCCCTGGCTACGATAg                                         >  1:7882632/1‑76 (MQ=255)
                                  gCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGGGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCt                                       >  1:6600647/1‑76 (MQ=255)
                                    cAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGGGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTcc                                     >  2:8040173/1‑75 (MQ=255)
                                     aGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGGGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGc                                    >  1:8350613/1‑76 (MQ=255)
                                        gCATGATTTGCAGCAACAGCGCAATGAGGAAAGGGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCAc                                 >  2:6108610/1‑76 (MQ=255)
                                        gCATGATTTGCAGCAACAGCGCAATGAGGAAAGGGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCAc                                 >  1:7180946/1‑76 (MQ=255)
                                         cATGATTTGCAGCAACAGCGCAATGAGGAAAGGGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACt                                >  2:4826802/1‑76 (MQ=255)
                                          aTGATTTGCAGCAACAGCGCAATGAGGAAAGGGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTTGCCACTa                               >  2:7187443/1‑76 (MQ=255)
                                          aTGATTTGCAGCAACAGCGCAATGAGGAAAGGGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCCCTa                               >  1:1773730/1‑76 (MQ=255)
                                          aTGATTTGCAGCAACAGCGCAATGAGGAAAGGGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTa                               >  2:4793232/1‑76 (MQ=255)
                                            gATTTGCAGCAACAGCGCAATGAGGAAAGGGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTAtt                             >  2:2731566/1‑76 (MQ=255)
                                            gATTTGCAGCAACAGCGCAATGAGGAAAGGGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTAtt                             <  2:3987467/76‑1 (MQ=255)
                                            gATTTGCAGCAACAGCGCAATGAGGAAAGGGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTAtt                             <  1:6827454/76‑1 (MQ=255)
                                             aTTTGCAGCAACAGCGCAATGAGGAAAGGGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTg                            <  2:2313583/76‑1 (MQ=255)
                                             aTTTGCAGCAACAGCGCAATGAGGAAAGGGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTg                            <  2:186/76‑1 (MQ=255)
                                             aTTTGCAGCAACAGCGCAATGAGGAAAGGGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTg                            >  1:6392165/1‑76 (MQ=255)
                                              tttGCAGAAACAGCGCAATGAGGAAAGGGAGCCAGATTACCCAGCGTCCCTGGCTACCATAGCTCGCCACTATTGc                           >  2:7479917/1‑76 (MQ=255)
                                                tGCAGCAACAGCGCAATGAGGAAAGGGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGccc                         >  2:3593974/1‑76 (MQ=255)
                                                tGCAGCAACAGCGCAATGAGGAAAGGGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGccc                         <  2:4362295/76‑1 (MQ=255)
                                                  cggcaACAGCGCAAGGAGGAAAGGGAGCCAGATTCCCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTc                       <  1:599764/74‑1 (MQ=255)
                                                  cagcaACAGCGCAATGAGGAAAGGGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTc                       <  2:3098944/76‑1 (MQ=255)
                                                    gcaACAGCGCAATGAGGAAAGGGAGCCAGATTCCCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTccc                     <  1:3581958/76‑1 (MQ=255)
                                                    gcaACAGCGCAATGAGGAAAGGGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTccc                     >  1:492003/1‑76 (MQ=255)
                                                     caACAGCGCAATGAGGAAAGGGAGCCAGATTCCCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCg                    <  1:1331775/76‑1 (MQ=255)
                                                     caACAGCGCAATGAGGAAAGGGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCg                    <  2:1468276/76‑1 (MQ=255)
                                                      aaCAGCGCAATGAGGAAAGGGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCATTATTGCCCTCCCgg                   >  1:4187216/1‑76 (MQ=255)
                                                      aaCAGCGCAATGAGGAAAGGGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCgg                   >  2:8420241/1‑76 (MQ=255)
                                                      aaCAGCGCAATGAGGAAAGGGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCgg                   <  1:1061663/76‑1 (MQ=255)
                                                       aCAGCGCAATGAGGAAAGGGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGgc                  <  1:5336512/76‑1 (MQ=255)
                                                        cAGCGCAATGAGGAAAGGGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGgcg                 >  1:6564168/1‑76 (MQ=255)
                                                         aGCGCAATGAGGAAAGGGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGgcgc                >  1:2726370/1‑76 (MQ=255)
                                                         aGCGCAATGAGGAAAGGGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGgcgc                <  2:34438/76‑1 (MQ=255)
                                                         aGCGCAATGAGGAAAGGGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCCgcgc                >  1:509027/1‑76 (MQ=255)
                                                         aGCGCAATGAGGAAAGGGAGCCAGATTACCCAGCGTCCATGGCTACGATAGCTCGCCACGATTGCCCTCCCGgcgc                >  1:8933470/1‑76 (MQ=255)
                                                          gcgcAATGAGGAAAGGGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCa               <  2:5226568/76‑1 (MQ=255)
                                                          gcgcAATGAGGAAAGGGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCa               >  1:6243546/1‑76 (MQ=255)
                                                           cgcAATGAGGAAAGGGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCAc              >  1:8191041/1‑76 (MQ=255)
                                                            gcAATGAGGAAAGGGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCAcg             >  1:7709948/1‑76 (MQ=255)
                                                             cAATGAGGAAAGGGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCAcgc            <  1:1285676/76‑1 (MQ=255)
                                                              aaTGAGGAAAGGGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCCTcccgcgcccgcg           >  2:1623649/1‑76 (MQ=255)
                                                               aggAGGAAGGGGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACgcgc          <  1:42758/74‑1 (MQ=255)
                                                               aTGAGGAAAGGGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACgcgc          >  1:4803383/1‑76 (MQ=255)
                                                               aTGAGGAAAGGGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCCCCCGGGCAcgccaa         >  1:3802591/1‑73 (MQ=255)
                                                                tGAGGAAAGGGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCa         <  1:5249776/76‑1 (MQ=255)
                                                                tGAGGAAAGGGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCa         <  1:7944484/76‑1 (MQ=255)
                                                                 gAGGAAAGGGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCAg        <  1:873430/76‑1 (MQ=255)
                                                                 gAGGAAAGGGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCAg        <  2:6670475/76‑1 (MQ=255)
                                                                  aGGAAAGGGACCCAGATTCCCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCAgg       <  2:2535304/76‑1 (MQ=255)
                                                                    gAAAGGGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCAGGCg     >  1:755942/1‑76 (MQ=255)
                                                                     aaaGGGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCAGGCgg    >  1:2626411/1‑76 (MQ=255)
                                                                       agagagCCAGCTTCCCCAGCGGCCCTGGCTACAATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCAGGCGGtt  <  1:6916355/76‑1 (MQ=255)
                                                                         |                                                                         
ACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGATTACCCAGCGTCCCTGGCTACGATAGCTCGCCACTATTGCCCTCCCGGCGCACGCGCAGGCGGTT  >  REL606/3326651‑3326797

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