Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,695,997 C→T *63* (TAG→TAA yhjR ← DUF2629 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,695,9970CT100.0% 235.3 / NA 73*63* (TAG→TAAyhjRDUF2629 family protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (41/32);  total (41/32)

CCGGGCAGGCATCGACCACCAGGACATTTTCTCCCAACATTTGTAATGACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATCAGGCAGTGAAAATGCCTGCTTTAACGCCACAATATCATTCTG  >  NC_000913/3695865‑3696131
                                                                                                                                    |                                                                                                                                      
ccGGGCAGGCATCGACCACCAGGACATTTTCTCCCAACATTTGTAATGACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACttt                                                                                                                                  >  2:816652/1‑139 (MQ=255)
   ggcaggcaGCGACCACCAGGACATTTTCTCCCAACATTTGTAATGACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCTCTTACTTttgt                                                                                                                               >  2:1048602/1‑139 (MQ=255)
                  ccaGGACATTTTCTCCCAACATTTGTAATGACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGcc                                                                                                                >  2:991935/1‑139 (MQ=255)
                      gACATTTTCTCCCAACATTTGTAATGACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCa                                                                                                            >  2:532741/1‑139 (MQ=255)
                           tttCTCCCAACATTTGTAATGACCAGGCTAATGCGGCGGTGATAGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGc                                                                                                       <  1:1048602/139‑1 (MQ=255)
                             tctcCCAACATTTGTAATGACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCa                                                                                                     <  1:1524171/139‑1 (MQ=255)
                                               gACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGc                                                                                   >  2:843481/1‑139 (MQ=255)
                                                aCCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGcc                                                                                  >  2:1011271/1‑139 (MQ=255)
                                                  cAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCaa                                                                                >  2:1240362/1‑139 (MQ=255)
                                                     gCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTg                                                                             >  1:525973/1‑139 (MQ=255)
                                                          tGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGc                                                                                                       <  2:1466111/108‑1 (MQ=255)
                                                          tGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGc                                                                                                       >  1:1466111/1‑108 (MQ=255)
                                                          tGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATg                                                                                          <  1:1098616/121‑1 (MQ=255)
                                                          tGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATg                                                                                          >  2:1098616/1‑121 (MQ=255)
                                                          tGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTcgc                                                                        >  2:1550093/1‑139 (MQ=255)
                                                               cggTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTttg                                                                                                                                <  1:1571090/78‑1 (MQ=255)
                                                               cggTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTttg                                                                                                                                >  2:1571090/1‑78 (MQ=255)
                                                               cggTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGc                                                                                                                 <  2:346794/93‑1 (MQ=255)
                                                               cggTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGc                                                                                                                 >  1:346794/1‑93 (MQ=255)
                                                                         tgttgtCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCgg                                                         >  1:836488/1‑139 (MQ=255)
                                                                          gttgtCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGc                                                        <  1:532741/139‑1 (MQ=255)
                                                                           ttgtCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCa                                                       <  2:413057/139‑1 (MQ=255)
                                                                           ttgtCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCa                                                       <  1:1240362/139‑1 (MQ=255)
                                                                              tCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATaa                                                    >  1:964681/1‑139 (MQ=255)
                                                                                cccACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATc                                                  >  1:1611949/1‑139 (MQ=255)
                                                                                cccACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATc                                                  >  1:990718/1‑139 (MQ=255)
                                                                                  cACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGc                                                        >  2:1395691/1‑131 (MQ=255)
                                                                                  cACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGc                                                        <  1:1395691/131‑1 (MQ=255)
                                                                                      ccTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATc                                            >  2:1404711/1‑139 (MQ=255)
                                                                                        tCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATCAg                                          >  1:1520681/1‑139 (MQ=255)
                                                                                         ccccGCACCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATCAGGc                                        >  2:797327/1‑139 (MQ=255)
                                                                                          cccGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGc                                                        >  2:1519470/1‑123 (MQ=255)
                                                                                          cccGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGc                                                        <  1:1519470/123‑1 (MQ=255)
                                                                                          cccGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATCAGGc                                        <  1:1266703/139‑1 (MQ=255)
                                                                                            cGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATCAGGCAg                                      >  2:1469641/1‑139 (MQ=255)
                                                                                                 ccccTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAAt                                                   >  2:330593/1‑121 (MQ=255)
                                                                                                 ccccTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAAt                                                   <  1:330593/121‑1 (MQ=255)
                                                                                                   ccTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTttgt                                                                                                                               <  1:814128/43‑1 (MQ=255)
                                                                                                   ccTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTttgt                                                                                                                               >  2:814128/1‑43 (MQ=255)
                                                                                                   ccTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTgcg                                                                                                                           >  1:132022/1‑47 (MQ=255)
                                                                                                   ccTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTgcg                                                                                                                           <  2:132022/47‑1 (MQ=255)
                                                                                                   ccTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGcc                                                                                                      >  2:852977/1‑68 (MQ=255)
                                                                                                   ccTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGcc                                                                                                      <  1:852977/68‑1 (MQ=255)
                                                                                                   ccTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCcgcg                                                                                     >  1:577166/1‑85 (MQ=255)
                                                                                                   ccTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCcgcg                                                                                     <  2:577166/85‑1 (MQ=255)
                                                                                                    cTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGtt                                                                           <  2:1275057/94‑1 (MQ=255)
                                                                                                    cTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGtt                                                                           >  1:1275057/1‑94 (MQ=255)
                                                                                                       cTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATCAGGCAGTGAAAATGCCt                           >  1:546307/1‑139 (MQ=255)
                                                                                                           cccAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCaaa                                                                                                                       <  1:1158779/43‑1 (MQ=255)
                                                                                                           cccAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCaaa                                                                                                                       >  2:1158779/1‑43 (MQ=255)
                                                                                                           cccAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATCAGGCAGTGAAAATGCCTGCtt                       >  2:1520735/1‑139 (MQ=255)
                                                                                                               gTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATCAGGCAGTGAAAATGCCTGCTTTAAc                   <  2:1611949/139‑1 (MQ=255)
                                                                                                                 aCGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACCCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATCAGGCAGTGAAAATGCCTGCTTTAACGc                 <  2:492478/139‑1 (MQ=255)
                                                                                                                   ggCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATCAGGCAGTGAAAATGCCTGCTTTAACGCca               <  2:1520681/139‑1 (MQ=255)
                                                                                                                   ggCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATCAGGCAGTGAAAATGCCTGCTTTAACGCca               <  1:843481/139‑1 (MQ=255)
                                                                                                                    gCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAata                                                             >  1:688057/1‑92 (MQ=255)
                                                                                                                    gCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAata                                                             <  2:688057/92‑1 (MQ=255)
                                                                                                                     ccATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCTGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATCAGGCAGTGAAAATGCCTGCTTTAACGCcaca             >  2:1539753/1‑139 (MQ=255)
                                                                                                                     ccATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGc                                                        <  2:1609243/96‑1 (MQ=255)
                                                                                                                     ccATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGc                                                        >  1:1609243/1‑96 (MQ=255)
                                                                                                                      cATTCACCACTCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCa                                                                                                     >  2:1082402/1‑50 (MQ=38)
                                                                                                                      cATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCa                                                                                                     <  1:1082402/50‑1 (MQ=255)
                                                                                                                       aTTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATCAGGCAGTGAAAATGCCTGCTTTAACGCCACAat           >  1:1014860/1‑139 (MQ=255)
                                                                                                                       aTTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATCAGGCAGTGAAAATGCCTGCTTTAACGCCACAat           <  1:1400058/139‑1 (MQ=255)
                                                                                                                        ttCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGc                                                                                                       <  2:698616/46‑1 (MQ=255)
                                                                                                                        ttCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGc                                                                                                       >  1:698616/1‑46 (MQ=255)
                                                                                                                         tCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATCAGGCAGTGAAAATGCCTGCTTTAACGCCACAatat         <  2:990718/139‑1 (MQ=255)
                                                                                                                               aTCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTggg                                                                  >  1:1215015/1‑76 (MQ=255)
                                                                                                                               aTCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTggg                                                                  <  2:1215015/76‑1 (MQ=255)
                                                                                                                               aTCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATCAGGCAGTGAAAATGCCTGCTTTAACGCCACAATATCATTCt   >  1:365623/1‑139 (MQ=255)
                                                                                                                                tCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGc                                                                                                       >  1:186162/1‑38 (MQ=37)
                                                                                                                                tCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGc                                                                                                       <  2:186162/38‑1 (MQ=37)
                                                                                                                                tCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATCAGGCAGTGAAAATGCCTGCTTTAACGCCACAATATCATTCTg  <  2:964681/139‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                      
CCGGGCAGGCATCGACCACCAGGACATTTTCTCCCAACATTTGTAATGACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATCAGGCAGTGAAAATGCCTGCTTTAACGCCACAATATCATTCTG  >  NC_000913/3695865‑3696131

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