Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 209,697 C→G 18.7% A577G (GCG→GGG)  dnaE → DNA polymerase III subunit alpha

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606209,6970C→G18.7% 2.4 ‑3.1 0.0251730dnaEDNA polymerase III subunit alpha
Reads supporting (aligned to +/- strand):  new base (5/0):  ref base (9/11):  total (14/11)
Fisher's exact test strand distribution p-value = 4.64e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

CAACTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGAGCCGCCGCTGGATATCGCCGC  >  REL606/209662‑209730
                                   |                                 
caACTGGGCGCTGGAGATGATCAACAAGCggcgggc                                   <  1:2551172/36‑1 (MQ=255)
 aaCTGGGCGCTGGAGATGATCAACAAGCggcgggcg                                  >  1:6252228/1‑36 (MQ=255)
   cTGGGCGCTGGAGATGATCAACAAGCGGCGGGCgaa                                <  1:2413108/36‑1 (MQ=255)
   cTGGGCGCTGGAGATGATCAACAAGCGGCGGGCgaa                                <  1:2478146/36‑1 (MQ=255)
   cTGGGCGCTGGAGATGATCAACAAGCGGCGGGCgaa                                <  1:7224356/36‑1 (MQ=255)
    tGGGCGCTGGAGATGATCAACAAGCGGCGGGCgaag                               <  1:7476327/36‑1 (MQ=255)
    tGGGCGCTGGAGATGATCAACAAGCGGCGGGCgaag                               <  1:485455/36‑1 (MQ=255)
       gcgcTGGAGATGATCAACAAGCGGCGGGGgaagaaa                            >  1:1023894/1‑35 (MQ=255)
       gcgcTGGAGATGATCAACAAGCGGCGGGCGAAGAAt                            <  1:7410989/36‑1 (MQ=255)
       gcgcTGGAGATGATCAACAAGCGGCGGGCGAAGAAt                            <  1:5632856/36‑1 (MQ=255)
       gcgcTGGAGATGATCAACAAGCGGCGGGCGAAGAAt                            <  1:4347422/36‑1 (MQ=255)
         gcTGGAGATGATCAACAAGCGGCGGGGGGAGAAAgg                          >  1:2846702/1‑36 (MQ=25)
            ggAGATGATCAACAAGCGGCGGGGGGAGAATTGCgg                       >  1:6760756/1‑35 (MQ=21)
            ggAGATGATCAACAAGCGGCGGGCGGAGGATGGCGa                       >  1:7775417/1‑36 (MQ=38)
             gagaTGATCAACAAGCGGCGGGCGAAGAATGGCGAg                      <  1:6269156/36‑1 (MQ=255)
             gagaTGATCAACAAGCGGCGGGCGAAGAATGGCGAg                      >  1:530865/1‑36 (MQ=255)
             gagaTGATCAACAAGCGGCGGGCGAAGAATGGCGAg                      <  1:7060626/36‑1 (MQ=255)
             gagaTGATCAACAAGCGCCGGGCGCAGAAAGGCGAg                      >  1:1109725/1‑36 (MQ=21)
              agaTGATCAACAAGCGGCGGGCGGAGAATGGCGAgc                     >  1:4247431/1‑36 (MQ=255)
               gatgatCAACAAGCGGCGGGCGGAGAATGGGGAggc                    >  1:4740122/1‑34 (MQ=25)
                atgatCAACAAGCGGCGGGGGCGGAAAGGCGGgccg                   >  1:5363834/1‑36 (MQ=16)
                  gatCAACAAGCGGCGGGCGGAGAATGGCGAggcgcc                 >  1:3669306/1‑36 (MQ=37)
                  gatCAACAAGCGGCGGGCGAAGAATGGCGAggcgcc                 >  1:1016095/1‑36 (MQ=255)
                  gatCAACAAGCGGCGGGCGAAGAATGGCGAgccgcc                 <  1:8190446/36‑1 (MQ=255)
                   atCAACAAGCGGCGGGGGAAGAATGGCGAGgcgccg                >  1:2858405/1‑36 (MQ=37)
                    tCAACAAGCGGCGGGCGAAGAATGGCGAGCcgccgc               >  1:6615196/1‑36 (MQ=255)
                      aacaaGCGGCGGGCGAAGAAAGGGGGGCcgccgccg             >  1:3162757/1‑34 (MQ=21)
                             ggcgggcgAAGAATGGCGAGCCGCCGCTGGATATcg      <  1:4125722/36‑1 (MQ=255)
                             ggcgggcgAAGAATGGCGAGCCGCCGCTGGATATcg      >  1:8354637/1‑36 (MQ=255)
                             ggcgggcgAAGAATGGCGAGCCGCCGCCGGATATcg      >  1:6508381/1‑36 (MQ=255)
                              gcgggagAAGAATGGCGAGCCGCCGCTGGATATcgc     <  1:7025265/36‑1 (MQ=255)
                               cgggcgAAGAATGGCGAGCCGCCGCTGGATATcgcc    >  1:6346460/1‑36 (MQ=255)
                                gggcgAAGAATGGCGAGCCGCCGCTGGATATcgccg   >  1:3577259/1‑36 (MQ=255)
                                 ggcgAAGAATGGCGAGCCGCCGCTGGATATcgccgc  >  1:4536242/1‑36 (MQ=255)
                                 ggcgAAGAATGCCGCCCCCCCGCGGGATATcgccgc  <  1:6611687/36‑1 (MQ=16)
                                   |                                 
CAACTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGAGCCGCCGCTGGATATCGCCGC  >  REL606/209662‑209730

Base quality scores:  ATCG  < 3 ≤  ATCG  < 4 ≤  ATCG  < 14 ≤  ATCG  < 27 ≤  ATCG  < 37 ≤  ATCG