Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 170,838 G→T 8.3% intergenic (+19/‑551) PP_0159 → / → PP_16SA family III CoA‑transferase/16S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947170,8380GT 8.3% 156.3 / 4.4 56intergenic (+19/‑551)PP_0159/PP_16SAfamily III CoA‑transferase/16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (24/27);  new base T (3/2);  total (27/29)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.30e-01

GAGACGCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAAGTAGTTA  >  NC_002947/170694‑170983
                                                                                                                                                 |                                                                                                                                                  
gagaCGCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGTGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAg                                                                                                                                                <  8:130148/149‑1 (MQ=255)
       cAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCTAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCTGACGGCCTccc                                                                                                                                         <  8:225715/149‑1 (MQ=255)
            gAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCtt                                                                                                                                    >  2:7704/1‑149 (MQ=255)
                  cGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGtttt                                                                                                                              >  8:282929/1‑149 (MQ=255)
                    ccGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCACGCCTACCCGc                                                                                                                                      >  6:266240/1‑139 (MQ=255)
                    ccGGGCGCCACCACTGCTGGGCGAGCATACCCAAGTGGTGCTGGGGTAAGTGTTGGGGCTGGATGCTGGCGAGGAACACCGTCTGCGTCATGCTGGGGTGCTCTATATTAAAAGCGGGCAGCCCGTCAGGCCTCCCCGc                                                                                                                                      <  5:266240/139‑1 (MQ=255)
                            caccacTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTa                                                                                                                     >  7:144045/1‑148 (MQ=255)
                                 atgctgGGCGAGAATACCTAAGTGGTGCTGGGGTAATTGTTGGGGCTGGTTGGTGGGGCGGTACAGCGGCTGCGTGCTGCTGTGGGGATTTAGTTGAAATGCGGTGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAg                                                                                                               <  6:107633/148‑1 (MQ=255)
                                  tgctgGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAgg                                                                                                              >  7:146013/1‑149 (MQ=255)
                                    ctgGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGt                                                                                                             >  6:258979/1‑148 (MQ=255)
                                         cGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTg                                                                                                        >  4:60134/1‑148 (MQ=255)
                                           aGCATACCGAAGTGGTGCTGGGGGAAGTGTGGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGaaa                                                                                                     <  7:282929/149‑1 (MQ=255)
                                               tACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGAGCTATGAAAGGTTATTGATATTa                                                                                                  >  4:101800/1‑148 (MQ=255)
                                               tACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGGCAGGCCAACCCGCATTTTTTTTTAGCGCAATGAAATGTAATTGAAATTa                                                                                                  >  6:116597/1‑148 (MQ=255)
                                                 ccGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTaaa                                                                                                <  8:144045/148‑1 (MQ=255)
                                                 ccGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTaaa                                                                                                <  5:258979/148‑1 (MQ=255)
                                                  cGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAgg                                                                                              >  4:260448/1‑149 (MQ=255)
                                                      gtggtgCTGGGGGAAGTGTTGGGGCTGTATGCTGGCTAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGa                                                                                          <  5:306303/149‑1 (MQ=255)
                                                            cTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTACTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCct                                                                                    >  4:246798/1‑149 (MQ=255)
                                                             tGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATCAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCctc                                                                                   <  2:243430/149‑1 (MQ=255)
                                                               ggggAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCtctc                                                                                 >  4:278340/1‑149 (MQ=255)
                                                                      tgtTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATccc                                                                          <  3:260448/149‑1 (MQ=255)
                                                                       gtTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTACTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATccc                                                                          <  3:246798/148‑1 (MQ=255)
                                                                          ggggCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCGGTGTTTTTCAGCGCGTTCAAAGGTAATTGCAAGTAAAGGTTGACGCCCTCTCGAATCGCCGta                                                                      >  4:297723/1‑149 (MQ=255)
                                                                           gggCTGGTTGGTTGCGACGTACAGCGTCGGCGTGATGCTGGGGTGCTCTAGATTAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTta                                                                      <  8:146013/148‑1 (MQ=39)
                                                                            ggatGGATGCTGGGGTTGTACAGCGTCGTCGTGATGCTGGGGTGCTCTTGTTGAAAAGCGGGGAGGCCTTCAGGCCGCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTTATTGAAATTAAAGGTTGACGCGCTCTCGAATCCCCTtata                                                                    <  7:228618/146‑1 (MQ=37)
                                                                                         cGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGAGCCCTATAAAAGGTAATTGAAATTAAAGGTTGACGCCCTATCGAATCCCCTTATAAAGCGCCacaca                                                        >  5:42636/1‑147 (MQ=38)
                                                                                              tACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAg                                                  <  3:60134/148‑1 (MQ=255)
                                                                                               aCAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCACCCCGCTTTGTTTTTAAGCGCAATCAAAGGTAATTTAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCCCCCCCCTTCCAGc                                                 >  2:76845/1‑148 (MQ=38)
                                                                                                cAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGa                                               >  5:248149/1‑149 (MQ=255)
                                                                                                 cgCGTCTGGTTTATGCTGGGGTTCTATAGATTAAAAGCGGGGAGGCCGTCTGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTc                                                      <  1:311009/141‑1 (MQ=37)
                                                                                                    gTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATc                                           <  6:248149/149‑1 (MQ=255)
                                                                                                     tCTGCGTGATGCTGGGGTGCTCTAGATTGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATc                                           >  1:174672/1‑149 (MQ=255)
                                                                                                           tGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATc                                                                            <  8:91218/110‑1 (MQ=255)
                                                                                                           tGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATc                                                                            >  7:91218/1‑110 (MQ=255)
                                                                                                                      tGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCTGCGACATCGGAACGAAAAACTCCTTg                         <  1:7704/149‑1 (MQ=255)
                                                                                                                            aGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCGTTGTTTTTGAGCGCTTTGAAAGGTAATTGATATTAAAGGTTGTCGCCCTGTCGAATACCCATATAATTCGCCCCACTTGCAGCGCGATCGGGACGAAAAACTCCTTGAGAtt                    >  8:54575/1‑148 (MQ=21)
                                                                                                                             gATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTc                   <  2:105214/148‑1 (MQ=255)
                                                                                                                             gATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCGTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTGTCGAATGCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTc                   >  5:297244/1‑148 (MQ=255)
                                                                                                                              aTGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCaa                 <  6:101320/149‑1 (MQ=255)
                                                                                                                                tGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCaa                 <  2:174672/148‑1 (MQ=255)
                                                                                                                                tGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAAt                <  7:54575/149‑1 (MQ=255)
                                                                                                                                   aaaGCGGGCAGGCCGTCAGGCCTCCCAGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCTAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAg             <  5:116597/149‑1 (MQ=255)
                                                                                                                                    aaGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCGCTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGt            >  5:240282/1‑149 (MQ=255)
                                                                                                                                     aGCGGGGAGGCCTGACGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTAACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGt            >  5:243377/1‑148 (MQ=18)
                                                                                                                                        ggggAGGCCTGACGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAag        >  2:93834/1‑149 (MQ=255)
                                                                                                                                        ggggAGGCCTGACGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAag        >  7:97588/1‑149 (MQ=255)
                                                                                                                                        ggggAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAag        <  3:278340/149‑1 (MQ=255)
                                                                                                                                         gggAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAag        >  3:119960/1‑148 (MQ=255)
                                                                                                                                         gggAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAag        <  4:119960/148‑1 (MQ=255)
                                                                                                                                          ggAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGaa                                       >  6:327243/1‑116 (MQ=255)
                                                                                                                                          ggAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGaa                                       <  5:327243/116‑1 (MQ=255)
                                                                                                                                            aGGCCGTCAGCCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACACCTTGAGATTCAATGAGTTAagtagt    <  6:42636/149‑1 (MQ=255)
                                                                                                                                             ggCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTgag                       >  2:82420/1‑129 (MQ=255)
                                                                                                                                             ggCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTgag                       <  1:82420/129‑1 (MQ=255)
                                                                                                                                               ccGTCAGGCCTCCCCGCTTTGTTTTTCATCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCATGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAAGTAGTTa  <  1:76845/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GAGACGCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAAGTAGTTA  >  NC_002947/170694‑170983

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: