Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,256,819 A→G 7.6% G73G (GGT→GGC rpoC ← RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,256,8190AG 7.6% 187.7 / 6.8 66G73G (GGT→GGCrpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base A (30/31);  new base G (3/2);  total (33/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.80e-01

CTTTAGTCTGGGTCACTTCAACGCCGCACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGTTCTGGT  >  minE/2256749‑2256888
                                                                      |                                                                     
cTTTAGTCTGGGTCACTTCAACGCCGCACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTa                                                                       >  1:956842/1‑71 (MQ=255)
   tAGTCTGGGTCACTTCAACGCCGCACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTttacttac                                                                      <  1:2455670/69‑1 (MQ=255)
   tAGTCTGGGTCACTTCAACGCCGCACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTa                                                                       >  1:3504351/1‑68 (MQ=255)
   tAGTCTGGGTCACTTCAACGCCGCACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTa                                                                       >  1:1893311/1‑68 (MQ=255)
   tAGTCTGGGTCACTTCAACGCCGCACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTa                                                                       >  1:2355965/1‑68 (MQ=255)
   tAGTCTGGGTCACTTCAACGCCGCACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTAc                                                                      <  1:2859699/69‑1 (MQ=255)
   tAGTCTGGGTCACTTCAACGCCGCACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCg                                                                    >  1:833099/1‑71 (MQ=255)
   tAGTCTGGGTCACTTCAACGCCGCACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCg                                                                    >  1:2964687/1‑71 (MQ=255)
   tAGTCTGGGTCACTTCAACGCCGCACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACATACCg                                                                    >  1:1261062/1‑71 (MQ=255)
     gTCTGGGTCACTTCAACGCCGCACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGca                                                                  <  1:4183438/71‑1 (MQ=255)
      tCTGGGTCACTTCAACGCCGCACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGcac                                                                 <  1:2688547/71‑1 (MQ=255)
        tGGGTCACTTCAACGCCGCACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAg                                                               >  1:3574906/1‑71 (MQ=255)
          ggTCACTTCAACGCCGCACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAg                                                               <  1:1342881/69‑1 (MQ=255)
             cACTTCAACGCCGCACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACt                                                          >  1:2773609/1‑71 (MQ=255)
             cACTTCAACGCCGCACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACt                                                          >  1:465325/1‑71 (MQ=255)
              aCTTCAACGCCGCACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTc                                                         >  1:4043579/1‑71 (MQ=255)
               cTTCGACGCCGCACTTTTCGCAGATCACGCCGCGGTGCTTCAGGCGCTTGTACTTGCCGCACAGGCATTCg                                                        >  1:3602722/1‑71 (MQ=255)
               cTTCGACGCCGCACTTTTCGCAGATCACGCCGCGGTGCTTCAGGCGCTTGTACTTGCCGCACAGGCATTCg                                                        >  1:4149514/1‑71 (MQ=255)
                ttCAACGCCGCACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGca                                                                  <  1:4205750/60‑1 (MQ=255)
                ttCAACGCCGCACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGca                                                                  <  1:292678/60‑1 (MQ=255)
                ttCAACGCCGCACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGca                                                                  <  1:2852610/60‑1 (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)
                             ttCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGcac                                                                 <  1:2212390/48‑1 (MQ=255)
                              tctcACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGc                                                                   >  1:1482223/1‑45 (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)
                                     gATGACGCCACGGTGTTTCAGGCGCTTGTACTTAc                                                                      >  1:3654770/1‑35 (MQ=255)
                                     gATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTACCGGCCCa                                      >  1:2587188/1‑67 (MQ=255)
                                     gATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTACCGGCCCa                                      >  1:2564452/1‑67 (MQ=255)
                                     gATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTACCGGCCCa                                      >  1:39560/1‑67 (MQ=255)
                                     gATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTACCGGCCCAAAGa                                  <  1:1145780/71‑1 (MQ=255)
                                     gATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTACCGGCCCAAAGa                                  <  1:2041196/71‑1 (MQ=255)
                                     gATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTACCGGCCCAAAGa                                  <  1:2046605/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:3884487/71‑1 (MQ=255)
                                                  tgtTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTACCGGCCCAAAGATACGGGCGCAGaa                     <  1:891174/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)
                                                                      |                                                                     
CTTTAGTCTGGGTCACTTCAACGCCGCACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGTTCTGGT  >  minE/2256749‑2256888

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: