Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 186,734 C→A 34.6% A928A (GCC→GCA dnaE → DNA polymerase III alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE186,7340CA34.6% 74.8 / 52.4 75A928A (GCC→GCAdnaEDNA polymerase III alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (38/11);  new base A (7/19);  total (45/30)
Fisher's exact test for biased strand distribution p-value = 4.43e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.89e-01

TAAAAGCGGCAGATCAACACGCGAAAGCGGAAGCTATCGGTCAGGCCGATATGTTCGGCGTGCTGGCCGAAGAGCCGGAACAAATTGAACAATCCTACGCCAGCTGCCAACCGTGGCCGGAGCAGGTGGTATTA  >  minE/186667‑186800
                                                                   |                                                                  
tAAAAGCGGCAGATCAACACGCGAAAGCGGAAGCTATCGGTCAGGCCGATATGTTCGGCGTGCTGGCCg                                                                   <  1:659784/69‑1 (MQ=255)
                  acGCGAAAGCTGAAGCTATCGGTCAGGCCGATATGTTCGGCGTGCTGGCCGAAGAGCCGGAACAAATTg                                                 <  1:1160730/69‑1 (MQ=255)
                  acGCGAAAGCGGAAGCTATCGGTCAGGCCGATATGTTCGGCGTGCTGGCCGAAGAGCCGGAACAAATTg                                                 <  1:1069809/69‑1 (MQ=255)
                  acGCGAAAGCGGAAGCTATCGGTCAGGCCGATATGTTCGGCGTGCTGGCCGAAGAGCCGGAACAAATTg                                                 <  1:560594/69‑1 (MQ=255)
                  acGCGAAAGCGGAAGCTATCGGTCAGGCCGATATGTTCGGCGTGCTGGCCGAAGAGCCGGAACAAATTg                                                 <  1:343719/69‑1 (MQ=255)
                  acGCGAAAGCGGAAGCTATCGGTCAGGCCGATATGTTCGGCGTGCTGGCCGAAGAGCCGGAACAAATTg                                                 <  1:2720474/69‑1 (MQ=255)
                      gAAAGCAGAAGCCATCGGTCAGGCCGATATGTTCGGCGTGCTGGCAGAAGAGCCGGAACaa                                                     <  1:1857985/61‑1 (MQ=255)
                      gAAAGCAGAAGCCATCGGTCAGGCCGATATGTTCGGCGTGCTGGCAGAAGAGCCGGAACaa                                                     <  1:2556696/61‑1 (MQ=255)
                      gAAAGCAGAAGCCATCGGTCAGGCCGATATGTTCGGCGTGCTGGCAGAAGAGCCGGAACaa                                                     <  1:1267831/61‑1 (MQ=255)
                      gAAAGCAGAAGCCATCGGTCAGGCCGATATGTTCGGCGTGCTGGCAGAAGAGCCGGAACaa                                                     <  1:2202044/61‑1 (MQ=255)
                      gAAAGCAGAAGCCATCGGTCAGGCCGATATGTTCGGCGTGCTGGCAGAAGAGCCGGAACaa                                                     <  1:2856232/61‑1 (MQ=255)
                      gAAAGCAGAAGCCATCGGTCAGGCCGATATGTTCGGCGTGCTGGCAGAAGAGCCGGAACaa                                                     <  1:2918423/61‑1 (MQ=255)
                      gAAAGCAGAAGCCATCGGTCAGGCCGATATGTTCGGCGTGCTGGCAGAAGAGCCGGAACaa                                                     <  1:2961428/61‑1 (MQ=255)
                      gAAAGCAGAAGCCATCGGTCAGGCCGATATGTTCGGCGTGCTGGCAGAAGAGCCGGAACaa                                                     <  1:2065326/61‑1 (MQ=255)
                      gAAAGCAGAAGCCATCGGTCAGGCCGATATGTTCGGCGTGCTGGCAGAAGAGCCGGAACaa                                                     <  1:2967210/61‑1 (MQ=255)
                      gAAAGCAGAAGCCATCGGTCAGGCCGATATGTTCGGCGTGCTGGCAGAAGAGCCGGAACaa                                                     <  1:1908746/61‑1 (MQ=255)
                      gAAAGCAGAAGCCATCGGTCAGGCCGATATGTTCGGCGTGCTGGCAGAAGAGCCGGAACaa                                                     <  1:3134202/61‑1 (MQ=255)
                      gAAAGCAGAAGCCATCGGTCAGGCCGATATGTTCGGCGTGCTGGCAGAAGAGCCGGAACaa                                                     <  1:3149931/61‑1 (MQ=255)
                      gAAAGCAGAAGCCATCGGTCAGGCCGATATGTTCGGCGTGCTGGCAGAAGAGCCGGAACaa                                                     <  1:3212905/61‑1 (MQ=255)
                      gAAAGCAGAAGCCATCGGTCAGGCCGATATGTTCGGCGTGCTGGCAGAAGAGCCGGAACaa                                                     <  1:1697520/61‑1 (MQ=255)
                      gAAAGCAGAAGCCATCGGTCAGGCCGATATGTTCGGCGTGCTGGCAGAAGAGCCGGAACaa                                                     <  1:3259176/61‑1 (MQ=255)
                      gAAAGCAGAAGCCATCGGTCAGGCCGATATGTTCGGCGTGCTGGCAGAAGAGCCGGAACaa                                                     <  1:3294557/61‑1 (MQ=255)
                      gAAAGCAGAAGCCATCGGTCAGGCCGATATGTTCGGCGTGCTGGCAGAAGAGCCGGAACaa                                                     <  1:710103/61‑1 (MQ=255)
                      gAAAGCAGAAGCCATCGGTCAGGCCGATATGTTCGGCGTGCTGGCAGAAGAGCCGGAACaa                                                     <  1:752365/61‑1 (MQ=255)
                       aaaGCAGAAGCCATCGGTCAGGCCGATATGTTCGGCGTGCTGGCAGAAGAGCCGGAACaa                                                     <  1:1233209/60‑1 (MQ=255)
                          gCGGAAGCTATCGGTCAGGCCGATATGTTCGGCGTGCTGGCCGAAGAGCCGGAACAAATTGAACAATcc                                         >  1:3123573/1‑69 (MQ=255)
                          gCGGAAGCTATCGGTCAGGCCGATATGTTCGGCGTGCTGGCCGAAGAGCCGGAACAAATTGAACAATcc                                         >  1:2948926/1‑69 (MQ=255)
                                     cGGTCAGGCCGATATGTTCGGCGTGCTGGCAGAAGAGCCGGAACAAATTGAACAAt                                           >  1:3137980/1‑56 (MQ=255)
                                     cGGTCAGGCCGATATGTTCGGCGTGCTGGCAGAAGAGCCGGAACAAATTGAACAATCCTACGCCAgtt                               >  1:2435829/1‑66 (MQ=255)
                                     cGGTCAGGCCGATATGTTCGGCGTGCTGGCAGAAGAGCCGGAACAAATTGAACAATCCTACGCCAgtt                               >  1:2264665/1‑66 (MQ=255)
                                     cGGTCAGGCCGATATGTTCGGCGTGCTGGCAGAAGAGCCGGAACAAATTGAACAATCCTACGCCAgtt                               >  1:2281890/1‑66 (MQ=255)
                                     cGGTCAGGCCGATATGTTCGGCGTGCTGGCAGAAGAGCCGGAACAAATTGAACAATCCTACGCCAgtt                               >  1:2508208/1‑66 (MQ=255)
                                     cGGTCAGGCCGATATGTTCGGCGTGCTGGCAGAAGAGCCGGAACAAATTGAACAATCCTACGCCAgtt                               >  1:143808/1‑66 (MQ=255)
                                     cGGTCAGGCCGATATGTTCGGCGTGCTGGCAGAAGAGCCGGAACAAATTGAACAATCCTACGCCAgtt                               >  1:134843/1‑66 (MQ=255)
                                      ggTCAGGCCGATATGTTCGGCGTGCTGGCCGAAGAGCCGGAACAAATTGAACAATCCTCCGCCAGCTg                              >  1:1366662/1‑68 (MQ=255)
                                      ggTCAGGCCGATATGTTCGGCGTGCTGGCCGAAGAGCCGGAACAAATTGAACAATCCTACGCCAGCTg                              >  1:958904/1‑68 (MQ=255)
                                      ggTCAGGCCGATATGTTCGGCGTGCTGGCCGAAGAGCCGGAACAAATTGAACAATCCTACGCCAGCTg                              >  1:1854601/1‑68 (MQ=255)
                                                     ttCGGCGTGCTGGCCGAAGCGCCGGAACAAATTGAACAATCCTACGCCAGCTGCCAACCGTGGc                   >  1:407098/1‑64 (MQ=255)
                                                     ttCGGCGTGCTGGCCGAAGAGCCGGAACAAATTGAACAATCCTACGCCAGCTGCCAACCGTgcc                   >  1:1574307/1‑62 (MQ=255)
                                                     ttCGGCGTGCTGGCCGAAGAGCCGGAACAAATTGAACAATCCTACGCCAGCTGCCAACCGTGGc                   >  1:2433796/1‑64 (MQ=255)
                                                     ttCGGCGTGCTGGCCGAAGAGCCGGAACAAATTGAACAATCCTACGCCAGCTGCCAACCGTGGc                   >  1:898809/1‑64 (MQ=255)
                                                     ttCGGCGTGCTGGCCGAAGAGCCGGAACAAATTGAACAATCCTACGCCAGCTGCCAACCGTGGc                   >  1:572156/1‑64 (MQ=255)
                                                     ttCGGCGTGCTGGCCGAAGAGCCGGAACAAATTGAACAATCCTACGCCAGCTGCCAACCGTGGc                   >  1:1050742/1‑64 (MQ=255)
                                                     ttCGGCGTGCTGGCCGAAGAGCCGGAACAAATTGAACAATCCTACGCCAGCTGCCAACCGTGGc                   >  1:715540/1‑64 (MQ=255)
                                                     ttCGGCGTGCTGGCCGAAGAGCCGGAACAAATTGAACAATCCTACGCCAGCTGCCAACCGTGGc                   >  1:2876299/1‑64 (MQ=255)
                                                     ttCGGCGTGCTGGCCGAAGAGCCGGAACAAATTGAACAATCCTACGCCAGCTGCCAACCGTGGCCgg                >  1:2162003/1‑67 (MQ=255)
                                                       cGGCGTGCTGGCCGAAGAGCCGGAACAAATTGAACAATCCTACGCCAGCTGCCAACCGTGGCCGGAGc             <  1:81495/68‑1 (MQ=255)
                                                       cGGCGTGCTGGCCGAAGAGCCGGAACAAATTGAACAATCCTACGCCAGCTGCCAACCGTGGCCGGAGc             <  1:2227217/68‑1 (MQ=255)
                                                       cGGCGTGCTGGCCGAAGAGCCGGAACAAATTGAACAATCCTACGCCAGCTGCCAACCGTGGCCGGAGc             <  1:1722144/68‑1 (MQ=255)
                                                            ggCTGGCCGAAGAGCCGGAAAAAATTGAACAATCCTACGCCAGCTGCCAACCGTGGCCGGAGCAggtgg       <  1:2186353/68‑1 (MQ=255)
                                                              cTGGCCGAAGAGCCGGAACAAATTGAACAATCCTACGCCAGCTGCCAACCGTGGCCGGAGCAggtggt      >  1:1518943/1‑68 (MQ=255)
                                                              cTGGCCGAAGAGCCGGAACAAATTGAACAATCCTACGCCAGCTGCCAACCGTGGCCGGAGCAggtggt      >  1:240231/1‑68 (MQ=255)
                                                              cTGGCCGAAGAGCCGGAACAAATTGAACAATCCTACGCCAGCTGCCAACCGTGGCCGGAGCAggtggt      >  1:812161/1‑68 (MQ=255)
                                                              cTGGCCGAAGAGCCGGAACAAATTGAACAATCCTACGCCAGCTGCCAACCGTGGCCGGAGCAggtggt      >  1:1203668/1‑68 (MQ=255)
                                                              cTGGCCGAAGAGCCGGAACAAATTGAACAATCCTACGCCAGCTGCCAACCGTGGCCGGAGCAggtggt      >  1:3278178/1‑68 (MQ=255)
                                                              cTGGCCGAAGAGCCGGAACAAATTGAACAATCCTACGCCAGCTGCCAACCGTGGCCGGAGCAggtggt      >  1:1683590/1‑68 (MQ=255)
                                                              cTGGCCGAAGAGCCGGAACAAATTGAACAATCCTACGCCAGCTGCCAACCGTGGCCGGAGCAggtggt      >  1:3252823/1‑68 (MQ=255)
                                                              cTGGCCGAAGAGCCGGAACAAATTGAACAATCCTACGCCAGCTGCCAACCGTGGCCGGAGCAggtggt      >  1:1429597/1‑68 (MQ=255)
                                                              cTGGCCGAAGAGCCGGAACAAATTGAACAATCCTACGCCAGCTGCCAACCGTGGCCGGAGCAggtggt      >  1:3154175/1‑68 (MQ=255)
                                                              cTGGCCGAAGAGCCGGAACAAATTGAACAATCCTACGCCAGCTGCCAACCGTGGCCGGAGCAggtggt      >  1:3033523/1‑68 (MQ=255)
                                                              cTGGCCGAAGAGCCGGAACAAATTGAACAATCCTACGCCAGCTGCCAACCGTGGCCGGAGCAggtggt      >  1:2257155/1‑68 (MQ=255)
                                                              cTGGCCGAAGAGCCGGAACAAATTGAACAATCCTACGCCAGCTGCCAACCGTGGCCGGAGCAggtggt      >  1:2603154/1‑68 (MQ=255)
                                                              cTGGCCGAAGAGCCGGAACAAATTGAACAATCCTACGCCAGCTGCCAACCGTGGCCGGAGCAggtggt      >  1:2760232/1‑68 (MQ=255)
                                                              cTGGCCGAAGAGCCGGAACAAATTGAACAATCCTACGCCAGCTGCCAACCGTGGCCGGAGCAggtggt      >  1:2858178/1‑68 (MQ=255)
                                                              cTGGCCGAAGAGCCGGAACAAATTGAACAATCCTACGCCAGCTGCCAACCGTGGCCGGAGCAggtggt      >  1:2862034/1‑68 (MQ=255)
                                                              cTGGCCGAAGAGCCGGAACAAATTGAACAATCCTACGCCAGCTGCCAACCGTGGCCGGAGCAggtggt      >  1:2172817/1‑68 (MQ=255)
                                                              cTGGCCGAAGAGCCGGAACAAATTGAACAATCCTACGCCAGCTGCCAACCGTGGCCGGAGCAggtggt      >  1:2087866/1‑68 (MQ=255)
                                                              cTGGCCGAAGAGCCGGAACAAATTGAACAATCCTACGCCAGCTGCCAACCGTGGCCGGAGCAggtggt      >  1:1968253/1‑68 (MQ=255)
                                                              cTGGCCGAAGAGCCGGAACAAATTGAACAATCCTACGCCAGCTGCCAACCGTGGCCGGAGCAggtggt      >  1:3033974/1‑68 (MQ=255)
                                                              cTGGCCGAAGAGCCGGAACAAATTGAACAATCCTACGCCAGCTGCCAACCGTGGCCGGAGCAggtggt      >  1:303686/1‑68 (MQ=255)
                                                              cTGGCCGAAGAGCCGGAACAAATTGAACAATCCTACGCCAGCTGCCAACCGTGGCCGGAGCAggtggt      >  1:3038092/1‑68 (MQ=255)
                                                              cTGGCCGAAGAGCCGGAACAAATTGAACAATCCTACGCCAGCTGCCAACCGTGGCCGGAGCAggtggt      >  1:3082137/1‑68 (MQ=255)
                                                              cTGGCCGAAGAGCCGGAACAAATTGAACAATCCTACGCCAGCTGCCAACCGTGGCCGGAGCAggtggt      >  1:1951091/1‑68 (MQ=255)
                                                              cTGGCCGAAGAGCCGGAACAAATTGAACAATCCTACGCCAGCTGCCAACCGTGGCCGGAGCAggtggt      >  1:1887841/1‑68 (MQ=255)
                                                                 gCCGAAGAGCCGGAACAAATTGAACAATCCTACGCCAGCTGCCAACCGTGGCCGGAGCAGGTGGTATTa  <  1:1848957/69‑1 (MQ=255)
                                                                   cGAAGAGCCGGAACAAATTGAACAATCCTACGCCAGCTGCCAACCGTGGCCGGAGCAGGTGGTATTa  <  1:1538384/67‑1 (MQ=255)
                                                                   cGAAGAGCCGGAACAAATTGAACAATCCTACGCCAGCTGCCAACCGTGGCCGGAGCAGGTGGTATTa  <  1:3405890/67‑1 (MQ=255)
                                                                   cGAAGAGCCGGAACAAATTGAACAATCCTACGCCAGCTGCCAACCGTGGCCGGAGCAGGTGGTATTa  <  1:3300119/67‑1 (MQ=255)
                                                                   cGAAGAGCCGGAACAAATTGAACAATCCTACGCCAGCTGCCAACCGTGGCCGGAGCAGGTGGTATTa  <  1:1497467/67‑1 (MQ=255)
                                                                   cGAAGAGCCGGAACAAATTGAACAATCCTACGCCAGCTGCCAACCGTGGCCGGAGCAGGTGGTATTa  <  1:317178/67‑1 (MQ=255)
                                                                   cGAAGAGCCGGAACAAATTGAACAATCCTACGCCAGCTGCCAACCGTGGCCGGAGCAGGTGGTACTa  <  1:763643/67‑1 (MQ=255)
                                                                   |                                                                  
TAAAAGCGGCAGATCAACACGCGAAAGCGGAAGCTATCGGTCAGGCCGATATGTTCGGCGTGCTGGCCGAAGAGCCGGAACAAATTGAACAATCCTACGCCAGCTGCCAACCGTGGCCGGAGCAGGTGGTATTA  >  minE/186667‑186800

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