Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,960,576 G→T A618E (GCG→GAG)  recC ← exodeoxyribonuclease V subunit RecC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,960,5760GT100.0% 208.3 / NA 59A618E (GCG→GAG) recCexodeoxyribonuclease V subunit RecC
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (36/23);  total (36/23)

GAGTACAAATGTTAACCGGTCCGGCGAGAAAACGCTGGCTGATACGTTCTTGATCCAGACGCTGTGCCAGTTCATCACGCAATAGTGACAGCGGCACCGCGTCGCCATACTGCGCACCTAAACCTTCGGCGATAATCGCCTGCCATTGTTGTTCGATCAGCGTCATCGCCGCTTCGGTTTCCGCATCCGGCAGGAAGAAGGCGTTGAGCATATCGCGACAAACCGGCAACCACTCTTCCAGCGGACGCTCCTGTGCCAGCCCGCGACGCC  >  NC_000913/2960439‑2960708
                                                                                                                                         |                                                                                                                                    
gagTACAAATGTTAACCGGTCCGGCGAGAAAACGCTGGCTGATACGTTCTTGATCCAGACGCTGTGCCAGTTCATCACGCAATAGTGACAGCGGCACCGCGTCGCCATACTGCGCACCTAAACCTTCGGCGATAATCTcc                                                                                                                                    <  2:359186/140‑1 (MQ=255)
       aaTGTTAACCGGTCCGGCGAGAAAACGCTGGCTGATACGTTCTTGATCCAGACGCTGTGCCAGTTCATCACGCAATAGTGACAGCGGCACCGCGTCGCCATACTGCGCACCTAAACCTTCGGCGATAATCTCCTGCCAtt                                                                                                                             <  1:784735/140‑1 (MQ=255)
           ttAACCGGTCCGGCGAGAAAACGCTGGCTGATACGTTCTTGATCCAGACGCTGTGCCAGTTCATCACGCAATAGTGACAGCGGCACCGCGTCGCCATACTGCGCACCTAAACCTTCGGCGATAATCTCCTGCCAttgttg                                                                                                                         >  1:580577/1‑140 (MQ=255)
           ttAACCGGTCCGGCGAGAAAACGCTGGCTGATACGTTCTTGATCCAGACGCTGTGCCAGTTCATCACGCAATAGTGACAGCGGCACCGCGTCGCCATACTGCGCACCTAAACCTTCGGCGATAATCTCCTGCCAttgttg                                                                                                                         >  2:615494/1‑140 (MQ=255)
           ttAACCGGTCCGGCGAGAAAACGCTGGCTGATACGTTCTTGATCCAGACGCTGTGCCAGTTCATCACGCAATAGTGACAGCGGCACCGCGTCGCCATACTGCGCACCTAAACCTTCGGCGATAATCTCCTGCCAttgttg                                                                                                                         >  2:118004/1‑140 (MQ=255)
                     cGGCGAGAAAACGCTGGCTGATACGTTCTTGATCCAGACGCTGTGCCAGTTCATCACGCAATAGTGACAGCGGCACCGCGTCGCCATACTGCGCACCTAAACCTTCGGCGATAATCTCCTGCCATTGTTGTTCGATCAGc                                                                                                               <  1:515851/140‑1 (MQ=255)
                                    ggctgATACGTTCTTGATCCAGACGCTGTGCCAGTTCATCACGCAATAGTGACAGCGGCACCGCGTCGCCATACTGCGCACCTAAACCTTCGGCGATAATCTCCTGCCATTGTTGTTCGATCAGCGTCATCGCCGCTTCg                                                                                                <  1:296601/140‑1 (MQ=255)
                                               tCTTGATCCAGACGCTGTGCCAGTTCATCACGCAATAGTGACAGCGGCACCGCGTCGCCATACTGCGCACCTAAACCTTCGGCGATAATCTCCTGCCATTGTTGTTCGATCAGCGTCATCGCCGCTTCGGTTTCCGCATc                                                                                     >  2:442166/1‑140 (MQ=255)
                                                   gATCCAGACGCTGTGCCAGTTCATCACGCAATAGTGACAGCGGCACCGCGTCGCCATACTGCGCACCTAAACCTTCGGCGATAATCTCCTGCCATTGTTGTTCGATCAGCGTCATCGCCGCTTCGGTTTCCGCATCCGGc                                                                                 <  1:118004/140‑1 (MQ=255)
                                                       cAGACGCTGTGCCAGTTCATCACGCAATAGTGACAGCGGCACCGCGTCGCCATACTGCGCACCTAAACCTTCGGCGATAATCTCCTGCCATTGTTGTTCGATCAGCGTCATCGCCGCTTCGGTTTCCGCATCCGGCAGga                                                                             >  1:730393/1‑140 (MQ=255)
                                                        aGACGCTGTGCCAGTTCATCACGCAATAGTGACAGCGGCACCGCGTCGCCATACTGCGCACCTAAACCTTCGGCGATAATCTCCTGCCATTGTTGTTCGATCAGCGTCATCGCCGCTTCGGTTTCCGAATCCGGCAGgaa                                                                            >  2:460204/1‑140 (MQ=255)
                                                         gACGCTGTGCCAGTTCATCACGCAATAGTGACAGCGGCACCGCGTCGCCATACTGCGCACCTAAACCTTCGGCGATAATCTCCTGCCATTGTTGTTCGATCAGCGTCATCGCCGCTTCGGTTTCCGCATCCGGCAGgaag                                                                           >  1:577520/1‑140 (MQ=255)
                                                            gCTGTGCCAGTTCATCACGCAATAGTGACAGCGGCACCGCGTCGCCATACTGCGCACCTAAACCTTCGGCGATAATCTCCTGCCATTGTTGTTCGATCAGCGTCATCGCCGCTTCGGTTTCCGCATCCGGCAGGaagaag                                                                        <  2:730393/140‑1 (MQ=255)
                                                            gCTGTGCCAGTTCATCACGCAATAGTGACAGCGGCACCGCGTCGCCATACTGCGCACCTAAACCTTCGGCGATAATCTCCTGCCATTGTTGTTCGATCAGCGTCATCGCCGCTTCGGTTTCCGCATCCGGCAGGaagaag                                                                        <  2:577520/140‑1 (MQ=255)
                                                             cTGTGCCAGTTCATCACGCAATAGTGACAGCGGCACCGCGTCGCCATACTGCGCACCTAAACCTTCGGCGATAATCTCCGGCCATTGTTGTTCGATCAGCGTCATCGCCGCTTCGGTTTCCGCATCCGGCAGGCAGAAgg                                                                       >  1:18933/1‑140 (MQ=255)
                                                              tgtgCCAGTTCATCACGCAATAGTGACAGCGGCACCGCGTCGCCATACTGCGCACCTAAACCTTCGGCGATAATCTCCTGCCATTGTTGTTCGATCAGCGTCATCGCCGCTTCGGTTTCCGCATCCGGCAGGAAGAAGGc                                                                      >  1:421799/1‑140 (MQ=255)
                                                                tgCCAGTTCATCACGCAATAGTGACAGCGGCACCGCGTCGCCATACTGCGCACCTAAACCTTCGGCGATAATCTCCTGCCATTGTTGTTCGATCAGCGTCATCGCCGCTTCGGTTTCCGCATCCGGCAGGAAGAAGGCGt                                                                    >  2:513269/1‑140 (MQ=255)
                                                                tgCCAGTTCATCACGCAATAGTGACAGCGGCACCGCGTCGCCATACTGCGCACCTAAACCTTCGGCGATAATCTCCTGCCATTGTTGTTCGATCAGCGTCATCGCCGCTTCGGTTTCCGCATCCGGCAGGAAGAAGGCGt                                                                    >  1:279059/1‑140 (MQ=255)
                                                                tgCCAGTACATCACGCAATAGTGACAGCGGCACCGCGTCGCCATACTGCGCACCTAAACCTTCGGCGATAATCTCCTGCCATTGTTGTTCGATCAGCGTCATCGCCGCTTCGGTTTCCGCATCCGGCAGGAAGAAGGCGt                                                                    >  2:95841/1‑140 (MQ=255)
                                                                    aGTTCATCACGCAATAGTGACAGCGGCACCGCGTCGCCATACTGCGCACCTAAACCTTCGGCGATAATCTCCTGCCATTGTTGTTCGATCAGCGTCATCGCCGCTTCGGTTTCCGCATCCGGCAGGAAGAAGGCGTTGAg                                                                <  2:580577/140‑1 (MQ=255)
                                                                          tcaCGCAATAGTGACAGCGGCACCGCGTCGCCATACTGCGCACCTAAACCTTCGGCGATAATCTCCTGCCATTGTTGTTCGATCAGCGTCATCGCCGCTTCGGTTTCCGCATCCGGCAGGAAGAAGGCGTTGAGCATATc                                                          >  1:396412/1‑140 (MQ=255)
                                                                          tcaCGCAATAGTGACAGCGGCACCGCGTCGCCATACTGCGCACCTAAACCTTCGGCGATAATCTCCTGCCATTGTTGTTCGATCAGCGTCATCGCCGCTTCGGTTTCCGCATCCGGCAGGAAGAAGGCGTTGAGCATATc                                                          >  2:387863/1‑140 (MQ=255)
                                                                           caCGCAATAGTGACAGCGGCACCGCGTCGCCATACTGCGCACCTAAACCTTCGGCGATAATCTCCTGCCATTGTTGTTCGATCAGCGTCATCGCCGCTTCGGTTTCCGCATCCGGCAGGAAGAAGGCGTTGAGCATATcg                                                         >  1:380801/1‑140 (MQ=255)
                                                                              gCAATAGTGACAGCGGCACCGCGTCGCCATACTGCGCACCTAAACCTTCGGCGATAATCTCCTGCCATTGTTGTTCGATCAGCGTCATCGCCGCTTCGGTTTCCGCATCCGGCAGGAAGAAGGCGTTGAGCATATCGCGa                                                      <  1:632321/140‑1 (MQ=255)
                                                                                      gACAGCGGCACCGCGTCGCCATACTGCGCACCTAAACCTTCGGCGATAATCTCCTGCCATTGTTGTTCGATCAGCGTCATCGCCGCTTCGGTTTCCGCATCCGGCAGGAAGAAGGCGTTGAGCATATCGCGACAAACCgg                                              <  1:615494/140‑1 (MQ=255)
                                                                                       aCAGCGGCACCGCGTCGCCATACTGCGCACCTAAACCTTCGGCGATAATCTCCTGCCATTGTTGTTCGATCAGCGTCATCGCCGCTTCGGTTTCCGCATCCGGCAGGAAGAAGGCGTTGAGCATATCGCGACAAACCGGc                                             >  2:24816/1‑140 (MQ=255)
                                                                                       aCAGCGGCACCGCGTCGCCATACTGCGCACCTAAACCTTCGGCGATAATCTCCTGCCATTGTTGTTCGATCAGCGTCATCGCCGCTTCGGTTTCCGCATCCGGCAGGAAGAAGGCGTTGAGCATATCGCGACAAACCGGc                                             <  1:241668/140‑1 (MQ=255)
                                                                                            ggCACCGCGTCGCCATACTGCGCACCTAAACCTTCGGCGATAATCTCCTGCCATTGTTGTTCGATCAGCGTCATCGCCGCTTCGGTTTCCGCATCCGGCAGGAAGAAGGCGTTGAGCATATCGCGACAAACCGGCAAc                                          >  1:130810/1‑138 (MQ=255)
                                                                                            ggCACCGCGTCGCCATACTGCGCACCTAAACCTTCGGCGATAATCTCCTGCCATTGTTGTTCGATCAGCGTCATCGCCGCTTCGGTTTCCGCATCCGGCAGGAAGAAGGCGTTGAGCATATCGCGACAAACCGGCAAc                                          <  2:130810/138‑1 (MQ=255)
                                                                                            ggCACCGCGTCGCCATACTGCGCACCTAAACCTTCGGCGATAATCTCCTGCCATTGTTGTTCGATCAGCGTCATCGCCGCTTCGGTTTCCGCATCCGGCAGGAAGAAGGCGTTGAGCATATCGCGACAAACCGGCAACCa                                        >  1:516538/1‑140 (MQ=255)
                                                                                            ggCACCGCGTCGCCATACTGCGCACCTAAACCTTCGGCGATAATCTCCTGCCATTGTTGTTCGATCAGCGTCATCGCCGCTTCGGTTTCCGCATCCGGCAGGAAGAAGGCGTTGAGCATATCGCGACAAACCGGCAACCa                                        <  1:387863/140‑1 (MQ=255)
                                                                                            ggCACCGCGTCGCCATACTGCGCACCTAAACCTTCGGCGATAATCTCCTGCCATTGTTGTTCGATCAGCGTCATCGCCGCTTCGGTTTCCGCATCCGGCAGGAAGAAGGCGTTGAGCATATCGCGACAAACCGGCAACCa                                        >  1:384919/1‑140 (MQ=255)
                                                                                            ggCACCGCGTCGCCATACTGCGCACCTAAACCTTCGGCGATAATCTCCTGCCATTGTTGTTCGATCAGCGTCATCGCCGCTTCGGTTTCCGCATCCGGCAGGAAGAAGGCGTTGAGCATATCGCGACAAACCGGCAACCa                                        >  2:87039/1‑140 (MQ=255)
                                                                                               aCCGCGTCGCCATACTGCGCTACTAAACCTTCCGCGATAATCTCCTGCCATTGTTGTTCGATCAGCGTCATCTCTGCCTCGGTTTCCGCATCCGGCAGGAAGAAGGCGTTGAGCATATCGTGACAAACCGGCAACCActc                                     <  2:708154/140‑1 (MQ=255)
                                                                                               aCCGCGTCGCCATACTGCGCACCTAAACCTTCGGCGATAATCTCCTGCCATTGTTGTTCGATCAGCGTCATCGCCGCTTCGGTTTCCGCATCCGGCAGGAAGAAGGCGTTGAGCATATCGCGACAAACCGGCAACCa                                        <  1:661084/137‑1 (MQ=255)
                                                                                               aCCGCGTCGCCATACTGCGCACCTAAACCTTCGGCGATAATCTCCTGCCATTGTTGTTCGATCAGCGTCATCGCCGCTTCGGTTTCCGCATCCGGCAGGAAGAAGGCGTTGAGCATATCGCGACAAACCGGCAACCa                                        >  2:661084/1‑137 (MQ=255)
                                                                                                   cgTCGCCATACTGCGCACCTAAACCTTCGGCGATAATCTCCTGCCATTGTTGTTCGATCAGCGTCATCGCCGCTTCGGTTTCCGCATCCGGCAGGAAGAAGGCGTTGAGCATATCGCGACAAACCGGCAACCACTCTTcc                                 >  1:393297/1‑140 (MQ=255)
                                                                                                    gTCGCCATACTGCGCACCTAAACCTTCGGCGATAATCTCCTGCCATTGTTGTTCGATCAGCGTCATCGCCGCTTCGGTTTCCGCATCCGGCAGGAAGAAGGCGTTGAGCATATCGCGACAAACCGGCAACCACTCTTCCa                                <  2:384919/140‑1 (MQ=255)
                                                                                                     tCGCCATACTGCGCACCTAAACCTTCGGCGATAATCTCCTGCCATTGTTGTTCGATCAGCGTCATCGCCGCTTCGGTTTCCGCATCCGGCAGGAAGAAGGCGTTGAGCATATCGCGACAAACCGGCAACCACTCTTCCAg                               >  2:309454/1‑140 (MQ=255)
                                                                                                                cgcACCTAAACCTTCGGCGATAATCTCCTGCCATTGTTGTTCGATCAGCGTCATCGCCGCTTCGGTTTCCGCATCCGGCAGGAAGAAGGCGTTGAGCATATCGCGACAAACCGGCAACCACTCTTCCAGCGGACGCTccg                    >  2:145125/1‑139 (MQ=255)
                                                                                                                cgcACCTAAACCTTCGGCGATAATCTCCTGCCATTGTTGTTCGATCAGCGTCATCGCCGCTTCGGTTTCCGCATCCGGCAGGAAGAAGGCGTTGAGCATATCGCGACAAACCGGCAACCACTCTTCCAGCGGACGCTc                      >  2:99276/1‑138 (MQ=255)
                                                                                                                cgcACCTAAACCTTCGGCGATAATCTCCTGCCATTGTTGTTCGATCAGCGTCATCGCCGCTTCGGTTTCCGCATCCGGCAGGAAGAAGGCGTTGAGCATATCGCGACAAACCGGCAACCACTCTTCCAGCGGACGCTc                      >  2:314608/1‑138 (MQ=255)
                                                                                                                cgcACCTAAACCTTCGGCGATAATCTCCTGCCATTGTTGTTCGATCAGCGTCATCGCCGCTTCGGTTTCCGCATCCGGCAGGAAGAAGGCGTTGAGCATATCGCGACAAACCGGCAACCACTCTTCCAGCGGACGCTc                      <  1:174445/138‑1 (MQ=255)
                                                                                                                cgcACCTAAACCTTCGGCGATAATCTCCTGCCATTGTTGTTCGATCAGCGTCATCGCCGCTTCGGTTTCCGCATCCGGCAGGAAGAAGGCGTTGAGCATATCGCGACAAACCGGCAACCACTCTTCCAGCGGACGCTc                      >  2:174445/1‑138 (MQ=255)
                                                                                                                cgcACCTAAACCTTCGGCGATAATCTCCTGCCATTGTTGTTCGATCAGCGTCATCGCCGCTTCGGTTTCCGCATCCGGCAGGAAGAAGGCGTTGAGCATATCGCGACAAACCGGCAACCACTCTTCCAGCGGACGCTc                      <  1:99276/138‑1 (MQ=255)
                                                                                                                cgcACCTAAACCTTCGGCGATAATCTCCTGCCATTGTTGTTCGATCAGCGTCATCGCCGCTTCGGTTTCCGCATCCGGCAGGAAGAAGGCGTTGAGCATATCGCGACAAACCGGCAACCACTCTTCCAGCGGACGCTc                      <  1:314608/138‑1 (MQ=255)
                                                                                                                cgcACCTAAACCTTCGGCGATAATCTCCTGCCATTGTTGTTCGATCAGCGTCATCGCCGCTTCGGTTTCCGCATCCGGCAGGAAGAAGGCGTTGAGCATATCGCGACAAACCGGCAACCACTCTTCCAGCGGACGCTCCt                    >  2:251420/1‑140 (MQ=255)
                                                                                                                cgcACCTAAACCTTCGGCGATAATCTCCTGCCATTGTTGTTCGATCAGCGTCATCGCCGCTTCGGTTTCCGCATCCGGCAGGAAGAAGGCGTTGAGCATATCGCGACAAACCGGCAACCACTCTTCCAGCGGACGCTCCt                    >  2:89205/1‑140 (MQ=255)
                                                                                                                cgcACCTAAACCTTCGGCGATAATCTCCTGCCATTGTTGTTCGATCAGCGTCATCGCCGCTTCGGTTTCCGCATCCGGCAGGAAGAAGGCGTTGAGCATATCGCGACAAACCGGCAACCACTCTTCCAGCGGACGCTCCt                    >  2:790460/1‑140 (MQ=255)
                                                                                                                cgcACCTAAACCTTCGGCGATAATCTCCTGCCATTGTTGTTCGATCAGCGTCATCGCCGCTTCGGTTTCCGCATCCGGCAGGAAGAAGGCGTTGAGCATATCGCGACAAACCGGCAACCACTCTTCCAGCGGACGCTCCt                    >  1:511518/1‑140 (MQ=255)
                                                                                                                cgcACCTAAACCTTCGGCGATAATCTCCTGCCATTGTTGTTCGATCAGCGTCATCGCCGCTTCGGTTTCCGCATCCGGCAGGAAGAAGGCGTTGAGCATATCGCGACAAACCGGCAACCACTCTTCCAGCGGACGCTCCt                    >  2:182936/1‑140 (MQ=255)
                                                                                                                  cACCTAAACCTTCGGCGATAATCTCCTGCCATTGTTGTTCGATCAGCGTCATCGCCGCTTCGGTTTCCGCATCCGGCAGGAAGAAGGCGTTGAGCATATCGCGACAAACCGGCAACCACTCTTCCAGCGGa                           >  1:418348/1‑131 (MQ=255)
                                                                                                                  cACCTAAACCTTCGGCGATAATCTCCTGCCATTGTTGTTCGATCAGCGTCATCGCCGCTTCGGTTTCCGCATCCGGCAGGAAGAAGGCGTTGAGCATATCGCGACAAACCGGCAACCACTCTTCCAGCGGa                           <  2:418348/131‑1 (MQ=255)
                                                                                                                  cACCTAAACCTTCGGCGATAATCTCCTGCCATTGTTGTTCGATCAGCGTCATCGCCGCTTCGGTTTCCGCATCCGGCAGGAAGAAGGCGTTGAGCATATCGCGACAAACCGGCAACCACTCTTCCAGCGGACGCTCCtgt                  >  1:623664/1‑140 (MQ=255)
                                                                                                                    ccTAAACCTTCGGCGATAATCTCCTGCCATTGTTGTTCGATCAGCGTCATCGCCGCTTCGGTTTCCGCATCCGGCAGGAAGAAGGCGTTGAGCATATCGCGACAAACCGGCAACCACTCTTCCAGCGGACGCTCCTGTGc                <  2:279059/140‑1 (MQ=255)
                                                                                                                       aaaCCTTCGGCGATAATCTCCTGCCATTGTTGTTCGATCAGCGTCATCGCCGCTTCGGTTTCCGCATCCGGCAGGAAGAAGGCGTTGAGCATATCGCGACAAACCGGCAACCACTCTTCCAGCGGACGCTCCTGTGCCAg             >  1:617769/1‑140 (MQ=255)
                                                                                                                        aaCCTTCGGCGATAATCTCCTGCCATTGTTGTTCGATCAGCGTCATCGCCGCTTCGGTTTCCGCATCCGGCAGGAAGAAGGCGTTGAGCATATCGCGACAAACCGGCAACCACTCTTCCAGCGGACGCTCCTGTGCCAGc            >  1:255344/1‑140 (MQ=255)
                                                                                                                               ggCGATAATCTCCTGCCATTGTTGTTCGATCAGCGTCATCGCCGCTTCGGTTTCCGCATCCGGCAGGAAGAAGGCGTTGAGCATATCGCGACAAACCGGCAACCACTCTTCCAGCGGACGCTCCTGTGCCAGCCCGCGAc     <  1:251420/140‑1 (MQ=255)
                                                                                                                               ggCGATAATCTCCTGCCATTGTTGTTCGATCAGCGTCATCGCCGCTTCGGTTTCCGCATCCGGCAGGAAGAAGGCGTTGAGCATATCGCGACAAACCGGCAACCACTCTTCCAGCGGACGCTCCTGTGCCAGCCCGCGAc     >  1:19055/1‑140 (MQ=255)
                                                                                                                                  gATAATCTCCTGCCATTGTTGTTCGATCAGCGTCATCGCCGCTTCGGTTTCCGCATCCGGCAGGAAGAAGGCGTTGAGCATATCGCGACAAACCGGCAACCACTCTTCCAGCGGACGCTCCTGTGCCAGCCCGCGACGcc  <  1:790460/140‑1 (MQ=255)
                                                                                                                                  gATAATCTCCTGCAATTGTTGTTCGATCAGCGTCATCGCCGCTTCGGTTTCCGCATCCGGCAGGAAGAAGGCGTTGAGCATATCGCGACAAACCGGCAACCACTCTTCCAGCGGACGCTCCTGTGCCAGCCCGCGACGcc  <  1:95841/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                    
GAGTACAAATGTTAACCGGTCCGGCGAGAAAACGCTGGCTGATACGTTCTTGATCCAGACGCTGTGCCAGTTCATCACGCAATAGTGACAGCGGCACCGCGTCGCCATACTGCGCACCTAAACCTTCGGCGATAATCGCCTGCCATTGTTGTTCGATCAGCGTCATCGCCGCTTCGGTTTCCGCATCCGGCAGGAAGAAGGCGTTGAGCATATCGCGACAAACCGGCAACCACTCTTCCAGCGGACGCTCCTGTGCCAGCCCGCGACGCC  >  NC_000913/2960439‑2960708

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: