Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 209,699 A→G 22.5% K578E (AAG→GAG)  dnaE → DNA polymerase III subunit alpha

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606209,6990A→G22.5% 5.3 ‑2.8 0.0311732dnaEDNA polymerase III subunit alpha
Reads supporting (aligned to +/- strand):  new base (7/0):  ref base (12/10):  total (20/11)
Fisher's exact test strand distribution p-value = 6.25e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

ACTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGAGCCGCCGCTGGATATCGCCGCGATC  >  REL606/209664‑209734
                                   |                                   
aCTGGGCGCTGGAGATGATCAACAAGCGGCGGGCga                                     <  1:7082495/36‑1 (MQ=255)
 cTGGGCGCTGGAGATGATCAACAAGCGGCGGGCgaa                                    <  1:1877095/36‑1 (MQ=255)
  tGGGCGCTGGAGATGATCAACAAGCGGCGGGCgaag                                   <  1:2131031/36‑1 (MQ=255)
  tGGGCGCTGGAGATGATCAACAAGCGGCGGGCgaag                                   <  1:4676734/36‑1 (MQ=255)
   gggCGCTGGAGATGATCAACAAGCGGCGGGCgaaga                                  <  1:195148/36‑1 (MQ=255)
         tGGAGATGATCAACAAGCGGCGGGGGGAGAAAgggg                            >  1:1891606/1‑34 (MQ=25)
          ggAGATGATCAACAAGCGGCGGGGGGAGAATggggg                           >  1:1086330/1‑33 (MQ=25)
          ggAGATGATCAACAAGCGGCGGGCGAAGAATGGCGa                           <  1:7463373/36‑1 (MQ=255)
           gagaTGAACAACAAGCGGCGGGCGAAGAATGGCGAg                          <  1:4926844/36‑1 (MQ=255)
            agaTGATCAACAAGCGGCGGGCGAAGAATGGCGAgc                         <  1:4360282/36‑1 (MQ=255)
               tgatCAACAAGCGGCGGGCGAAGAATGGCGAgccgc                      >  1:2436327/1‑36 (MQ=255)
                gatCAACAAGCGGCGGGGGGAGGAAGGCGAgacggc                     >  1:4396990/1‑34 (MQ=12)
                gatCAACAAGCGGCGGGCGGAGAAAGGGGAgcggcc                     >  1:4813057/1‑36 (MQ=16)
                gatCAACAAGCGGCGGGCGCAGAATGGCGAgccgcc                     <  1:1863538/36‑1 (MQ=255)
                 atCAACAAGCCGCGGGAGGACAATTGCCAgccaccc                    >  1:442696/1‑35 (MQ=0)
                  tCAACAAGCGGCGGGCGAAGAATGGCGAGCcgccgc                   >  1:2268331/1‑36 (MQ=255)
                   caacaaGCGGCGGGCGGAGAAAGGCGAGgcgccggg                  >  1:4765826/1‑34 (MQ=25)
                   caacaaGCGGCGGGCGCAGAATGGCGAGGcgccgcc                  >  1:7547777/1‑35 (MQ=25)
                   caacaaGCGGCGGGCGAAGAATGGGGAGGAGCCGCt                  >  1:2742657/1‑36 (MQ=37)
                    aacaaGCGGCGGGGGGAGAATGGCGAGCCGCCGCTg                 >  1:5866200/1‑36 (MQ=25)
                     acaaGCGGCGGGCGAAGAATGGCGGGCCGCCGCTgg                >  1:3316020/1‑36 (MQ=255)
                      caaGCGGCGGGCGAAGAATGGTGCGCCGCCGCTGGa               <  1:5263316/36‑1 (MQ=37)
                       aaGCGGCGGGCGAAGAATGGCGAGCCGCCGCTgggt              >  1:3163138/1‑34 (MQ=255)
                       aaGCGGCGGGCGAAGAATGGCGAGCCGCCGCTGGat              <  1:3166815/36‑1 (MQ=255)
                         gcggcgGGCGAAGAATGGCGAGCCGCCGCTGGatat            <  1:6048149/36‑1 (MQ=255)
                          cggcggGCGAAGAATGGCGAGCCGCCGCGGGGTATc           >  1:1011690/1‑36 (MQ=37)
                          cggaggGCAAAGAATGGCGAGCCGCCGCTGGATATc           <  1:4988342/36‑1 (MQ=25)
                           ggcgggcgAAGAATGGCGAGCCGCCGCTGGATATcg          <  1:5939050/36‑1 (MQ=255)
                            gcgggcgAAGAATGGCGAGCCGCCGCTGGATATcgc         <  1:3945851/36‑1 (MQ=255)
                            gcgggcgAAGAATGGCGAGCCGCCGCTGGATATcgc         <  1:6478028/36‑1 (MQ=255)
                              gggcgAAGAATGGCGAGCCGCCGCTGGATATcgccg       >  1:470436/1‑36 (MQ=255)
                              gggcgAAGAATGGCGAGCCGCCGCTGGATATcgcag       >  1:6975673/1‑34 (MQ=255)
                               ggcgAAGAATGGCGAGCCGCCGCTGGATATcgccgc      >  1:6479942/1‑36 (MQ=255)
                                gcgAAGAATGGCGAGCCGCCGCTGGATATCGCcgcg     >  1:4503375/1‑36 (MQ=255)
                                gcgAAGAATGGCGAGCCGCCGCTGGATATCGCcgcg     >  1:4231002/1‑36 (MQ=255)
                                 cgAAGAATGGCGAGCCGCCGCTGGATATCGCCGCGa    >  1:197464/1‑36 (MQ=255)
                                  gaagaaTGGCGAGCCGCCGCCGGATATCGCCGCGAt   >  1:2018725/1‑36 (MQ=255)
                                  gaagaaTGGCGAGCCGACGCTGGATATCGCCGAGa  >  1:2559290/1‑35 (MQ=25)
                                  gaagaaTGGCGAGCAGCCGCTGGATATCGACGCGAt   >  1:524135/1‑36 (MQ=37)
                                   aagagTGGCGACCCGCCGCTGGATATCGCCGCGATc  <  1:387579/36‑1 (MQ=21)
                                   |                                   
ACTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGAGCCGCCGCTGGATATCGCCGCGATC  >  REL606/209664‑209734

Base quality scores:  ATCG  < 3 ≤  ATCG  < 6 ≤  ATCG  < 12 ≤  ATCG  < 25 ≤  ATCG  < 37 ≤  ATCG