Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,253,725 C→T 15.1% A1105T (GCG→ACG)  rpoC ← RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,253,7250CT15.1% 232.0 / 29.1 86A1105T (GCG→ACG) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base C (38/35);  new base T (5/8);  total (43/43)
Fisher's exact test for biased strand distribution p-value = 5.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.67e-01

GATTCCTGCGGAATACGCGCCAGGGTGTCACCAGAGCTGATCTGTACGCCATCTTCCAGCTGAACAATCGCTTTACCCGGCAGGAAGTACTGCGCTGGCATATCGGTACCTGGGATCAGAACGTCGTTACCCTGAGCATCA  >  minE/2253655‑2253795
                                                                      |                                                                      
gATTCCTGCGGAATACGCGCCAGGGTGTCACCAGAGCTGATCTGTACGCCATCTTCCAGCTGAACAATCGc                                                                        >  1:1180354/1‑71 (MQ=255)
gATTCCTGCGGAATACGCGCCAGGGTGTCACCAGAGCTGATCTGTACGCCATCTTCCAGCTGAACAATCGc                                                                        >  1:2178440/1‑71 (MQ=255)
  ttCCTGCGGAATACGCGCCAGGGTGTCACCAGAGCTGATCTGTACGCCATCTTTCAGCTGAACAATCGCtt                                                                      >  1:195302/1‑71 (MQ=255)
  ttCCTGCGGAATACGCGCCAGGGTGTCACCAGAGCTGATCTGTACGCCATCTTCCAGCTGAACAATCGttt                                                                      >  1:3509953/1‑71 (MQ=255)
     cTGCGGAATACGCGCCAGGGTGTCACCAGATCTGATCTGTACGCCATCTTCCAGCTGAACAATCGCTTTAc                                                                   >  1:2299246/1‑71 (MQ=255)
       gCGGAATACGCGCCAGGGTGTCACCAGAGCTGATCTGTACGCCATCTTCCAGCTGAACAATCGCTTTa                                                                    >  1:1268353/1‑68 (MQ=255)
        cGGAATACGCGCCAGGGTGTCACCAGAGCTGATCTGTACGCCATCTTCCAGCTGAACAATCGCTTTACCCg                                                                >  1:706410/1‑71 (MQ=255)
         ggAATACGCGCCAGGGTGTCACCAGAGCTGATCTGTACGCCATCTTCCAGCTGAACAATCGCTTTACCCgg                                                               <  1:2111944/71‑1 (MQ=255)
         ggAATACGCGCCAGGGTGTCACCAGAGCTGATCTGTACGCCATCTTCCAGCTGAACAATCGCTTTACCCgg                                                               >  1:2078582/1‑71 (MQ=255)
            aTACGCGCCAGGGTGTCACCAGAGCTGATCTGTACGCCATCTTCCAGCTGAACAATCGCTTTACCCGGCa                                                             <  1:480087/70‑1 (MQ=255)
             tACGCGCCAGGGTGTCACCAGAGCTGATCTGTACGCCATCTTCCAGCTGAACAATCGCTTTACCCGGCAgg                                                           <  1:3830278/71‑1 (MQ=255)
             tACGCGCCAGGGTGTCACCAGAGCTGATCTGTACGCCATCTTCCAGCTGAACAATCGCTTTACCCGGCAgg                                                           <  1:3830283/71‑1 (MQ=255)
               cgcgCCAGGGTGTCACCAGAGCTGATCTGTACGCCATCTTCCAGCTGAACAATCGCTTTACCCGGCAgg                                                           >  1:1663411/1‑69 (MQ=255)
               cgcgCCAGGGTGTCACCAGAGCTGATCTGTACGCCATCTTCCAGCTGAACAATCGCTTTACCCGGCAgg                                                           >  1:436468/1‑69 (MQ=255)
               cgcgCCAGGGTGTCACCAGAGCTGATCTGTACGCCATCTTCCAGCTGAACAATCGCTTTACCCGGCAgg                                                           >  1:744639/1‑69 (MQ=255)
               cgcgCCAGGGTGTCACCAGAGCTGATCTGTACGCCATCTTCCAGCTGAACAATCGCTTTACCCGGCAgg                                                           >  1:3889197/1‑69 (MQ=255)
               cgcgCCAGGGTGTCACCAGAGCTGATCTGTACGCCATCTTCCAGCTGAACAATCGCTTTACCCGGCAGGa                                                          <  1:3471915/70‑1 (MQ=255)
               cgcgCCAGGGTGTCACCAGAGCTGATCTGTACGCCATCTTCCAGCTGAACAATCGCTTTACCCGGCAGGa                                                          <  1:1165963/70‑1 (MQ=255)
               agcgcCAGGGTGTCACCAGAGCTGATCTGTACGCCATCTTCCAGCTGAACAATCGCTTTACCCGGCAgg                                                           >  1:3159559/2‑69 (MQ=255)
                gcgcCAGGGTGTCACCAGAGCTGATCTGTACGCCATCTTCCAGCTGAACAATCGCTTTACCCGGCAGGAAg                                                        <  1:3521941/71‑1 (MQ=255)
                gcgcCAGGGTGTCACCAGAGCTGATCTGTACGCCATCTTCCAGCTGAACAATCGCTTTACCCGGCAGGAAg                                                        >  1:2450595/1‑71 (MQ=255)
                gcgcCAGGGTGTCACCAGAGCTGATCTGTACGCCATCTTCCAGCTGAACAATCGCTTTACCCGGCAGGAAg                                                        >  1:459451/1‑71 (MQ=255)
                gcgcCAGGGTGTCACCAGAGCTGATCTGTACGCCATCTTCCAGCTGAACAATCGCTTTACCCGGCAGGAAg                                                        >  1:464357/1‑71 (MQ=255)
                 cgcCAGGGTGTCACCAGAGCTGATCTGTACGCCATCTTCCAGCTGAACAATCGCTTTACCCGGCAGGAAg                                                        >  1:3469598/1‑70 (MQ=255)
                 cgcCAGGGTGTCACCAGAGCTGATCTGTACGCCATCTTCCAGCTGAACAATCGCTTTACCCGGCAGGAAGt                                                       <  1:2020902/71‑1 (MQ=255)
                  gcCAGGGTGTCACCAGAGCTGATCTGTACGCCATCTTCCAGCTGAACAATCGCTTTACCCGGCAGGaa                                                         >  1:1871369/1‑68 (MQ=255)
                          gtCACCAGAGCTGATCTGTACGCCATCTTCCAGCTGAACAATCGCTTTACCCGGCAGGAAg                                                        <  1:796267/61‑1 (MQ=255)
                          gtCACCAGAGCTGATCTGTACGCCATCTTCCAGCTGAACAATCGCTTTACCCGGCAGGAAg                                                        <  1:118371/61‑1 (MQ=255)
                          gtCACAAGAGCTGATCTGTACGCCATCTTCCAGCTGAACAATCGTTTTACCCGGCAGGAAGTACTGCGCTg                                              >  1:4083460/1‑71 (MQ=255)
                           tCACCAGAGCTGATCTGTACGCCATCTTCCAGCTGAACAATCGTTTTa                                                                    <  1:3613206/48‑1 (MQ=255)
                           tCACCAGAGCTGATCTGTACGCCATCTTCCAGCTGAACAATCGTTTTa                                                                    <  1:4056580/48‑1 (MQ=255)
                           tCACCAGAGCTGATCTGTACGCCATCTTCCAGCTGAACAATCGTTTTa                                                                    <  1:1544809/48‑1 (MQ=255)
                           tCACCAGAGCTGATCTGTACGCCATCTTCCAGCTGAACAATCGTTTTa                                                                    <  1:263315/48‑1 (MQ=255)
                           tCACCAGAGCTGATCTGTACGCCATCTTCCAGCTGAACAATCGCTTTACCCGGCAgg                                                           <  1:1967683/57‑1 (MQ=255)
                           tCACCAGAGCTGATCTGTACGCCATCTTCCAGCTGAACAATCGCTTTACCCGGCAGGAAGt                                                       <  1:1412254/61‑1 (MQ=255)
                             aCCAGAGCTGATCTGTACGCCATCTTCCAGCTGAACAATCGTTTTAcc                                                                  >  1:2964414/1‑48 (MQ=255)
                             aCCAGAGCTGATCTGTACGCCATCTTCCAGCTGAACAATCGTTTTAcc                                                                  >  1:1583440/1‑48 (MQ=255)
                             aCCAGAGCTGATCTGTACGCCATCTTCCAGCTGAACAATCGCTTTACCCCGCAGGAAGt                                                       <  1:1506028/59‑1 (MQ=255)
                              ccAGAGCTGATCTGTACGCCATCTTCCAGCTGAACAATCGCTTTACCCGGCAGGAAGTACTGCGCTGGCa                                           >  1:202006/1‑70 (MQ=255)
                               cAGAGCTGATCTGTACGCCATCTTCCAGCTGAACAATCGCTTTACCCGGCAGGAAGTACTGCGCTGGCata                                         >  1:3714008/1‑71 (MQ=255)
                                agagCTGATCTGTACGCCATCTTCCAGCTGAACAATCGCTTTACCCGGCAGGAAGTACTGCGCTGGCatat                                        <  1:2597266/71‑1 (MQ=255)
                                 gagCTGATCTGTACGCCATCTTCCAGCTGAACAATCGCTTTACCCGGCAGGaa                                                         >  1:204650/1‑53 (MQ=255)
                                 gagCTGATCTGTACGCCATCTTCCAGCTGAACAATCGCTTTACCCGGCAGGAAGTACTGCGCTgg                                             <  1:395736/65‑1 (MQ=255)
                                 gagCTGATCTGTACGCCATCTTCCAGCTGAACAATCGCTTTACCCGGCAGGAAGTACTGCGCTgg                                             <  1:2688190/65‑1 (MQ=255)
                                 gagCTGATCTGTACGCCATCTTCCAGCTGAACAATCGCTTTACCCGGCAGGAAGTACTGCGATGGCATATc                                       >  1:3373879/1‑71 (MQ=255)
                                  agCTGATCTGTACGCCATCTTCCAGCTGAACAATCGCTTTACCCGGCAGGAAGTACTGCGCTgg                                             <  1:3397865/64‑1 (MQ=255)
                                  agCTGATCTGTACGCCATCTTCCAGCTGAACAATCGCTTTACCCGGCAGGAAGTACTGCGCTGGCATATCg                                      <  1:1974857/71‑1 (MQ=255)
                                    cTGATCTGTACGCCATCTTCCAGCTGAACAATCGCTTTACCCGGCAg                                                            <  1:4113501/47‑1 (MQ=255)
                                    cTGATCTGTACGCCATCTTCCAGCTGAACAATCGCTTTACCCGGCAGGAAGTACTGCGCTGGCATATCgg                                     >  1:3363968/1‑70 (MQ=255)
                                    cTGATCTGTACGCCATCTTCCAGCTGAACAATCGCTTTACCCGGCAGGAAGTACTGCGCTGGCATATCGGt                                    <  1:1609867/71‑1 (MQ=255)
                                     tGATCTGTACGCCATCTTCCAGCTGAACAATCGCTTTACCCGGCAGGAAGTACTGCGCTGGCATATCGGTa                                   >  1:2692015/1‑71 (MQ=255)
                                      gATCTGTACGCCATCTTCCAGCTGAACAATCGCTTTACCCGGCAGGAAGTACTGCGCTGGCATATCGGTAc                                  >  1:3474260/1‑71 (MQ=255)
                                      gATCTGTACGCCATCTTCCAGCTGAACAATAGCTTTACCCGGCAGGAAGTACTGCGCTgg                                             <  1:2683027/60‑1 (MQ=255)
                                       aTCTGTACGCCATCTTCCAGCTGAACAATCGCTTTACCCGGCAGGAAGTACTGCGCTGGCATATCGGTa                                   >  1:276730/1‑69 (MQ=255)
                                          tGTACGCCATCTTCCAGCTGAACAATCGCTTTACCCGGCAGGAAGTACTGCGCTGGCATATCGGTa                                   <  1:3194099/66‑1 (MQ=255)
                                             aCGCCATCTTCCAGCTGAACAATCGCTTTACCCGGCa                                                             >  1:1743743/1‑37 (MQ=255)
                                              cGCCATCTTCCAGCTGAACAATCGCTTTACCCGGCAGGAAGTAc                                                     >  1:3651838/1‑44 (MQ=255)
                                               gCCATCTTCCAGCTGAACAATCGCTTTACCCGGCAGGAAGTACTGCGCTGGCATATCGGTACCTGGGATCa                         <  1:3400295/71‑1 (MQ=255)
                                                ccATCTTCCAGCTGAACAATCGTTTTACCCGGCAGGAAGTACTGCGCTGGCATATCGGTACCTGGGATCAg                        <  1:2557460/71‑1 (MQ=255)
                                                ccATCTTCCAGCTGAACAATCGCTTTACCCGGCAgg                                                           >  1:3514044/1‑36 (MQ=255)
                                                 cATCTTCCAGCTGAACAATCGCTTTGCCCGGCAGGAAGTACTGCGCTGGCATATCGGTACCTGGGATCAGa                       <  1:2012734/71‑1 (MQ=255)
                                                 cATCTTCCAGCTGAACAATCGCTTTACCCGGCAGGAAGTACTGCGCTGGCATATCGGTACCTGGGATCAGa                       <  1:3058345/71‑1 (MQ=255)
                                                  aTCTTCCAGCTGAACAATCGCTTTACCCGGCAGGAAGTACTGCGCTgg                                             <  1:1930532/48‑1 (MQ=255)
                                                  aTCTTCCAGCTGAACAATCGCTTTACCCGGCAGGAAGTACTGCGCTGGCATATCGGTACCTGGGATCAg                        <  1:1389032/69‑1 (MQ=255)
                                                   tCTTCCAGCTGAACAATCGCTTTACCCGGCAGGAAGTACTGCGCTGGCATATCGGTACCTGGGATCAGAAc                     >  1:1918618/1‑71 (MQ=255)
                                                    cTTCCAGCTGAACAATCGCTTTACCCGGCAGGAAGTACTGCGCTGGCATATCGGTACCTGGGATCAGAAcg                    >  1:1682085/1‑71 (MQ=255)
                                                      tCCAGCTGAACAATCGTTTTACCCGGCAGGAAGTACTGCGCTgg                                             >  1:534677/1‑44 (MQ=255)
                                                      tCCAGCTGAACAATCGCTTTACCCGGCAGGAAGTACTGCGCt                                               >  1:2165551/1‑42 (MQ=255)
                                                       ccAGCTGAAGAATCGCTTTACCCGGCAGGAAGTACTGCGCTGGCATATCGGTACCTGGGATCAGAAcgtcg                 >  1:3585223/1‑71 (MQ=255)
                                                       ccAGCTGAACAATCGCTTTACCCGGCAGGAAGTACTGCGCTGGCATATCGGTACCTGGGATCAGAAcgtcg                 >  1:4063877/1‑71 (MQ=255)
                                                       ccAGCTGAACAATCGCTTTACCCGGCAGGAAGTACTGCGCTGGCATATCGGTACCTGGGATCAGAAcgtc                  <  1:697059/70‑1 (MQ=255)
                                                       ccAGCTGAACAATCGCTTTACCCGGCAGGAAGTACTGCGCTGGCATATCGGTACCTGGGATCAGAAcgtc                  <  1:1868320/70‑1 (MQ=255)
                                                       ccAGCTGAACAATCGCTTTACCCGGCAGGAAGTACTGCGCTGGCATATCGGTACCTGGGATCAGAAcg                    <  1:244123/68‑1 (MQ=255)
                                                       ccAGCTGAACAATCGCTTTACCCGGCAGGAAGTACTGCGCTGGCATATCGGTACCTGGGATCAGAAc                     <  1:1817067/67‑1 (MQ=255)
                                                        cAGCTGAACAATCGCTTTACCCGGCAGGAAGTACTGCGCTGGCATATCg                                      >  1:3594654/1‑49 (MQ=255)
                                                        cAGCTGAACAATCGCTTTACCCGGCAGGAAGTACTGCGCTGGCATATCGGTAc                                  >  1:3522503/1‑53 (MQ=255)
                                                        cAGCTGAACAATCGCTTTACCCGGCAGGAAGTACTGCGCTGGCATATCGGTAc                                  >  1:2766134/1‑53 (MQ=255)
                                                        cAGCTGAACAATCGCTTTACCCGGCAGGAAGTACTGCGCTGGCATATCGGTAc                                  >  1:900495/1‑53 (MQ=255)
                                                         aGCTGAACAATCGTTTTACCCGGCAGGAAGTACTGCGCTGGCATATCGGTACCTGGGATCAGAAcgtcgt                >  1:620272/1‑70 (MQ=255)
                                                          gCTGAACAATCGCTTTACCCGGCAGGAAGTACTGCGCTGGCATATCGGTACCTGGGATCAGAAAGTCGTTa              <  1:1355955/71‑1 (MQ=255)
                                                           cTGAACAATCGCTTTACCCGGCAGGAAGTACTGCGCTGGCATATCGGTACCTGGGATCAGAACGTCGTTa              <  1:3761809/70‑1 (MQ=255)
                                                              aaCAATCGTTTTACCCGGCAGGAAGTACTGCGCTGGCATATCGGTACCTGGGATCAGAACGTCGTTAcc            <  1:3607064/69‑1 (MQ=255)
                                                              aaCAATCGCTTTACCCGGCAGGAAGTACTGCGCTGGCATATCGGTACCTGGGATCAGAACGTCGTTACCCt          >  1:623088/1‑71 (MQ=255)
                                                                cAATCGTTTTACCCGGCAGGAAGTACTGCGCTGGCATATCGGTACCTGGGATCAGAACGTCGTTACCCTGa        <  1:474892/71‑1 (MQ=255)
                                                                  aTCGTTTTACCCGGCAGGAAGTACTGCGCTGGCATATCGGTACCTGGGATCAg                        <  1:1058404/53‑1 (MQ=255)
                                                                  aTCGCTTTACCCGGCAGGAAGTACTGCGCTGGCATATCGGTACCTGGGATCAGAACGTCGTTACCCTGAg       >  1:2952229/1‑70 (MQ=255)
                                                                   tCGCTTTACCCGGCAGGAAGTACTGCGCTGGCATATCGGTACCTGGGATCAGAACGTCGTTACCCTGAGCa     <  1:1532491/71‑1 (MQ=255)
                                                                     gCTTTACCCGGCAGGAAGTACTGCGCTGGCATATCGGTACCTGGGATCAGAACGTCGTTACCCTGAGCAt    >  1:3350966/1‑70 (MQ=255)
                                                                     gCTTTACCCGGCAGGAAGTACTGCGCTGGCATATCGGTACCTGGGATCAGAACGTCGTTACCCTGAGCATc   <  1:312922/71‑1 (MQ=255)
                                                                      cTTTACCCGGCAGGAAGTACTGCGCTGGCATATCGGTACCTGGGATCAGAACGTCGTTACCCTGAGCATCa  <  1:2205306/71‑1 (MQ=255)
                                                                      |                                                                      
GATTCCTGCGGAATACGCGCCAGGGTGTCACCAGAGCTGATCTGTACGCCATCTTCCAGCTGAACAATCGCTTTACCCGGCAGGAAGTACTGCGCTGGCATATCGGTACCTGGGATCAGAACGTCGTTACCCTGAGCATCA  >  minE/2253655‑2253795

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: