Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,518,356 A→G intergenic (‑121/‑16) fecI ← / → ykiE RNA polymerase sigma factor FecI/protein YkiE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,518,3560AG100.0% 138.2 / NA 37intergenic (‑121/‑16)fecI/ykiERNA polymerase sigma factor FecI/protein YkiE
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (18/19);  total (18/19)

CAATGTCATCTGCATCAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTATAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTATTATTGACATCCTCACTGCCCTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCAC  >  NC_000913/4518112‑4518597
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                 
cAATGTCATCTGCATCAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAaggt                                                                                                                                                                                                                                             <  1:216355/251‑1 (MQ=255)
       aTCTGCATAAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGtat                                                                                                                                                                                                                                       >  2:135056/1‑250 (MQ=255)
               cAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCACAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGGTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATggg                                                                                                                                                                                                                               >  1:127360/1‑250 (MQ=255)
                                    gcgTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCAt                                                                                                                                                                                                         <  1:378044/251‑1 (MQ=255)
                                                            cGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCatgat                                                                                                                                                                                 <  1:77485/251‑1 (MQ=255)
                                                               gATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGGAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGAtgtg                                                                                                                                                                              >  2:70727/1‑251 (MQ=255)
                                                                    gtgCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCt                                                                                                                                                                          <  1:169854/250‑1 (MQ=255)
                                                                      gCCATATAGCGCGTCGAACGTTAAGGAAGCTGTGGGAGTGGCGAGGTGAGACATGCGGAGTGCAGCAAAAGTTAATGATCACATCGACATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGGTTCCAATAGACATTTTTAACATGTTGTTTTGCGAAGTGTGGTAAGGTAGGTATAAAATGGGATGGTGTTGCTGCTTCTGGCATGTGTCGGTCAGAATGGTTCATGATCTGGgctgct                                                                                                                                                                       <  2:410981/251‑1 (MQ=255)
                                                                                gACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGGCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGACATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtg                                                                                                                                                             <  1:288827‑M1/251‑10 (MQ=255)
                                                                                           tAAGGATGCTGTGGTAGTGGCGCGGTCAGAACTGCGGAGTGGATCAAAAGATAATTATCACGTAGTCCTATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGAGTTGTAAGGTAGGTATTAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGAATCATGATGTGGTCTGCTAatcaaggaggaattgtgtcg                                                                                                                                                  >  1:109660‑M1/1‑231 (MQ=255)
                                                                                           tAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtc                                                                                                                                                   <  2:152033‑M1/250‑20 (MQ=255)
                                                                                               ggagcagtggtagtgTCGAGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGGAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAGTTGTGTCCAATAGGCATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGGATAAAATGGGATGGAGCCTCTGCTACTGGCATGTGACGGTCAGAATGACACATGATGTGGTCTGCTActcaaggtggaatcgtgtcgcagg                                                                                                                                              <  2:37071‑M1/245‑25 (MQ=255)
                                                                                                     gtgGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtatacc                                                                                                                                         <  2:86705‑M1/250‑30 (MQ=255)
                                                                                                     gtgGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtatac                                                                                                                                          <  1:105967‑M1/249‑29 (MQ=255)
                                                                                                        gTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgcca                                                                                                                                     <  1:135056‑M1/251‑34 (MQ=255)
                                                                                                           gTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGTTTCGGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGTTActcaaggtggaattgtgtcgcagtataccgccagaa                                                                                                                                  <  2:323620‑M1/251‑37 (MQ=255)
                                                                                                                        aCATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaat                                                                                                                      >  2:306406‑M1/1‑202 (MQ=255)
                                                                                                                            gCGGAGTGCATCAAAAGTTAATTATCACGTAGCCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactct                                                                                                                 <  2:69257‑M1/251‑54 (MQ=255)
                                                                                                                             cGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctt                                                                                                                <  1:269452‑M1/251‑55 (MQ=255)
                                                                                                                                       cAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttac                                                                                                      <  1:105255‑M1/251‑65 (MQ=255)
                                                                                                                                         aaaGTGAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAACATTGTTTCCAATAGACAGTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATGAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggtgattaatactctttacctgttaccc                                                                                                    <  2:109660‑M1/251‑67 (MQ=255)
                                                                                                                                          aaGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccg                                                                                                   <  1:370056‑M1/251‑68 (MQ=255)
                                                                                                                                                 ttATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacgg                                                                                            >  1:162808‑M1/1‑177 (MQ=255)
                                                                                                                                                   aTCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggt                                                                                           >  1:337674‑M1/1‑175 (MQ=255)
                                                                                                                                                          aGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcc                                                                                   <  1:22527‑M1/251‑84 (MQ=255)
                                                                                                                                                               tatTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactcgttacctgttacccggctacggtatctatcccagt                                                                               <  1:120729‑M1/250‑88 (MQ=255)
                                                                                                                                                                  tAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttga                                                                            >  2:145327‑M1/1‑160 (MQ=255)
                                                                                                                                                                      atGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGCGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgc                                                                            >  2:271532‑M1/1‑156 (MQ=255)
                                                                                                                                                                      atGAGAAAGGTTATCATTACAATTGGAAATAAAATTGTTTCAAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTATGGTAGGTATAAAATGGGATGTCGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGCTGTGGTCTGCTActcaaggtggaattgtgtcgcagtatacccccagaaacggggattaatactcattacctgttacccggctacgctagctatcccagttga                                                                            >  2:387980‑M1/1‑156 (MQ=255)
                                                                                                                                                                           aaTGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcaccc                                                                   >  2:108667‑M1/1‑151 (MQ=255)
                                                                                                                                                                                  aTCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTTAGAGTTGTAAGGAAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTAAGAATGACTCATGATGTGACCTGCTActcacggcggaattgtgtcgcagtataccgccagacacggggattaatactctttacctgttacccggctacggtagctatcccagttgacgatcacccacaaaacg                                                           >  2:365023‑M1/1‑144 (MQ=255)
                                                                                                                                                                                            ttGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggc                                                 >  2:397059‑M1/1‑134 (MQ=255)
                                                                                                                                                                                                         tGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccat                                    >  1:107971‑M1/1‑121 (MQ=255)
                                                                                                                                                                                                                           gtaACATGTTCTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcgcccagca                   >  2:291539‑M1/2‑103 (MQ=255)
                                                                                                                                                                                                                             aaCATGTTGTTTTTCTAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcgcccagcaact                >  1:140179‑M1/1‑101 (MQ=255)
                                                                                                                                                                                                                                            tAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcgcccagcaacttggtcatctcttcc  >  2:214688‑M1/1‑86 (MQ=255)
                                                                                                                                                                                                                                            tAAGTGTTGTAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTActcaaggtggaattgtgtcgcagtataccgccagaaacggggattaatactctttacctgttacccggctacggtagctatcccagttgaagatcacccacaaactgttacccaggcgcatacggtccatcacgcgttcgcccagcaacttggtcatctcttcc  >  1:62281‑M1/1‑86 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                 
CAATGTCATCTGCATCAAAAGCAGACTGGAGTTTGCGCGTCAGCCAGCTTTTCAACCAGCCGTGATGTGTGCCATAAAGCGACTCGAACGTTAAGGAAGCTGTGGTAGTGGCGCGGTCAGACATGCGGAGTGCATCAAAAGTTAATTATCACGTAGTCATATTAATATGAGAATGGTTATCATTACAATTGGAAATAAAATTGTTTCCAATAGACATTTTTAACATGTTGTTTTTCTAAGTGTTATAAGGTAGGTATAAAATGGGATGGAGCCTCTGCTTCTGGCATGTGTCGGTCAGAATGACTCATGATGTGGTCTGCTATTATTGACATCCTCACTGCCCTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCAC  >  NC_000913/4518112‑4518597

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