Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,271,674:1 +T 9.6% intergenic (+58/‑42) PP_t32 → / → PP_t33 tRNA‑Asp/tRNA‑Asp

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,271,6741.T 9.6% 189.3 / 16.7 52intergenic (+58/‑42)PP_t32/PP_t33tRNA‑Asp/tRNA‑Asp
Reads supporting (aligned to +/- strand):  ref base . (23/24);  new base T (3/2);  total (26/26)
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.94e-01

AGGAACGCACTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTCAAGGTCCCTAGAACACCTT  >  NC_002947/2271527‑2271821
                                                                                                                                                    |                                                                                                                                                    
agGAACGCACTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCATGTGCCTTgt                                                                                                                                                      <  5:230138/149‑2 (MQ=255)
     cGCACTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTGGGCTGa                                                                                                                                                <  1:233835/149‑1 (MQ=255)
       cACTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTGGGCTGATc                                                                                                                                              >  4:325292/1‑149 (MQ=255)
        aCTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCATTAACAGATGGGTGCGTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTGGGCTGATc                                                                                                                                              >  1:2869/1‑148 (MQ=255)
          tGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTggg                                                                                                                                                    >  8:47899/1‑140 (MQ=255)
          tGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTggg                                                                                                                                                    >  8:41886/1‑140 (MQ=255)
          tGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTggg                                                                                                                                                    <  7:47899/140‑1 (MQ=255)
          tGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTggg                                                                                                                                                    <  7:41886/140‑1 (MQ=255)
            cgcAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGGCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTGGGCTGATCCCAGt                                                                                                                                         <  4:122902/149‑1 (MQ=255)
              cAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTGCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTGGGCTGATCCCAGttt                                                                                                                                       <  7:94030/149‑1 (MQ=255)
               aGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTGGGCTGATCCCAGTTTc                                                                                                                                      <  2:90158/149‑1 (MQ=255)
                      gTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCgg                                                                                                                                <  3:325292/148‑1 (MQ=255)
                        tCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCGGgc                                                                                                                              <  6:75075/148‑1 (MQ=255)
                               gTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCt                                                                                                                      >  5:255240/1‑149 (MQ=255)
                                             cTGGCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGc                                                                                                        <  6:312999/149‑1 (MQ=255)
                                                tCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCg                                                                                                       <  6:255240/147‑1 (MQ=255)
                                                  aCGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTa                                                                                                    <  3:334411/148‑1 (MQ=255)
                                                    gCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGtt                                                                                                 <  5:256689/149‑1 (MQ=255)
                                                            tCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCg                                                                                                       <  7:308388/135‑1 (MQ=255)
                                                            tCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCg                                                                                                       >  8:308388/1‑135 (MQ=255)
                                                              gcgGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCATTAAAAGATCCTTGCAAACCAACCGTCTAAAACCACAAAACACATCCCTTTATCTGATAAAATTTTCAATCCGTCACAAGCATGAACGATACGCAGCGGTAGTTCAATCGGTTa                                                                                       >  6:137753/1‑149 (MQ=255)
                                                                      gAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATa                                                                                                              <  3:172186/118‑1 (MQ=255)
                                                                      gAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATa                                                                                                              >  4:172186/1‑118 (MQ=255)
                                                                            cGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGt                                                                         >  2:332439/1‑149 (MQ=255)
                                                                                    gcgcCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCgg                                                                 >  3:347392/1‑149 (MQ=255)
                                                                                     cgcCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCggg                                                               <  4:347392/149‑1 (MQ=255)
                                                                                     cgcCATTCTTCAGTATCACCTGGGCGCCTGGCACCTGTCTGAAGCCACAAGGCACATGCCTTGGGCTGATCCCAGTTTCATTCGGGCTGAATCCTGATCTATTCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCggg                                                               <  6:135979/149‑1 (MQ=25)
                                                                                        cATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGt                                                            >  5:215784/1‑149 (MQ=255)
                                                                                           ttttCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACAGGCCTGTCACGCCGGGGGTcgc                                                         >  1:228514/1‑149 (MQ=255)
                                                                                                cgcaaCAGATCGGTGCCTGGCACCGGGCTGAAGCCACAAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAATGATACGCAGCGGTAGTTCAGTCGGTTa                                                                                       <  6:348940/112‑1 (MQ=25)
                                                                                                     cAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGt                                               >  1:281581/1‑149 (MQ=255)
                                                                                                     cAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGt                                               >  1:286440/1‑149 (MQ=255)
                                                                                                        aTGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTccc                                            >  2:174392/1‑149 (MQ=255)
                                                                                                          gggTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGt                                          >  5:100718/1‑149 (MQ=255)
                                                                                                           ggTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTc                                         <  1:332439/149‑1 (MQ=255)
                                                                                                                       cGGTCTGAAGCCACAAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCgg                                                                              <  6:255250/101‑1 (MQ=255)
                                                                                                                               aGCCACAAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTc                     <  2:281581/149‑1 (MQ=255)
                                                                                                                               aGCCACAAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTc                     <  2:286440/149‑1 (MQ=255)
                                                                                                                                agcccaAGGCATGTGCCTTGTGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGttt                      >  3:151036/5‑148 (MQ=14)
                                                                                                                                 gcccaAGGCATGTGCCTTGTGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGtt                                                    <  5:260048/114‑1 (MQ=17)
                                                                                                                                 gcccaAGGCATGTGCCTTGTGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGtt                                                    >  6:260048/4‑117 (MQ=17)
                                                                                                                                 ccACAAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGCGGGTCGCGGGTTCGAGTCCCGTCCGCTgcg                                  <  1:225757/135‑1 (MQ=255)
                                                                                                                                 ccACAAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTgcg                                  >  2:225757/1‑134 (MQ=255)
                                                                                                                                  cccaAGGCATGTGCCTTGTGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTc                     >  6:144136/3‑147 (MQ=255)
                                                                                                                                  cacaAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTCAAg                  >  5:53560/1‑149 (MQ=255)
                                                                                                                                   acaAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTCAAgg                 >  3:46372/1‑149 (MQ=255)
                                                                                                                                   acaAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTCAAgg                 >  4:212397/1‑149 (MQ=255)
                                                                                                                                       ggcatgtGCCTTGTGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTCAAGGTcc              <  4:151036/143‑1 (MQ=255)
                                                                                                                                       ggCACATTCCTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACTATACGCTGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTCAAGGTccc             <  2:2869/149‑1 (MQ=255)
                                                                                                                                          acaTGCCTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTCAAGGTccc             >  2:233904/1‑146 (MQ=255)
                                                                                                                                          acaTGCCTTGCGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTAAGAATACCGGCCTGTCACGCCGGGGGTCGCGCGTACGAGTCCCGTCCGCTGCGCCATCTTCGTATCAAGGTccc             >  6:250803/1‑146 (MQ=17)
                                                                                                                                               ccTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTCAAGGTCCCTAGAacac     <  6:53560/149‑1 (MQ=255)
                                                                                                                                                cTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGGCTGAACGATACGGAGCGGTAGTTCAGTGGTTTAGAATACGGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCGGTCCGCTGCGCCATCCTCGTTTCAAGGTCCCTAGAACAcc    <  5:250803/149‑1 (MQ=17)
                                                                                                                                                cTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGa                                                 >  6:232842/1‑104 (MQ=255)
                                                                                                                                                cTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGa                                                 <  5:232842/104‑1 (MQ=255)
                                                                                                                                                 ttGGGCTGATCCCAGTTTCAATCGAGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTCAAGGTCCCTAGAACACCt   >  1:230300/1‑149 (MQ=255)
                                                                                                                                                  tGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTCAAGGTCCCTAGAACACCtt  >  7:301660/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AGGAACGCACTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTCAAGGTCCCTAGAACACCTT  >  NC_002947/2271527‑2271821

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