Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,470,951 A→C 11.9% intergenic (‑38/+35) ridA ← / ← pyrI enamine/imine deaminase, redox regulated chaperone/aspartate carbamoyltransferase, PyrI subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,470,9510AC11.9% 119.5 / 6.9 58intergenic (‑38/+35)ridA/pyrIenamine/imine deaminase, redox regulated chaperone/aspartate carbamoyltransferase, PyrI subunit
Reads supporting (aligned to +/- strand):  ref base A (21/28);  new base C (2/5);  total (25/33)
Fisher's exact test for biased strand distribution p-value = 6.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.32e-01

CGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTC  >  NC_000913/4470812‑4471093
                                                                                                                                           |                                                                                                                                              
cGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGTTAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCGTta                                                                                                                                             >  1:482839/1‑143 (MQ=255)
  ggATCTGCGCGGAGGTGATGATCAGGTTGCCGAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATtata                                                                                                                                           <  2:767259/143‑1 (MQ=255)
  ggATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAACAGTCTGGCTaataatc                                                                                                                                           >  2:491978/1‑142 (MQ=255)
  ggATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATtata                                                                                                                                           <  2:938650/143‑1 (MQ=255)
    aTCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATCATAgg                                                                                                                                         >  2:677661/1‑143 (MQ=255)
     tCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCAGTAAAAAATATGGCACATTATAggg                                                                                                                                        >  1:613620/1‑143 (MQ=255)
        gcccGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGc                                                                                                                                     <  2:137594/141‑1 (MQ=255)
         aCCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATGTCTCCTTTATTACAGCGGTAAGAAGTCTCGCTCATTAGAGGGAGcg                                                                                                                                    >  2:344195/1‑142 (MQ=255)
           cGGAGGTGATGGTCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCTGTCGCGATAGTTTTGCTAATGATTTCTCCTTTATTACACCGGTAAAAACTCTGGCTCATTATACGTAGGCAg                                                                                                                                  >  2:477178/1‑143 (MQ=255)
            ggAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCCCTATAATGAGCCAGa                                                                                                                                 <  1:401960/143‑1 (MQ=255)
             gAGGTTATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCCCTATAATGAGCCAGAc                                                                                                                                <  1:227206/143‑1 (MQ=255)
                 tgatgaTCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACtttt                                                                                                                            <  1:621660/143‑1 (MQ=255)
                   atgatCATATTGCCCAGTGCTACGCCCTGTACGTAAGGACCGAGAGTTGGCGTGGCATTTTCCGTCGGTATAGGTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGTTCCCTATAATGAGCCAGACttttta                                                                                                                          <  2:559495/143‑3 (MQ=255)
                   aggctCATATTGCCCGGATCACCGCCTTGTACGTAAGGCCCGCTAGCTGCGGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATGTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCCCTATAATGAGCCAGACttttta                                                                                                                          <  2:370539/139‑3 (MQ=25)
                     gatCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTCTAGGGAGCCAGACTTttatta                                                                                                                        >  2:427560/1‑143 (MQ=255)
                     gatCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTttatta                                                                                                                        <  1:701521/143‑1 (MQ=255)
                         gtatTGCCCACGTCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGCGCATTTTCCGTCGCGATAGTTTCGCTCATGATTTCTCCTTTATTACCGCGGTAATTAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCaa                                                                                                                    <  2:202085/142‑1 (MQ=255)
                                     tCAGCGCCCTGTACGTAGGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAACGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCcattaatg                                                                                                        <  2:880088/143‑2 (MQ=255)
                                      cAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCaattaatta                                                                                                       >  1:129597/1‑142 (MQ=255)
                                       atcGCCCTGTACGTTAGGTCCGTTAGCTGCCGGTGCATTTTCCGGCGTGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTgg                                                                                                      <  2:16260/141‑1 (MQ=255)
                                       aaCGCCCTGTGCGTAATGACCGTTAGCTGCCGGTGCTTTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACTGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTgg                                                                                                      <  2:235816/143‑1 (MQ=255)
                                         cGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGcc                                                                                                    <  2:7999/143‑1 (MQ=255)
                                           cccTGTACGTAAGGACCGATAGCCGCCGGTCCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAg                                                                                                  <  2:46714/143‑1 (MQ=255)
                                            ccTGTACGTCAGGACCAATAGCTCCCCGTGCCCTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCCTTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGc                                                                                                 <  1:902137/143‑1 (MQ=255)
                                                tACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCGTTATAGGGAGCCAGACCTTTATTACTAAACGGAATCAACTGGCCAGaacca                                                                                             >  2:679305/1‑143 (MQ=255)
                                                   gTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACAt                                                                                          >  2:308595/1‑143 (MQ=255)
                                                        gACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATgg                                                                                     <  1:317433/143‑1 (MQ=255)
                                                            gagtgCTGTCGGCGCAGTTTCCGGCGCGATCGTGTTGCTGTTGGGTGCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGaaa                                                                                 <  1:298945/139‑1 (MQ=255)
                                                                gCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAActc                                                                             >  1:748121/1‑143 (MQ=255)
                                                                  cgCCGGTGCATTTGCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCtt                                                                           <  1:786983/142‑1 (MQ=255)
                                                                          cATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCGAGCACCACATTATGGGAAAAATCTTTTTccaag                                                                   >  1:722148/1‑139 (MQ=255)
                                                                            ttttCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACATTTATCACCAACCGCAATTACTTGGCCAGCACCACATGATGGGAACACTATTATACACAgca                                                                 >  2:225025/1‑141 (MQ=255)
                                                                                  gTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGCGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCAGTTTCACAGTATTTACa                                                           >  2:20429/1‑143 (MQ=255)
                                                                                   tCGCGATAGGTTTGCTCATGATTTGTCCTTTATTACGAAGCTGAAAAGTCTGGCTCATTATAGGGAGCCGGACGCTTAAGACCAACCGCAATCAATTTGCCACCACCACATGATGGGAAAACTCTTTTTCACAGTATTGGCAt                                                          >  2:239494/1‑143 (MQ=255)
                                                                                     gcgATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGTCTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCAttt                                                        <  1:308595/143‑1 (MQ=255)
                                                                                           gTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAgcgc                                                  <  1:849950/143‑1 (MQ=255)
                                                                                            ttttGCTCATGATTTCTCCTTTATTACCGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGcgcg                                                 <  1:859295/143‑1 (MQ=255)
                                                                                             tttGCTCATGATTTCTCCTTTATTACAGCGTTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTCCCCACTATAATTAGTTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGa                                                <  2:871844/143‑1 (MQ=255)
                                                                                              ttGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCATTTTCACAGTATTTGCATTTGAGCGCCat                                               >  2:639470/1‑143 (MQ=255)
                                                                                                   cttGATTTCGCCTTCCTTACAGCGTTAGAAAGTCTGGCTCATGATAGGGAGCCTCACTTTCATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCAt                                          <  2:280997/141‑1 (MQ=255)
                                                                                                        tttCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGgcg                                     >  2:578142/1‑143 (MQ=255)
                                                                                                            gccgtttttaCAGCGGGAAAAGGCCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGtt                                 <  1:499254/136‑1 (MQ=255)
                                                                                                            gccTTTATTACAGCGGTAAAAAGGCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAATCTCGTGTTCTTAGTATTTGCATTTGAGCGCGATATCATTGGCGCGtt                                 <  1:179021/142‑1 (MQ=255)
                                                                                                             ccTTTATTACATCTGTAAAAAGTCTGCTGCATTCTAGGGAGCAAGGCTTTAGCCTCCATGCGGAATTCGTTCGCAAGCACCACATTATGGGAATACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGttt                                <  2:918897/143‑1 (MQ=255)
                                                                                                               tttATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTc                              <  2:209768/143‑1 (MQ=255)
                                                                                                                 tattaCAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGTATTTGCGCGCGATAGCCTTGGCGCGGTTTCgg                            >  1:585057/1‑142 (MQ=255)
                                                                                                                       aGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAACACTCTTTTTCCCAGTATTTGCATTTGAGCGCGCTCTCCTTGGCGCCTTTTCGCACGGCa                      >  2:96522/1‑143 (MQ=255)
                                                                                                                             aaaaaGTCTGGATCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTg                <  1:221997/143‑1 (MQ=255)
                                                                                                                                 aGTCTGGCTCATTATAGGGAGCCAGACTTTTTTTATCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATg            >  2:334758/1‑143 (MQ=255)
                                                                                                                                  gTCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCGTTGGCGCGTCTTCGCACGTCAACACTGGATGa           >  1:802069/1‑143 (MQ=255)
                                                                                                                                  gTCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGCTTTCGCACGGCAAAGCTGGATGa           >  1:942628/1‑143 (MQ=255)
                                                                                                                                    cTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCCAAGCTGGCTGaaa         >  2:485391/1‑143 (MQ=255)
                                                                                                                                     tGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTTTTTTTCACAGTATTTGCGTTTGAGCGCGCTATCATTGGCGCGGTTTCGCACGGCAAAGCTGGGTGAAAc        >  1:362515/1‑143 (MQ=255)
                                                                                                                                     tGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCaaa                    <  2:708011/131‑1 (MQ=255)
                                                                                                                                     tGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCaaa                    >  1:708011/1‑131 (MQ=255)
                                                                                                                                      ggCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAAcc       <  1:862894/143‑1 (MQ=255)
                                                                                                                                       gCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCg      <  1:907164/143‑1 (MQ=255)
                                                                                                                                         tCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATAAAACCGGt    >  2:911572/1‑143 (MQ=255)
                                                                                                                                          cATTATAGGGAGCCAGACTGTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCGTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGtt   <  1:6363/143‑1 (MQ=255)
                                                                                                                                           aTTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTc  <  1:400153/143‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                              
CGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTC  >  NC_000913/4470812‑4471093

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