Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,134,986 T→A 8.8% intergenic (‑47/+18) PP_4520 ← / ← PP_4521 hypothetical protein/aerotaxis receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,134,9860TA 8.8% 150.0 / 3.1 56intergenic (‑47/+18)PP_4520/PP_4521hypothetical protein/aerotaxis receptor
Reads supporting (aligned to +/- strand):  ref base T (30/21);  new base A (3/2);  total (33/23)
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

AGCAGGATCACACGGGTGCGGCTGTCGAGCGCGGGGATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTCCTGGCTGATTTCCGCGCCTTGCCTGGCTTGCCCGGCACTCTGGTCGCACAGCTGGGCAATTCTCACGATCTGC  >  NC_002947/5134842‑5135133
                                                                                                                                                |                                                                                                                                                    
agcagGATCACACGGGTGCGGCTGTCGAGCGCGGGGATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCGCTTTGc                                                                                                                                                  <  3:55715/149‑1 (MQ=255)
  cagGATCACACGGGTGCGGCTGTCGAGCGCGGGGATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCgg                                                                                                                                                <  6:326700/149‑1 (MQ=255)
  cagGATCACACGGGTGCGGCTGTCGAGCGCGGGGATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCgg                                                                                                                                                >  8:314581/1‑149 (MQ=255)
        cacaCGGGTGCGGCTGTCGAGCGCGGGGATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCGGGCCCcg                                                                                                                                          >  4:152828/1‑149 (MQ=255)
          cacGGGTGCGGCTGTCGAGCGCGGGGATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCGGGCCCcgc                                                                                                                                         >  4:251177/1‑148 (MQ=255)
                    gCTGTCGAGCGCGGGGATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCgg                                                                                                                             >  1:120808/1‑149 (MQ=255)
                            gcgcGGGGATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCt                                                                                                                     >  7:52483/1‑149 (MQ=255)
                               cgGGGATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTc                                                                                                                  <  5:94402/149‑1 (MQ=255)
                                ggggATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTcc                                                                                                                 >  7:62407/1‑149 (MQ=255)
                                         tGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTa                                                                                                        >  3:311404/1‑149 (MQ=255)
                                           ggCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTCTCCGCCAGGCTATg                                                                                                      <  7:314581/149‑1 (MQ=255)
                                              aTAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAg                                                                                                   >  2:359205/1‑149 (MQ=255)
                                                      tGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCg                                                                                           >  7:17755/1‑149 (MQ=255)
                                                         gcGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGTCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCa                                                                                        <  6:373495/148‑1 (MQ=255)
                                                           gAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGGACCTTTCCGCCAGGCTATGCAGGTACTCCGCCAt                                                                                       <  6:329603/148‑1 (MQ=255)
                                                           gAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGGACCTTTCCGCCAGGCTATGCAGGTACTCCGCCAt                                                                                       >  5:329603/1‑148 (MQ=255)
                                                           gAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATg                                                                                      >  4:104701/1‑149 (MQ=255)
                                                            aGCGGGCGCGGGCAGCGCGCGCAGGGGGGTCTGGGCAGGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATgc                                                                                     <  8:52483/149‑1 (MQ=255)
                                                                ggCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTc                                                                                 <  3:251177/149‑1 (MQ=255)
                                                                ggCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTc                                                                                 >  3:58564/1‑149 (MQ=255)
                                                                 gcgcTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCg                                                                                <  4:311404/149‑1 (MQ=255)
                                                                     tGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGt                                                                            >  7:89516/1‑149 (MQ=255)
                                                                         aGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGc                                                                                                     <  2:206911/120‑1 (MQ=255)
                                                                         aGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGc                                                                                                     >  1:206911/1‑120 (MQ=255)
                                                                         aGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCACGCGCTCGAGGTCCTg                                                                        >  6:77325/1‑149 (MQ=255)
                                                                              tCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTCCTGGCTGa                                                                   >  1:326798/1‑149 (MQ=255)
                                                                              tCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTCCTGGCTGa                                                                   >  7:280627/1‑149 (MQ=255)
                                                                                         tCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCg                                                                                <  1:146127/125‑1 (MQ=255)
                                                                                           tGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTCCTGGCTGATTTCCGCGCCTTg                                                      >  8:35313/1‑149 (MQ=255)
                                                                                                      gTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTCCTGGCTGATTTCCGCGCCTTGCCTGGCTTGcc                                           >  4:251515/1‑149 (MQ=255)
                                                                                                               gCGGGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTCCTGGCTGATTTCCGCGCCTTGCCTGGCTTGCCCGGCActct                                  >  8:34146/1‑149 (MQ=255)
                                                                                                               gCGGGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTCCTGGCTGATTTCCGCGCCTTGCCTGGCTTGCCCGGCActc                                   <  7:35313/148‑1 (MQ=255)
                                                                                                               gCGGGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTCCTGGCTGATTTCCGCGCCTTGCCTGGCTTGCCCGGCActc                                   <  2:120808/148‑1 (MQ=255)
                                                                                                                cGGGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTCCTGGCTGATTTCCGCGCCTTGCCTGGCTTGCCCGGCACTCTg                                 <  2:320657/149‑1 (MQ=255)
                                                                                                                 ggggACAGTTTGCCACATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTCCTGGCTGATTTCCGCGCCTTGCCTGGCTTGCCCGGCACTCTgg                                <  3:152828/149‑1 (MQ=255)
                                                                                                                 ggggACAGTTTGCCACATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTCCTGGCTGATTTACGCTCCTTTCCTGGCTTGCCCGGCACTCTgg                                >  4:216147/1‑149 (MQ=255)
                                                                                                                           tGCCACATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTCCTGGCTGATTTCCGCGCCTTGCCTGGCTTGCCCGGCACTCTGGTCGCACAGCt                      >  7:9557/1‑149 (MQ=255)
                                                                                                                             ccACATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTCCTGGCTGATTTCCGCGCCTTGCCTGGCTTGCCCGGCACTCTGGTCGCACAGCTgg                    <  8:280627/149‑1 (MQ=255)
                                                                                                                             ccACATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTCCTGGCTGATTTCCGCGCCTTGCCTGGCTTGCCCGGCACTCTGGTCGCACAGCTgg                    <  8:17755/149‑1 (MQ=255)
                                                                                                                              cacaTAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTAAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTCCTGGCTGATTTCCGCGCCTTGCCTGGCTTGCCCGGCACTCTGGTCGCACAGCTggg                   >  8:150369/1‑149 (MQ=255)
                                                                                                                                caTAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCCGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTCCTGGCTGATTTCCGCGCCTTGCCTGGCTTGCCCGGCACTCTGGTCGCACAGCTGGGc                  >  2:323414/1‑148 (MQ=255)
                                                                                                                                caTAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTCCTGGCTGATTTCCGCGCCTTGCCTGGCTTGCCCGGCACTCTGGTCGCACAGCt                      <  8:9925/144‑1 (MQ=255)
                                                                                                                                caTAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTCCTGGCTGATTTCCGCGCCTTGCCTGGCTTGCCCGGCACTCTGGTCGCACAGCt                      >  7:9925/1‑144 (MQ=255)
                                                                                                                                caTAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTCCTGGCTGATTTCCGCGCCTTGCCTGGCTTGCCCGGCACTCTGGTCGCACAGCTGGGc                  >  1:211709/1‑148 (MQ=255)
                                                                                                                                caTAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTACTGGCTGATTTCCGCGCCTTGCCTGGCTTGCCCGGCACTCTGGTCGCACAGCTGGGc                  >  8:286082/1‑148 (MQ=255)
                                                                                                                                   aaCCGGGGCCCGCTTTGCGGGCCCCGCGGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGGCCTGGCGGATTTCCGCGCCTTGCCTGGCTTGCCCGGCACTCTGGTCGCACAGCTGGGCAAt               >  1:209552/1‑149 (MQ=255)
                                                                                                                                      cGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTCCTGGCTGATTTCCGCGCCTTGCCTGGCTTGCCCGGCACTCTGGTCGCACAGCTGGGCAATtctc           <  3:104701/149‑1 (MQ=255)
                                                                                                                                      cGGGGCCCGCAAAGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTCCTGGCTGATTTCCGCGCCTTGCCTGGCTTGCCCGGCACTCTGGTCGCACAGCTGGGCAATtct            >  8:30640/1‑148 (MQ=255)
                                                                                                                                        gggCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTCCTGGCTGATTTCCGCGCCTTGCCTGGCTTGCCCGGCACTCTGGTCGCACAGCTGGGCAATTCTCAc         >  2:52756/1‑149 (MQ=255)
                                                                                                                                        gggCCCGCAACGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTCCTGGCTGATTTCCGCGCCTTGCCTGGc                                                <  8:323114/110‑1 (MQ=255)
                                                                                                                                        gggCCCGCAAAGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGCGGTCCTGGCTGATTTCCGCGCCTTGCCTGGc                                                >  7:327986/1‑110 (MQ=255)
                                                                                                                                        gggCCCGCAAAGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTCCTGGCTGATTTCCGCGCCTTGCCTGGc                                                <  8:327986/110‑1 (MQ=255)
                                                                                                                                        gggCCCGCAAAGCGGGCCCCGCGAGTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCGCGAGGTCCTGGCTGATTTCCGCGCCTTGCCTGGc                                                >  7:323114/1‑110 (MQ=255)
                                                                                                                                           cccGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTCCTGGCTGATTTCCGCGCCTTGCCTGGCTTGCCCGGCACTCTg                                 <  8:205212/122‑1 (MQ=255)
                                                                                                                                           cccGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTCCTGGCTGATTTCCGCGCCTTGCCTGGCTTGCCCGGCACTCTg                                 >  7:205212/1‑122 (MQ=255)
                                                                                                                                               cTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTCCTGGCTGATTTCCGCGCCTTGCCTGGCTTGCCCGGCACTCTGGTCGCACAGCTGGGCAATTCTCACGATctgc  <  5:308414/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
AGCAGGATCACACGGGTGCGGCTGTCGAGCGCGGGGATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCGCTTTGCGGGCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTCCTGGCTGATTTCCGCGCCTTGCCTGGCTTGCCCGGCACTCTGGTCGCACAGCTGGGCAATTCTCACGATCTGC  >  NC_002947/5134842‑5135133

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: