Predicted mutation
evidence seq id position mutation annotation gene description
RA USA300TCH1516_ALE 2,010,363 C→T C146Y (TGT→TAT)  USA300HOU_RS10125 ← Putative multidrug export ATP‑binding/permease protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,010,3630CT98.4% 193.7 / NA 62C146Y (TGT→TAT) USA300HOU_RS10125Putative multidrug export ATP‑binding/permease protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (25/36);  minor base G (0/1);  total (25/37)
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

TTTCACGTGTCAATTTTCTTAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATACAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTAAAGCTTGTAAATGGTTGTAT  >  USA300TCH1516_ALE/2010226‑2010500
                                                                                                                                         |                                                                                                                                         
tttCACGTGTCAATTTTCTTAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGACCATTATGGATAGTGCAATAATAATTGTTATATAAt                                                                                                                                        >  2:259838/1‑141 (MQ=255)
tttCACGTGTCAATTTTCTTAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAAt                                                                                                                                        >  1:649622/1‑141 (MQ=255)
     cGTGTCAATTTTCTTAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAAc                                                                                                                                   <  1:138844/141‑1 (MQ=255)
      gtgtCAATTTTCTTAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAAcc                                                                                                                                  >  2:36133/1‑141 (MQ=255)
      gtgtCAATTTTCTTAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAAcc                                                                                                                                  >  1:281956/1‑141 (MQ=255)
             ttttCTTAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATGAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTCTATAATCTAACCAAATAtt                                                                                                                           >  1:613991/1‑141 (MQ=255)
               ttCTTAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAACTATTCa                                                                                                                         >  2:661601/1‑141 (MQ=255)
                              aGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATCATCATTGTTTTATAATCTAACCAAATATTCATTAACCCGGTTaaaa                                                                                                          >  1:251742/1‑141 (MQ=255)
                                 aaaCGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTCAGATGa                                                                                                       >  2:497059/1‑141 (MQ=255)
                                 aaaCGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGa                                                                                                       <  2:391823/141‑1 (MQ=255)
                                 aaaCGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGa                                                                                                       <  1:892933/141‑1 (MQ=255)
                                   aCGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGaaa                                                                                                     >  2:167381/1‑141 (MQ=255)
                                    cGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAAt                                                                                                    <  1:917269/141‑1 (MQ=255)
                                      tACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAAATTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCt                                                                                                  >  1:423390/1‑141 (MQ=255)
                                       acacCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCtt                                                                                                 <  2:514251/141‑1 (MQ=255)
                                             ttAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTtttgttt                                                                                           <  1:267142/141‑1 (MQ=255)
                                             ttAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTtttgttt                                                                                           <  1:850406/141‑1 (MQ=255)
                                                 aatgtaaaatgGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTc                                                                                       >  2:805075/1‑141 (MQ=255)
                                                     taaaatgGGAAGATAAACAGTGCTGCTAAATTCAATTTCACATATAAAAAGAACATTATGGATAGTGGAATAATAATTGTTATATAATCTAGCCAAATATTCGTTAACCCGGTTAAAATGAAATCTTGTGTTTGTTCAAca                                                                                   <  2:371918/141‑1 (MQ=255)
                                                     taaaatgGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAAca                                                                                   <  2:620866/141‑1 (MQ=255)
                                                     gaaaatgGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAAca                                                                                   <  2:371704/140‑1 (MQ=255)
                                                         atgGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAAcatcat                                                                               <  2:304466/141‑1 (MQ=255)
                                                           gggAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTa                                                                             <  1:36133/141‑1 (MQ=255)
                                                           gggAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTa                                                                             <  2:281956/141‑1 (MQ=255)
                                                              aaGATAAACAGTGCTGCTAAAGTCAATTCCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCACATATTCATTAACCTGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATc                                                                          <  2:377078/141‑1 (MQ=255)
                                                                gATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCAc                                                                        <  1:279891/141‑1 (MQ=255)
                                                                gATAAACAGTGCTGCTAAAGTCAATTCCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACGAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCAc                                                                        <  1:556867/141‑1 (MQ=255)
                                                                aaTAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCAc                                                                        <  2:523722/140‑1 (MQ=255)
                                                                aaTAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCAc                                                                        <  1:903359/140‑1 (MQ=255)
                                                                 aTAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAAACTTTTGTTTGTTCAACATCATTAATCAct                                                                       >  2:633833/1‑141 (MQ=255)
                                                                   aaaCAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCActct                                                                     >  2:800439/1‑141 (MQ=255)
                                                                    aaCAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTa                                                                    <  1:259415/141‑1 (MQ=255)
                                                                      cAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGa                                                                  <  1:728248/141‑1 (MQ=255)
                                                                         ggctgctAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGatat                                                               <  1:434725/140‑1 (MQ=255)
                                                                          gctgctAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATAtt                                                              <  2:77309/141‑1 (MQ=255)
                                                                          gctgctAAAGTCAATTTCACATCCAAAAATAACATTATAGACAGTGCAATAATAATTGGTATATAATCTAACCAAATATTCATTAACCCGGCTAAAATGAAATCTTTTTTTTGTTCAACATCATTAATCACTCTAGATAtt                                                              <  2:703823/141‑1 (MQ=255)
                                                                             gctAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTCAATATTACt                                                           >  2:155121/1‑141 (MQ=255)
                                                                               tAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTg                                                         <  2:528836/141‑1 (MQ=255)
                                                                                aaaGTCAACTTCACATCTAAAAAGAACATTATCGATAGTGCAATAAGAATTGTTATATAAGCTAACCAAATATTCAGTAACCCGGGGAAAATGAAATCTTTTGTTTGTTCAACATCATGAATCACGCTAGATATTACTTGa                                                        <  1:413394/141‑1 (MQ=255)
                                                                                  aGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGAcc                                                      <  1:326173/141‑1 (MQ=255)
                                                                                  aGTCAATTTAACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGAcc                                                      <  2:154419/141‑1 (MQ=255)
                                                                                   gTCAATTTCCCATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCt                                                     <  1:610079/141‑1 (MQ=255)
                                                                                   gTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCt                                                     <  1:668462/141‑1 (MQ=255)
                                                                                   gTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAAACAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTCATCACTCTCTATATTACTTGTCCt                                                     >  1:896810/1‑141 (MQ=255)
                                                                                    tCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCtt                                                    >  1:768243/1‑140 (MQ=255)
                                                                                    tCAATCCCACATCTAAAAAGAACATTATGGATAGTGCACTAATAATTGTTATATAATCTAACGAAATATTCATTAACCCGGCTAAAATGACAGCTTTTGTTGGTTCAACATCATTAATCACTCTAGATATTACTTGACCTa                                                    >  2:405776/1‑141 (MQ=255)
                                                                                      aaTTTCACATCTAAAAAGAACATTATGGATAGTTCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACt                                                  <  2:401303/141‑1 (MQ=255)
                                                                                        tttCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTg                                                <  1:224571/141‑1 (MQ=255)
                                                                                         ttCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGa                                               >  1:510894/1‑141 (MQ=255)
                                                                                          tCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGat                                              >  1:236255/1‑141 (MQ=255)
                                                                                          tCACATCTAAAAAGAACACTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGat                                              <  1:330451/141‑1 (MQ=255)
                                                                                                   aaaaGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCAt                                     <  1:805075/141‑1 (MQ=255)
                                                                                                          cATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAAtct                              <  1:874906/141‑1 (MQ=255)
                                                                                                          cATTAATGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAAtct                              <  2:737078/141‑1 (MQ=255)
                                                                                                           aTTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAAtctc                             >  1:847062/1‑141 (MQ=255)
                                                                                                           aTTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAAtctc                             >  1:299286/1‑141 (MQ=255)
                                                                                                               tGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCAc                         >  2:89921/1‑141 (MQ=255)
                                                                                                                ggATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACt                        >  1:57821/1‑141 (MQ=255)
                                                                                                                 gATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACtt                       >  1:550783/1‑141 (MQ=255)
                                                                                                                   tAGGGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTaa                     <  2:341138/141‑1 (MQ=255)
                                                                                                                                  tGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTAAAGCTTGTAAATGGtt      >  1:231830/1‑141 (MQ=255)
                                                                                                                                      atatAATGTAACCAAATAGTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATGACTTGACCTACTCGATTATTAGCATAAAATCTCGCACTTAAAGCTTGTAACTGGTTGtat  <  2:592602/141‑1 (MQ=255)
                                                                                                                                      atatAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTAAAGCTTGTAAATGGTTGtat  <  1:56917/141‑1 (MQ=255)
                                                                                                                                      ataGAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTCGTTTGTTCAACATCATTAATCACTCTAGATACTACTTGACCTACCTGATTATTAGCATAAAATCTCGCACTTAAAGCTTGTAAATGGTTGtat  <  1:889567/141‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                         
TTTCACGTGTCAATTTTCTTAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATACAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTAAAGCTTGTAAATGGTTGTAT  >  USA300TCH1516_ALE/2010226‑2010500

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