Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 5,908,523 G→A 6.4% A291V (GCG→GTG)  argA ← Amino‑acid acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24405,908,5230GA 6.4% 181.2 / 2.3 62A291V (GCG→GTG) argAAmino‑acid acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (27/31);  new base A (2/2);  total (29/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.29e-01

GCCCTCGCGCTCGACCACGCTGAACTGCTCGATCTCTCGCTCCAGCACTTCGCGCGAGCGGCGCACCAGGATGCCCTGCTCTTCCAGCGGGCTGATCAATTCCAGCAGGCCACCGACGTCCTCGATGGTCGCCTCGCGCACCACTTCGAACTGCTCCTGCGACACCAGCGTACCGCCACCGCCCCGTGTGAACAGCTCGGTGAGCAGCGCGCCGTCCTCGGCATAGCTGACGATATGGCTGCGCGCCACACCGCCTTT  >  KT2440/5908393‑5908650
                                                                                                                                  |                                                                                                                               
gCCCTCGCGCTCGACCACGCTGAACTGCTCGATCTCTCGCTCCATCACTTCGCGCGAGCGGCGCACCAGGATGCCCTGCTCTTCCAGCGGGCTGATCAATTCGAGCAGGCCACCGACGTCCTCGATGGTCGCCTc                                                                                                                             <  5:82405/135‑1 (MQ=255)
  ccTCGCGCTCGACCACGCTGAACTGCTCGATCTCTCGCTCCAGCACTTCGCGCGAGCGGCGCACCAGGATGCCCTGCTCTTCCAGCGGGCTGATCAATTCCAGCCGGCCACCGACGTCCTCGATGGTCGCCTcgc                                                                                                                           >  5:49866/1‑135 (MQ=255)
  ccTCGCGCTCGACCACGCTGAACTGCTCGATCTCTCGCTCCAGCACTTCGCGCGAGCGGCGCACCAGGATGCCCTGCTCTTCCAGCGGGCTGATCAATTCCAGCAGGCCACCGACGTCCTCGATGGTCGCCTcgcg                                                                                                                          >  8:143984/1‑136 (MQ=255)
  ccTCGCGCTCGACCACGCTGAACTGCTCGATCTCTCGCTCCAGCACTTCGCGCGAGCGGCGCACCAGGATGCCCTGCTCTTCCAGCGGGCTGATCAATTCCAGCAGGCCACCGACGTCCTCGATGGTCGCCTcgc                                                                                                                           >  2:35069/1‑135 (MQ=255)
  ccTCGCGCTCGACCACGCTGAACTGCTCGATCTCTCGCTCCAGCACTTCGCGCGAGCGGCGCACCAGGATGCCCTGCTCTTCCAGCGGGCTGATCAATTCCAGCAGGCCACCGACGTCCTCGATGGTCGCCTcgc                                                                                                                           >  3:50081/1‑135 (MQ=255)
  ccTCGCGCTCGACCACGCTGAACTGCTCGATCTCTCGCTCCAGCACTTCGCGCGAGCGGCGCACCAGGATGCCCTGCTCTTCCAGCGGGCTGATCAATTCCAGCAGGCCACCGACGTCCTCGATGGTCGCCTcgc                                                                                                                           <  3:101029/135‑1 (MQ=255)
  ccTCGCGCTCGACCACGCTGAACTGCTCGATCTCTCGCTCCAGCACTTCGCGCGAGCGGCGCACCAGGATGCCCTGCTCTTCCAGCGGGCTGATCAATTCCAGCAGGCCACCGACGTCCTCGATGGTCACCTcgcg                                                                                                                          >  7:16042/1‑136 (MQ=255)
  ccTCGCGCTCGACCACGCTGAACTGCTCGATCTCTCGCTCCAGCACTTCGCGCGAGCGGCGCACCAGGATGCCCTGCTCTTCCAGCGGGCTGATCAATTCCAGCAGGCCACCGACGTCCTCGATGGTCACCTcgcg                                                                                                                          >  7:11737/1‑136 (MQ=255)
  ccTCGCGCTCGACCACGCTGAACTGCTCGATCTCTCGCTCCAGCACTTCGCGCGAGCGGCGCACCAGGATGCCCTGCTCTTCCAGCGGGCTGATCAATTCCAGCAGGCCACCGACGTCCTCGATGGGCGCCTcgc                                                                                                                           <  5:9097/135‑1 (MQ=255)
      gcgcTCGACCACGCTGAACTGCTCGATCTCTCGCTCCAGCACTTCGCGCGAGCGGCGCACCAGGATGCCCTGCTCTTCCAGCGGGCTGATCAATTCCAGCAGGCCACCGACGTCCTCGATGGTCGCCTCGCGcacc                                                                                                                      <  7:241375/136‑1 (MQ=255)
       cgcTCGACCACGCTGAACTGCTCGATCTCTCGCTCCAGCACTTCGCGCGAGCGGCGCACCAGGATGCCCTGCTCTTCCAGCGGGCTGATCAATTCCAGCAGGCCACCGACGTCCTCGATGGTCGCCTCGCGcacc                                                                                                                      <  6:141088/135‑1 (MQ=255)
        gcTCGACCACGCTGAACTGCTCGATCTCTCGCTCCAGCACTTCGCGCGAGCGGCGCACCAGGATGCCCTGCTCTTCCAGCGGGCTGATCAATTCCAGCAGGCCACCGACGTCCTCGATGGTCGCCTCGCGcacca                                                                                                                     >  1:266296/1‑135 (MQ=255)
               cACGCTGAACTGCTCGATCTCTCGCTCCAGCACTTCGCGCGAGCGGCGCACCAGGATGCCCTGCTCTTCCAGCGGGCTGATCAATTCCAGCAGGCCACCGACGTCCTCGATGGTCGCCTCGCGCACCACTTCGaa                                                                                                              <  6:179876/135‑1 (MQ=255)
                        cTGCTCGATCTCTCGCTCCAGCACTTCGCGCGAGCGGCGCACCAGGATGCCCTGCTCTTCCAGCGGGCTGATCAATTCCAGCAGGCCACCGACGTCCTCGATGGTCGCCTCGCGCACCACTTCGAACTGCTCCTg                                                                                                     <  6:244396/135‑1 (MQ=255)
                        cTGCTCGATCTCTCGCTCCAGCACATCGCGCGAGCGGCGCACCAGGATGCCCTGCTCTTCCAGCGGGCTGATCAAATCCAGCAGGCCACCGACGTCCTCGATGGTCGCCTCGCGCACCACTTCGAACTGCTCCTg                                                                                                     <  1:138407/135‑1 (MQ=255)
                         tGCTCGATCTCTCGCTCCAGCACTTCGCGCGAGCGGCGCACCAGGATGCCCTGCTCTTCCAGCGGGCTGATCAATTCCAGCAGGCCACCGACGTCCTCGATGGTCGCCTCGCGCACCACTTc                                                                                                                 >  6:126751/1‑122 (MQ=255)
                         tGCTCGATCTCTCGCTCCAGCACTTCGCGCGAGCGGCGCACCAGGATGCCCTGCTCTTCCAGCGGGCTGATCAATTCCAGCAGGCCACCGACGTCCTCGATGGTCGCCTCGCGCACCACTTc                                                                                                                 <  5:126751/122‑1 (MQ=255)
                         tGCTCGATCTCTCGCTCCAGCACTTCGCGCGAGCGGCGCACCAGGATGCCCTGCTCTTCCAGCGGGCTGATCAATTCCAGCAGGCCACCGACGTCCTCGATGGTCGCCTCGCGCACCACTTCGAACTGCTCCTg                                                                                                     >  6:118663/1‑134 (MQ=255)
                              gATCTCTCGCTCCAGCACTTCGCGCGAGCGGCGCACCAGGATGCCCTGCTCTTCCAGCGGGCTGATCAATTCCAGCAGGCCACCGACGTCCTCGATGGTCGCCTCGCGCACCACTTCGAACTGCTCCTGCGacac                                                                                               >  3:111521/1‑135 (MQ=255)
                              gATCTCTCGCTCCAGCACTTCGCGCGAGCGGCGCACCAGGATGCCCTGCTCTTCCAGCGGGCTGATCAATTCCAGCAGGCCACCGACGTCCTCGATGGTCGCCTCGCGCACCACTTCGAACTGCTCCTGCGacac                                                                                               >  6:142020/1‑135 (MQ=255)
                              gATCTCTCGCTCCAGCACTTCGCGCGAGCGGCGCACCAGGATGCCCTGCTCTTCCAGCGGGCTGATCAATTCCAGCAGGCCACCGACGTCCTCGATGGTCGCCTCGCGCACCACTTCGAACTGCTCCTGCGacac                                                                                               >  6:150191/1‑135 (MQ=255)
                              gATCTCTCGCTCCAGCACTTCGCGCGAGCGGCGCACCAGGATGCCCTGCTCTTCCAGCGGGCTGATCAATTCCAGCAGGCCACCGACGTCCTCGATGGTCGCCTCGCGCACCACTTCGAACTGCTCCTGCGaca                                                                                                >  2:86267/1‑134 (MQ=255)
                              gATCTCTCGCTCCAGCACTTCGCGCGAGCGGCGCACCAGGATGCCCTGCTCTTCCAGCGGGCTGATCAATTCCAGCAGGCCACCGACGTCCTCGATGGTCGCCTCGCGCACCACTTCGAACTGCTCCTGCGaca                                                                                                >  6:117515/1‑134 (MQ=255)
                              gATCTCTCGCTCCAGCACTTCGCGCGAGCGGCGCACCAGGATGCCCTGCTCTTCCAGCGGGCTGATCAATTCCAGCAGGCCACCGAAGTCCGCGAAGGTCGCCTCGCGCACCACTTCGAACTGCTCCTTCGACAcc                                                                                              >  8:105409/1‑136 (MQ=255)
                                   ctcGCTCCAGCACTTCGCGCGAGCGGCGCACCAGGATGCCCTGCTCTTCCAGCGGGCTGATCAATTCCAGCAGGCCACCGACGTCCTCGATGGTCGCCTCGCGCACCACTTCGAACTGCTCCTGCGACACCAGc                                                                                           >  3:40925/1‑134 (MQ=255)
                                     cGCTCCAGCACTTCGCGCGAGCGGCGCACCATGATGCCCTGCTCTTGCAGCGGGCTGATCAATTCCAGCAGGACACTGAAGTCCTCGATGGTCGCCTCGCTCACTACTTCGAACTGCTCCTGCGACTCCAGCg                                                                                          <  3:47846/133‑1 (MQ=255)
                                     cGCTCCAGCACTTCGCGCGAGCGGCGCACCAGGATGCCCTGCTCTTCCAGCGGGCTGATCAATTCCAGCAGGCCACCGACGTCCTCGATGGTCGCCTCGCGCACCACTTCGAACTGCTCCTGCGACACCAGCg                                                                                          <  1:159118/133‑1 (MQ=255)
                                     cGCTCCAGCACTTCGCGCGAGCGGCGCACCAGGATGCCCTGCTCTTCCAGCGGGCTGATCAATTCCAGCAGGCCACCGACGTCCTCGATGGTCGCCTCGCGCACCACTTCGAACTGCTCCTGCGACACCAGCg                                                                                          <  5:4254/133‑1 (MQ=255)
                                     cGCTCCAGCACTTCGCGCGAGCGGCGCACCAGGATGCCCTGCTCTTCCAGCGGGCTGATCAATTCCAGCAGGCCACCGACGTCCTCGATGGTCGCCTCGCGCACCACTTCGAACTGCTCCTGCGACACCAGCg                                                                                          <  5:238294/133‑1 (MQ=255)
                                     cGCTCCAGCACTTCGCGCGAGCGGCGCACCAGGATGCCCTGCTCTTCCAGCGGGCTGATCAATTCCAGCAGGCCACCGACGTCCTCGATGGTCGCCTCGCGCACCACTTCGAACTGCTCCTGCGACACCAGCg                                                                                          <  3:224314/133‑1 (MQ=255)
                                        tCCAGCACTTCGCGCGAGCGGCGCACCAGGATGCCCTGCTCTTCCAGCGGGCTGATCAATTCCAGCAGGCCACCGACGTCCTCGATGGTCGCCTCGCGCACCACTTCGAACTGCTCCTGCGACACCAGCGTaccg                                                                                     >  5:62768/1‑135 (MQ=255)
                                           aGCACTTCGCGCGAGCGGCGCACCAGGATGCCCTGCTCTTCCAGCGGGCTGATCAATTCCAGCAGGCCACCGACGTCCTCGATGGTCGCCTCGCGCACCACTTCGAACTGCTCCTGCGACACCAGCGTaccgccac                                                                                 <  7:31115/136‑1 (MQ=255)
                                            gCACTTCGCGCGAGCGGCGCACCAGGATGCCCTGCTCTTCCAGCGGGCTGATCAATTCCAGCAGGCCACCGACGTCCTCGATGGTCGCCTCGCGCACCTCTTCGAACTGCTCCTGCGACACCAGCGTaccgccac                                                                                 <  6:204755/135‑1 (MQ=255)
                                            gCACTTCGCGCGAGCGGCGCACCAGGATGCCCTGCGCTTCCAGCGGGCTGATCAATTCCAGCAGGCCACCGACGTCCTCGATGGTCGCCTCGCGCACCACTTCGAACTGCTCCTGCGACACCAGCGTaccgccac                                                                                 <  1:72229/135‑1 (MQ=255)
                                                   gcgcgAGCGGCGCACCAGGATGCCCTGCTCTTCCAGCGGGCTGATCAATTCCAGCAGGCCACCGACGTCCTCGATGGTCGCCTCGCGCACCACTTCGAACTGCTCCTGCGACACCAGCGTACCGCCACCGccc                                                                            <  1:198382/133‑1 (MQ=255)
                                                     gcgAGCGGCGCACCAGGATGCCCTGCTCTTCCAGCGGGCTGATCAATTCCAGCAGGCCACCGACGTCCTCGATGGTCGCCTCGCGCACCCCTTCGAACTGCTCCTGCGACCCCAGCGTACCACCACCGCCCCgtg                                                                        >  2:119828/1‑135 (MQ=255)
                                                     gcgAGCGGCGCACCAGGATGCCCTGCTCTTCCAGCGGGCTGATCAATTCCAGCAGGCCACCGACGTCCTCGATGGTCGCCTCGCGCACCACTTCGAACTGCTCCTGCGACACCAGCGTACCGCCCCCGCCCCgtg                                                                        >  4:18321/1‑135 (MQ=255)
                                                     gcgAGCGGCGCACCAGGATGCCCTGCTCTTCCAGCGGGCTGATCAATTCCAGCAGGCCACCGACGTCCTCGATGGACGCCTCGCGCACCACTTCGAACTGCTCCTGCGACACCAGCGTACCGCCACCGCCCCgtgt                                                                       >  8:148421/1‑136 (MQ=255)
                                                     gcgAGCGGCGCACCAGGATGCCCTGCTCTTCCAGCGGGCTGATCAATTCCAGCAGGCCACCGACGTCCTCGATCGTCGCCTCGCGCACCACTTCGAACTGCTCCTGCGACACCAGCCTCCCGCCACCCCCCCgag                                                                        >  6:27230/1‑133 (MQ=255)
                                                     gcgAGCGGCGCACCAGGATGCACTGCTCTTCCAGCGGGCTGATCAATTCCAGCAGGCAACCGACGTCCTCGATGGTCGCCTCGCGCCCCACTTCGAACTGCTCCTGCGCCACCAGCGTACCGCCACCGCCCCgtgt                                                                       >  8:221902/1‑136 (MQ=255)
                                                       gAGCGGCGCACCAGGATGCCCTGCTCTTCCAGCGGGCTGATCAATTCCAGCAGGCCACCGACGTCCTCGATGGTCGCCTCGCGCACCACTTCGAACTGCTCCTGCGACACCAGCGTACCGCCACCGCCCCGtgtg                                                                      >  8:216297/1‑135 (MQ=255)
                                                        aGCGGCGCACCAGGATGCCCTGCTCTTCCAGCGGGCTGATCAATTCCAGCAGGCCACCGACGTCCTCGATGGTCGCCTCGCGCACCACTTCGAACTGCTCCTGCGACACCAGCGTACCGCCACCGCCCCGTGTGa                                                                     <  4:244976/135‑1 (MQ=255)
                                                                         cccTGCTCTTCCAGCGGGCTGATCAATTCCAGCAGGCCACCGACGTCCTCGATGGTCGCCTCGCGCACCACTTCGAACTGCTCCTGCGACACCAGCGTACCGCCACCGCCCCGTGTGAACAGCTCGGTGagcagc                                                    <  8:129539/135‑1 (MQ=255)
                                                                          ccTGCTCTTCCAGGGGGCTGATCAATTCCAGCAGGCCACCGACGTCCTCGATGGTCGCCTCGCGCACCACTTCGAACTGCTCCTGCGACACCAGCGTACCGCCACCGCCCCGTGTGAACAGCTCGGTGagcagc                                                    <  6:196268/134‑1 (MQ=255)
                                                                          ccTGCTCTTCCAGCGGGCTGATCAATTCCAGCAGGCCACCGACGTCCTCGATGGTCGCCTCGCGCACCACTTCGAACTGCTCCTGCGACACCAGCGTACCGCCACCGCCCCGTGTGAACAGCTCGGTGagcagc                                                    <  5:112576/134‑1 (MQ=255)
                                                                          ccTGCTCTTCCAGCGGGCTGATCAATTCCAGCAGGCCACCGACGTCCTCGATGGTCGCCTCGCGCACCACTTCGAACTGCTCCTGCGACACCAGCGTACCGCCACCGCCCCGTGTGAACAGCTCGGTGagcagc                                                    <  3:32695/134‑1 (MQ=255)
                                                                          ccTGCTCTTCCAGCGGGCTGATCAATTCCAGCAGGCCACCGACGTCCTCGATGGTCGCCTCGCGCACCACTTCGAACTGCTCCTGCGACACCAGCGTACCGCCACCGCCCCGTGTGAACAGCTCGGTGagcagc                                                    <  1:64309/134‑1 (MQ=255)
                                                                          ccTGCTCTTCCAGCGGGCTGATCAATTCCAGCAGGCCACCGACGTCCTCGATGGTCGCCTCGCGCACCACTTCGAACTGCTCCTGCGACACCAGCGTACCGCCACCGCCCCGTGTGAACAGCTCGGTGagcagc                                                    <  5:27230/134‑1 (MQ=255)
                                                                           cTGCTCTTCCAGCGGGCTGATCAATTCCAGCAGGCCACCGACGTCCTCGATGGTCGCCTCGCGCACCACTTCGAACTGCTCCTGCGACACCAGCGTACCGCCACCGCCCCGTGTGAACAGCTCGGTGAGCAgcgc                                                  <  2:163370/135‑1 (MQ=255)
                                                                               tctTCCAGCGGGCTGATCAATTCCAGCAGGCCACCGACGTCCTCGATGGTCGCCTCGCGCACCACTTCGAACTGCTCCTGCGACACCAGCGTACCGCCACCGCCCCGTGTGAACAGCTCGGTGAGCAGCGCGcc                                               <  2:266296/134‑1 (MQ=255)
                                                                                      gCGGGCTGATCAATTCCAGCAGGCCCCCGACGTCCTCGATGGTCGCCTCGCGCACCACTTCGAACTGCTCCTGCGACACCAGCGTACCGCCACCGCCCCGTGTGAACAGCTCGGTGAGCAGCGCCCCGTCCTCg                                        >  1:171700/1‑134 (MQ=255)
                                                                                         ggCTGATCAATTCCAGCAGGCCACCGACGTCCTCGATGGTCACCTCGCGCACCACTTCGAACTGCTCCTGCGACACCAGCGTACCGCCACCGCCCCGTGTGAACAGCTCGGTGAGCAGCGCGCCGTCCTCGGCATa                                   <  8:16042/136‑1 (MQ=255)
                                                                                         ggCTGATCAATTCCAGCAGGCCACCGACGTCCTCGATGGTCACCTCGCGCACCACTTCGAACTGCGCCTGCGACACCAGCGTACCGCCACCGCCCCGTGTGACCAGCTCGGTGAGCAGCGCGCCGTCCTCGGCATa                                   <  8:11737/136‑1 (MQ=255)
                                                                                          gCTGATCAATTCCAGCAGGCCACCGACGTCCTCGATGGTCGCCTCGCGCACCACTTCGAACTGCTCCTGCGACACCAGCGTACCGCCACCGCGCCGTGTGAACAGCTCGGTGAGCAGCGCGCCGTCCTCGGCATa                                   <  2:158275/135‑1 (MQ=255)
                                                                                          gCTGATCAATTCCAGCAGGCCACCGACGTCCTCGATGGTCGCCTCGCGCACCACTTCGAACTGCTCCTGCGACACCAGCGTACCGCCACCGCCCCGTGTGAACAGCTCGGTGAGCAGCGCGCCGTCCTCGGCATa                                   <  5:150191/135‑1 (MQ=255)
                                                                                          gCTGATCAATTCCAGCAGGCCACCGACGTCCTCGATGGTCGCCTCGCGCACCACTTCGAACTGCTCCTGCGACACCAGCGTACCGCCACCGCCCCGTGTGAACAGCTCGGTGAGCAGCGCGCCGTCCTCGGCATa                                   <  2:41673/135‑1 (MQ=255)
                                                                                              aTCAATTCCAGCAGGCCACCGACGTCCTCGATGGTCGCCTCGCGCACCACTTCGAACTGCTCCTGCGCCACCAGCGTACCGCCACCGCCCCGTGTGAACAGCTCGGTGAGCATCGCGCCGTCCTCGGCATAGCTg                               >  6:207237/1‑135 (MQ=255)
                                                                                                      cagcagGCCACCGACGTCCTCGATGGTCGCCTCGCGCACCACTTCGAACTGCTCCTGCGACACCAGCGTACCGCCACCGCCCCGTGTGAACAGCTCGGTGAGCAGCGCGCCGTCCTCGGCATAGCTGACGATATg                       >  3:43940/1‑135 (MQ=255)
                                                                                                              cACCGACGTCCTCGATGGTCGCCTCGCGCACCACTTCGAACTGCTCCTGCGACACCAGCGTACCGCCACCGCCCCGTGTGAACAGCTCGGTGAGCAGCGCGCCGTCCTCGGCATAGCTGACGATATGGCTgcgcg               >  4:169771/1‑135 (MQ=255)
                                                                                                                ccGACGTCCTCGATGGTCGCCTCGCGCACCACTTCGAACTGCTCCTGCGACACCAGCGTACCGCCACCGCCCCGTGTGAACAGCTCGGTGAGCAGCGCGCCGTCCTCGGCATAGCTGACGATATGGCTGCGCGcc             <  5:142020/135‑1 (MQ=255)
                                                                                                                         tCGATGGTCGCCTCGCGCACCACTTCGAACTGCTCCTGCGACCCCAGCGTACCGCCACCGCCCCGTGTGAACAGCTCGGTGAGCAGCGCGCCGTCCTCGGCATAGCTGACGATATGGCTGCGCGCCACACCGCCt    >  3:121506/1‑135 (MQ=255)
                                                                                                                           gATGGTCGCCTCGCGCACCACTTCGAACTGCTCCTGCGACCCCAGCGTACCGCCACCGCCCCGTGTGAACAGCTCGGTGAGCAGCGCGCCGTCCTCGGCATAGCTGACGATATGGCTGCGCGCCACACCGCCttt  >  4:32182/1‑135 (MQ=255)
                                                                                                                                  |                                                                                                                               
GCCCTCGCGCTCGACCACGCTGAACTGCTCGATCTCTCGCTCCAGCACTTCGCGCGAGCGGCGCACCAGGATGCCCTGCTCTTCCAGCGGGCTGATCAATTCCAGCAGGCCACCGACGTCCTCGATGGTCGCCTCGCGCACCACTTCGAACTGCTCCTGCGACACCAGCGTACCGCCACCGCCCCGTGTGAACAGCTCGGTGAGCAGCGCGCCGTCCTCGGCATAGCTGACGATATGGCTGCGCGCCACACCGCCTTT  >  KT2440/5908393‑5908650

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