Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,256,831 C→T 6.0% E69E (GAG→GAA rpoC ← RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,256,8310CT 6.0% 219.7 / 3.9 66E69E (GAG→GAArpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base C (29/33);  new base T (2/2);  total (31/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

CACTTCAACGCCGCACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGTTCTGGTTTGAACGTACGG  >  minE/2256762‑2256900
                                                                     |                                                                     
cACTTCAACGCCGCACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACt                                                                      >  1:465325/1‑71 (MQ=255)
cACTTCAACGCCGCACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACt                                                                      >  1:2773609/1‑71 (MQ=255)
 aCTTCAACGCCGCACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTc                                                                     >  1:4043579/1‑71 (MQ=255)
  cTTCGACGCCGCACTTTTCGCAGATCACGCCGCGGTGCTTCAGGCGCTTGTACTTGCCGCACAGGCATTCg                                                                    >  1:4149514/1‑71 (MQ=255)
  cTTCGACGCCGCACTTTTCGCAGATCACGCCGCGGTGCTTCAGGCGCTTGTACTTGCCGCACAGGCATTCg                                                                    >  1:3602722/1‑71 (MQ=255)
         gccgcACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCtt                                                             <  1:3641170/71‑1 (MQ=255)
            gcACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTa                                                          <  1:721655/71‑1 (MQ=255)
            gcACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTa                                                          <  1:267855/71‑1 (MQ=255)
            gcACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTa                                                          <  1:2852369/71‑1 (MQ=255)
             cACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTAc                                                         <  1:155797/71‑1 (MQ=255)
              aCTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTAcc                                                        <  1:3956610/71‑1 (MQ=255)
               cTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCttt                                                            <  1:4324115/66‑1 (MQ=255)
                 tctcACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTa                                                          >  1:634881/1‑66 (MQ=255)
                 tctcACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTa                                                          >  1:3535335/1‑66 (MQ=255)
                  ctcACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTACCgga                                                     >  1:2215632/1‑69 (MQ=255)
                  ctcACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTACCgg                                                      <  1:3844695/69‑1 (MQ=255)
                     acaGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTAc                                                         >  1:2744608/1‑63 (MQ=255)
                     acaGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTACCGGCCCa                                                  >  1:1594818/1‑70 (MQ=255)
                      caGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTACCGGccc                                                   <  1:235584/68‑1 (MQ=255)
                      caGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTACCGGCCCaaa                                                <  1:4395260/71‑1 (MQ=255)
                       aGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTa                                                          >  1:2945820/1‑60 (MQ=255)
                       aGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTa                                                          >  1:240790/1‑60 (MQ=255)
                       aGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTACCGGCCCaaa                                                <  1:1589896/70‑1 (MQ=255)
                        gATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTACCGGCCCa                                                  >  1:39560/1‑67 (MQ=255)
                        gATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTACCGGCCCa                                                  >  1:2564452/1‑67 (MQ=255)
                        gATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTACCGGCCCa                                                  >  1:2587188/1‑67 (MQ=255)
                        gATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTACCGGCCCAAAGa                                              <  1:1145780/71‑1 (MQ=255)
                        gATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTACCGGCCCAAAGa                                              <  1:2046605/71‑1 (MQ=255)
                        gATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTACCGGCCCAAAGa                                              <  1:2041196/71‑1 (MQ=255)
                         aTGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTa                                                          >  1:963957/1‑58 (MQ=255)
                         aTGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTa                                                          >  1:2376058/1‑58 (MQ=255)
                                cACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCTTAATCTTTTACCGGCCCAAAGATACGGgcg                                      <  1:1673729/71‑1 (MQ=255)
                                cACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTACCGGCCCAAAGATACGGgcg                                      <  1:1521080/71‑1 (MQ=255)
                                 aCGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTa                                                                  <  1:3575386/42‑1 (MQ=255)
                                  cGGTGCTTCAGGCGCTTGTACTTGCCGCACAGGCATTCGTAGTCCTTGATCGGGCCAAAGATCTTGGCGCa                                    <  1:1482452/71‑1 (MQ=255)
                                   ggTGCTTCAGGCGCTTGTACTTGCCGCACAGGCATTCGTAGTCCTTGATCGGGCCAAAGATCTTGGCGCa                                    <  1:3895856/70‑1 (MQ=255)
                                     tgtTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTACCGGCCCAAAGATACGGGCGCAGaa                                 <  1:891174/71‑1 (MQ=255)
                                     tgtTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTACCGGCCCAAAGATACGGGCGCAGaa                                 <  1:3884487/71‑1 (MQ=255)
                                       tttCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTACCGGCCCAAAGATACGGGCGCAGaaaa                               <  1:2673864/71‑1 (MQ=255)
                                       tttCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTACCGGCCCAAAGATACGGGCGCAGaaa                                <  1:4409524/70‑1 (MQ=255)
                                          cAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTACCGGCCCAAAGATa                                            <  1:4054929/55‑1 (MQ=255)
                                          cAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTACCGGCCCAAAGATa                                            <  1:1226017/55‑1 (MQ=255)
                                            ggCGCTTGTACTTGCCGCACAGGCACTCGTAGTCCTTCACCGGGCCGAAGATCTTGGCGCAGAACAGGCCg                          >  1:547731/1‑71 (MQ=255)
                                            ggCGCTTGTACTTACCGCACAGGCTCTCGTAATCTTTTACCGGCCCAAAGATAc                                           >  1:298424/1‑54 (MQ=255)
                                                  tGTACTTACCGCACAGGCACTCGTAATCTTTTACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGt                    >  1:1622404/1‑71 (MQ=255)
                                                   gTACTTACCGCACAGGCACTCGTAATCTTTTACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCa                       >  1:258386/1‑67 (MQ=255)
                                                   gTACTTACCGCACAGGCACTCGTAATCTTTTACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGtt                   >  1:739846/1‑71 (MQ=255)
                                                   gTACTTACCGCACAGGCACTCGTAATCTTTTACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGt                    >  1:3443924/1‑70 (MQ=255)
                                                     aCTTACCGCACAGGCACTCGTAATCTTTTACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCAc                      >  1:3716180/1‑66 (MQ=255)
                                                     aCTTACCGCACAGGCACTCGTAATCTTTTACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCAc                      >  1:1136026/1‑66 (MQ=255)
                                                       ttACCGCACAGGCACTCGTAATCTTTTACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGTTCTgg               >  1:3562054/1‑71 (MQ=255)
                                                        tACCGCACAGGCACTCGTAATCTTTTACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGTTCTGGt              >  1:2436464/1‑71 (MQ=255)
                                                          ccGCACAGGCACTCGTAATCTTTTACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACTTTCTGGtt             <  1:2982129/70‑1 (MQ=255)
                                                          ccGCACAGGCACTCGTAATCTTTTACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGTTCTGGtt             >  1:1429869/1‑70 (MQ=255)
                                                            gCACAGGCACTCGTAATCTTTTACCGGCCCAAAGATAc                                           >  1:4226562/1‑38 (MQ=255)
                                                            gCACAGGCACTCGTAATCTTTTACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGtt                   >  1:265228/1‑62 (MQ=255)
                                                            gCACAGGCACTCGTAATCTTTTACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGTTCTGGtt             >  1:2468014/1‑68 (MQ=255)
                                                             cacaGGCACTCGTAATCTTTTACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGTTCTGGt              <  1:445381/66‑1 (MQ=255)
                                                             cacaGGCACTCGTAATCTTTTACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGTTCTGGt              <  1:1637287/66‑1 (MQ=255)
                                                             cacaGGCACTCGTAATCTTTTACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGTTCTGGTTTGaa         <  1:123186/71‑1 (MQ=255)
                                                                 ggCACTCGTAATCTTTTACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGTTCTGGTTTGAACGTa     <  1:1505071/71‑1 (MQ=255)
                                                                 ggCACTCGTAATCTTTTACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGTTCTGGTTTGAACGTa     <  1:395793/71‑1 (MQ=255)
                                                                  gCACTCGTAATCTTTTACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGTTCTGGTTTGAACGTAc    <  1:3306271/71‑1 (MQ=255)
                                                                   cACTCGTAATCTTTTACCGGCCCAAAGATACGGGCGCAGAAAAGGcc                           <  1:1704308/47‑1 (MQ=255)
                                                                    aCTCGTAATCTTTTACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGTTCTGGTTTGAACGTACgg  <  1:1487578/71‑1 (MQ=255)
                                                                    aCTCGTAATCTTTTACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGTTCTGGTTTGAACGTACg   <  1:4181078/70‑1 (MQ=255)
                                                                     |                                                                     
CACTTCAACGCCGCACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGTTCTGGTTTGAACGTACGG  >  minE/2256762‑2256900

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: