Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,310,322 A→T 9.9% V28V (GTT→GTA ada ← fused DNA‑binding transcriptional dual regulator/O6‑methylguanine‑DNA methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,310,3220AT 9.9% 131.7 / 2.6 54V28V (GTT→GTAadafused DNA‑binding transcriptional dual regulator/O6‑methylguanine‑DNA methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (20/28);  new base T (3/3);  total (23/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTTTTTCTGGCTGACAACGTTTGCAGGGGCGAAAGCCAGCGGCGAGTGCCTCGCTGGCATTTGCGTAGAAGGAGACGTTTTCCCGCAAAGCATGTCTGGCGCGGCAAGACGGACGGCAAAAGATGCCTGTGGTACGCACGGCGAAAACGAATTCGCCGTCGGCATTCGGGTCGCGGGCTAAGACAGATTGCCAGCGTTGATCGTCAGTTAAGCATGTGGCTTTTTTCATAATCAGCTCCCTGGTTAAGGATAGCCTTTAGGCTGCCCGGTCACCATCACGCAAAAACCAACAAT  >  NC_000913/2310176‑2310470
                                                                                                                                                  |                                                                                                                                                    
gCTTTTTCTGGCTGACAACGTTTGCAGGGGCGAAAGCCAGCGGCGAGTGCCTCGCTGGCATTTGCGTAGAAGGAGACGTTTTCCCGCAAAGCATGTCAGGCGCGGCAAGACGGACGGCAAAAGATGCCTGTGGTACGCACGGCGAAAAc                                                                                                                                                    >  2:873489/1‑149 (MQ=255)
  tttttCTGGCTGACAACGTTTGCAGGGGCGAAAGCCAGCGGCGAGTGCCTCGCTGGCATTTGCGTAGAAGGAGACGTTTTCCCGCAAAGCATGTCTGGCGCGGCAAGACGGACGGCAAAAGATGCCTGTGGTACGCACGGCGAAAACGa                                                                                                                                                  >  2:341502/1‑149 (MQ=255)
      tCTGGCTGACAACGTTTGCAGGGGCGAAAGCCAGCGGCGAGTGCCTCGCTGGCATTTGCGTAGAAGGAGACGTTTTCCCGCAAAGCATGTCTGGCGCGGCAAGACGGACGGCAAAAGATGCCTGTGGTACGCACGGCGAAAACGAATTc                                                                                                                                              >  1:948368/1‑149 (MQ=255)
        tGGCTGACAACGTTTGCAGGGGCGAAAGCCAGCGGCGAGTGCCTCGCTGGCATTTGCGTAGAAGGAGACGTTTTCCCGCAAAGCATGTCTGGCGCGGCAAGACGGACGGCAAAAGATGCCTGTGGTACGCACGGCGAAAACGAATTCGc                                                                                                                                            >  2:272472/1‑149 (MQ=255)
        tGGCTGACAACGTTTGCAGGGGCGAAAGCCAGCGGCGAGTGCCTCGCTGGCATTTGCGTAGAAGGAGACGTTTTCCCGCAAAGCATGTCTGGCGCGGCAAGACGGACGGCAAAAGATGCCTGTGGTACGCACGGCGAAAACGAATTCGc                                                                                                                                            >  1:4711/1‑149 (MQ=255)
         ggCTGACAACGTTTGCAGGGGCGAAAGCCAGCGGCGAGTGCCTCGCTGGCATTTGCGTAGAAGGAGACGTTTTCCCGCAAAGCATGTCTGGCGCGGCAAGACGGACGGCAAAAGATGCCTGTGGTACGCACGGCGAAAACGAATTCGcc                                                                                                                                           >  2:718370/1‑149 (MQ=255)
           cTGACATCGTTTGCTTGGCCGAAAGCCAGCGGCGCGTGCCTCGCTGGCATTTGCGTAGAAGGAGACGTTTTCCCCCAAAGCATGTCTGGCGCGGCAAGACGGACGGCAAAAGATGCCTGTTGTACGCTCGGCGAATTCGTTTTCGCCGt                                                                                                                                         <  1:117169/149‑1 (MQ=255)
           cTGACACCGTTTGCAGGGGCGAAAGCCAGCGGCGAGTGCCGCGCTGGCATTTGCGTAGAAGGAGACGTTTTCCCGCAAAGCATGTCTGGCGCGGCAAGACGGACGGCAAAAGATGCCTGTGGTACGCACGGCGAATTCGTTTTCGCCGt                                                                                                                                         <  2:364493/149‑1 (MQ=255)
             gACCACGTTTGCAGGGGCGAACGCCAGCGGCGAGTGCCTCGCTGGCATTTGCGTAGAAGGAGACGGTTTCCCGCAAAGCATGTCTGGCGCGGCAAGACGGACGGCAAAAGATGCCTGTGGTACGCACGGCGAAAACGAATTCGCCGTCg                                                                                                                                       <  2:52910/149‑1 (MQ=255)
             gACAACGTTTGCAGGGGCGAAAGCCAGCGGCGAGTGCCTCGCTGGCATTTGCGTAGAAGGAGACGTTTTCCCGCAAAGCATGTCTGGCGCGGCAAGACGGACGGCAAAAGCTGCCTGTGGTACGCACGGCGAAAACGAATTCGCCGTCg                                                                                                                                       >  1:263511/1‑149 (MQ=255)
             gACAACGTTTGCAGGGGCGAAAGCCAGCGGCGAGTGCCTCGCTGGCATTTGCGTAGAAGGAGACGTTTTCCCGCAAAGCATGTCTGGCGCGGCAAGACGGACGGCAAAAGATGCCTGTGGTACGCACGGCGAAAACGAATTCGCCGTCg                                                                                                                                       >  2:20104/1‑149 (MQ=255)
              aCAACGTTTGCAGGGGCGAAAGCCAGCGGCGAGTGCCTCGCTGGCATTTGCGTAGAAGGAGACGTTTTCCCGCAAAGCATGTCTGTCGCGGCAAGACGGACGGCAAAAGATGCCTGTGGTAAGCACGGCGAAAACGAATTCGCCGTCgg                                                                                                                                      >  2:248215/1‑149 (MQ=255)
              aCAACGTTTGCAGGGGCGAAAGCCAGCGGCGAGTGCCTCGCTGGCATTTGCGTAGAAGGAGACGTTTTCCCGCAAAGCATGTCTGGCGCGGCAAGACGGACGGCAAAAGATGCCTGTGGTACGCACGGCGAAAACGAATTCGCCGTCgg                                                                                                                                      <  1:27816/149‑1 (MQ=255)
              aCAACGTTTGCAGGGGCGAAAGCCAGCGGCGAGTGCCTCGCTGGCATTTGCGTAGAAGGAGACGTTTTCCCGCAAAGCATGTCTGGCGCGGCAAGACGGACGGCAAAAGATGCCTGTGGTACGCACGGCGAAAACGAATTCGCCGTCgg                                                                                                                                      <  1:393043/149‑1 (MQ=255)
                     ttGCAGGGGCGAAAGCCTGCGGCGAGTGCCTCGCTGGCATTTGCGTAGAAGGAGACGTTTTCCCGCAAAGCATGTCTGGCGCGGCAAGACGGACGGCAAAAGATGCCTGTGGTACGCACGGCGAAAACGAATTCGCCGTCGGCATTCgg                                                                                                                               <  2:232859/149‑1 (MQ=255)
                       gCAGGGGCGAAAGCCAGCGGCGAGTGCCTCGCTGGCATTTGCGTAGAAGGAGACGTTTTCCCGCAAAGCATGTCTGGCGCGGCAAGACGGACGGCAAAAGATGCCTGTGGTACGCACGGCGAAAACGAATTCGCCGTCGGCATTCGGGt                                                                                                                             >  1:140844/1‑149 (MQ=255)
                             gCGAAAGCCAGCGGCGAGTGCCTCGCTGGCATTTGCGTTGAAGGAGACGTTTTCCCGCAAAGCATGTCTGGCGCGGCAAGACGGACGGCAAAAGATGCCTGTGGTACGCACGGCGAAAACGAATTCGCCGTCGGCATTCGGGTCGCggg                                                                                                                       >  1:216460/1‑149 (MQ=255)
                             gCGAAAGCCAGCGGCGAGTGCCTCGCTGGCATTTGCGTAGAAGGAGACGTTTTCCCGCAAAGCATGTCTGGCGCGGCATGACGGACGGCAAAAGATGCCTGTGGTACGCACGGCGAATACGAATTCGCCGTCGGCATTCGGGTCGCggg                                                                                                                       >  2:671880/1‑149 (MQ=255)
                                   gCCAGCGGCGAGTGCCTCGCTGGCATTTGCGTAGAAGGAGACGTTTCCCGCAAAGCATGTCTGGCGCGGCAAGACGGACGGCAAAAGATGCCTGTGGTACGCACGGCGAAAACGAATTCGCCGTCGGCATTCGGGTCGCGGGCTAAGAc                                                                                                                <  1:635039/149‑1 (MQ=255)
                                    ccAGCGGCGAGTGCCGCGCTGGCATTTGCGTAGAAGGAGACGGTTTCCCGCAAAGCATGTCTGGCGCGGCAAGACGGACGGCAAAAGATGCCTGTGGTACGCACGGCGAAAACGAATTCGCCGTCGGCATTCGGGGCGCGGGCTAAGAc                                                                                                                <  1:910511/149‑1 (MQ=255)
                                     cAGCGGCGAGTGCCTCGCTGGCATTTGCGTAGAAGGAGACGTTTTCCCGCAAAGCATGTCTGGCGCGGCAAGACGGACGGCAAAAGATGCCTGTGGTACGCACGGCGAAAACGAATTCGCCGTCGGCATTCGGGTCGCGGGCTAAGACa                                                                                                               <  2:307537/149‑1 (MQ=255)
                                     cAGCGGCGAGTGCCTCGCTGGCATTTGAGTAGAAGGAGACGTTTTCCCGCAAAGCATGTCTGGCGCGGCAAGACGGACGGCAAAAGATGAATGTGGTACGCACGGCGAAAACGACTTCGCCGTCGGCATTCGGGTCGCGGCCTAAGCCa                                                                                                               >  2:211963/1‑149 (MQ=255)
                                      aGCGGCGAGTGCCTCGCTGGCATTTGCGTAGAAGGAGACGTTTTCCCGCAAAGCATGTCTGGCGCGGCAAGACGGACGGCAAAAGATGCCTGTGGTACGCACGGCGAAAACGAATTCGCCGTCGGCATTCGGGTCGCGGGCTAAGACAg                                                                                                              >  1:218739/1‑149 (MQ=255)
                                           cgAGTGCCTCGCTGGCATTTGCGTAGAAGGAGACGTTTTCCCGCAAAGCATGTCTGGCGCGGCAAGACGGACGGCAAAAGATGCCTGTGGTACGCACGGCGAAAACGAATTCGCCGTCGGCATTCGGGTCGCGGGCTAAGACAGATTGc                                                                                                         <  2:139442/149‑1 (MQ=255)
                                            gAGTGCCTCGCTGGCATTTGCGTAGAAGGAGACGTTTTCCCGCAAAGCATGTCTGGCGCGGCAAGACGGACGGCAAAAGATGCCTGTGGTACGCACGGCGAAAACGAATTCGCCGTCGGCATTCGGGTCGCGGGCTAAGACAGATTGcc                                                                                                        >  1:53336/1‑149 (MQ=255)
                                             aGTGCCTCGCTGGCATTTGCGTAGAAGGAGACGTTTTCCCGCAAAGCATGTCTGGCGCGGCAAGACGGACGGCAAAAGATGCCTGTGGTACGCACGGCGAAAACGAATTCGCCGTCGGCATTCGGGTCGCGGGCTAAGACAGATTGCCa                                                                                                       <  2:643739/149‑1 (MQ=255)
                                              atGCCGCGCTGGCATTTGCGTAGAAGGAGACGTTTTCCCGCTACGGATGTCTCGCGCGGCCCGCCGGAGGGCAAAAGCTGCCTGTCGTACGCACGGCGAATACGAATTCGCCGTCGTCATTCGGGGCGCGGGCTCAGACAGATTGCCAg                                                                                                      <  2:510582/148‑1 (MQ=255)
                                                 ccTCGCTGGCATTTGCGTAGAAGGAGACGTTTTCCCGCAAAGCATGTCTGGCGCGGCAAGACGGACGGCAAAAGATGCCTGTGGTACGCACGGCGAAAACGAATTCGCCGTCGGCATTCGGGTCGCGGGCTAAGACAGAATGCCAGCgg                                                                                                   >  2:51155/1‑148 (MQ=255)
                                                  cTCGCTGGCATTTGCGTAGAAGGAGACGTTTTCCCGCAAAGCATGTCTGGCGCGGCAAGACGGACGGCAAAAGATGCCTGTGGTACGCACGGCGAAAACGAATTCGCCGTCGGCATTCGGGTCGCGGGCTAAGACAGATTGCCAGCGtt                                                                                                  <  1:898194/149‑1 (MQ=255)
                                                    cGCTGGCATTTGCGTAGAAGGAGACGTTTTCCCGCAAAGCATGTCTGGCGCGGCAAGACGGACGGCAAAAGATGCCTGTGGTACGCACGGCGAAAACGAATTCGCCGTCGGCATTCGGGTCGCGGGCTAAGACAGATTGCCAGCGTTGa                                                                                                <  1:206274/149‑1 (MQ=255)
                                                      cgggCATTTGCGTAGAAGTAGCGGTTTTCCCGCAAAGCAGGTCTGGCGCGGCAAGACGGACGGCTAAAGATGCCTGTGGTACGCACGGCGATAACGAATTCGCCGGCGGCATGCGGGTCGCGGGCTAAGACAGATTGCCAGCGTTGATc                                                                                              <  2:419519/147‑1 (MQ=255)
                                                      cgggCATTTGCGTAGAAGGAGACGGTTTCCCGCAAAGCATGTCTGGCGCGGCAAGACGGACGGCAAAAGATGCCTGTGGTACGCACGGCGAAAACGAATTCGCCGTCGGCATTCGGGTCGCGGGCTAAGACAGATTGCCAGCGTTGATc                                                                                              <  2:931914/147‑1 (MQ=255)
                                                         gCATTTGCGTAGAAGGAGACGTTTTCCCGCAAAGCATGTCTGGCGCGGCAAGACGGACGGCAAAAGATGCCTGTGGTACGCACGGCGAAAACGAATTCGCCGTCGGCATTCGGGTCGCGGGCTAAGACAGATTGCCAGCGTTGATCGTc                                                                                           >  2:348011/1‑149 (MQ=255)
                                                          cATTTGCGTAGAAGGAGACGTTTTCCCGCAAAGCATGTCTGGCGCGGCAAGACGGACGGCAAAAGATGCCTGTGGTACGCACGGCGAAAACGAATTCGCCGTCGGCATTCGGGTCGCGGGCTAAGACAGATTGCCAGCGTTGATCGTCa                                                                                          <  2:454854/149‑1 (MQ=255)
                                                          cATTTGCGTAGAAGGAGACGTTTTCCCGCAAAGCATGTCTGGCGCGGCAAGACGGACGGCAAAAGATGCCTGTGGTACGCACGGCGAAAACGAATTCGCCGTCGGCATTCGGGTCGCGGGCTAAGACAGATTGCCAGCGTTGATCGTCa                                                                                          <  2:560371/149‑1 (MQ=255)
                                                               gCGTAGAAGGAGACGTGTTCCCGCAAAGCATGTCTGGCGCGTCAAGACGGAGGGCAATAGATGCCTGTGGTACGCACGGCGAAAACGAATTCGCCGTCGGCATTCGGGTCGCGGGCTAAGACAGATTGCCAGCGTTGATCGTCAGTTaa                                                                                     <  1:74371/149‑1 (MQ=255)
                                                                cGTAGAAGGAGACGTTTTCCCGCAAAGCATGTATGGCGCGGCAAGACGGACGGCAAAAGATGCCTGTGGTACGCACGGCGAAAACGAATTCGCCGTCGGAATTCGGGTCGCGGGCTAAGACAGATTGACAGAGTTGATCGTCTGTTAAg                                                                                    >  1:573622/1‑149 (MQ=255)
                                                                 gTAGAAGGAGACGTTTTCCCGCAAAGCATGTCTGGCTCGGCAAGACGGACGGCAAAAGATGCCTGTGGTACGCACGGCGAAAACGAATTCGCCGTCGGCATTCGGGTCGCGGGCTAAGACAGATTGCCAGCGTTGATCGTCAGTTAAGc                                                                                   <  2:762476/149‑1 (MQ=255)
                                                                     aaGGAGACGTTTTCACGCAAAGCATGTCTGGCGCGGCCAGACGGACGGCAAAAGATGCCTGTGGTACGCACGGCGATAACGAATTCGACGTCGGCATTCGGGTCTCGGGCGATGACAGATTGGCAGCGTTGATCGTCAGTTTTGCAtgt                                                                               >  2:950692/1‑149 (MQ=255)
                                                                       ggAGACGTTTTCCCGCAAAGCATGTCTGGCGCGGCAAGACGGACGGCAAAAGATGCCTGTGGTACGCACGGCGAAAACGAATTCGCCGTCGGCATTCGGGTCGCGGGCTAAGACAGATTGCCAGCGTTGATCGTCAGTTAAGCATGTgg                                                                             <  1:341502/149‑1 (MQ=255)
                                                                            cGTTTTCCCGCAAAGCATGTCTGGCGCGGCAAGACGGACGGCAAAAGATGCCTGTGGTACGCACGGCGAAAACGAATTCGCCGTCGGCATTCGGGTCGCGGGCTAAGACAGATTGCCAGCGTTGATCGTCAGTTAAGCATGTGGCtttt                                                                        <  1:535579/149‑1 (MQ=255)
                                                                             gTTTTCCCGCAAAGCATGTCTGGCGCGGCAAGACGGACGGCAAAAGATGCCTGTGGTACGCACGGCGAAAACGAATTCGCCGTCGGCCTTCGGGTCGCGGGCTAAGACAGATTGCCAGCGTTGATCGTCAGTTAAGCATGTGGCttttt                                                                       >  1:1014480/1‑149 (MQ=255)
                                                                               gttCCCGCAAAGCAGGTCTGGCGCGGCAAGACGGACGGCAAAAGCTGCCTGTGGTACGCACGGCGAAAACGAATTCGCCGTCGGCATTCGGGGCGCGGGCTAAGACAGATTGCCAGCGTTGATCGTCAGTTAAGCATGTGGCTTTTTTc                                                                     <  1:744910/148‑1 (MQ=255)
                                                                                  cccGCAAAGCAAGGCTGGCGCGGCAAGACGGACGGCAAAAGATGCCTGTGGTACGCACGGCGAAAACGAAATCGCCGTCGGCATTCGGGGCGCGGGCTAAGACAGATTGCCAGCGTTGATCGTCAGTTAAGCATGTGGCTTTTTTCATa                                                                  <  2:934615/149‑1 (MQ=255)
                                                                                                 tGGCGCGGCAAGACGGACGGCAAAAGATGCTGTGGTACGCACGGCGAAAACGAATTCGCCGTCGGCATTCGGGTCGCGGGCTACGACAGATTGCCAGCGTTGATCGGCAGTTAAGCATGGGGCTTTGTTCATAATCAGCTCCCTGGTTa                                                  >  1:823219/1‑149 (MQ=255)
                                                                                                              cggacggCAAAAGATGCCTGTGGTACGCACGGCGAAAACGAATTCGCCGTCGGCATTCGGGTCGCGGGCTAAGACAGATTGCCAGCGTTGATCGTCAGTTAAGCATGTGGCTTTTTTCATAATCAGCTCCCTGGTTAAGGATAGCCttt                                      <  1:233904/149‑1 (MQ=255)
                                                                                                               ggacggCAAAAGATGCCTGTGGTACGCACGGCGAAAACGAATTCGCCGTCGGCATTCGGGTCGCGGGCTAAGACAGAGTGCCAGCGTTGATCGTCAGTTAAGCATGTGGCTTTTTTCATAATCAGCTCCCTGGTTAAGGATAGCCTTTa                                     <  2:216460/149‑1 (MQ=255)
                                                                                                                gacggCAAAAGATGCCTGTGGTACGCACGGCGAAAACGAATTCGCCGTCGGCATTCGGGTCGCGGGCTAAGACAGATTGCCAGCGTTGATCGTCAGTTAAGCATGTGGCTTTTTTCATAATCAGCTCCCTGGTTAAGGATAGCCTTTAg                                    <  1:348011/149‑1 (MQ=255)
                                                                                                                 acgtCAAAAGATGCCTGTGGTACGCACGGCGAAAACGAATTCGCCGTCGGCATTCGGGTCGCGGGCTAAGACAGATTGCCAGCGTTGATCGTCAGTTAAGCATGTGGCTATTTTCATACTCAGCTCCCTGGTTAAGGATAGCCTTTAgg                                   <  1:912879/149‑1 (MQ=255)
                                                                                                                        aaGATGCCTGTGGTACGCACGGCGAAAACGAATTCGCCGTCGGCATTCGGGTCGCGGGCTAAGACAGATTGCCAGCGTTGATCGTCAGTTAAGCATGTGGCTTTTTTCATAATCAGCTCCCTGGTTAAGGATAGCCTTTAGGCTGCCCg                            <  1:305407/149‑1 (MQ=255)
                                                                                                                          gATGCCTGTGGTACGCACGGCGAAAACGAATTCGCCGTCGGCATTCGGGTCGCGGGCTAAGACAGATTGCCAGCGTTGATCGTCAGTTAAGCATGTGGCTTTTTTCATAATCAGCTCCCTGGTTAAGGATAGCCTTTAGGCTGCCCGGt                          <  2:400495/149‑1 (MQ=255)
                                                                                                                                        cgaCGGCGAATTCGTTTTCGCCGTCGGCATTCGGGTCGCGGGCTAAGACAGATTGCCAGCGTTGATCGTCAGTTAAGCATGTGGCTTTTTTCATAATCAGCTCCCTGGTTAAGGATAGCCTTTAGGCTGCCCGGTCACCATCCCGCaaa            >  1:627525/3‑149 (MQ=255)
                                                                                                                                        cgaCGGCGAATTCGTTTTCGCCGTCGGCATTCGGGTCGCGGGCTAAGACAGATTGCCAGCGTTGATCGTCAGTTAAGAATGTGGCTTTTTTCATAATCAGCTCACTGGTTAATGATATCCTTTAGGGTGCCCGGTCACCATCACGCaaa            >  1:689044/3‑149 (MQ=255)
                                                                                                                                           cGGCGAAAACGAATTCGCCGTCGGCATTCGGGTCGCGGGCTAAGACAGATTGCCAGCGTTGATCGTCAGTTAAGCATGTGGCTTTTTTCATAATCAGCTCCCTGGTTAAGGATAGCCTTTAGGCTGCCCGGTCACCATCACGCAAAAAc         <  1:369800/149‑1 (MQ=255)
                                                                                                                                                  aaCGAATTCGCCGTCGGCATTCGGGTCGCGGGCTAAGACAGATTGCCAGCGTTGATCGTCAGTTAAGCATGTGGCTTTTTTCATAATCAGCTCCCTGGTTAAGGATAGCCTTTAGGCTGCCCGGTCACCATCACGCAAAAACCAACAAt  <  2:945165/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GCTTTTTCTGGCTGACAACGTTTGCAGGGGCGAAAGCCAGCGGCGAGTGCCTCGCTGGCATTTGCGTAGAAGGAGACGTTTTCCCGCAAAGCATGTCTGGCGCGGCAAGACGGACGGCAAAAGATGCCTGTGGTACGCACGGCGAAAACGAATTCGCCGTCGGCATTCGGGTCGCGGGCTAAGACAGATTGCCAGCGTTGATCGTCAGTTAAGCATGTGGCTTTTTTCATAATCAGCTCCCTGGTTAAGGATAGCCTTTAGGCTGCCCGGTCACCATCACGCAAAAACCAACAAT  >  NC_000913/2310176‑2310470

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