Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,254,767 C→G 13.5% T757T (ACG→ACC rpoC ← RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,254,7670CG13.5% 194.0 / 19.8 74T757T (ACG→ACCrpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base C (23/41);  new base G (6/4);  total (29/45)
Fisher's exact test for biased strand distribution p-value = 1.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

CGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAATCTGTGC  >  minE/2254699‑2254835
                                                                    |                                                                    
cGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGt                                                                    <  1:3602964/71‑1 (MQ=255)
    cACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCg                                                                <  1:3674579/71‑1 (MQ=255)
    cACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCg                                                                <  1:3680419/71‑1 (MQ=255)
    cACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCg                                                                <  1:2490041/71‑1 (MQ=255)
    cACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCg                                                                <  1:1036823/71‑1 (MQ=255)
    cACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCg                                                                <  1:866263/71‑1 (MQ=255)
      ccGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCgat                                                              <  1:3583096/71‑1 (MQ=255)
        gtgggtggAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCgatga                                                            >  1:1140902/1‑71 (MQ=255)
        gtgggtggAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCgatga                                                            >  1:2141417/1‑71 (MQ=255)
        gtgggtggAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCgatga                                                            <  1:2328737/71‑1 (MQ=255)
        gtgggtggAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCgatga                                                            >  1:4459948/1‑71 (MQ=255)
        gtgggtggAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCgatga                                                            >  1:2753853/1‑71 (MQ=255)
         tgggtggAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCgatgat                                                           >  1:4466266/1‑71 (MQ=255)
         tgggtggAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCgatgat                                                           >  1:3934677/1‑71 (MQ=255)
           ggtggAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGatgatg                                                          <  1:250970/70‑1 (MQ=255)
               gagaTGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGc                                                      >  1:1973328/1‑70 (MQ=255)
                agaTGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCa                                                    <  1:3638102/71‑1 (MQ=255)
                 gaTGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCAt                                                   <  1:3370327/71‑1 (MQ=255)
                 gaTGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCAt                                                   <  1:1619030/71‑1 (MQ=255)
                     aaGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGttt                                                                  <  1:740292/52‑1 (MQ=255)
                     aaGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGttt                                                                  <  1:4158060/52‑1 (MQ=255)
                     aaGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGttt                                                                  <  1:2658214/52‑1 (MQ=255)
                     aaGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGttt                                                                  <  1:4338288/52‑1 (MQ=255)
                     aaGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGttt                                                                  <  1:2083262/52‑1 (MQ=255)
                     aaGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGttt                                                                  <  1:5078207/52‑1 (MQ=255)
                     aaGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGttt                                                                  <  1:593299/52‑1 (MQ=255)
                     aaGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCgg                                               >  1:3747987/1‑71 (MQ=255)
                      aGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGc                                              <  1:2367583/71‑1 (MQ=255)
                      aGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGc                                              <  1:3297121/71‑1 (MQ=255)
                            ggAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCg                                          >  1:3662170/1‑69 (MQ=255)
                            ggAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGc                                         >  1:2802122/1‑70 (MQ=255)
                            ggAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGc                                         >  1:4010258/1‑70 (MQ=255)
                            ggAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGc                                         >  1:3499390/1‑70 (MQ=255)
                                tACGTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCgg                                               <  1:4690708/60‑1 (MQ=255)
                                tACGTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCgg                                               <  1:3641668/60‑1 (MQ=255)
                                tACGTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCgg                                               <  1:3306006/60‑1 (MQ=255)
                                 aCGTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCgg                                               <  1:1160143/59‑1 (MQ=255)
                                  cGTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCGGCt                                             >  1:4047195/1‑60 (MQ=255)
                                  cGTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCGGCt                                             >  1:294702/1‑60 (MQ=255)
                                  cGTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCGGCt                                             >  1:2760828/1‑60 (MQ=255)
                                  cGTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCGGCt                                             >  1:1737997/1‑60 (MQ=255)
                                  cGTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCGGCt                                             >  1:1581500/1‑60 (MQ=255)
                                  agTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCGGCt                                             >  1:461485/2‑60 (MQ=255)
                                    ttCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGAcc                                <  1:1159017/71‑1 (MQ=255)
                                    ttCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGAc                                 >  1:1283729/1‑70 (MQ=255)
                                    ttCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGAc                                 >  1:4213329/1‑70 (MQ=255)
                                      cAGAGCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCAc                              >  1:1040516/1‑71 (MQ=255)
                                      cAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCAc                              >  1:3658936/1‑71 (MQ=255)
                                       aGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCAc                              <  1:3303827/70‑1 (MQ=255)
                                       aGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCAc                              <  1:4016088/70‑1 (MQ=255)
                                       aGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACg                             <  1:3205649/71‑1 (MQ=255)
                                       aGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATACGGCTTCGCCATCAGACCAc                              <  1:3090693/70‑1 (MQ=255)
                                           cTTCACGGAAGTTCGCGGTGATTGGCGTTTCGAAGATGGAGCCATCCgg                                               <  1:2967508/49‑1 (MQ=255)
                                            ttCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCgg                                               <  1:4973791/48‑1 (MQ=255)
                                             tCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCa                           <  1:886499/67‑1 (MQ=255)
                                              cACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCa                           <  1:4818177/66‑1 (MQ=255)
                                              cACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATAc                        >  1:3823471/1‑69 (MQ=255)
                                              cACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATAc                        >  1:1726566/1‑69 (MQ=255)
                                                cGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGc                    <  1:2769271/71‑1 (MQ=255)
                                                cGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGc                    <  1:1202514/71‑1 (MQ=255)
                                                cGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGc                    <  1:7006/71‑1 (MQ=255)
                                                cGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGc                    <  1:4048458/71‑1 (MQ=255)
                                                cGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGc                    <  1:4027435/71‑1 (MQ=255)
                                                 ggAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGc                    >  1:1557053/1‑70 (MQ=255)
                                                    aGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCgg                                               <  1:2906056/40‑1 (MQ=255)
                                                    aGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGc                >  1:3509438/1‑71 (MQ=255)
                                                          cgGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGa          <  1:1402842/71‑1 (MQ=255)
                                                            gTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAAt        <  1:4762503/71‑1 (MQ=255)
                                                            gTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAAt        >  1:3248259/1‑71 (MQ=255)
                                                            gTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAAt        >  1:4995036/1‑71 (MQ=255)
                                                            gTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAAt        <  1:5072192/71‑1 (MQ=255)
                                                            gTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAAt        <  1:4459727/71‑1 (MQ=255)
                                                             tGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAATc       >  1:1870126/1‑71 (MQ=255)
                                                                   gCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAATCTGTGc  <  1:5027828/70‑1 (MQ=255)
                                                                    |                                                                    
CGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAATCTGTGC  >  minE/2254699‑2254835

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: