Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,713,371 C→T 6.8% A159V (GCA→GTA)  srmB → ATP‑dependent RNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,713,3710CT 6.8% 186.0 / 3.8 59A159V (GCA→GTA) srmBATP‑dependent RNA helicase
Reads supporting (aligned to +/- strand):  ref base C (28/27);  new base T (2/2);  total (30/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.01e-01

ATCACCGGCGGCGTAGCCTATATGAACCACGCGGAAGTGTTCAGCGAAAATCAGGACATCGTGGTCGCCACGACCGGACGTCTGCTGCAATACATAAAAGAAGAGAACTTCGATTGCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTAAACAGACCCTGCTCTTTTCGGCAACGCTGGAAGGCGATGCGATTCAGGACTTTGCCGAGCGTCTGCTGGAAGATC  >  NC_000913/2713223‑2713514
                                                                                                                                                    |                                                                                                                                               
aTCACCGGCGGCGTAGCCTATATGAACCACGCGGAAGTGTTCAGCGAAAATCAGGACATCGTGGTCGCCACGACCGGACGTCTGCTGCAATACATAAAAGAAGAGAACTTCGATTGCCGCGCGGTTGAAACGCTGATCCTCGACGAAGc                                                                                                                                                 >  1:155592/1‑149 (MQ=255)
               gCCTATATGAACCACGCGGAAGTGTTCAGCGAAAATCAGGACATCGTGGTCGCCACGACCGGACGTCTGCTGCAATACATAAAAGAAGAGAACTTCGATTGCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGa                                                                                                                                  >  2:433549/1‑149 (MQ=255)
                   atatGAACCACGCGGAAGTGTTCAGCGAAAATCAGGACATCGTGGTCGCCACGACCGGACGTCTGCTGCAATACATAAAAGAAGAGAACTTCGATTGCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATg                                                                                                                              <  4:388614/149‑1 (MQ=255)
                    tatGAACCACGCGGAAGTGTTCAGCGAAAATCAGGACATCGTGGTCGCCACGACCGGACGTCTGCTGCAATACATAAAAGAAGAGAACTTCGATTGCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATgg                                                                                                                             <  3:265764/149‑1 (MQ=255)
                    tatGAACCACGCGGAAGTGTTCAGCGAAAATCAGGACATCGTGGTCGCCACGACCGGACGTCTGCTGCAATACATAAAAGAAGAGAACTTCGATTGCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATgg                                                                                                                             >  1:281017/1‑149 (MQ=255)
                      cgAACCACGCGGAAGTGTTCAGCGAAAATCAGGACATCGTGGTCGCCACGACCGGACGTCTGCTGCAATACATAAAAGAAGAGAACTTCGATTGCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGc                                                                                                                           <  1:3119/148‑1 (MQ=255)
                       gAACCACGCGGAAGTGTTCAGCGAAAATCAGGACATCGTGGTCGCCACGACCGGACGTCTGCTGCAATACATAAAAGAAGAGAACTTCGATTGCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCt                                                                                                                          >  1:290311/1‑149 (MQ=255)
                              gcgGAAGTGTTCAGCGAAAATCAGGACATCGTGGTCGCCACGACCGGACGTCTGCTGCAATACATAAAAGAAGAGAACTTCGATTGCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCa                                                                                                                   <  1:194794/149‑1 (MQ=255)
                                ggAAGTGTTCAGCGAAAATCAGGACATCGTGGTCGCCACGACCGGACGTCTGCTGCAATACATAAAAGAAGAGAACTTCGATTGCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAgg                                                                                                                 <  1:368540/149‑1 (MQ=255)
                                      gtTCAGCGAAAATCAGGACATCGTGGTCGCCACGACCGGACGTCTGCTGCAATACATAAAAGAAGAGAACTTCGATTGCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCg                                                                                                           <  3:7724/149‑1 (MQ=255)
                                            cGAAAATCAGGACATCGTGGTCGCCACGACCGGACGTCTGCTGCAATACATAAAAGAAGAGAACTTCGATTGCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACata                                                                                                     <  2:515850/149‑1 (MQ=255)
                                             gAAAATCAGGACATCGTGGTCGCCACGACCGGACGTCTGCTGCAATACATAAAAGAAGAGAACTTCGATTGCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACatat                                                                                                    >  1:218492/1‑149 (MQ=255)
                                                 aTCAGGACATCGTGGTCGCCACGACCGGACGTCTGCTGCAATACATAAAAGAAGAGAACTTCGATTGCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCt                                                                                                >  3:20607/1‑149 (MQ=255)
                                                      gACATCGTGGTCGCCACGACCGGACGTCTGCTGCAATACATAAAAGAAGAGAACTTCGATTGCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGa                                                                                           >  4:345645/1‑149 (MQ=255)
                                                        cATCGTGGTCGCCACGACCGGACGTCTGCTGCAATACATAAAAGAAGAGAACTTCGATTGCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGaaa                                                                                         <  2:290311/149‑1 (MQ=255)
                                                         aTCGTGGTCGCCACGACCGGACGTCTGCTGCAATACATAAAAGAAGAGAACTTCGATTGCCGCGCGGTTGAAACGATGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAAc                                                                                        >  4:405028/1‑149 (MQ=255)
                                                          tCGTGGTCGCCACGACCGGACGTCTGCTGCAATACATAAAAGAAGAGAACTTCGATTGCCGCGCGGTTGAAACGCTGATCCTCGACGAAGTAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAAcg                                                                                       >  3:545806/1‑149 (MQ=255)
                                                          tCGTGGTCGCCACGACCGGACGTCTGCTGCAATACATAAAAGAAGAGAACTTCGATTGCCGCGCGGTTGAAACGCTGATCCTCGACGAAGTAGACCGTATGCAGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAAcg                                                                                       >  1:184695/1‑149 (MQ=255)
                                                            gTGGTCGCCACGACCGGACGTCTGCTGCAATACATAAAAGAAGAGAACTTCGATTGCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAAcgcg                                                                                     >  2:72831/1‑149 (MQ=255)
                                                                 cGCCACGACCGGACGTCTGCTGCAATACATAAAAGAAGAGAACTTCGATTGCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGc                                                                                <  1:514568/149‑1 (MQ=255)
                                                                   ccACGACCGGACGTCTGCTGCAATACATAAAAGAAGAGAACTTCGATTGCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGt                                                                              >  2:239567/1‑149 (MQ=255)
                                                                          cGGACGTCTGCTGCAATACATAAAAGAAGAGAACTTCGATTGCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTAAACAGa                                                                       <  4:171443/149‑1 (MQ=255)
                                                                          cGGACGTCTGCTGCAATACATAAAAGAAGAGAACTTCGATTGCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTAAACAGa                                                                       >  3:588001/1‑149 (MQ=255)
                                                                               gTCTGCTGCAATACATAAAAGAAGAGAACTTCGATTGCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTAAACAGACCCTg                                                                  >  4:129305/1‑149 (MQ=255)
                                                                                tCTGCTGCAATACATAAAAGAAGAGAACTTCGATTGCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTAAACAGACCCTGc                                                                 <  3:450556/149‑1 (MQ=255)
                                                                                    ctgcAATACATAAAAGAAGAGAACTTCGATTGCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTAAACa                                                                         >  3:272569/1‑137 (MQ=255)
                                                                                    ctgcAATACATAAAAGAAGAGAACTTCGATTGCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTAAACa                                                                         <  4:272569/137‑1 (MQ=255)
                                                                                    ctgcAATACATAAAAGAAGAGAACTTCGATTGCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTAAACa                                                                         <  2:623136/137‑1 (MQ=255)
                                                                                    ctgcAATACATAAAAGAAGAGAACTTCGATTGCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTAAACa                                                                         >  1:623136/1‑137 (MQ=255)
                                                                                         aaaCATAAAAGAAGAAAACTTCGATTGCCGCGCGGTTGAAACGCTGATCCTCGACGGAGCAGACCGTATGCTGGACATCGGCTTCGCTCAGGCAAGCGAACATATTGCTCACGAAACGCGCTGGCGTAAACAGACCGTGCTCTTTTCgg                                                        >  4:446958/3‑149 (MQ=255)
                                                                                          tACATAAAAGAAGAGAACTTCGATTGCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTAAACAGACCCTGCTCTTTTCGGc                                                       >  3:250704/1‑149 (MQ=255)
                                                                                               aaaaGAAGAGAACTTCGATTGCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGc                                                                                >  2:624250/1‑119 (MQ=255)
                                                                                               aaaaGAAGAGAACTTCGATTGCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGc                                                                                <  1:624250/119‑1 (MQ=255)
                                                                                                 aagaagAGAACTTCGATTGCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTAAACAGACCCTGCTCTTTTCGGCAACGCTg                                                >  1:135147/1‑149 (MQ=255)
                                                                                                  agaagaGAACTTCGATTGCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTAAACAGACCCTGCTCTTTTCGGCAACGCTgg                                               >  2:586327/1‑149 (MQ=255)
                                                                                                     agagAACTTCGATTGCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTAAACAGACCCTGCTCTTTTCGGCAACGCTGGAAg                                            <  4:20607/149‑1 (MQ=255)
                                                                                                       agaACTTCGATTGCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTAAACAGACCCTGCTCTTTTCGGCAACGCTGGAAGgc                                          >  1:25227/1‑149 (MQ=255)
                                                                                                            ttCGATTGCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATc                                                                                                            >  2:95577/1‑78 (MQ=255)
                                                                                                            ttCGATTGCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATc                                                                                                            <  1:95577/78‑1 (MQ=255)
                                                                                                            ttCGATTGCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGaa                                                                                                         >  1:551412/1‑81 (MQ=255)
                                                                                                            ttCGATTGCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGaa                                                                                                         <  2:551412/81‑1 (MQ=255)
                                                                                                            ttCGATTGCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTAAACAGACCCTGCTCTTTTCGGCAACGCTGGAAGgcgatgc                                     >  1:38792/1‑149 (MQ=255)
                                                                                                              cGATTGCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGATGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTAAACAGACCCTGCTCTTTTCGGCAACGCTGGAAGgcgatgcga                                   <  2:187240/149‑1 (MQ=255)
                                                                                                                aTTGCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTAAACAGACCCTGCTCTTTTCGGCAACGCTGGAAGGCGATGCGAtt                                 <  2:155592/149‑1 (MQ=255)
                                                                                                                 ttGCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTAAACAGACCCTGCTCTTTTCGGCAACTCTGGAAGGCGATGCGATTc                                >  4:422605/1‑149 (MQ=255)
                                                                                                                 ttGCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTAAACAGACCCTGCTCTTTTCGGCAACGCTGGAAGGCGATGCGATTc                                >  4:80885/1‑149 (MQ=255)
                                                                                                                   gCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTAAACAGACCCTGCTCTTTTCGGCAACGCTGGAAGGCGATGCGATTCAg                              <  2:218492/149‑1 (MQ=255)
                                                                                                                       cgcgGTTGAAACGCTGATCCTCGACGAAGTAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTAAACAGACCCTGCTCTTTTCGGCAACGCTGGAAGGCGATGCGATTCAGGACt                          <  2:184695/149‑1 (MQ=255)
                                                                                                                       cgcgGTTGAAACGCTGATCCTCGACGAAGTAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTAAACAGACCCTGCTCTTTTCGGCAACGCTGGAAGGCGATGCGATTCAGGACt                          <  4:545806/149‑1 (MQ=255)
                                                                                                                         cgGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTAAACAGACCCTGCTCTTTTCGGCAACGCTGGAAGGCGATGCGATTCAGGACttt                        <  2:135147/149‑1 (MQ=255)
                                                                                                                                aaCGCTGATCCTCGACGAAGCAGAGCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTAAACAGACCCTGCTCTTTTCGGCAACGCTGGAAGGCGATGCGATTCAGGACTTTGCCGAGc                 <  2:25227/149‑1 (MQ=255)
                                                                                                                                aaCGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTAAACAGACCCTGCTCTTTTCGGCAACGCTGGAAGgcgat                                       >  1:414896/1‑127 (MQ=255)
                                                                                                                                aaCGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTAAACAGACCCTGCTCTTTTCGGCAACGCTGGAAGgcgat                                       <  2:414896/127‑1 (MQ=255)
                                                                                                                                   gCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTAAACAGACCCTGCTCTTTTCGGCAACGCTGGAAGGCGATGCGATTCAGGACTTTGCCGAGCGTc              >  4:305856/1‑149 (MQ=255)
                                                                                                                                     tGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTAAACAGACCCTGCTCTTTTCGGCAACGCTGGAAGGCGATGCGATTCAGGACTTTGCCGAGCGTctg            >  3:306815/1‑149 (MQ=255)
                                                                                                                                     tGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTAAACAGACCCTGCTCTTTTCGGCAACGCTGGAAGGCGATGCGATTCAGGACTTTGCCGAGCGTctg            >  4:329799/1‑149 (MQ=255)
                                                                                                                                       aTCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTAAACAGACCCTGCTCTTTTCGGCAACGCTGGAAGGCGATGCGATTCAGGACTTTGCCGAGCGTctgct          <  4:514473/149‑1 (MQ=255)
                                                                                                                                           tCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTAAACAGACCCTGCTCTTTTCGGCAACGCTGGAAGGCGATGCGATTCAGGACTTTGCCGAGCGTCTGCTGGaa      <  3:305856/149‑1 (MQ=255)
                                                                                                                                             gacgaAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTAAACAGACCCTGCTCTTTTCGGCAACGCTGGAAGGCGATGCGATTCAGGACTTTGCCGAGCGTCTGCTGGAAGa    <  3:129305/149‑1 (MQ=255)
                                                                                                                                               cgaAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTAAACAGACCCTGCTCTTTTCGGCAACGCGGGAAGGCGATGCGATTCAGGACTTTGCCGAGCGTCTGCTGGAAGATc  <  3:422605/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
ATCACCGGCGGCGTAGCCTATATGAACCACGCGGAAGTGTTCAGCGAAAATCAGGACATCGTGGTCGCCACGACCGGACGTCTGCTGCAATACATAAAAGAAGAGAACTTCGATTGCCGCGCGGTTGAAACGCTGATCCTCGACGAAGCAGACCGTATGCTGGATATGGGCTTCGCTCAGGATATCGAACATATTGCTGGCGAAACGCGCTGGCGTAAACAGACCCTGCTCTTTTCGGCAACGCTGGAAGGCGATGCGATTCAGGACTTTGCCGAGCGTCTGCTGGAAGATC  >  NC_000913/2713223‑2713514

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