Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,702,600 T→G K263Q (AAA→CAA)  era ← membrane‑associated, 16S rRNA‑binding GTPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,702,6000TG100.0% 216.5 / NA 57K263Q (AAA→CAA) eramembrane‑associated, 16S rRNA‑binding GTPase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (28/29);  total (28/29)

GCGCTGCCAGCCTTCCATCGGAGTTACTCTTAAAGATCGTCAACGTAACCGAGACTGCGCAGTGCGCGTTCGTCGTCGGCCCAACCGGATTTCACTTTTACCCACAGTTCAAGGTGAACAGGCGCTTCGAACATTTCCTGCATGTCTTTACGCGCTTCAATCCCGATGGTTTTGATCTTGGCCCCTTTGTTGCCAATGACCATCTTCTTCTGCCCTTCACGCTCAACGAGAATCAAACCGTTGATGTCATAACCACCGCGTTCGTTAGAGACGAAACGTTCGAT  >  NC_000913/2702452‑2702735
                                                                                                                                                    |                                                                                                                                       
gcgcTGCCAGCCTTCCATCGGAGTTACTCTTAAAGATCGTCAACGTAACCGAGACTGCGCAGTGCGCGTTCGTCGTCGGCCCAACCGGATTTCACTTTTACCCACAGTTCAAGGTGAACAGGCGCTTCGAACATTTCCTGCATGTCTTGAc                                                                                                                                       >  2:23285/1‑151 (MQ=255)
  gcTGCCAGCCTTCCATCGGAGTTACTCTTAAAGATCGTCAACGTAACCGAGACTGCGCAGTGCGCGTTCGTCGTCGGCCCAACCGGATTTCACTTTTACCCACAGTTCAAGGTGAACAGGCGCTTCGAACATTTCCTGCATGTCTTGAcgc                                                                                                                                     >  2:641018/1‑151 (MQ=255)
  gcTGCCAGCCTTCCATCGGAGTTACTCTTAAAGATCGTCAACGTAACCGAGACTGCGCAGTGCGCGTTCGTCGTCGGCCCAACCGGATTTCACTTTTACCCACAGTTCAAGGTGAACAGGCGCTTCGAACATTTCCTGCATGTCTTGAcgc                                                                                                                                     <  2:498195/151‑1 (MQ=255)
   cTGCCAGCCTTCCATCGGAGTTACTCTTAAAGATCGTCAACGTAACCGAGACTGCGCAGTGCGCGTTCGTCGTCGGCCCAACCGGATTTCACTTTTACCCACAGTTCAAGGTGAACAGGCGCTTCGAACATTTCCTGCATGTCTTGAcgcg                                                                                                                                    >  1:271698/1‑151 (MQ=255)
     gCCAGCCTTCCATCGGAGTTACTCTTAAAGATCGTCAACGTAACCGAGACTGCGCAGTGCGCGTTCGTCGTCGGCCCAACCGGATTTCACTTTTACCCACAGTTCAAGGTGAACAGGCGCTTCGAACATTTCCTGCATGTCTTGAcgcg                                                                                                                                    <  2:192438/149‑1 (MQ=255)
     gCCAGCCTTCCATCGGAGTTACTCTTAAAGATCGTCAACGTAACCGAGACTGCGCAGTGCGCGTTCGTCGTCGGCCCAACCGGATTTCACTTTTACCCACAGTTCAAGGTGAACAGGCGCTTCGAACATTTCCTGCATGTCTTGACgcgc                                                                                                                                   <  2:621201/150‑1 (MQ=255)
      ccAGCCTTCCATCGGAGTTACTCTTAAAGATCGTCAACGTAACCGAGACTGCGCAGTGCGCGTTCGTCGTCGGCCCAACCGGATTTCACTTTTACCCACAGTTCAAGGTGAACAGGCGCTTCGAACATTTCCTGCATGTCTTGACGCGCtt                                                                                                                                 <  1:589430/151‑1 (MQ=255)
              ccATCGGAGTTACTCTTAAAGATCGTCAACGTAACCGAGACTGCGCAGTGCGCGTTCGTCGTCGGCCCAACCGGATTTCACTTTTACCCACAGTTCAAGGTGAACAGGCGCTTCGAACATTTCCTGCATGTCTTGACGCGCTTCAATCCCg                                                                                                                         >  2:757318/1‑151 (MQ=255)
                aTCGGAGTTACTCTTAAAGATCGTCAACGTAACCGAGACTGCGCAGTGCGCGTTCGTCGTCGGCCCAACCGGATTTCACTTTTACCCACAGTTCAAGGTGAACAGGCGCTTCGAACATTTCCTGCATGTCTTGACGCGCTTCAATCCCGAt                                                                                                                       <  1:402036/151‑1 (MQ=255)
                 tCGGAGTTACTCTTAAAGATCGTCAACGTAACCGAGACTGCGCAGTGCGCGTTCGTCGTCGGCCCAACCGGATTTCACTTTTACCCACAGTTCAAGGTGAACAGGCGCTTCGAACATTTCCTGCATGTCTTGACGCGCTTCAATCCCGATg                                                                                                                      >  2:650890/1‑151 (MQ=255)
                    gAGTTACTCTTAAAGATCGTCAACGTAACCGAGACTGCGCAGTGCGCGTTCGTCGTCGGCCCAACCGGATTTCACTTTTACCCACAGTTCAAGGTGAACAGGCGCTTCGAACATTTCCTGCATGTCTTGACGCGCTTCAATCCCGATGGtt                                                                                                                   >  1:87049/1‑151 (MQ=255)
                      gTTACTCTTAAAGATCGTCAACGTAACCGAGACTGCGCAGTGCGCGTTCGTCGTCGGCCCAACCGGATTTCACTTTTACCCACAGTTCAAGGTGAACAGGCGCTTCGAACATTTCCTGCATGTCTTGACGCGCTTCAATCCCGATGGtttt                                                                                                                 <  2:576764/151‑1 (MQ=255)
                         aCTCTTAAAGATCGTCAACGTAACCGAGACTGCGCAGTGCGCGTTCGTCGTCGGCCCAACCGGATTTCACTTTTACCCACAGTTCAAGGTGAACAGGCGCTTCGAACATTTCCTGCATGTCTTGACGCGCTTCAATCCCGATGGTTTTGa                                                                                                               <  1:854913/150‑1 (MQ=255)
                           tctTAAAGATCGTCAACGTAACCGAGACTGCGCAGTGCGCGTTCGTCGTCGGCCCAACCGGATTTCACTTTTACCCACAGTTCAAGGTGAACAGGCGCTTCGAACATTTCCTGCATGTCTTGACGCGCTTCAATCCCGATGGTTTTGATCt                                                                                                            >  1:26099/1‑151 (MQ=255)
                                aaGATCGTCAACGTAACCGAGACTGCGCAGTGCGCGTTCGTCGTCGGCCCAACCGGATTTCACTTTTACCCACAGTTCAAGGTGAACAGGCGCTTCGAACATTTCCTGCATGTCTTGACGCGCTTCAATCCCGATGGTTTTGATCTTGGcc                                                                                                       >  1:318482/1‑151 (MQ=255)
                                   aTCGTCAACGTAACCGAGACTGCGCAGTGCGCGTTCGTCGTCGGCCCAACCGGATTTCACTTTTACCCACAGTTCAAGGTGAACAGGCGCTTCGAACATTTCCTGCATGTCTTGACGCGCTTCAATCCCGATGGTTTTGATCTTGGCCCCt                                                                                                    <  1:881509/151‑1 (MQ=255)
                                       tCAACGTAACCGAGACTGCGCAGTGCGCGTTCGTCGTCGGCCCAACCGGATTTCACTTTTACCCACAGTTCAAGGTGAACAGGCGCTTCGAACATTTCCTGCATGTCTTGACGCGCTTCAATCCCGATGGTTTTGATCTTGGCCCCTttgt                                                                                                >  1:250904/1‑151 (MQ=255)
                                        cAACGTAACCGAGACTGCGCAGTGCGCGTTCGTCGTCGGCCCAACCGGATTTCACTTTTACCCACAGTTCAAGGTGAACAGGCGCTTCGAACATTTCCTGCATGTCTTGACGCGCTTCAATCCCGATGGTTTTGATCTTGGCCCCTttgtt                                                                                               <  2:638822/151‑1 (MQ=255)
                                           cGTAACCGAGACTGCGCAGTGCGCGTTCGTCGTCGGCCCAACCGGATTTCACTTTTACCCACAGTTCAAGGTGAACAGGCGCTTCGAACATTTCCTGCATGTCTTGACGCGCTTCAATCCCGATGGTTTTGATCTTGGCCCCTTTGTTGcc                                                                                            >  1:558864/1‑151 (MQ=255)
                                           cGTAACCGAGACTGCGCAGTGCGCGTTCGTCGTCGGCCCAACCGGATTTCACTTTTACCCACAGTTCAAGGTGAACAGGCGCTTCGAACATTTCCTGCATGTCTTGACGCGCTTCAATCCCGATGGTTTTGATCTTGGCCCCTTTGTTGcc                                                                                            >  2:382927/1‑151 (MQ=255)
                                                ccGAGACTGCGCAGTGCGCGTTCGTCGTCGGCCCAACCGGATTTCACTTTTACCCACAGTTCAAGGTGAACAGGCGCTTCGAACATTTCCTGCATGTCTTGACGCGCTTCAATCCCGATGGTTTTGATCTTGGCCCCTTTGTTGCCAATGa                                                                                       <  1:134050/151‑1 (MQ=255)
                                                   agaCTGCGCAGTGCGCGTTCGTCGTCGGCCCAACCGGATTTCACTTTTACCCACAGTTCAAGGTGAACAGGCGCTTCGAACATTTCCTGCATGTCTTGACGCGCTTCAATCCCGATGGTTTTGATCTTGGCCCCTTTGTTGCCAATGACCa                                                                                    <  2:878024/151‑1 (MQ=255)
                                                     aCTGCGCAGTGCGCGTTCGTCGTCGGCCCAACCGGATTTCACTTTTACCCACAGTTCAAGGTGAACAGGCGCTTCGAACATTTCCTGCATGTCTTGACGCGCTTCAATCCCGATGGTTTTGATCTTGGCCCCTTTGTTGCCAATGACCAtc                                                                                  <  1:114307/151‑1 (MQ=255)
                                                       tGCGCAGTGCGCGTTCGTCGTCGGCCCAACCGGATTTCACTTTTACCCACAGTTCAAGGTGAACAGGCGCTTCGAACATTTCCTGCATGTCTTGACGCGCTTCAATCCCGATGGTTTTGATCTTGGCCCCTTTGTTGCCAATGACCAtctt                                                                                >  1:423983/1‑151 (MQ=255)
                                                          gcAGTGCGCGTTCGTCGTCGGCCCAACCGGATTTCACTTTTACCCACAGTTCAAGGTGAACAGGCGCTTCGAACATTTCCTGCATGTCTTGACGCGCTTCAATCCCGATGGTTTTGATCTTGGCCCCTTTGTTGCCAATGACCAtcttctt                                                                             <  1:405113/151‑1 (MQ=255)
                                                           cAGTGCGCGTTCGTCGTCGGCCCAACCGGATTTCACTTTTACCCACAGTTCAAGGTGAACAGGCGCTTCGAACATTTCCTGCATGTCTTGACGCGCTTCAATCCCGATGGTTTTGATCTTGGCCCCTTTGTTGCCAATGACCAtcttcttc                                                                            <  1:170370/151‑1 (MQ=255)
                                                            aGTGCGCGTTCGTCGTCGGCCCAACCGGATTTCACTTTTACCCACAGTTCAAGGTGAACAGGCGCTTCGAACATTTCCTGCATGTCTTGACGCGCTTCAATCCCGATGGTTTTGATCTTGGCCCCTTTGTTGCCAATGACCAtcttcttct                                                                           >  1:582053/1‑151 (MQ=255)
                                                              tGCGCGTTCGTCGTCGGCCCAACCGGATTTCACTTTTACCCACAGTTCAAGGTGAACAGGCGCTTCGAACATTTCCTGCATGTCTTGACGCGCTTCAATCCCGATGGTTTTGATCTTGGCCCCTTTGTTGCCAATGACCATCTTCTTCTGc                                                                         >  1:887959/1‑151 (MQ=255)
                                                                   gTTCGTCGTCGGCCCAACCGGATTTCACTTTTACCCACAGTTCAAGGTGAACAGGCGCTTCGAACATTTCCTGCATGTCTTGACGCGCTTCAATCCCGATGGTTTTGATCTTGGCCCCTTTGTTGCCAATGACCATCTTCTTCTGCCCTTc                                                                    >  1:732219/1‑151 (MQ=255)
                                                                    ttCGTCGTCGGCCCAACCGGATTTCACTTTTACCCACAGTTCAAGGTGAACAGGCGCTTCGAACATTTCCTGCATGTCTTGACGCGCTTCAATCCCGATGGTTTTGATCTTGGCCCCTTTGTTGCCAATGACCATCTTCTTCTGCCCTTc                                                                    >  2:694136/1‑150 (MQ=255)
                                                                    ttCGTCGTCGGCCCAACCGGATTTCACTTTTACCCACAGTTCAAGGTGAACAGGCGCTTCGAACATTTCCTGCATGTCTTGACGCGCTTCAATCCCGATGGTTTTGATCTTGGCCCCTTTGTTGCCAATGACCATCTTCTTCTGCCCTTc                                                                    >  2:85754/1‑150 (MQ=255)
                                                                         cgtcgGCCCAACCGGATTTCACTTTTACCCACAGTTCAAGGTGAACAGGCGCTTCGAACATTTCCTGCATGTCTTGACGCGCTTCAATCCCGATGGTTTTGATCTTGGCCCCTTTGTTGCCAATGACCATCTTCTTCTGCCCTTCACGCt                                                               <  2:263388/150‑1 (MQ=255)
                                                                            cgGCCCAACCGGATTTCACTTTTACCAACAGTTCAAGGTGAACAGGCGCTTCGAACATTTCCTGCATGTCTTGACGCGCTTCAATCCCGATGGTTTTGATCTTGGCCCCTTTGTTGCCAATGACCATCTTCTTCTGCCCTTCACGCTCaa                                                            <  2:773305/150‑1 (MQ=255)
                                                                             ggCCCAACCGGATTTCACTTTTACCCACAGTTCAAGGTGAACAGGCGCTTCGAACATTTCCTGCATGTCTTGACGCGCTTCAATCCCGATGGTTTTGATCTTGGCCCCTTTGTTGCCAATGACCATCTTCTTCTGCCCTTCACGCTCAACg                                                          <  1:172434/151‑1 (MQ=255)
                                                                                ccAACCGGATTTCACTTTTACCCACAGTTCAAGGTGAACAGGCGCTTCGAACATTTCCTGCATGTCTTGACGCGCTTCAATCCCGATGGTTTTGATCTTGGCCCCTTTGTTGCCAATGACCATCTTCTTCTGCCCTTCACGCTCAACg                                                          >  1:100928/1‑148 (MQ=255)
                                                                                ccAACCGGATTTCACTTTTACCCACAGTTCAAGGTGAACAGGCGCTTCGAACATTTCCTGCATGTCTTGACGCGCTTCAATCCCGATGGTTTTGATCTTGGCCCCTTTGTTGCCAATGACCATCTTCTTCTGCCCTTCACGCTCAACg                                                          >  2:905784/1‑148 (MQ=255)
                                                                                ccAACCGGATTTCACTTTTACCCACAGTTCAAGGTGAACAGGCGCTTCGAACATTTCCTGCATGTCTTGACGCGCTTCAATCCCGATGGTTTTGATCTTGGCCCCTTTGTTGCCAATGACCATCTTCTTCTGCCCTTCACGCTCAACg                                                          >  1:40456/1‑148 (MQ=255)
                                                                                  aaCCGGATTTCACTTTTACCCACAGTTCAAGGTGAACAGGCGCTTCGAACATTTCCTGCATGTCTTGACGCGCTTCAATCCCGATGGTTTTGATCTTGGCCCCTTTGTTGCCAATGACCATCTTCTTCTGCCCTTCACGCTCAACGAGAAt                                                     <  2:797835/151‑1 (MQ=255)
                                                                                     cGGATTTCACTTTTACCCACAGTTCAAGGTGAACAGGCGCTTCGAACATTTCCTGCATGTCTTGACGCGCTTCAATCCCGATGGTTTTGATCTTGGCCCCTTTGTTGCCAATGACCATCTTCTTCTGCCCTTCACGCTCAACGAGAATCaa                                                  <  1:679748/151‑1 (MQ=255)
                                                                                     cGGATTTCACTTTTACCCACAGTTCAAGGTGAACAGGCGCTTCGAACATTTCCTGCATGTCTTGACGCGCTTCAATCCCGATGGTTTTGATCTTGGCCCCTTTGTTGCCAATGACCATCTTCTTCTGCCCTTCACGCTCAACGAGAATCa                                                   >  2:729745/1‑150 (MQ=255)
                                                                                       gATTTCACTTTTACCCACAGTTCAAGGTGAACAGGCGCTTCGAACATTTCCTGCATGTCTTGACGCGCTTCAATCCCGATGGTTTTGATCTTGGCCCCTTTGTTGCCAATGACCATCTTCTTCTGCCCTTCACGCTCAACGAGAATCAAAc                                                >  1:896829/1‑151 (MQ=255)
                                                                                           tCACTTTTACCCACAGTTCAAGGTGAACAGGCGCTTCGAACATTTCCTGCATGTCTTGACGCGCTTCAATCCCGATGGTTTTGATCTTGGCCCCTTTGTTGCCAATGACCATCTTCTTCTGCCCTTCACGCTCAACGAGAATCAAACCGt                                             <  1:624940/150‑1 (MQ=255)
                                                                                                tttACCCACAGTTCAAGGTGAACAGGCGCTTCGAACATTTCCTGCATGTCTTGACGCGCTTCAATCCCGATGGTTTTGATCTTGGCCCCTTTGTTGCCAATGACCATCTTCTTCTGCCCTTCACGCTCAACGAGAATCAAACCGTTGATg                                        <  2:366028/150‑1 (MQ=255)
                                                                                                tttACCCACAGTTCAAGGTGAACAGGCGCTTCGAACATTTCCTGCATGTCTTGACGCGCTTCAATCCCGATGGTTTTGATCTTGGCCCCTTTGTTGCCAATGACCATCTTCTTCTGCCCTTCACGCTCAACGAGAATCAAACCGTTGATGt                                       >  2:784748/1‑151 (MQ=255)
                                                                                                 ttACCCACAGTTCAAGGTGAACAGGCGCTTCGAACATTTCCTGCATGTCTTGACGCGCTTCAACCCCGATGTTTTTGATCTTGGCCCCTTTGTTGCCAATGACCATTTTCTTCTGCCCTTCACGCTCAACGAGAATCAAACCGTTGATGt                                       <  2:866343/150‑1 (MQ=255)
                                                                                                   aCCCACAGTTCAAGGTGAACAGGCGCTTCGAACATTTCCTGCATGTCTTGACGCGCTTCAATCCCGATGGTTTTGATCTTGGCCCCTTTGTTGCCAATGACCATCTTCTTCTGCCCTTCACGCTCAACGAGAATCAAACCGTTGATGTCa                                     <  1:654314/150‑1 (MQ=255)
                                                                                                      cacaGTTCAAGGTGAACAGGCGCTTCGAACATTTCCTGCATGTCTTGACGCGCTTCAATCCCGATGGTTTTGATCTTGGCCCCTTTGTTGCCAATGACCATCTTCTTCTGCCCTTCACGCTCAACGAGAATCAAACCGTTGATGTCATAac                                 <  2:632722/151‑1 (MQ=255)
                                                                                                           ttCAAGGTGAACAGGCGCTTCGAACATTTCCTGCATGTCTTGACGCGCTTCAATCCCGATGGTTTTGATCTTGGCCCCTTTGTTGCCAATGACCATCTTCTTCTGCCCTTCACGCTCAACGAGAATCAAACCGTTGATGTCATAACCACcg                            >  1:770490/1‑151 (MQ=255)
                                                                                                                    aaCAGGCGCTTCGAACATTTCCTGCATGTCTTGACGCGCTTCAATCCCGATGGTTTTGATCTTGGCCCCTTTGTTGCCAATGACCATCTTCTTCTGCCCTTCACGCTCAACGAGAATCAAACCGTTGATGTCATAACCACCGcgttcgtt                    >  2:797909/1‑150 (MQ=255)
                                                                                                                              tCGAACATTTCCTGCATGTCTTGACGCGCTTCAATCCCGATGGTTTTGATCTTGGCCCCTTTGTTGCCAATGACCATCTTCTTCTGCCCTTCACGCTCAACGAGAATCAAACCGTTGATGTCATAACCACCGCGTTCGTTAGAGACGaaa          <  1:904032/150‑1 (MQ=255)
                                                                                                                              tCGAACATTTCCTGCATGTCTTGACGCGCTTCAATCCCGATGGTTTTGATCTTGGCCCCTTTGTTGCCAATGACCATCTTCTTCTGCCCTTCACGCTCAACGAGAATCAAACCGTTGATGTCATAACCACCGCGTTCGTTAGAGACGaa           <  2:835950/149‑1 (MQ=255)
                                                                                                                              tCGAACATTTCCTGCATGTCTTGACGCGCTTCAATCCCGATGGTTTTGATCTTGGCCCCTTTGTTGCCAATGACCATCTTCTTCTGCCCTTCACGCTCAACGAGAATCAAACCGTTGATGTCATAACCACCGCGTTCGTTAGAGACGAAAc         >  1:254639/1‑151 (MQ=255)
                                                                                                                                 aaCATTTCCTGCATGTCTTGACGCGCTTCAATCCCGATGGTTTTGATCTTGGCCCCTTTGTTGCCAATGACCATCTTCTTCTGCCCTTCACGCTCAACGAGAATCAAACCGTTGATGTCATAACCACCGCGTTCGTTAGAGACGAAACGtt      >  1:833999/1‑151 (MQ=255)
                                                                                                                                 aaCATTTCCTGCATGTCTTGACGCGCTTCAATCCCGATGGTTTTGATCTTGGCCCCTTTGTTGCCAATGACCATCTTCTTCTGCCCTTCACGCTCAACGAGAATCAAACCGTTGATGTCATAACCACCGCGTTCGTTAGAGACGAAACGtt      <  2:465337/151‑1 (MQ=255)
                                                                                                                                 aaCATTTCCTGCATGTCTTGACGCGCTTCAATCCCGATGGTTTTGATCTTGGCCCCTTTGTTGCCAATGACCATCTTCTTCTGCCCTTCACGCTCAACGAGAATCAAACCGTTGATGTCATAACCACCGCGTTCGTTAGAGACGAAACGtt      <  1:12074/151‑1 (MQ=255)
                                                                                                                                   cATTTCCTGCATGTCTTGACGCGCTTCAATCCCGATGGTTTTGATCTTGGCCCCTTTGTTGCCAATGACCATCTTCTTCTGCCCTTCACGCTCAACGAGAATCAAACCGTTGATGTCATAACCACCGCGTTCGTTAGAGACGAAACGTTCg    >  2:594638/1‑151 (MQ=255)
                                                                                                                                      ttCCTGCATGTCTTGACGCGCTTCAATCCCGATGGTTTTGATCTTGGCCCCTTTGTTGCCAATGACCATCTTCTTCTGCCCTTCACGCTCAACGAGAATCAAACCGTTGATGTCATAACCACCGCGTTCGTTAGAGACGAAACGTTCGAt  <  2:271693/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                       
GCGCTGCCAGCCTTCCATCGGAGTTACTCTTAAAGATCGTCAACGTAACCGAGACTGCGCAGTGCGCGTTCGTCGTCGGCCCAACCGGATTTCACTTTTACCCACAGTTCAAGGTGAACAGGCGCTTCGAACATTTCCTGCATGTCTTTACGCGCTTCAATCCCGATGGTTTTGATCTTGGCCCCTTTGTTGCCAATGACCATCTTCTTCTGCCCTTCACGCTCAACGAGAATCAAACCGTTGATGTCATAACCACCGCGTTCGTTAGAGACGAAACGTTCGAT  >  NC_000913/2702452‑2702735

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