Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,254,740 A→G 10.8% G766G (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,254,7400AG10.8% 256.4 / 18.0 95G766G (GGT→GGCrpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base A (43/41);  new base G (6/4);  total (50/45)
Fisher's exact test for biased strand distribution p-value = 7.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.81e-01

TTCAGTGCGGTATCCGCCAGACCTTTACGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACG  >  minE/2254672‑2254808
                                                                    |                                                                    
ttCAGTGCGGTATCCGCCAGACCTTTACGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGAcc                                                                    >  1:115373/1‑71 (MQ=255)
  cAGCGCCGTATCCGCCAGACCTTTACGAGCGCCGTGGGTCGAGATGAAGTACTGCAGAACGTTCAGTCCtt                                                                  >  1:750925/1‑71 (MQ=255)
       cGGTATCCGCCAGACCTTTACGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGa                                                                      >  1:5107191/1‑62 (MQ=255)
       cGGTATCCGCCAGACCTTTACGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGAc                                                                     >  1:3612339/1‑63 (MQ=255)
       cGGTATCCGCCAGACCTTTACGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGAc                                                                     >  1:4386182/1‑63 (MQ=255)
       cGGTATCCGCCAGACCTTTACGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGAc                                                                     >  1:4585495/1‑63 (MQ=255)
       cGGTATCCGCCAGACCTTTACGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGAc                                                                     >  1:2289451/1‑63 (MQ=255)
       cGGTATCCGCCAGACCTTTACGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGAc                                                                     >  1:4220552/1‑63 (MQ=255)
       cGGTATCCGCCAGACCTTTACGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGAc                                                                     >  1:1644566/1‑63 (MQ=255)
       cGGTATCCGCCAGACCTTTACGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGAc                                                                     >  1:5072010/1‑63 (MQ=255)
       cGGTATCCGCCAGACCTTTACGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGAc                                                                     >  1:1410621/1‑63 (MQ=255)
       cGGTATCCGCCAGACCTTTACGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGAc                                                                     >  1:1301247/1‑63 (MQ=255)
       cGGTATCCGCCAGACCTTTACGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGAc                                                                     >  1:649130/1‑63 (MQ=255)
       cGGTATCCGCCAGACCTTTACGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGAc                                                                     >  1:753149/1‑63 (MQ=255)
       cGGTATCCGCCAGACCTTTACGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGAc                                                                     >  1:812465/1‑63 (MQ=255)
       cGGTATCCGCCAGACCTTTACGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGAc                                                                     >  1:1070557/1‑63 (MQ=255)
       cGGTATCCGCCAGACCTTTACGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGAc                                                                     >  1:2804398/1‑63 (MQ=255)
          tATCCGCCAGACCTTTACGAGCACCGTGGGTGGAAATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAg                                                          >  1:4768096/1‑71 (MQ=255)
          tATCCGCCAGACCTTTACGAGCACCGTCGGTGGAAATGAAGTACTGGAGTACGTTCAGACCTTCACGGaa                                                           >  1:3912804/1‑70 (MQ=255)
          tATCCGCCAGACCTTTACGAGCACCGAGGGTGGAAATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAg                                                          >  1:4922957/1‑71 (MQ=255)
              cGCCAGACCTTTACGAGCACCGTGGGTGGAGATGAAGTACTGGAGTGCGTTCAGACCTTCACGGAAGtt                                                        >  1:4246128/1‑69 (MQ=255)
              cGCCAGACCTTTACGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGaa                                                           >  1:1829973/1‑66 (MQ=255)
              cGCCAGACCTTTACGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGtt                                                        >  1:3861999/1‑69 (MQ=255)
              cGCCAGACCTTTACGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGtt                                                        >  1:18555/1‑69 (MQ=255)
              cGCCAGACCTTTACGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTc                                                       <  1:712658/70‑1 (MQ=255)
              cGCCAGACCTTTACGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTc                                                       <  1:1048281/70‑1 (MQ=255)
              cGCCAGACCTTTACGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTc                                                       <  1:1835123/70‑1 (MQ=255)
              cGCCAGACCTTTACGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTc                                                       <  1:17428/70‑1 (MQ=255)
              cGCCAGACCTTTACGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTc                                                       <  1:1205087/70‑1 (MQ=255)
               gCCAGACCTTTACGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTc                                                       <  1:3495751/69‑1 (MQ=255)
                   gACCTTTACGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCa                                                                >  1:4848455/1‑56 (MQ=255)
                   gACCTTTACGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGt                                                  >  1:3964406/1‑70 (MQ=255)
                   gACCTTTACGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGt                                                  >  1:1590413/1‑70 (MQ=255)
                      cTTTACGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGt                                                         <  1:482234/60‑1 (MQ=255)
                      cTTTACGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGt                                                         <  1:4914160/60‑1 (MQ=255)
                      cTTTACGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGt                                                         <  1:1133119/60‑1 (MQ=255)
                      cTTTACGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGt                                                         <  1:3699838/60‑1 (MQ=255)
                         tACGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGAca                                                                    >  1:3035947/1‑45 (MQ=255)
                           cGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGt                                         <  1:3602964/71‑1 (MQ=255)
                               cACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCg                                     <  1:3680419/71‑1 (MQ=255)
                               cACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCg                                     <  1:3674579/71‑1 (MQ=255)
                               cACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCg                                     <  1:1036823/71‑1 (MQ=255)
                               cACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCg                                     <  1:866263/71‑1 (MQ=255)
                               cACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCg                                     <  1:2490041/71‑1 (MQ=255)
                                 ccGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCgat                                   <  1:3583096/71‑1 (MQ=255)
                                   gtgggtggAGATGAAGTACTGGAGTACGTTCAGACCtt                                                                  >  1:4664979/1‑38 (MQ=255)
                                   gtgggtggAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCgatga                                 >  1:1140902/1‑71 (MQ=255)
                                   gtgggtggAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCgatga                                 >  1:2753853/1‑71 (MQ=255)
                                   gtgggtggAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCgatga                                 >  1:4459948/1‑71 (MQ=255)
                                   gtgggtggAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCgatga                                 <  1:2328737/71‑1 (MQ=255)
                                   gtgggtggAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCgatga                                 >  1:2141417/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)
                                        tggAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGt                                                         <  1:5082925/42‑1 (MQ=255)
                                          gagaTGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGc                           >  1:1973328/1‑70 (MQ=255)
                                           agaTGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCa                         <  1:3638102/71‑1 (MQ=255)
                                            gaTGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTcg                                                      <  1:504150/41‑1 (MQ=255)
                                            gaTGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTcg                                                      <  1:3058502/41‑1 (MQ=255)
                                            gaTGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTcg                                                      <  1:1291735/41‑1 (MQ=255)
                                            gaTGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTcg                                                      <  1:4113205/41‑1 (MQ=255)
                                            gaTGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCAt                        <  1:1619030/71‑1 (MQ=255)
                                            gaTGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCAt                        <  1:3370327/71‑1 (MQ=255)
                                                aaGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGttt                                       <  1:740292/52‑1 (MQ=255)
                                                aaGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGttt                                       <  1:2083262/52‑1 (MQ=255)
                                                aaGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGttt                                       <  1:4158060/52‑1 (MQ=255)
                                                aaGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGttt                                       <  1:593299/52‑1 (MQ=255)
                                                aaGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGttt                                       <  1:2658214/52‑1 (MQ=255)
                                                aaGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGttt                                       <  1:4338288/52‑1 (MQ=255)
                                                aaGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGttt                                       <  1:5078207/52‑1 (MQ=255)
                                                aaGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCgg                    >  1:3747987/1‑71 (MQ=255)
                                                 aGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGc                   <  1:3297121/71‑1 (MQ=255)
                                                 aGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGc                   <  1:2367583/71‑1 (MQ=255)
                                                       ggAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCg               >  1:3662170/1‑69 (MQ=255)
                                                       ggAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGc              >  1:3499390/1‑70 (MQ=255)
                                                       ggAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGc              >  1:2802122/1‑70 (MQ=255)
                                                       ggAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGc              >  1:4010258/1‑70 (MQ=255)
                                                           tACGTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCgg                    <  1:3306006/60‑1 (MQ=255)
                                                           tACGTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCgg                    <  1:4690708/60‑1 (MQ=255)
                                                           tACGTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCgg                    <  1:3641668/60‑1 (MQ=255)
                                                            aCGTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCgg                    <  1:1160143/59‑1 (MQ=255)
                                                             cGTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCGGCt                  >  1:4047195/1‑60 (MQ=255)
                                                             cGTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCGGCt                  >  1:1581500/1‑60 (MQ=255)
                                                             cGTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCGGCt                  >  1:1737997/1‑60 (MQ=255)
                                                             cGTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCGGCt                  >  1:294702/1‑60 (MQ=255)
                                                             cGTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCGGCt                  >  1:2760828/1‑60 (MQ=255)
                                                             agTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCGGCt                  >  1:461485/2‑60 (MQ=255)
                                                               ttCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGAcc     <  1:1159017/71‑1 (MQ=255)
                                                               ttCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGAc      >  1:4213329/1‑70 (MQ=255)
                                                               ttCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGAc      >  1:1283729/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)
                                                                    |                                                                    
TTCAGTGCGGTATCCGCCAGACCTTTACGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACG  >  minE/2254672‑2254808

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: