Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,201,428:1 +C 7.0% coding (756/1368 nt) PP_2808 → oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,201,4281.C 7.0% 218.0 / 11.5 57coding (756/1368 nt)PP_2808oxidoreductase
Reads supporting (aligned to +/- strand):  ref base . (23/30);  new base C (2/2);  total (25/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.39e-01

CCATCGAGGCGGGGGTGAAGCTCTACGAACGCACGCCGTTGCTGAGCTTCAGCGAAGGGCCGGTCATTACCTGCAAGACTGCCCATGGCCGTGCAACCGCCCCGGTGATGGTGCTGGCAACCAACGCCTATACCCCGCAGCTGGGGGTGTTGCGCAACAAGGTCGCGGCGATTCGGGTGTCGGCAATCGAAACCCAAACGCTGTCGACGCAGCAGTTGGCTGCCCTCGGCTGGCAAGGGCGGGAAGGGATCATCACGCCGCACTACACCATGGAGAGCCACCGCCTGAC  >  NC_002947/3201286‑3201574
                                                                                                                                               |                                                                                                                                                  
ccATCGAGGCGGGGGTGAAGCTCTACGAACGCACGCCGTTGCTGAGCTTCAGCGAAGGGCCGGTCATTACCTGCAAGACTGCCCATGGCCGTGCAACCGCCCCGGTGATGGTGCTGGCAACCAACGCCTATACCCCGCAGCTGGGGgtg                                                                                                                                              >  3:200105/1‑149 (MQ=255)
   tCGAGGCGGGGGTGAAGCTCTACGAACGCACGCCGTTGCTGAGCTTCAGCGAAGGGCCGGTCATTACCTGCAAGACTGCCCATGGCCGTGCAACCGCCCCGGTGATGGTGCTGGCAACCAACGCCTATACCCCGCAGCTGGGGGTGTTg                                                                                                                                           <  4:200105/149‑1 (MQ=255)
   tCGAGGCGGGGGGGAAGCTCTACGAACGCACGCCGTTGCTGAGCTTCAGCGAAGGGCCGGTCATTACCTGCAAGACTGCCCATGGCCGTGCAACCGCCCCGGTGATGGTGCTGGCAACCAACGCCTATACCCCGCAGCTGGGGGTGTTg                                                                                                                                           <  5:96244/149‑1 (MQ=255)
        gCGGGGGGGAAGCTCTACGAACGCACGCCGTTGCTGAGCTTCAGCGAAGGGCCGGTCATTACCTGCAAGACTGCCCATGGCCGGGCAACCGCCCCGGTGATGGTGCTGGCAACCAACGCCTATACCCCGCAGCTGGGGGTGTTGCGcaa                                                                                                                                      <  8:202659/149‑1 (MQ=255)
            gggTGAAGCTCTACGAACGCACGCCGTTGCTGAGCTTCAGCGAAGGGCCGGTCATTACCTGCAAGACTGCCCATGGCCGTGCAACCGCCCCGGTGATGGTGCTGGCAACCAACGCCTATACCCCGCAGCTGGGGGTGTTGCGCAACAAg                                                                                                                                  <  8:162451/149‑1 (MQ=255)
            gggTGAAGCTCTACGAACGCACGCCGTTGCTGAGCTTCAGCGAAGGGCCGGTCATTACCTGCAAGACTGCCCATGGCCGTGCAACCGCCCCGGTGATGGTGCTGGCAACCAACGCCTATACCCCGCAGCTGGGGGTGTTGCGCAACAAg                                                                                                                                  <  5:165251/149‑1 (MQ=255)
                 aaGCTCTACGAACGCACGCCGTTGCTGAGCTTCAGCGAAGGGCCGGTCATTACCTGCAAGACTGCCCATGGCCGTGCAACCGCCCCGGTGATGGTGCTGGCAACCAACGCCTATACCCCGCAGCTGGGGGTGTTGCGCAACAAGGTcgc                                                                                                                             <  2:14524/149‑1 (MQ=255)
                   gCTCTACGAACGCACGCCGTTGCTGAGCTTCAGCGAAGGGCCGGTCATTACCTGCAAGACTGCCCATGGCCGTGCAACCGCCCCGGTGATGGTGCTGGCAACCAACGCCTATACCCCGCAGCTGGGGGTGTTGCGCAACAAGGTCgcgg                                                                                                                           >  3:224502/1‑149 (MQ=255)
                     tctACGAACGCACGCCGTTGCTGAGCTTCAGCGAAGGGCCGGTCATTACCTGCAAGACTGCCCATGGCCGTGCAACCGCCCCGGTGATGGTGCTGGCAACCAACGCCTATACCCCGCAGCTGGGGGTGTTGCGCAACAAGGTCgcggc                                                                                                                          >  4:316977/1‑148 (MQ=255)
                           aaCGCACGCCGTTGCTGAGCTTCAGCGAAGGGCCGGTCATTACCTGCAAGACTGCCCATGGCCGTGCAACCGCCCCGGTGATGGTGCTGGCAACCAACGCCTATACCCCGCAGCTGGGGGTGTTGCGCAACAAGGTCGCGGCGATTCgg                                                                                                                   <  2:16659/149‑1 (MQ=255)
                             cgcacgcCGTTGCTGAGCTTCAGCGAAGGGCCGGTCATTACCTGCAAGACTGCCCATGGCCGTGCAACCGCCCCGGTGATGGTGCTGGCAACCAACGCCTATACCCCGCAGCTGGGGGTGTTGCGCAACAAGGTCGCGGCGATTCGGgt                                                                                                                 <  8:105532/149‑1 (MQ=255)
                                acgcCGTTGCTGAGCTTCAGCGAAGGGCCGGTCATTACCTGCAAGACTGCCCATGGCCGTGCAACCGCCCCGGTGATGGTGCTGGCAACCAACGCCTATACCCCGCAGCTgg                                                                                                                                                   >  4:307998/1‑112 (MQ=255)
                                acgcCGTTGCTGAGCTTCAGCGAAGGGCCGGTCATTACCTGCAAGACTGCCCATGGCCGTGCAACCGCCCCGGTGATGGTGCTGGCAACCAACGCCTATACCCCGCAGCTgg                                                                                                                                                   <  3:307998/112‑1 (MQ=255)
                                     gTTGCTGAGGTTCAGCGAAGGGCCGGTCATTACCTGCAAGACTGCCCATGGCCGTTCAACCGCCCCGGTGATGGTGCTGTCAACCAACGCCTATACCCCGGAGCTGGGGGTGTTGCGCAACAAGGTCGCGGCGATTCGGGTGTCGGCaa                                                                                                         <  7:49722/149‑1 (MQ=255)
                                        gCTGAGCTTCAGCGAAGGGCCGGTCATTACCTGCAAGACTGCCCATGGCCGTGCAACCGCCCCGGTGATGGTGCTGGCAACCAACGCCTATACCCCGCAGCTGGGGGTGTTGCGCAACAAGGTCGCGGCGATTCGGGTGTCGGCAATc                                                                                                       <  4:224502/148‑1 (MQ=255)
                                                  ggCGAAGGGCCGGTCATTACCTGCAAGACTGCCCATGGCCGTGCAACCGCCCCGGTGATGGTGCTGGCAACCAACGCCTATACCCCGCAGCTGGGGGTGTTGCGCAACAAGGTCGCGGCGATTCGGGTGTCgg                                                                                                            >  7:164608/2‑133 (MQ=255)
                                                  ggCGAAGGGCCGGTCATTACCTGCAAGACTGCCCATGGCCGTGCAACCGCCCCGGTGATGGTGCTGGCAACCAACGCCTATACCCCGCAGCTGGGGGTGTTGCGCAACAAGGTCGCGGCGATTCGGGTGTCgg                                                                                                            <  8:164608/132‑1 (MQ=255)
                                                  aGCGACGGGCCGGTCATTACCTGCAAGACTGCCCATGGCCGTGCAACCGCCCCGGTGATGGTGCTGGCAACCAACGCCTATACCCCTCAGCTGGGGGTGTTGCACAACAAGGGCGCGGCGATTCGGGTGTCGGCAATCGAAACCCaaa                                                                                             <  3:257122/148‑1 (MQ=255)
                                                  aGCGAAGGGCCGGTCATTACCTGCAAGACTGCCCATGGCCGTGCAACCGCCCCGGTGATGGTGCTGGCAACCAACGCCTATACCCCGCAGCTGGGGGTGTTGCGCAACAAGGTCACGGCGATTCGGGTGTCGGCAATCGAAACCCAAAc                                                                                            <  3:316977/149‑1 (MQ=255)
                                                       aGGGCCGGTCATTACCTGCAAGACTGCCATGGCCGTGCAACCGCCCCGGTGATGTTGCTGGCAACCAACGCCTATACCCCGCAGCTGGGGGTGTTGCGCAACAAGGTCGCGGCGATTCGGGTGTCGGCAATCGAAACCCAAACGCTGTc                                                                                      <  2:283945/149‑1 (MQ=255)
                                                         ggCCGGTCATTACCTGCAAGACTGCCCATGGCCGTGCAACCGCCCCGGTGATGGTGCTGGCAACCAACGCCTATACCCCGCAGCTGGGGGTGTTGCGCAACAAGGTCGCGGCGATTCGGGTGTCGGCAATCGAAAc                                                                                                  <  3:172792/136‑1 (MQ=255)
                                                         ggCCGGTCATTACCTGCAAGACTGCCCATGGCCGTGCAACCGCCCCGGTGATGGTGCTGGCAACCAACGCCTATACCCCGCAGCTGGGGGTGTTGCGCAACAAGGTCGCGGCGATTCGGGTGTCGGCAATCGAAAc                                                                                                  >  4:172792/1‑136 (MQ=255)
                                                              gTCATTACCTGCAAGACTGCCCATGGCCGTGCAACCGCCCCGGTGATGGTGCTGGCAACCAACGCCTATACCCCGCAGCTGGGGGTGTTGCGCAACAAGGTCGCGGCGATTCGGGTGTCGGCAATCGAAACCCAAACGCTGTCGACgca                                                                                <  8:21570/149‑1 (MQ=255)
                                                                       tGCAAGACTGCCCATGGCCGTGCAACCGCCCCGGTGATGGTGCTGGCAACCAACGCCTATACCCCGCAGCTGGGGGTGTTGCGCAACAAGGTCGCGGCGATTCGGGTGTCGGCAAt                                                                                                        <  4:286628/116‑1 (MQ=255)
                                                                       tGCAAGACTGCCCATGGCCGTGCAACCGCCCCGGTGATGGTGCTGGCAACCAACGCCTATACCCCGCAGCTGGGGGTGTTGCGCAACAAGGTCGCGGCGATTCGGGTGTCGGCAAt                                                                                                        >  3:286628/1‑116 (MQ=255)
                                                                        gCAAGACTGCCCATGGCCGTGCAACCGCCCCGGTGATGGTGCTGGCAACCAACGCCTATACCCCGCAGCTGGGGGTGTTGCGCAACAAGGTCGCGGCGATTCGGGTGTCGGCAATCGAAACCCAAACGCTGTCGACGCAGCAGTTGGCt                                                                      >  7:74005/1‑149 (MQ=255)
                                                                            gACTGCCCATGGCCGTGCAACCGCCCCGGTGATGGTGCTGGCAACCAACGCCTATACCCCGCAGCTGGGGGTGTTGCGCAACAAGGTCGCGGCGATTCGGGTGTCGGCAATCGAAACCCAAACGCTGTCGACGCAGCAGTTGGCTGccc                                                                  >  8:111750/1‑149 (MQ=255)
                                                                                 cccATGGCCGTGCAACCGCCCCGGTGATGGTGCTGGCAACCAACGCCTATACCCCGCAGCTGGGGGTGTTGCGCAACAAGGTCGCGGCGATTCGGGTGTCGGCAATCGAAACCCAAACGCTGTCGACGCAGCAGTTGGCTGCCCTCGGc                                                             >  3:33365/1‑149 (MQ=255)
                                                                                  ccATGGCCGTGCAACCGCCCCGGTGATGGTGCTGGCAACCAACGCCTATACCCCGCAGCTGGGGGTGTTGCGCAACAAGGTCGCGGCGATTCGGGTGTCGGCAATCGAAACCCAAACGCTGTCGACg                                                                                  >  4:169017/1‑127 (MQ=255)
                                                                                  ccATGGCCGTGCAACCGCCCCGGTGATGGTGCTGGCAACCAACGCCTATACCCCGCAGCTGGGGGTGTTGCGCAACAAGGTCGCGGCGATTCGGGTGTCGGCAATCGAAACCCAAACGCTGTCGACg                                                                                  <  3:169017/127‑1 (MQ=255)
                                                                                  ccATGGCCGTGCAACCGCCCCGGTGATGGTGCTGGCAACCAACGCCTATACCCCGCAGCTGGGGGTGTTGCGCAACAAGGTCGCGGCGATTCGGGTGTCGGCAATCGAAACCCAAACGCTGTCGACGCAGCAGTTGGCTGCCCTCGGCt                                                            >  5:253095/1‑149 (MQ=255)
                                                                                  ccATGGCCGTGCAACCGCCCCGGTGATGGTGCTGGCAACCAACGCCTATACCCCGCAGCTGGGGGTGTTGCGCAACAAGGTCGCGGCGATTCGGGTGTCGGCAATCGAAACCCAAACGCTGTCGACGCAGCAGTTGGCTGCCCTCGGCt                                                            >  5:369288/1‑149 (MQ=255)
                                                                                   cATGGCCGTGCAACCGCCCCGGTGATGGTGCTGGCAACCAACGCCTATACCCCGCAGCTGGGGGTGTTGCGCAACAAGGTCGCGGCGATTCGGGTGTCGGCAATCGAAACCCAAACGCTGTCGACGCAGCAGTTGGCTGCCCTCGGCTg                                                           >  2:251257/1‑149 (MQ=255)
                                                                                    aTGGCCGTGCAACCGCCCCGGTGATGGTGCTGGCAACCAACGCCTATACCCCGCAGCTGGGGGTGTTGCGCAACAAGGTCGCGGCGATTCGGGTGTCGGCAATCGAAACCCAAACGCTGTCGACGCAGCAGTTGGCTGCCCTCGGCTgg                                                          >  5:364026/1‑149 (MQ=255)
                                                                                         cGTGCAACCGCCCCGGTGATGGTGCTGGCAACCAACGCCTATACCCCGCAGCTGGGGGTGTTGCGCAACAAGGTCGCGGCGATTCGGGTGTCGGCAATCGAAACCCAAACGCTGTCGACGCAGCAGTTGGCTGCCCTCGGCTGGCAAgg                                                     <  7:111750/149‑1 (MQ=255)
                                                                                                  gCCCCGGTGATGGTGCTGGCAACCAACGCCTATACCCCGCAGCTGGGGGTGTTGCGCAACAAGGTCGCGGCGATTCGGGTGTCGGCAATCGAAACCCAAACGCTGTCGACGCAGCAGTTGGCTGCCCTCGGc                                                             <  2:50315/132‑1 (MQ=255)
                                                                                                  gCCCCGGTGATGGTGCTGGCAACCAACGCCTATACCCCGCAGCTGGGGGTGTTGCGCAACAAGGTCGCGGCGATTCGGGTGTCGGCAATCGAAACCCAAACGCTGTCGACGCAGCAGTTGGCTGCCCTCGGc                                                             >  1:50315/1‑132 (MQ=255)
                                                                                                   ccccGGTGATGGTGCTGGCAACCAACGCCTATACCCCGCAGCTGGGGGTGTTGCGCAACAAGGTCGCGGCGATTCGGGTGTCGGCAATCGAAACCCAAACGCTGTCGACGCAGCAGTTGGCTGCCCTCGGCTGGCAAGGGCGGGAAggg                                           <  8:74005/149‑1 (MQ=255)
                                                                                                   ccccGGTGATGGTGCTGGCAACCAACGCCTATACCCCGCAGCTGGGGGTGTTGCGCAACAAGGTCGCGGCGATTCGGGTGTCGGCAATCGAAACCCAAACGCTGTCGACGCAGCAGTTGGCTGCCCTCGGCTGGCAAGGGCGGGAAggg                                           <  4:33365/149‑1 (MQ=255)
                                                                                                    cccGGTGATGGTGCTGGCAACCAACGCCTATACCCCGCAGCTGGGGGTGTTGCGCAACAAGGTCGCGGCGATTCGGGTGTCGGCAATCGAAACCCAAACGCTGTCGACGCAGCAGTTGGCTGCCCTCGGCTGGCAAGGGCGGGAAGGGa                                          <  6:364026/149‑1 (MQ=255)
                                                                                                              gTGCTGGCAACCAACGCCTATACCCCGCAGCTGGGGGTGTTGCGCAACAAGGTCGCGGCGATTCGGGTGTCGGCAATCGAAACCCAAACGCTGTCGACGCAGCAGTTGGCTGCCCTCGGCTGGCAAGGGCGGGAAGGGATCATCAcgc                                 <  1:251257/148‑1 (MQ=255)
                                                                                                                gCTGGCAACCAACGCCTATACCCCGCAGCTGGGGGTGTTGCGCAACAAGGTCGCGGCGATTCGGGTGTCGGCAATCGAAACCCAAACGCTGTCGACGCAGCAGTTGGCTGCCCTCGGCTGGCAAGGGCGGGAAGGGATCATCAcgccg                               >  1:33527/1‑148 (MQ=255)
                                                                                                                     caaccaacGCCTATACCCCGCAGCTGGGGGTGTTGCGCAACAAGGTCGCGGCGATTCGGGTGTCGGCAATCGAAACCCAAACGCTGTCGACGCAGCAGTTGGCTGCCCTCGGCTGGCAAGGGCGGGAAGGGATCATCACGCCGCACTa                          >  7:316061/1‑148 (MQ=255)
                                                                                                                                 ataCCCCGCAGCTGGGGGTGTTGCGCAACAAGGTCGCGGCGATTCGGGTGTCGGCAATCGAAACCCAAACGCTGTCGACGCAGCAGTTGGCTGCCCTCGGCTGGc                                                         >  6:270753/1‑105 (MQ=255)
                                                                                                                                 ataCCCCGCAGCTGGGGGTGTTGCGCAACAAGGTCGCGGCGATTCGGGTGTCGGCAATCGAAACCCAAACGCTGTCGACGCAGCAGTTGGCTGCCCTCGGCTGGc                                                         <  5:270753/105‑1 (MQ=255)
                                                                                                                                 ataCCCCGCAGCTGGGGGTGTTGCGCAACAAGGTCGCGGCGATTCGGGTGTCGGCAATCGAAACCCAAACGCTGTCGACGCAGCAGTTGGCTGCCCTCGGCTGGCAAGGGCGGGAAGGGATCATCAcgccg                               <  3:141309/131‑1 (MQ=255)
                                                                                                                                 ataCCCCGCAGCTGGGGGTGTTGCGCAACAAGGTCGCGGCGATTCGGGTGTCGGCAATCGAAACCCAAACGCTGTCGACGCAGCAGTTGGCTGCCCTCGGCTGGCAAGGGCGGGAAGGGATCATCAcgccg                               >  4:141309/1‑131 (MQ=255)
                                                                                                                                    accccCAGCTGCGGGGTGTTGCGCAACAAGGTCGCGGCGATTCGGGTGTCGGCAATCGAAACCCAAACGCTGTCGACGCAGCAGTTGGCTGCCCTCGGCTgg                                                          >  7:191289/2‑102 (MQ=255)
                                                                                                                                    accccCAGCTGCGGGGTGTTGCGCAACAAGGTCGCGGCGATTCGGGTGTCGGCAATCGAAACCCAAACGCTGTCGACGCAGCAGTTGGCTGCCCTCGGCTgg                                                          <  8:191289/101‑1 (MQ=255)
                                                                                                                                     ccccCAGCTGCGGGGTGTTGCGCAACAAGGTCGCGGCGATTCGGGTGTCGGCAATCGAAACCCAAACGCTGTCGACGCAGCAGTTGGCTGCCCTCGGCTGGCAAGGGCGGGAAGGGatc                                        <  7:261571/119‑1 (MQ=255)
                                                                                                                                     ccccCAGCTGCGGGGTGTTGCGCAACAAGGTCGCGGCGATTCGGGTGTCGGCAATCGAAACCCAAACGCTGTCGACGCAGCAGTTGGCTGCCCTCGGCTGGCAAGGGCGGGAAGGGatc                                        >  8:261571/1‑119 (MQ=255)
                                                                                                                                     cccGCAGCTGGGGGTGTTGCGCAACAAGGTCGCGGCGATTCGGGTGTCGGCAATCGAAACCCAAACGCTGTCGACGCAGCAGTTGGCTGCCCTCGGCTGGCAAg                                                      <  7:68866/104‑1 (MQ=255)
                                                                                                                                     cccGCAGCTGGGGGTGTTGCGCAACAAGGTCGCGGCGATTCGGGTGTCGGCAATCGAAACCCAAACGCTGTCGACGCAGCAGTTGGCTGCCCTCGGCTGGCAAg                                                      >  8:68866/1‑104 (MQ=255)
                                                                                                                                     cccGCAGCTGGGGGTGTTGCGCAACAAGGTCGCGGCGATTCGGGTGTCGGCAATCGAAACCCAAACGCTGTCGACGCAGCAGTTGGCTGCCCTCGGCTGGCAAGGGCGGGAAGGGATCATCACGCCGCACTACACCATGGAGAGCCAcc         >  8:120649/1‑149 (MQ=255)
                                                                                                                                     cccGCAGCTGGGGGTGTTGCGCAACAAGGTCGCGGCGATTCGGGTGTCGGCAATCGAAACCCAAACGCTGTCGACGCAGCAGTTGGCTGCCCTCGGCTGGCAAGGGCGGGAAGGGATCATCACGCCGCACTACACCATGGAGAGCCAcc         >  5:282621/1‑149 (MQ=255)
                                                                                                                                       cGCAGCTGGGGGTGTTGCGCAACAAGGTCGCGGCGATTCGGGTGTCGGCAATCGAAACCCAAACGCTGTCGACGCAGCAGTTGGCTGCCCTCGGCTGGCAAGGGCGGGAAGGGATCATCACGCCGCACTACACCATGGAGAGCCACCGc       <  2:33527/149‑1 (MQ=255)
                                                                                                                                       cGCAGCTGGGGGTGTTGCGCAACAAGGTCGCGGCGATTCGGGTGTCGGCAATCGAAACCCAAACGCTGTCGACGCAGCAGTTGGCTGCCCTCGGCTGGCAAGGGCGGGAAGGGATCATCACGCCGCACTACACCATGGAGAGCCACCGc       <  8:316061/149‑1 (MQ=255)
                                                                                                                                       cGCAGCTGGGGGTGTTGCGCAACAAGGTCGCGGCGATTCGGGTGTCGGCAATCGAAACCCAAACGCTGTCGACGCAGCAGTTGGCTGCCCTCGGCTGGCAAGGGCGGGAAGGGATCATCACGCCGCACTACACCATGGAGAGCCACCGc       <  6:369288/149‑1 (MQ=255)
                                                                                                                                            cTGGGGGTGTTGCGCAACAAGGTCGCGGCGATTCGGGTGTCGGCAATCGAAACCCAAACGCTGTCGACGCAGCAGTTGGCTGCCCTCGGCTGGCAAGGGCGGGAAGGGATCATCACGCCGCACTACACCATGGAGAGCCACCGCCTGac  >  8:222774/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
CCATCGAGGCGGGGGTGAAGCTCTACGAACGCACGCCGTTGCTGAGCTTCAGCGAAGGGCCGGTCATTACCTGCAAGACTGCCCATGGCCGTGCAACCGCCCCGGTGATGGTGCTGGCAACCAACGCCTATACCCCGCAGCTGGGGGTGTTGCGCAACAAGGTCGCGGCGATTCGGGTGTCGGCAATCGAAACCCAAACGCTGTCGACGCAGCAGTTGGCTGCCCTCGGCTGGCAAGGGCGGGAAGGGATCATCACGCCGCACTACACCATGGAGAGCCACCGCCTGAC  >  NC_002947/3201286‑3201574

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: