Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,222,944 4 bp→AGCA coding (2671‑2674/3093 nt) b3076 → evolved beta‑D‑galactosidase alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,222,9440CA98.4% 362.5 / ‑4.8 127Q891S (CAA→AGC) b3076evolved beta‑D‑galactosidase alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base A (66/58);  total (66/61)
Fisher's exact test for biased strand distribution p-value = 2.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.64e-01
*NC_0009133,222,9450AG100.0% 338.0 / NA 126Q891S (CAA→AGC) b3076evolved beta‑D‑galactosidase alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (66/60);  total (66/60)
*NC_0009133,222,9460AC99.2% 340.5 / NA 126Q891S (CAA→AGC) b3076evolved beta‑D‑galactosidase alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (64/60);  total (66/60)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009133,222,9470GA99.2% 360.3 / NA 125A892T (GCC→ACC) b3076evolved beta‑D‑galactosidase alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base A (66/57);  minor base C (0/1);  total (66/59)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 128 total aligned reads displayed.

GGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCCAAGCCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCA  >  NC_000913/3222881‑3223021
                                                               ||||                                                                          
ggTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGGACCGTCGATg                                                                   >  1:1193096/1‑76 (MQ=255)
ggTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATg                                                                   >  1:1674647/1‑76 (MQ=255)
 ggcGTGGACCGGGCAAAAACACCCCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCCGCACCGTCGATGc                                                                  <  2:354758/74‑1 (MQ=255)
    ggggACCGGGCGAAAACTCCGCCGCCAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCAt                                                               <  2:3142373/74‑1 (MQ=255)
    gTGGACCGGGCGAAAACTACNNCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCAt                                                               >  1:1819250/1‑76 (MQ=255)
    gTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCAt                                                               >  2:4718170/1‑76 (MQ=255)
    gTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGATAACATCATCGATATCTGGCGCAGCACCGTCGATGCCAt                                                               >  1:4668056/1‑76 (MQ=255)
    gTGGACCGGGAGAAAACTCCGCCGCCAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCAt                                                               <  1:1158340/76‑1 (MQ=255)
     tGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATg                                                              >  2:4125191/1‑76 (MQ=255)
     tGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATg                                                              >  2:3612322/1‑76 (MQ=255)
     gggACCGGGCGAAAACTCCGCCGCCAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATg                                                              <  2:1676231/75‑1 (MQ=255)
       gcccGGGCAAAACCCCCGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGtt                                                            <  1:795627/74‑1 (MQ=255)
       gACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGtt                                                            >  1:2165273/1‑76 (MQ=255)
       gACCGGCCAAAAACTCCCCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGtt                                                            <  2:5168195/76‑1 (MQ=255)
         ccGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCg                                                          >  1:2543316/1‑76 (MQ=255)
         ccGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGATCg                                                          >  2:5327252/1‑76 (MQ=255)
          cGGGCGAAACCTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCga                                                         <  1:3514331/76‑1 (MQ=255)
          cGGGCGAAAACTACGCCGCCAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCCCCGTCGATGCCATGTTCga                                                         <  2:4754805/76‑1 (MQ=255)
           gggCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATCTCTGGCGCAGCACCGGCGATGCCATGTTCgag                                                        >  2:4798714/1‑76 (MQ=255)
           gggCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCgag                                                        >  1:5052028/1‑76 (MQ=255)
           gggCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCgag                                                        >  2:2298265/1‑76 (MQ=255)
            ggCGAAACCCGCGCCTCCAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCgaga                                                       <  2:5322740/76‑1 (MQ=255)
             gCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGaa                                                      >  1:1094342/1‑76 (MQ=255)
              cGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAAc                                                     >  2:3730929/1‑76 (MQ=255)
               gAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACt                                                    >  1:5059301/1‑76 (MQ=255)
                aaaaCTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGCTGCCATGTTCGAGAACTa                                                   >  2:5676050/1‑76 (MQ=255)
                aaaaCTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGCTCGAGAACTa                                                   >  2:2130803/1‑76 (MQ=255)
                  aaCTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATc                                                 >  2:3643554/1‑76 (MQ=255)
                   aCTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATcc                                                <  1:3437259/76‑1 (MQ=255)
                    ctgcGCCGTGAGCCAGCAGGCTACCATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATccc                                               <  1:5808266/73‑1 (MQ=255)
                    cTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATccc                                               >  2:1619365/1‑76 (MQ=255)
                    cTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATccc                                               >  2:5382976/1‑76 (MQ=255)
                    cTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATccc                                               >  2:3891129/1‑76 (MQ=255)
                     tACGCCGCCAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCt                                              <  1:2939859/76‑1 (MQ=255)
                     tACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCt                                              >  2:2943707/1‑76 (MQ=255)
                     tACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCt                                              >  1:5371313/1‑76 (MQ=255)
                     tACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCACGTTCGAGAACTATCCCt                                              >  2:584094/1‑76 (MQ=255)
                      aCGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCtt                                             >  2:4122393/1‑76 (MQ=255)
                       cgccgcCAGCCAGCAGGCTAACATCATCGGTATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTc                                            <  2:3699259/76‑1 (MQ=255)
                       cgccgcCAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTc                                            <  1:3427520/76‑1 (MQ=255)
                       ccccGCCAGCCAGCAGGCTAGCATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGATAACTATCCCTTc                                            <  2:3135929/74‑1 (MQ=255)
                        gCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGACGCCATGTTCGAGAACTATCCCTTcc                                           >  2:4305023/1‑76 (MQ=255)
                         ccGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTccc                                          <  2:1393185/76‑1 (MQ=255)
                          cgccagccagcAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCg                                         <  2:1190672/73‑1 (MQ=255)
                          cGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCg                                         >  1:3077015/1‑76 (MQ=255)
                          cGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCg                                         >  1:3701997/1‑76 (MQ=255)
                           gACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGc                                        >  1:1731797/1‑76 (MQ=255)
                           gACAGCCAGCAGGCGAACATCATCGATAACAGGCCCAGCACCGTAGATGCCACCGTCGAGAACGAGCCCTGCCCga                                        >  2:1760970/1‑75 (MQ=255)
                            ccagccagcAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTAt                                                  <  2:3740256/64‑1 (MQ=255)
                            ccagccagcAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCa                                       <  1:157877/75‑1 (MQ=255)
                            aCAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGCTGCCATGTTCGAGAACTATCCCTTCCCGCa                                       >  2:343735/1‑76 (MQ=255)
                            aCAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTAt                                                  >  1:3740256/1‑65 (MQ=255)
                            aCAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCa                                       >  1:3213831/1‑76 (MQ=255)
                            aCAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCa                                       >  1:1449297/1‑76 (MQ=255)
                             cagccagcAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAg                                      >  1:1832685/1‑76 (MQ=255)
                              cgccagcAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGa                                     <  2:2792532/75‑1 (MQ=255)
                              agccagcAGGCTAACCTCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGa                                     >  2:3385947/1‑76 (MQ=255)
                              agccagcAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGa                                     <  1:3519264/76‑1 (MQ=255)
                              agccagcAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGa                                     >  2:4243722/1‑76 (MQ=255)
                               gccagcAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaa                                    <  2:3742442/76‑1 (MQ=255)
                                ccagcAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaac                                   >  1:4804708/1‑76 (MQ=255)
                                ccagcAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaac                                   >  1:55030/1‑76 (MQ=255)
                                ccagcAGGCTAACATCATCGATATCGGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaac                                   <  2:5675858/76‑1 (MQ=255)
                                 cagcagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaaca                                  <  1:3291564/76‑1 (MQ=255)
                                 cagcagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaaca                                  <  1:1582330/76‑1 (MQ=255)
                                 cagcagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaaca                                  <  1:2708822/76‑1 (MQ=255)
                                  agcagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaacaa                                 <  1:2885059/76‑1 (MQ=255)
                                  agcagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaacaa                                 <  2:4281987/76‑1 (MQ=255)
                                  agcagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaacaa                                 <  2:2081701/76‑1 (MQ=255)
                                    cgtgCTAACAGCGTCGATATCTGGCGCGGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACg                               <  2:5654061/73‑1 (MQ=255)
                                    cagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaacaact                               <  2:4711764/76‑2 (MQ=255)
                                    cagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACg                               <  1:118356/76‑1 (MQ=255)
                                    cagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACg                               >  2:5159657/1‑76 (MQ=255)
                                    cagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACg                               <  1:4755028/76‑1 (MQ=255)
                                    cagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACg                               <  1:5157005/76‑1 (MQ=255)
                                     agGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACgg                              <  1:1289237/76‑1 (MQ=255)
                                     agGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACgg                              >  1:797726/1‑76 (MQ=255)
                                       gCTAACATCATCTATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTa                            <  1:4857502/76‑1 (MQ=255)
                                       gCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTa                            >  1:4945944/1‑76 (MQ=255)
                                       gCTAACATCAGCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTa                            <  2:3260071/76‑1 (MQ=255)
                                        gtAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTaa                           <  2:4724811/75‑1 (MQ=255)
                                        cTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTaa                           <  2:1731730/76‑1 (MQ=255)
                                         tAACTTCATCGATGTCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAAc                          <  2:5269871/76‑1 (MQ=255)
                                         tAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAAc                          >  1:2174705/1‑76 (MQ=255)
                                          aaCATCATCGATATCTGGCGCAGCACTGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAAcc                         >  2:2306622/1‑76 (MQ=255)
                                          aaCATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAAcc                         >  1:1744775/1‑76 (MQ=255)
                                           aCATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCg                        <  2:391813/76‑1 (MQ=255)
                                           aCATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCg                        <  2:438999/76‑1 (MQ=255)
                                            catcatCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGt                       <  2:3206894/76‑1 (MQ=255)
                                             atcatcGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTc                      <  1:1504182/76‑1 (MQ=255)
                                              tcatcGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCa                     >  1:4733076/1‑76 (MQ=255)
                                               catcGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACCGTAACCGTCac                    >  2:29160/1‑75 (MQ=255)
                                                atcGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGc                   >  2:5366918/1‑76 (MQ=255)
                                                atcGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGc                   >  2:5509221/1‑76 (MQ=255)
                                                atcGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGc                   >  1:4204802/1‑76 (MQ=255)
                                                 tcGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAAcc                         >  1:566605/1‑69 (MQ=255)
                                                   gATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATg                <  2:1313953/76‑1 (MQ=255)
                                                    atatCTGGCGCAGCGCCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGt               <  1:1257398/76‑1 (MQ=255)
                                                      atCTGGCGCAGCACCGTCGATGCCATGTTAGAGAACTATGCCTTCCCGCAGAACAACGGTAACCGTCAGCGTGTcc             <  1:139456/76‑1 (MQ=255)
                                                                 ggCCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCa  <  2:3416628/75‑1 (MQ=255)
                                                               ||||                                                                          
GGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCCAAGCCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCA  >  NC_000913/3222881‑3223021

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