Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,799,177 C→T 16.9% intergenic (‑191/+25) PP_2452 ← / ← ansB hypothetical protein/glutaminase‑asparaginase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,799,1770CT16.9% 116.2 / 14.5 53intergenic (‑191/+25)PP_2452/ansBhypothetical protein/glutaminase‑asparaginase
Reads supporting (aligned to +/- strand):  ref base C (20/24);  new base T (4/5);  total (24/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 = 5.93e-01

CTGTAAGTAAATTCGCAACAGCTTGTCAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTC  >  NC_002947/2799029‑2799324
                                                                                                                                                    |                                                                                                                                                   
cTGTAATTAAATTCGCTACTGCTTGTCAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGTCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGgc                                                                                                                                                     <  4:86837/149‑1 (MQ=255)
  gtaagtaaATTCGCAACAGCTTGTCAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGgcg                                                                                                                                                    >  3:122160/1‑148 (MQ=255)
   taagtaaATTCGCAACAGCTTGTCAGGAGAAATCCAGCGTTGGCCAGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCa                                                                                                                                                  >  1:243277/1‑149 (MQ=255)
       taaATTCGCAACAGCTTGTCAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCgg                                                                                                                                              >  5:340292/1‑149 (MQ=255)
           ttCGCAACAGCTTGTCAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCt                                                                                                                                          >  5:14636/1‑149 (MQ=255)
              gCAACAGCTTGTCAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGACGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAg                                                                                                                                       >  2:122955/1‑149 (MQ=255)
                aaCAGCTTGTCAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGGCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGGGCCGGTGCGCCGGCTCTGAggg                                                                                                                                     <  2:105165/149‑1 (MQ=255)
                 actgCTTGTCAGGTGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAATGGTCAACGATGTCGCCGGTTTGTAACCTCTGTGGCTACGCGGTCGAGGAATTGCCAACGAGGCCCTCAGGGCCGGTGCGCCGGCTCTGAgggc                                                                                                                                    <  1:229374/146‑2 (MQ=255)
                    gCTTGTCAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGaa                                                                                                                                 <  6:340292/149‑1 (MQ=255)
                    gCTTGTCAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGaa                                                                                                                                 <  2:243277/149‑1 (MQ=255)
                    gCTTGTCAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGaa                                                                                                                                 <  7:266532/149‑1 (MQ=255)
                     cTTGTCAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGa                                                                                                                                        >  7:24617/1‑141 (MQ=255)
                     cTTGTCAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGa                                                                                                                                        <  8:24617/141‑1 (MQ=255)
                     cTTGTCAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAg                                                                                                                                >  4:128158/1‑149 (MQ=255)
                      ttGTCAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAg                                                                                                                                >  5:105737/1‑148 (MQ=255)
                       tGTCAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCa                                                                                                                              >  7:54033/1‑149 (MQ=255)
                         tCAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGACGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACt                                                                                                                            <  1:122955/149‑1 (MQ=255)
                           aGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAAGCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCa                                                                                                                          <  1:109678/149‑1 (MQ=255)
                           aGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCa                                                                                                                          <  1:105315/149‑1 (MQ=255)
                                     cAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCag                                                                                                                <  8:54033/149‑1 (MQ=255)
                                     cAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCa                                                                                                                 <  6:14636/148‑1 (MQ=255)
                                            ggCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATc                                                                                                          <  4:122160/148‑1 (MQ=255)
                                                  tAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTg                                                                                                     >  3:215854/1‑147 (MQ=255)
                                                    aaCCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAg                                                                                                  <  6:105737/148‑1 (MQ=255)
                                                         tGAGGCCGGACAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCttt                                                                                            >  1:286722/1‑149 (MQ=255)
                                                           aGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGATTTGTCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTg                                                                                           <  1:76577/148‑1 (MQ=255)
                                                                ggTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTcc                                                                                     >  6:141643/1‑149 (MQ=255)
                                                                  tCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTg                                                                                           <  6:120347/141‑1 (MQ=255)
                                                                  tCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTg                                                                                           >  5:120347/1‑141 (MQ=255)
                                                                    aGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTgg                                                                                  <  8:182919/148‑1 (MQ=255)
                                                                        aGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCtt                                                                             <  2:286722/149‑1 (MQ=255)
                                                                        aGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCtt                                                                             >  6:278595/1‑149 (MQ=255)
                                                                           tAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTgg                                                                           <  4:215854/148‑1 (MQ=255)
                                                                            aaaCGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGt                                                                          >  4:159985/1‑148 (MQ=255)
                                                                                      gCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCa                                                               >  5:235261/1‑149 (MQ=255)
                                                                                      gCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCa                                                               >  7:260908/1‑149 (MQ=255)
                                                                                                  ccTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTc                                                                                      <  4:228149/114‑1 (MQ=255)
                                                                                                  ccTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTc                                                                                      >  3:228149/1‑114 (MQ=255)
                                                                                                              cgcgGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCg                                       <  3:128158/149‑1 (MQ=255)
                                                                                                                        gAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTc                             <  7:285808/149‑1 (MQ=255)
                                                                                                                                aaCGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGt                           >  2:209108/1‑143 (MQ=255)
                                                                                                                                aaCGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGt                           <  1:209108/143‑1 (MQ=255)
                                                                                                                                     ggCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGGCTCACGAcc                <  6:235261/149‑1 (MQ=255)
                                                                                                                                        ccTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAg             >  8:107195/1‑149 (MQ=255)
                                                                                                                                          tCAGGGCCGGTGCGCCGGCTCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGt                                                                          >  7:104097/1‑86 (MQ=38)
                                                                                                                                          tCAGGGCCGGTGCGCCGGCTCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGt                                                                          <  8:104097/86‑1 (MQ=38)
                                                                                                                                          tCAGGGCCGGTGCGCCGGCTCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGa                         <  2:173954/135‑1 (MQ=255)
                                                                                                                                          tCAGGGCCGGTGCGCCGGCTCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGa                         >  1:173954/1‑135 (MQ=255)
                                                                                                                                          tCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGc                                        >  2:211676/1‑120 (MQ=255)
                                                                                                                                          tCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGc                                        <  1:211676/120‑1 (MQ=255)
                                                                                                                                            agggCCGGTGCGCCGGCTCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGt         >  1:2451/4‑149 (MQ=255)
                                                                                                                                            agagCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTc           <  5:141643/147‑1 (MQ=255)
                                                                                                                                            agagCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCg          <  5:278595/148‑1 (MQ=255)
                                                                                                                                              ggCCGGTGCGCCGGCTCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGt            >  1:169845/2‑144 (MQ=255)
                                                                                                                                              ggCCGGTGCGCCGGCTCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGt            <  2:169845/143‑1 (MQ=255)
                                                                                                                                                   gcgcACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTc  <  3:159985/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CTGTAAGTAAATTCGCAACAGCTTGTCAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTC  >  NC_002947/2799029‑2799324

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