Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,256,893 +C coding (405/423 nt) yhaQ ← putative L‑serine dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,256,8911.C97.4% 298.8 / 2.6 115coding (407/423 nt)yhaQputative L‑serine dehydratase
Reads supporting (aligned to +/- strand):  ref base . (2/1);  new base C (56/56);  total (58/57)
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.35e-01

Only 100 of 122 total aligned reads displayed.

CAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTT  >  NC_000913/3256816‑3256961
                                                                             |                                                                       
cAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGt                                                                          >  1:4569987/1‑76 (MQ=255)
cAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGt                                                                          >  1:1836498/1‑76 (MQ=255)
 aGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTc                                                                        >  2:635033/1‑76 (MQ=255)
  ggCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTcc                                                                       >  2:3071397/1‑76 (MQ=255)
   gCCGAAGTGTTCTTCTTCGACAATAAGCCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTccc                                                                      >  1:1374553/1‑75 (MQ=255)
   gCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTccc                                                                      >  2:4191055/1‑75 (MQ=255)
    ccGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTccct                                                                     <  2:2859294/76‑3 (MQ=255)
         gtgtTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAg                                                                <  2:1912314/76‑1 (MQ=255)
         gtgtTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAg                                                                >  2:3611270/1‑76 (MQ=255)
          tgtTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAgg                                                               <  1:2897694/76‑1 (MQ=255)
           gtTCTTCTTCGCCAAGAAACACGCCGCCGACAGAGTAATAGGTTTTACTTACAAGCTCGTCCTGGCCCTTCCAGGc                                                              >  1:4514014/1‑76 (MQ=255)
           gtTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGc                                                              <  1:3548270/76‑1 (MQ=255)
           gtTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGc                                                              >  1:3347788/1‑76 (MQ=255)
            ttcttcttcGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCa                                                             <  1:4018596/76‑1 (MQ=255)
            ttcttcttcGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCa                                                             >  1:4704630/1‑76 (MQ=255)
              cttcttcGACAGTAACCCCTCCGCCGACAGAGTGATAGGTTTTACTTGATGGCTCTTCCTGTCCGTTCCAGGCAGt                                                           <  2:860986/76‑1 (MQ=255)
              cttcttcGACAATAAACCCGCCGCCGCCAGAGTAATAGGTTTTACTTAATACCTCTTCCTTTCCCTTCCCGGcccg                                                           >  2:58796/1‑73 (MQ=255)
               ttcttcGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTg                                                          <  2:1851636/76‑1 (MQ=255)
               ttcttcGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTg                                                          >  1:3223267/1‑76 (MQ=255)
                tcttcGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGa                                                         <  1:266640/76‑1 (MQ=255)
                 cttcGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGAt                                                        >  1:3587318/1‑76 (MQ=255)
                  ttcGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGacc                                                       <  1:2321941/76‑3 (MQ=255)
                  ttcGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATc                                                       >  1:2544044/1‑76 (MQ=255)
                  ttcGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATc                                                       <  1:2962644/76‑1 (MQ=255)
                   tcGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATcc                                                      >  1:707155/1‑76 (MQ=255)
                     gACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGc                                                    >  1:3202351/1‑76 (MQ=255)
                     gACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGc                                                    >  2:4013063/1‑76 (MQ=255)
                        aaTAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCAtt                                                 >  2:3158995/1‑76 (MQ=255)
                        aaTAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCAtt                                                 <  1:1237442/76‑1 (MQ=255)
                         aTAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCCTTc                                                >  2:233837/1‑76 (MQ=255)
                         aTAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTc                                                <  2:1631194/76‑1 (MQ=255)
                          tAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTcc                                               >  1:2688780/1‑76 (MQ=255)
                             aCCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGtt                                            <  2:1930137/76‑1 (MQ=255)
                              cccGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTtc                                           <  1:3503436/76‑1 (MQ=255)
                               ccgccgccgACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTtct                                          >  2:4055961/1‑76 (MQ=255)
                                cgccgccgACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTtctc                                         >  1:3529156/1‑76 (MQ=255)
                                 gccgccgACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCa                                        >  1:2417090/1‑76 (MQ=255)
                                 gccgccgACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCa                                        >  1:1942806/1‑76 (MQ=255)
                                 gccgccgACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCa                                        <  1:1863437/76‑1 (MQ=255)
                                  ccgccgACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCAt                                       >  1:4572724/1‑76 (MQ=255)
                                   cgccgACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATg                                      <  2:1740479/76‑1 (MQ=255)
                                    gccgACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCTGTGATCCGCATTCCGTTCTCATgg                                     >  1:1174782/1‑76 (MQ=255)
                                    gccgACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATgg                                     >  1:1029243/1‑76 (MQ=255)
                                     ccgACATAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCCGTGATCCGCATTCCGTTCTCATGGc                                    <  2:3671562/76‑1 (MQ=255)
                                     ccgACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGc                                    >  2:1802555/1‑76 (MQ=255)
                                      cgACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTACCAGGCAGTGATCCGCATTCCGTTCTCATGGCg                                   >  2:3574147/1‑76 (MQ=255)
                                       gACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGCGa                                  <  2:4756259/76‑1 (MQ=255)
                                       gACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGCGa                                  <  1:3326438/76‑1 (MQ=255)
                                        aCAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGCGAg                                 >  1:2190671/1‑76 (MQ=255)
                                          agagTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGCGAGGc                               >  1:3357989/1‑76 (MQ=255)
                                          agagTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGCGAGGc                               >  1:4292714/1‑76 (MQ=255)
                                           gagTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGCGAGGCa                              >  1:2420363/1‑76 (MQ=255)
                                           gagTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGCGAGGCa                              >  2:2670726/1‑76 (MQ=255)
                                           gagTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGCGAGGCa                              <  1:277965/76‑1 (MQ=255)
                                           gagTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGCGAGACa                              >  2:1206780/1‑76 (MQ=255)
                                            agTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGCGAGGCaa                             >  2:2374430/1‑76 (MQ=255)
                                            agTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGCGAGGCaa                             <  2:3072184/76‑1 (MQ=255)
                                              tAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGCGAGGCAACa                           <  2:4616759/76‑1 (MQ=255)
                                              tAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGCGAGGCAACa                           >  1:3275388/1‑76 (MQ=255)
                                                aTAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGCGAGGCAACAtt                         <  1:74294/76‑1 (MQ=255)
                                                  aGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGCGAGGCAACATTTc                       <  1:2620250/76‑1 (MQ=255)
                                                  aGGTTGTGCTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGCGAGGCAACATTTc                       <  2:1646019/76‑1 (MQ=255)
                                                  aGGGTTTACTTAATAGCTCTTCCTGGCCCTTCCAGGCAGTGATCCGCATTCCCTGCGCATGGCGAGGCAACATTTc                       >  2:1348554/1‑76 (MQ=255)
                                                   gggtttACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGCGAGGCAACCTTTCg                      >  2:4661470/4‑76 (MQ=255)
                                                   ggTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGCGAGGCAACATTTCg                      <  1:525880/76‑1 (MQ=255)
                                                   ggTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGCGAGGCAACATTTCg                      >  2:2304918/1‑76 (MQ=255)
                                                    gTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGCGAGGCAACATTTCgg                     <  2:824747/76‑1 (MQ=255)
                                                    gTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGCGAGGCAACATTTCgg                     <  2:1437430/76‑1 (MQ=255)
                                                     ttttACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGCGAGGCAACATTTCggg                    <  2:3816911/76‑1 (MQ=255)
                                                      tttACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGCGAGGCAACATTTCGGGa                   <  1:4852677/76‑1 (MQ=255)
                                                       ttACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGCGAGGCAACATTTCGGGAt                  <  2:3389083/76‑1 (MQ=255)
                                                       ttACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTCATCATGGCGAGGCAACATTTCGGGAt                  >  1:3965618/1‑76 (MQ=255)
                                                         aCTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGCGAGGCAACCTTTCGGGATgg                >  2:4173431/1‑76 (MQ=255)
                                                         aCTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGCGAGGCAACATTTCGGGATgg                >  2:2145691/1‑76 (MQ=255)
                                                          cTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGCGAGGCAACATTTCGGGATGGa               >  1:2653977/1‑76 (MQ=255)
                                                          cTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGCGAGGCAACATTTCGGGATGGa               <  1:1971462/76‑1 (MQ=255)
                                                          cTTAATAGATCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGCGAGGCAACATTTCGGGATGGa               <  1:1984225/76‑1 (MQ=255)
                                                           ttAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGCGAGGCAACATTTCGGGATGGaa              >  1:3232065/1‑76 (MQ=255)
                                                           ttAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGCGAGGCAACATTTCGGGATGGaa              >  1:1770375/1‑76 (MQ=255)
                                                           ttAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGCGAGGCAACATTTCGGGATGGaa              <  2:1033431/76‑1 (MQ=255)
                                                            tAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGCGAGGCAACATTTCGGGATGGAAg             <  2:168802/76‑1 (MQ=255)
                                                             aaTAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGCGAGGCAACATTTCGGGATGGAAGa            >  1:1237567/1‑76 (MQ=255)
                                                             aaTAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGCGAGGCAACATTTCGGGATGGAAGa            <  2:3648507/76‑1 (MQ=255)
                                                              aTAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGCGAGGCTACATTTCGGGATGGAAGAt           >  1:1942305/1‑76 (MQ=255)
                                                              aTAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGCGAGGCAACATTTCGGGATGGAAGAt           >  1:546392/1‑76 (MQ=255)
                                                               tAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGCGAGGCAACATTTCGGGATGGAAGATa          >  1:667814/1‑76 (MQ=255)
                                                               tAGCTCTTCCTGCCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGTGAGGCAACATTTCGGGATGGAAGATa          >  2:4300907/1‑76 (MQ=255)
                                                                aGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGCGAGGCAACATTTCGGGATGGAAGATaa         >  2:1192042/1‑76 (MQ=255)
                                                                aGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGCGAGGCAACATTTCGGGATGGAAGATaa         <  1:229281/76‑1 (MQ=255)
                                                                 gCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGCGAGGCAACATTTTGGGATGGAAGATAAt        <  2:3928980/76‑1 (MQ=255)
                                                                 gCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGCGAGGCAACATTTCGGGATGGAAGATAAt        <  2:3685527/76‑1 (MQ=255)
                                                                 gCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGCGAGGCAACATTTCGGGATGGAAGATAAt        <  1:3362607/76‑1 (MQ=255)
                                                                  ctctTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGCGAGGCAACATTTCGGGATGGAAGATAATg       <  1:4819738/76‑1 (MQ=255)
                                                                   tctTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGCGAGGCAACATTTCGGGATGGAAGCTAATGt      >  1:4061943/1‑76 (MQ=255)
                                                                   tctTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGCGAGGCAACATTTCGGGATGGAAGATAATGt      <  2:4016319/76‑1 (MQ=255)
                                                                   tctTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGCGAGGCAACATTTCGGGATGGAAGATAATGt      <  2:1718725/76‑1 (MQ=255)
                                                                    ctTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGCGAGGCAACATTTCGGGATGGAAGATAATGtt     <  2:1033184/76‑1 (MQ=255)
                                                                     ttCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGCGAGGCAACATTTCGGGATGGAAGATAATGTTc    <  2:3783031/76‑1 (MQ=255)
                                                                     ttCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGCGAGGCAACATTTCGGGATGGAAGATAATGTTc    >  2:1796927/1‑76 (MQ=255)
                                                                       ccTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGCGAGGCAACATTTCGGGATGGAAGATAATGTTCtt  <  1:3846112/76‑1 (MQ=255)
                                                                             |                                                                       
CAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTT  >  NC_000913/3256816‑3256961

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