Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,489,241 T→C 10.7% R80R (CGT→CGC dnaN → DNA polymerase III, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,489,2410TC10.7% 143.9 / 10.2 56R80R (CGT→CGCdnaNDNA polymerase III, beta subunit
Reads supporting (aligned to +/- strand):  ref base T (26/24);  new base C (2/4);  total (28/28)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.11e-01

GCGCTGGTTCAGCCACACGAGCCAGGAGCGACGACCGTTCCGGCGCGCAAATTCTTTGATATCTGCCGTGGTCTGCCTGAAGGCGCGGAAATTGCCGTGCAGCTGGAAGGTGAACGGATGCTGGTACGCTCCGGGCG  >  minE/2489173‑2489309
                                                                    |                                                                    
gcgcTGGTTCAGCCACACGAGCCAGGAGCGACGACCGTTCCGGCGCGCAAATTCTTTGATATCTGCCGTgg                                                                    >  1:4159789/1‑71 (MQ=255)
   cTGGTTCAGCCACACGAGCCAGGAGCGACGACCGTTCCGGCGCGCAAATTCTTTGATATCTGCCGTGGTc                                                                  <  1:43094/70‑1 (MQ=255)
        tCAGCCACACGAGCCAGGAGCGACGACCGTTCCGGCGCGCAAATTCTTTGATATCTGCCGTGGTCTGcc                                                              >  1:3715827/1‑69 (MQ=255)
         cAGCCACACGAGCCAGGAGCGACGACCGTTCCGGCGCGCAAATTCTTTGATATCTGCCGTGGTCTGCCTg                                                            <  1:779182/70‑1 (MQ=255)
         cAGCCACACGAGCCAGGAGCGACGACCGTTCCGGCGCGCAAATTCTTTGATATCTGCCGTGGTCTGCCTg                                                            <  1:646194/70‑1 (MQ=255)
         cAGCCACACGAGCCAGGAGCGACGACCGTTCCGGCGCGCAAATTCTTTGATATCTGCCGTGGTCTGCCTg                                                            <  1:1347543/70‑1 (MQ=255)
          aGCCACACGAGCCAGGAGCGACGACCGTTCCGGCGCGCAAATTCTTTGATATCTGCCGTGGTCTGCCTg                                                            <  1:2248596/69‑1 (MQ=255)
           gCCACACGAGCCAGGAGCGACGACCGTTCCGGCGCGCAAATTCTTTGATATCTGCCGTGGTCTGCCTGAAg                                                         >  1:151926/1‑71 (MQ=255)
           gCCACACGAGCCAGGAGCGACGACCGTTCCGGCGCGCAAATTCTTTGATATCTGCCGTGGTCTGCCTGAAg                                                         >  1:418294/1‑71 (MQ=255)
                 cGAGCCAGGAGCGACGACCGTTCCGGCGCGCAAATTCTTTGATATCTGCCGTGGTCTGCCTGAAGGcgcg                                                    >  1:1891825/1‑70 (MQ=255)
                  gAGCCAGGAGCGACGACCGTTCCGGCGCGCAAATTCTTTGATATCTGCCGTGGTCTGCCTGAAGGCGCGGa                                                  <  1:2490903/71‑1 (MQ=255)
                     ccAGGAGCGACGACCGTTCCGGCGCGCAAATTCTTTGATATCTGCCGTGGTCTGCCTGAAGGCGCGGAAAt                                               >  1:2366033/1‑71 (MQ=255)
                      cAGGAGCGACGACCGTTCCGGCGCGCAAATTCTTTGATATCTGCCGTGGTCTGCCTGAAGGCGCGGAAAtt                                              >  1:1426775/1‑71 (MQ=255)
                      cAGGAGCGACGACCGTTCCGGCGCGCAAATTCTTTGATATCTGCCGTGGTCTGCCTGAAGGCGCGGAAAtt                                              >  1:303584/1‑71 (MQ=255)
                         gAGCGACGACCGTTGCGGCGCGCAAATTCTTTGATATCTGCCGTGGTCTGCCTGa                                                           <  1:648713/55‑1 (MQ=255)
                         gAGCGACGACCGTTCCGGCGCGCAAATTCTTTGATATCTGCCGTGGTCTGCCTGa                                                           <  1:3667946/55‑1 (MQ=255)
                         gAGCGACGACCGTTCCGGCGCGCAAATTCTTTGATATCTGCCGTGGTCTGCCTGa                                                           <  1:319162/55‑1 (MQ=255)
                         gAGCGACGACCGTTCCGGCGCGCAAATTCTTTGATATCTGCCGTGGTCTGACTg                                                            <  1:3009964/54‑1 (MQ=255)
                          aGCGACGACCGTTCCGGCGCGCAAATTCTTTGATATCTGCCGTGGTCTGCCTGAAGGCGCGGAAATTGCCg                                          >  1:3724855/1‑71 (MQ=255)
                          aGCGACGACCGTTCCGGCGCGCAAATTCTTTGATATCTGCCGTGGTCTGCCTGAAGGCGCGGAAATTGCCg                                          >  1:1252572/1‑71 (MQ=255)
                          aGCGACGACCGTTCCGGCGCGCAAATTCTTTGATATCTGCCGTGGTCTGCCTGAAGGCGCGGAAATTGCCg                                          >  1:3907471/1‑71 (MQ=255)
                          aGCGACGACCGTTCCGGCGCGCAAATTCTTTGATATCTGCCGTGGTCTGCCTGAAGGCGCGGAAATTGCCg                                          >  1:3926920/1‑71 (MQ=255)
                           gCGACGACCGTTCCGGCGCGCAAATTCTTTGATATCTGCCGTGGTCTGCCTGAAg                                                         >  1:2060697/1‑55 (MQ=255)
                             gacgacCGTTCCGGCGCGCAAATTCTTTGATATCTGCCGTGGTCTGCCTGAAGGCGCGGAAATTGCCGTGc                                       <  1:231390/71‑1 (MQ=255)
                                     ttCCGGCGCGCAAATTCTTTGATATCTGCCGTGGTCTGCCTg                                                            <  1:1980147/42‑1 (MQ=255)
                                     ttCCGGCGCGCAAATTCTTTGATATCTGCCGTGGTCTGCCTg                                                            <  1:1886695/42‑1 (MQ=255)
                                       ccGGCGCGCAAATTCTTTGATATCTGCCGTGGTCTGCCTGAAGGCGCGGa                                                  <  1:1197965/50‑1 (MQ=255)
                                       ccGGCGCGCAAATTCTTTGATATCTGCCGTGGTCTGCCTGAAGGCGCGGAAATTGCCGTGCAGCTGGAAgg                             >  1:844761/1‑71 (MQ=255)
                                       ccGGCGCGCAAATTCTTTGATATCTGCCGTGGTCTGCCTGAAGGCGCGGAAATTGCCGTGCAGCTGGAAg                              >  1:1587671/1‑70 (MQ=255)
                                         ggCGCGCAAATTCTTTGATATCTGCCGTGGTCTGCCTGAAGGCGCGGAAATTGCCGTGCAGCTGGAAGGTg                           <  1:2392815/71‑1 (MQ=255)
                                              gcAAATTCTTTGATATCTGCCGTGGTCTGCCTGAAGGCGCGGAAATTGCCGTGCAGCTGGAAGGTGAACg                       <  1:3906644/70‑1 (MQ=255)
                                              gcAAATTCTTTGATATCTGCCGTGGTCTGCCTGAAGGCGCGGAAATTGCCGTGCAGCTGGAAGGTGAACg                       >  1:3880914/1‑70 (MQ=255)
                                                aaaTTCTTTGATATCTGCCGTGGTCTGCCTGAAGGCGCGGAAATTGCCGTGCAGCTGGAAGGTGAACGGAt                    <  1:620191/71‑1 (MQ=255)
                                                aaaTTCTTTGATATCTGCCGTGGTCTGCCTGAAGGCGCGGAAATTGCCGTGCAGCTGGAAGGTGAACGGAt                    <  1:3211840/71‑1 (MQ=255)
                                                   ttCTTTGATATTTGCCGTGGTCTGCCTGAAGGCGCGGAAATTGCCGTGCAGCTGGAAGGTGAACGGATGCt                 >  1:119195/1‑71 (MQ=255)
                                                   ttCTTTGATATCTGCCGTGGTCTGCCTGAAGGCGCGGAAATTGCCGTGCAGCTGGAAGGTGaa                         >  1:1761695/1‑63 (MQ=255)
                                                    tCTTTGATATCTGCCGTGGTCTGCCTGAAGGCGCGGAAATTGCCGTGCAGCTGGAAGGTGAACGGa                     >  1:988331/1‑66 (MQ=255)
                                                      tttGATATCTGCCGTGGTCTGCCTGAAGGCGCGGAAATTGCCGTGCAGCTGGAAGGTGAACGGATGCTgg               >  1:3309828/1‑70 (MQ=255)
                                                      tttGATATCTGCCGTGGTCTGCCTGAAGGCGCGGAAATTGCCGTGCAGCTGGAAGGTGAACGGATGCTgg               >  1:3938104/1‑70 (MQ=255)
                                                        tGATATCTGCCGTGGTCTGCCTGAAGGCGCGGAAATTGCCGTGCAGCTGGAAGGTGAACGGATGCTgg               <  1:2733676/68‑1 (MQ=255)
                                                        tGATATCTGCCGCGGTCTGCCTGAAGGCGCGGAAATTGCCGTGCAGCTGGAAGGTGAACGGATGCTGGt              <  1:385185/69‑1 (MQ=255)
                                                        tGATATCTGCCGCGGTCTGCCTGAAGGCGCGGAAATTGCCGTGCAGCTGGAAGGTGAACGGATGCTGGt              <  1:3295806/69‑1 (MQ=255)
                                                        tGATATCTGCCGCGGTCTGCCTGAAGGCGCGGAAATTGCCGTGCAGCTGGAAGGTGAACGGATGCTGGt              <  1:2015781/69‑1 (MQ=255)
                                                        tGATATCTGCCGCGGTCTGCCTGAAGGCGCGGAAATTGCCGTGCAGCTGGAAGGTGAACGGATGCTGGt              <  1:1456042/69‑1 (MQ=255)
                                                          atatCTGCCGTGGTCTGCCTGAAGGCGCGGAAATTGCCGTGCAGCTGGAAGGTGAACgg                      <  1:862757/59‑1 (MQ=255)
                                                               tGCCGTGGTCTGCCTGAAGGCGCGGAAATTGCCGTGCAGCTGGAAGGTGAACGGATGCTGGTACGCTCCgg     >  1:687061/1‑71 (MQ=255)
                                                               tGCCGTGGTCTGCCTGAAGGCGCGGAAATTGCCGTGCAGCTGGAAGGTGAACGGATGCTGGTACGCTCCgg     >  1:3918990/1‑71 (MQ=255)
                                                               tGCCGTGGTCTGCCTGAAGGCGCGGAAATTGCCGTGCAGCTGGAAGGTGAACGGATGCTGGTACGCTCCgg     >  1:2530513/1‑71 (MQ=255)
                                                               tGCCGTGGTCTGCCTGAAGGCGCGGAAATTGCCGTGCAGCTGGAAGGTGAACGGATGCTGGTACGCTCCgg     >  1:2073631/1‑71 (MQ=255)
                                                               tGCCGTGGTCTGCCTGAAGGCGCGGAAATTGCCGTGCAGCTGGAAGGTGAACGGATGCTGGTACGCTCCgg     >  1:1580315/1‑71 (MQ=255)
                                                               tGCCGCGGTCTGCCTGAAGGCGCGGAAATTGCCGTGCAGCTGGAAgg                             >  1:1032726/1‑47 (MQ=255)
                                                               tGCCGCGGTCTGCCTGAAGGCGCGGAAATTGCCGTGCAGCTGGAAgg                             >  1:1027654/1‑47 (MQ=255)
                                                                gCCGTGGTCTGCCTGAAGGCGCGGAAATTGCCGTGCAGCTGGa                                <  1:3833435/43‑1 (MQ=255)
                                                                 ccGTGGTCTGCCTGAAGGCGCGGAAATTGCCGTGCAGCTGGAAGGTGAACGGATGCTGGTACGCTCCGGGc   <  1:3774532/71‑1 (MQ=255)
                                                                  cGTGGTCTGCCTGAAGGCGCGGAAATTGCCGTGCAGCTGGAAGGTGAACGGATGCTGGTACGCTCCGGGCg  <  1:3091598/71‑1 (MQ=255)
                                                                  cGTGGTCTGCCTGAAGGCGCGGAAATTGCCGTGCAGCTGGAAGGTGAACGGATGCTGGTACGCTCCGGGCg  <  1:1474917/71‑1 (MQ=255)
                                                                    |                                                                    
GCGCTGGTTCAGCCACACGAGCCAGGAGCGACGACCGTTCCGGCGCGCAAATTCTTTGATATCTGCCGTGGTCTGCCTGAAGGCGCGGAAATTGCCGTGCAGCTGGAAGGTGAACGGATGCTGGTACGCTCCGGGCG  >  minE/2489173‑2489309

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