Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 4,108,360 T→C N48S (AAC→AGC)  metJ ← transcriptional repressor protein MetJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,108,3600TC99.3% 1146.5 / ‑4.0 278N48S (AAC→AGC) metJtranscriptional repressor protein MetJ
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (123/153);  total (124/154)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 291 total aligned reads displayed.

CGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCT  >  REL606/4108287‑4108432
                                                                         |                                                                        
cGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGgctg                                                                        >  1:4194667/1‑76 (MQ=255)
 ggCAAAGGATGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGgctgg                                                                       >  2:3382863/1‑75 (MQ=255)
   cAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGCTGTTc                                                                     <  2:3110197/76‑1 (MQ=255)
   cAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGCTGTTc                                                                     <  2:1998143/76‑1 (MQ=255)
       ggTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGCTGTTCACCt                                                                 <  2:4136356/76‑1 (MQ=255)
        gTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGCTGTTCACCTg                                                                <  2:4650592/76‑1 (MQ=255)
         ttGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGCTGTTCACCTgg                                                               <  2:5690760/76‑2 (MQ=255)
         ttGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGCTGTTCACCTga                                                               >  2:5950157/1‑76 (MQ=255)
         ttGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGCTGTTCACCTga                                                               <  2:5206931/76‑1 (MQ=255)
          tGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGCTGTTCa                                                                    <  2:2377576/70‑1 (MQ=255)
          tGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGCTGTTCACCTgac                                                              >  1:4411473/1‑76 (MQ=255)
          tGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGCTGTTCACCTgac                                                              <  1:2893310/76‑1 (MQ=255)
          tGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGCTGTTCACCTgac                                                              <  1:7155967/76‑1 (MQ=255)
          tGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGCTGTTCACCTgac                                                              >  1:7947410/1‑76 (MQ=255)
           gCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGCTGTTCACCTgacg                                                             <  1:7652723/76‑1 (MQ=255)
           gCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGCTGTTCACCTgacg                                                             <  1:3625227/76‑1 (MQ=255)
           gCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGCTGTTCACCTgacg                                                             <  2:4604680/76‑1 (MQ=255)
            cccGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGCTGTTCATCTgtcga                                                            >  2:6307667/1‑76 (MQ=255)
             ccGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGCTGTTCACCTgacgac                                                           <  2:3220347/76‑1 (MQ=255)
              cGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGCTGTTCACCTGacgacg                                                          >  1:5646261/1‑76 (MQ=255)
               ggTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGCTGTTCACCTGACGAcgc                                                         >  1:4760779/1‑76 (MQ=255)
                 tAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGCTGTTCACCTGACGACgcgt                                                       >  2:8308517/1‑76 (MQ=255)
                 tAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGCTGTTCACCTGACGACgcgt                                                       >  1:3936993/1‑76 (MQ=255)
                  aaaGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGCTGTTCACCTGACGAcgcggg                                                      >  1:5921784/1‑74 (MQ=255)
                   aaGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGCTGTTCACCTGACGACgcgtgc                                                     >  1:1972551/1‑76 (MQ=255)
                    aGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGCTGTTCACCTGACGACgcgtgcg                                                    <  1:5056477/76‑1 (MQ=255)
                    aGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGCTGTTCACCTGACGACgcgtgcg                                                    >  2:7130111/1‑76 (MQ=255)
                       cATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGCTGTTCACCTGACGACGCGTGCGTTc                                                 <  2:7619364/76‑1 (MQ=255)
                       cATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGCTGTTCACCTGACGACGCGTGCGTTc                                                 >  1:3111774/1‑76 (MQ=255)
                        aTGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGCTGTTCACCTGACGACGCGTGCGTTCa                                                >  2:6833194/1‑76 (MQ=255)
                        aTGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGCTGTTCACCTGACGACGCGGGCGTTCa                                                >  2:6653883/1‑76 (MQ=255)
                         tGCAGAAACGCTTCGCGCAGCAGCTCGCTGTTGGTAGCGTGACGCAGGCTGTTCACCTGACGACGCGTGAGTTCAt                                               <  1:6094194/76‑1 (MQ=39)
                         tGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGCTGTTCACCTGACGACGCGTGCGTTCAt                                               <  1:1758621/76‑1 (MQ=255)
                          gCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGCTGTTCACCTGACGACGCGTGCGTTCATc                                              <  2:6021795/76‑1 (MQ=255)
                            aGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGCTGTTCACCTGACGACGCGTGCGTTCATCgg                                            >  1:3736349/1‑76 (MQ=255)
                             gAAACGCTTCGCACAGCAGCTTGCTGTTGGTAGCGTGACGCAGGCTGTTCACCTGACGACGCGTGCGTTCATCGGt                                           <  2:2067287/76‑1 (MQ=255)
                             gAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGCTGTTCACCTGACGACGCGTGCGTTCATCGGt                                           <  2:2169366/76‑1 (MQ=255)
                              aaaCGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGCTGTTCACCTGACGACGCGTGCGTTCATCGGTg                                          <  2:2321409/76‑1 (MQ=255)
                               aaCGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGCTGTTCACCTGACGACGCGTGCGTTCATCGGTGa                                         >  2:5831610/1‑76 (MQ=255)
                                aCGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGCTGTTCACCTGACGACGCGTGCGTTCATCGGTGAg                                        <  2:1672527/76‑1 (MQ=255)
                                  gCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGCTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGa                                      <  1:6414480/76‑1 (MQ=255)
                                  gCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGCTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGa                                      >  2:6577863/1‑76 (MQ=255)
                                   cTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGCTGTTCACCTGACGACGCGTGCGTTCATCGGTGAgggt                                     >  1:6725971/1‑74 (MQ=255)
                                     tCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGCTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGAttt                                   <  1:511781/76‑1 (MQ=255)
                                      cGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGCTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGAtttt                                  >  1:7759760/1‑76 (MQ=255)
                                      cGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGCTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGAtttt                                  <  2:124956/76‑1 (MQ=255)
                                      cGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGCTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGAtttt                                  >  1:3529361/1‑76 (MQ=255)
                                       gCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGCTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGAttttt                                 <  1:1161307/76‑1 (MQ=255)
                                       gCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGCTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGAttttt                                 <  2:4935170/76‑1 (MQ=255)
                                        cacaGCAGCTCGCTGTTGGTAGCGTGACGCAGGCTGTTCACCTGACGACGCGTGCGTTCATCGGGGAGGATTTTTa                                >  2:6939169/1‑76 (MQ=255)
                                          cagcagCTCGCTGTTGGTAGCGTGACGCAGGCTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAac                              <  1:7669841/76‑1 (MQ=255)
                                           agcagcTCGCTGTTGGTAGCGTGACGCAGGCTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAaca                             <  1:1673541/76‑1 (MQ=255)
                                             cagcTCGCTGTTGGTAGCGTGACGCAGGCTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACAcc                           >  2:2560370/1‑76 (MQ=255)
                                             cagcTCGCTGTTGGTAGCGTGACGCAGGCTGTTCACCTGACGACGCGGGCGTTCATCGGTGAGGATTTTTAAAAcc                           >  2:1214858/1‑76 (MQ=255)
                                               gcTCGCTGTTGGTAGCGTGACGCAGGCTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCtt                         <  2:2745798/76‑1 (MQ=255)
                                               gcTCGCTGTTGGTAGCGTGACGCAGGCTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCtt                         <  2:8031523/76‑1 (MQ=255)
                                               gcTCGCTGTTGGTAGCGTGACGCAGGCTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCtt                         <  1:1513652/76‑1 (MQ=255)
                                                cTCGCTGTTGGTAGCGTGACGCAGGCTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTa                        <  2:6287262/76‑1 (MQ=255)
                                                 tCGCTGTTGGTAGCGTGACGCAGGCTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTaa                       <  2:4397747/76‑1 (MQ=255)
                                                  cGCTGTTGGTAGCGTGACGCAGGCTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAag                      <  1:2380311/76‑1 (MQ=255)
                                                  cGCTGTTGGTAGCGTGACGCAGGCTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAag                      <  2:5539373/76‑1 (MQ=255)
                                                   gCTGTTGGTAGCGTGACGCAGGCTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAaga                     <  1:2695806/76‑1 (MQ=255)
                                                     tGTTGGTAGCGTGACGCAGGCTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAgg                   >  1:5208764/1‑76 (MQ=255)
                                                     tGTTGGTAGCGTGACGCAGGCTGTTCACCTGACGACGCGTGCGTTCATCGGGGAGGATTTTTAACACCTTAAgggg                   >  1:5392781/1‑76 (MQ=255)
                                                      gTTGGTAGCGTGACGCAGGCTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGa                  <  2:2671207/76‑1 (MQ=255)
                                                       ttGGTAGCGTGACGCAGGCTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGaa                 <  2:3705491/76‑1 (MQ=255)
                                                       ttGGTAGCGTGACGCAGGCTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGaa                 <  1:872449/76‑1 (MQ=255)
                                                       ttGGTAGCGTGACGCAGGCTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGaa                 <  1:1604379/76‑1 (MQ=255)
                                                       ttGGTAGCGTGACGCAGGCTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGaa                 <  1:6664178/76‑1 (MQ=255)
                                                        tGGTAGCGTGACGCAGGCTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTCAGAGGAAt                >  2:1305558/1‑76 (MQ=255)
                                                          gTAGCGTGACGCAGGCTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATgg              <  1:2502782/76‑1 (MQ=255)
                                                          gTAGCGTGACGCAGGCTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATgg              <  1:752011/76‑1 (MQ=255)
                                                           tAGCGTGACGCAGGCTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGa             <  2:3070394/76‑1 (MQ=255)
                                                           tAGCGTGACGCAGGCTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGa             >  1:4941791/1‑76 (MQ=255)
                                                           tAGCGTGACGCAGGCTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGa             <  1:1208950/76‑1 (MQ=255)
                                                           tAGCGTGACGCAGGCTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGa             <  1:8158783/76‑1 (MQ=255)
                                                           tAGCGTGACGCAGGCTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGa             <  2:7824753/76‑1 (MQ=255)
                                                           tAGCGTGACGCAGGCTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGa             <  1:7222621/76‑1 (MQ=255)
                                                            aGCGTGACGCAGGCTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGaa            <  1:2213939/76‑1 (MQ=255)
                                                            aGCGTGACGCAGGCTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGaa            >  1:3809239/1‑76 (MQ=255)
                                                            aGCGTGACGCAGGCTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGaa            >  1:6172538/1‑76 (MQ=255)
                                                              cGTGACGCAGGCTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGGTTTTTAACACCTTAAGAGGAATGGAAAc          >  1:5861957/1‑76 (MQ=255)
                                                              cGTGACGCAGGCTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAAc          >  1:6251508/1‑76 (MQ=255)
                                                               gTGACGCAGGTTGTTCACCTGACGACGCGTGCGGTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAAcc         <  1:3966562/76‑1 (MQ=255)
                                                               gTGACGCAGGCTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAAcc         >  1:1295451/1‑76 (MQ=255)
                                                               gTGACGCAGGCTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAAcc         <  2:3526240/76‑1 (MQ=255)
                                                               gTGACGCAGGCTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAAcc         <  2:3250869/76‑1 (MQ=255)
                                                                tGACGCAGGCTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCg        >  2:5430398/1‑76 (MQ=255)
                                                                tGACGCAGGCTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCg        <  1:4644784/76‑1 (MQ=255)
                                                                 gACGCAGGCTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGt       >  1:996421/1‑76 (MQ=255)
                                                                 gACGCAGGCTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGt       <  1:4515807/76‑1 (MQ=255)
                                                                 gACGCAGGCTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGt       >  2:174794/1‑76 (MQ=255)
                                                                  aCGCAGGCTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTa      <  2:4729139/76‑1 (MQ=255)
                                                                   cGCAGGCTGTTCACCTGACGACGCGTTCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTaa     <  1:8364799/76‑1 (MQ=255)
                                                                   cGCAGGCTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTaa     <  2:2956394/76‑1 (MQ=255)
                                                                   cGCAGGCTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTaa     <  2:7433413/76‑1 (MQ=255)
                                                                   cGCAGGCTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTaa     <  2:6059/76‑1 (MQ=255)
                                                                    gCAGGCTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAAt    >  1:6503787/1‑76 (MQ=255)
                                                                    gCAGGCTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAAt    <  1:695004/76‑1 (MQ=255)
                                                                    gCAGGCTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAAt    <  2:3904172/76‑1 (MQ=255)
                                                                         |                                                                        
CGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCT  >  REL606/4108287‑4108432

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