Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,253,075 G→C 31.3% S1321S (TCC→TCG rpoC ← RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,253,0750GC31.3% 88.0 / 57.7 67S1321S (TCC→TCGrpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base G (29/17);  new base C (13/8);  total (42/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.21e-01

TCGTCGCGTTTGCCCGCAACGGCTGCTTCGGTCAGCACGCGAGTGGTCTCCTGGAACGATGCCGCGGAGATGAAGGACTCGGTTGCCAGAGACGCTTTGGTGATACCCAGCAGATCGCGGGAGTAAGTTG  >  minE/2253010‑2253139
                                                                 |                                                                
tcgtcgCGTTTGCCCGCAACGGCTGCTTCGGTCAGCACGCGAGTGGTCTCCTGGAACGATGCCGCGGAGAt                                                             >  1:446644/1‑71 (MQ=255)
  gtcgCGTTTGCCCGCAACGGCTGCTTCGGTCAGCACGCGAGTGGTCTCCTGGAACGATGCCGCGGAGATGa                                                           <  1:4677691/71‑1 (MQ=255)
  gtcgCGTTTGCCCGCAACGGCTGCTTCGGTCAGCACGCGAGTGGTCTCCTGGAACGATGCCGCGGAGATGa                                                           <  1:3792987/71‑1 (MQ=255)
          tggccGCAACGGCTGCTTCGGTCAGCACGCGAGTGGTGTCCTGGAACGATGCCGCGGAGATGAAGGACTc                                                    <  1:4065730/67‑1 (MQ=255)
           gCCCGCAACGGCTGCTTCGGTCAGCACGCGAGTGGTCTCCTGGAACGATGCCGCGGAGATGAAGGACTCg                                                   <  1:1987193/70‑1 (MQ=255)
             ccGCAACGGCTGCTTCGGTCAGCACGCGAGTGGTCTCCTGGAACGATGCCGCGGAGATGAAGGACTCGGtt                                                >  1:1041079/1‑71 (MQ=255)
             ccGCAACGGCTGCTTCGGTCAGCACGCGAGTGGTCTCCTGGAACGATGCCGCGGAGATGAAGGACTCGGtt                                                >  1:4323849/1‑71 (MQ=255)
                           tCGGTCAGCACGCGAGTGGTCTCCTGGAACGATGCCGCGGAGATGAAGGACTCGGTTGCCAGAGACGCttt                                  >  1:3491352/1‑71 (MQ=255)
                            cGGTCAGCACGCGAGTGGTCTCCTGGAACGATGCCGCGGAGATGAAggt                                                       >  1:3801863/1‑48 (MQ=255)
                            cGGTCAGCACGCGAGTGGTCTCCTGGAACGATGCCGCGGAGATGAAGGACTCGGTTGCCAGAGACGCtt                                   <  1:1634887/69‑1 (MQ=255)
                             ggTCAGCACGCGAGTGGTCTCCTGGAACGATGCCGCGGAGATGAAGGACTCGGTTGCCAGAGACGCTTTgg                                <  1:1251468/71‑1 (MQ=255)
                              gTCAGCACGCGAGTGGTCTCCTGGAACGATGCCGCGGAGATGAAGGACTCGGTTGCCAGAGACGCTTTGGt                               >  1:2425062/1‑71 (MQ=255)
                                 aGCACGGGGGTGGTTTCCTGGAACGATGCGGCCGAGATGaac                                                         <  1:2737307/42‑2 (MQ=255)
                                 aGCACGCGGGTGGTTTCCTGGAACGATGCGGCCGAGATGaac                                                         <  1:2295718/42‑2 (MQ=255)
                                 aGCACGCGGGTGGTTTCCTGGAACGATGCGGCCGAGATGaac                                                         <  1:2681639/42‑2 (MQ=255)
                                 aGCACGCGGGTGGTTTCCTGGAACGATGCGGCCGAGATGaac                                                         <  1:4302292/42‑2 (MQ=255)
                                 aGCACGCGGGTGGTTTCCTGGAACGATGCGGCCGAGATGaac                                                         <  1:4396902/42‑2 (MQ=255)
                                 aGCACGCGGGTGGTTTCCTGGAACGATGCGGCCGAGATGaac                                                         <  1:1256645/42‑2 (MQ=255)
                                 aGCACGCGGGTGGTTTCCTGGAACGATGCGGCCGAGATGaac                                                         <  1:1253965/42‑2 (MQ=255)
                                 aGCACGCGGGTGGTTTCCTGGAACGATGCGGCCGAGATGaac                                                         <  1:1224086/42‑2 (MQ=255)
                                 aGCACGCGAGTGGTCTCCTGGAACGATGCCGCGGAGATGAAGGACTCGGTTGCCAGAGACGCTTTGGTGa                             >  1:3678840/1‑70 (MQ=255)
                                 aGCACGCGAGTGGTCTCCTGGAACGATGCCGCGGAGATGAAGGACTCGGTTGCCAGAGACGCTTTGGTGa                             >  1:2187203/1‑70 (MQ=255)
                                 aGCACGCGAGTGGTCTCCTGGAACGATGCCGCGGAGATGAAGGACTCGGTTGCCAGAGACGCTTTGGTGa                             >  1:4706076/1‑70 (MQ=255)
                                 aGCACGCGAGTGGTCTCCTGGAACGATGCCGCGGAGATGAACGACTCGGTTGCCa                                            >  1:2363949/1‑55 (MQ=255)
                                  gCACGCGAGTGGTCTCCTGGAACGATGCCGCGGAGATGAAGGACTCGGTTGCCAGAGACGCTTTGGTGATa                           <  1:638836/71‑1 (MQ=255)
                                  gCACGCGAGTGGTCTCCTGGAACGATGCCGCGGAGATGAAGGACTCGGTTGCCAGAGACGCTTTGGTGATa                           <  1:1362392/71‑1 (MQ=255)
                                  gCACGCGAGTGGTCTCCTGGAACGATGCCGCGGAGATGAAGGACTCGGTTGCCAGAGACGCTTTGGTGATa                           <  1:3874895/71‑1 (MQ=255)
                                  gCACGCGAGTGGTCTCCTGGAACGATGCCGCGGAGATGAAGGACTCGGTTGCCAGAGACGCTTTGGTGATa                           <  1:845628/71‑1 (MQ=255)
                                  gCACGCGAGTGGTCTCCTGGAACGATGCCGCGGAGATGAAGGACTCGGTTGCCAGAGACGCTTTGGTGATa                           <  1:3052781/71‑1 (MQ=255)
                                   cACGCGGGTGGTTTCCTGGAACGATGCGGCCGAGATGAACGa                                                       >  1:838693/1‑42 (MQ=255)
                                   cACGCGGGTGGTTTCCTGGAACGATGCGGCCGAGATGAACGa                                                       >  1:3768434/1‑42 (MQ=255)
                                   cACGCGGGTGGTTTCCTGGAACGATGCGGCCGAGATGAACGa                                                       >  1:860139/1‑42 (MQ=255)
                                   cACGCGGGTGGTTTCCTGGAACGATGCGGCCGAGATGAACGa                                                       >  1:3801056/1‑42 (MQ=255)
                                   cACGCGGGTGGTTTCCTGGAACGATGCGGCCGAGATGAACGa                                                       >  1:58809/1‑42 (MQ=255)
                                   cACGCGGGTGGTTTCCTGGAACGATGCGGCCGAGATGAACGa                                                       >  1:4950113/1‑42 (MQ=255)
                                   cACGCGGGTGGTTTCCTGGAACGATGCGGCCGAGATGAACGa                                                       >  1:1059048/1‑42 (MQ=255)
                                   cACGCGGGTGGTTTCCTGGAACGATGCGGCCGAGATGAACGa                                                       >  1:357429/1‑42 (MQ=255)
                                   cACGCGGGTGGTTTCCTGGAACGATGCGGCCGAGATGAACGa                                                       >  1:2832008/1‑42 (MQ=255)
                                   cACGCGGGTGGTTTCCTGGAACGATGCGGCCGAGATGAACGa                                                       >  1:1504887/1‑42 (MQ=255)
                                   cACGCGGGTGGTTTCCTGGAACGATGCGGCCGAGATGAACGa                                                       >  1:1985820/1‑42 (MQ=255)
                                   cACGCGGGTGGTTTCCTGGAACGATGCGGCCGAGATGAACGa                                                       >  1:2244842/1‑42 (MQ=255)
                                   cACGCGGGTGGTTTCCTGGAACGACGCGGCCGAGATGAACGa                                                       >  1:4894285/1‑42 (MQ=16)
                                   cACGCGAGTGGTCTCCTGGAACGATGCCGCGGAGATGAAGGACTCGGTTGCCAGAGACGCTTTGGTGATa                           <  1:2065545/70‑1 (MQ=255)
                                   cACGCGAGTGGTCTCCTGGAACGATGCCGCGGAGATGAAGGACTCGGTTGCCAGAGACGCTTTGGTGATAc                          >  1:4366853/1‑71 (MQ=255)
                                               ctcCTGGAACGATGCCGCGGAGATGAAGGACTCGGTTGCCAGAGACGCTTTGGTGATACCCAGCAGATcgc              >  1:1149570/1‑71 (MQ=255)
                                               ctcCTGGAACGATGCCGCGGAGATGAAGGACTCGGTTGCCAGAGACGCTTTGGTGATACCCAGCAGATcgc              >  1:127303/1‑71 (MQ=255)
                                                  cTGGAACGATGCCGCGGAGATGAAGGACTCGGTTGCCAGAGACGCtt                                   <  1:3973221/47‑1 (MQ=255)
                                                  cTGGAACGATGCCGCGGAGATGAAGGACTCGGTTGCCAGAGACGCt                                    >  1:2139130/1‑46 (MQ=255)
                                                  cTGGAACGATGCCGCGGAGATGAAGGACTCGGTTGCCAGAGACGCTTTGGt                               >  1:877257/1‑51 (MQ=255)
                                                  cTGGAACGATGCCGCGGAGATGAAGGACTCGGTTGCCAGAGACGCTTTGGt                               >  1:4074258/1‑51 (MQ=255)
                                                  cTGGAACGATGCCGCGGAGATGAAGGACTCGGTTGCCAGAGACGCTTTGGt                               >  1:4051475/1‑51 (MQ=255)
                                                  cTGGAACGATGCCGCGGAGATGAAGGACTCGGTTGCCAGAGACGCTTTGGTGATACCCAGCAGATCGCgg            <  1:1870147/70‑1 (MQ=255)
                                                  cTGGAACGATGCCGCGGAGATGAAGGACTCGGTTGCCAGAGACGCTTTGGTGATACCCAGCAGATCGCgg            >  1:3199377/1‑70 (MQ=255)
                                                  cTGGAACGATGCCGCGGAGATGAAGGACTCGGTTGCCAGAGACGCTTTGGTGATACCCAGCAGATCGCgg            <  1:4253263/70‑1 (MQ=255)
                                                     gAACGATGCCGCGGAGATGAAGGACTCGGTTGCCAGAGACGCTTTGGTGAt                            >  1:3636370/1‑51 (MQ=255)
                                                     gAACGATGCCGCGGAGATGAAGGACTCGGTTGCCAGAGACGCTTTGGTGAt                            >  1:2321255/1‑51 (MQ=255)
                                                     gAACGATGCCGCGGAGATGAAGGACTCGGTTGCCAGAGACGCTTTGGTGAt                            >  1:1778829/1‑51 (MQ=255)
                                                     gAACGATGCCGCGGAGATGAAGGACTCGGTTGCCAGAGACGCTTTGGTGAt                            >  1:3847572/1‑51 (MQ=255)
                                                     gAACGATGCCGCGGAGATGAAGGACTCGGTTGCCAGAGACGCTTTGGTGATACCCAGCAGATCGCGGGAg         >  1:1242992/1‑70 (MQ=255)
                                                       aCGCTGCCGCGGAGATGAAGGACTCGGTTGCCAgaga                                        >  1:2061427/1‑37 (MQ=255)
                                                       aCGATGCCGCGGAGATGAAGGACTCGGTTGCCaga                                          >  1:1988304/1‑35 (MQ=255)
                                                       aCGATGCCGCGGAGATGAAGGACTCGGTTGCCAgaga                                        >  1:5085138/1‑37 (MQ=255)
                                                       aCGATGCCGCGGAGATGAAGGACTCGGTTGCCAgaga                                        >  1:705834/1‑37 (MQ=255)
                                                       aCGATGCCGCGGAGATGAAGGACTCGGTTGCCAgaga                                        >  1:5100485/1‑37 (MQ=255)
                                                       aCGATGCCGCGGAGATGAAGGACTCGGTTGCCAgaga                                        >  1:4865944/1‑37 (MQ=255)
                                                           tGCCGCGGAGATGAAGGACTCGGTTGCCAGAGATGCTTTGGTGATACCCAGCAGATCGCGGGAGTAAGTTg  <  1:2661028/71‑1 (MQ=255)
                                                           tGCCGCGGAGATGAAGGACTCCGTTGCCAGAGACGCTTTGGTGATACCCAGCAGATCGCGGGAGTAAGTTg  <  1:408744/71‑1 (MQ=255)
                                                               gcgGAGATGAAGGACTCGGTTGCCAGAGACGCTTTGGt                               <  1:1682570/38‑1 (MQ=255)
                                                                cgGAGATGAAGGACTCGGTTGCCAGAGACGCTTTGGt                               <  1:5145576/37‑1 (MQ=255)
                                                                 |                                                                
TCGTCGCGTTTGCCCGCAACGGCTGCTTCGGTCAGCACGCGAGTGGTCTCCTGGAACGATGCCGCGGAGATGAAGGACTCGGTTGCCAGAGACGCTTTGGTGATACCCAGCAGATCGCGGGAGTAAGTTG  >  minE/2253010‑2253139

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: