Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,254,647 A→G 13.5% T797T (ACT→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,6470AG13.5% 256.6 / 25.2 89T797T (ACT→ACCrpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base A (35/42);  new base G (3/9);  total (38/51)
Fisher's exact test for biased strand distribution p-value = 2.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATACCTTCATGGGTACCACAATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGTTTTCAGTGCGGTATCCGCCAGACCTTTACGAGCACCGTGGGTG  >  minE/2254579‑2254713
                                                                    |                                                                  
aTACCTTCATGGGTACCACAATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAg                                                                   >  1:1895339/1‑70 (MQ=255)
aTACCTTCATGGGTACCACAATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGt                                                                  >  1:1057772/1‑71 (MQ=255)
aTACCTTCATGGGTACCACAATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGt                                                                  >  1:368078/1‑71 (MQ=255)
 tACCTACATGGGTACCACAATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGt                                                                  >  1:1361052/1‑70 (MQ=255)
      tCATGGGTACCACAATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGAcgacga                                                                    >  1:1817848/1‑63 (MQ=255)
      tCATGGGTACCACAATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGAcgacga                                                                    >  1:3122633/1‑63 (MQ=255)
      tCATGGGTACCACAATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGAcgacga                                                                    >  1:3084721/1‑63 (MQ=255)
      tCATGGGTACCACAATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGt                                                             >  1:2967108/1‑70 (MQ=255)
      tCATGGGTACCACAATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGt                                                             >  1:1600744/1‑70 (MQ=255)
       cATGGGTACCACAATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGt                                                             <  1:3270563/69‑1 (MQ=255)
       cATGGGTACCACAATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGAAGAGTCAGGTa                                                            <  1:714221/70‑1 (MQ=255)
        aTGGGTACCACAATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAAc                                                          <  1:1723953/71‑1 (MQ=255)
              aCCACAATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCGGAGt                                                    >  1:1814300/1‑71 (MQ=255)
               ccACAATCGTCTTCGGTAACCACCAGGTTCTGCGCCACGGCAACCAGACGACGAGTCAGGTAACCGGAGt                                                    <  1:478107/70‑1 (MQ=255)
                cacaATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCGGAGTTc                                                  >  1:2811659/1‑71 (MQ=255)
                  caATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCg                                                 >  1:2517912/1‑70 (MQ=255)
                   aaTCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTc                                                                 <  1:2761326/53‑1 (MQ=255)
                    aTCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTa                                                            <  1:411812/57‑1 (MQ=255)
                    aTCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTa                                                            <  1:1391619/57‑1 (MQ=255)
                    aTCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAg                                              <  1:3231330/71‑1 (MQ=255)
                    aTCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAg                                              <  1:746673/71‑1 (MQ=255)
                     tCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGt                                             >  1:2539298/1‑71 (MQ=255)
                     tCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGt                                             >  1:1741901/1‑71 (MQ=255)
                        tCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTa                                                            <  1:145764/53‑1 (MQ=255)
                        tCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGtttt                                          <  1:3479414/71‑1 (MQ=255)
                        tCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGtttt                                          <  1:1208562/71‑1 (MQ=255)
                          ttCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAAc                                                          >  1:281974/1‑53 (MQ=255)
                            cGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAg                                                               >  1:1319297/1‑46 (MQ=255)
                              gTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGTTTTCAGTg                                     <  1:2128214/70‑1 (MQ=255)
                              gTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGTTTTCAGTg                                     <  1:846624/70‑1 (MQ=255)
                               tcaccaccAGATCTTGCGTCACGTCCACCAGACGACGGGTc                                                                 <  1:606162/39‑1 (MQ=16)
                               tcaccaccAGATCTTGCGTCACGTCCACCAGACGACGGGTc                                                                 <  1:196598/39‑1 (MQ=16)
                               tcaccaccAGATCTTGCGTCACGTCCACCAGACGACGGGTc                                                                 <  1:691904/39‑1 (MQ=16)
                               tcaccaccAGATCTTGCGTCACGTCCACCAGACGACGGGTc                                                                 <  1:3434053/39‑1 (MQ=16)
                               tcaccaccAGATCTTGCGTCACGTCCACCAGACGACGGGTc                                                                 <  1:801822/39‑1 (MQ=16)
                               tcaccaccAGATCTTGCGTCACGTCCACCAGACGACGGGTc                                                                 <  1:3329271/39‑1 (MQ=16)
                               tAACCACCAGGTCCTGCGCCACGTCAACCAGAcgacga                                                                    <  1:428827/38‑1 (MQ=255)
                               tAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCgg                                                       <  1:2749328/51‑1 (MQ=255)
                               tAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCgg                                                       <  1:2366973/51‑1 (MQ=255)
                               tAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCgg                                                       <  1:2330950/51‑1 (MQ=255)
                               tAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCgg                                                       <  1:1947422/51‑1 (MQ=255)
                               tAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCgg                                                       <  1:2208723/51‑1 (MQ=255)
                               tAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCgg                                                       <  1:1200365/51‑1 (MQ=255)
                               tAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGTTTTCAGTGCg                                   >  1:762294/1‑71 (MQ=255)
                                 accaccAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCgg                                                       <  1:660541/49‑1 (MQ=255)
                                 accaccAGATCTTGCGTCACGTCCACCAGACGACGGGTCAg                                                               >  1:1191944/1‑41 (MQ=19)
                                 accaccAGATCTTGCGTCACGTCCACCAGACGACGGGTCAg                                                               >  1:2836454/1‑41 (MQ=19)
                                 accaccAGATCTTGCGTCACGTCCACCAGACGACGGGTCAg                                                               >  1:2577199/1‑41 (MQ=19)
                                         gTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGtttt                                          <  1:3193477/54‑1 (MQ=255)
                                         gTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGTTTTCAGTGCGGTATCCGCCa                         <  1:1632490/71‑1 (MQ=255)
                                         gTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGTTTTCAGTGCGGTATCCGCCa                         <  1:1529348/71‑1 (MQ=255)
                                         gTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGTTTTCAGTGCGGTATCCGCCa                         <  1:2112265/71‑1 (MQ=255)
                                               cgcCACGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGTTTTCAGTGCGGTATCCGCCAGa                       >  1:1208169/1‑67 (MQ=255)
                                               cgcCACGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGTTTTCAGTGCGGTATCCGCCAGAc                      >  1:1407179/1‑68 (MQ=255)
                                               cgcCACGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGTTTTCAGTGCGGTATCCGCCAGAc                      >  1:2030528/1‑68 (MQ=255)
                                               cgcCACGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGTTTTCAGTGCGGTATCCGCCAGAc                      >  1:2045586/1‑68 (MQ=255)
                                               cgcCACGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGTTTTCAGTGCGGTATCCGCCAGACCtt                   <  1:2442844/71‑1 (MQ=255)
                                                 ccACGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGtttt                                          <  1:2868052/46‑1 (MQ=255)
                                                 ccACGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGtttt                                          <  1:2858556/46‑1 (MQ=255)
                                                 ccACGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGttt                                           <  1:3006562/45‑1 (MQ=255)
                                                 ccACGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGttt                                           <  1:3264060/45‑1 (MQ=255)
                                                  cACGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGTTTTCAGTGCGGTATCCGCCAGACCTTTAc                <  1:3055801/71‑1 (MQ=255)
                                                   aCGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGTTTTc                                         >  1:943406/1‑45 (MQ=255)
                                                   aCGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGTTTTc                                         >  1:1354788/1‑45 (MQ=255)
                                                   aCGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGTTTTc                                         >  1:1068827/1‑45 (MQ=255)
                                                   aCGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGTTTTCa                                        >  1:279901/1‑46 (MQ=255)
                                                   aCGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGTTTTCAGTGCGGTATCCGCCAGACCTTTa                 >  1:218428/1‑69 (MQ=255)
                                                    cGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGTTTTCAGTGCGGTATCCGCCAGACCttt                  <  1:267509/67‑1 (MQ=255)
                                                    cGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGTTTTCAGTGCGGTATCCGCCAGACCttt                  <  1:2648786/67‑1 (MQ=255)
                                                    cGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGTTTTCAGTGCGGTATCCGCCAGACCttt                  <  1:2718938/67‑1 (MQ=255)
                                                    cGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGTTTTCAGTGCGGTATCCGCCAGACCttt                  <  1:1171308/67‑1 (MQ=255)
                                                    cGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGTTTTCAGTGCGGTATCCGCCAGACCTTTACg               >  1:1078566/1‑70 (MQ=255)
                                                     gTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGTTTTCAGTGCGGTATCCGCCAGACCttt                  <  1:1302728/66‑1 (MQ=255)
                                                      tCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGttt                                           <  1:1169482/40‑1 (MQ=255)
                                                      tCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGTTTTc                                         <  1:163596/42‑1 (MQ=255)
                                                        aaCCAGACGACGAGTCAGGTAACCGGAGTTCGCAGTTTTCAGTGCGGTATCCGCCAGACCTTTACGAGCAc          >  1:1027300/1‑71 (MQ=255)
                                                         aCCAGACGACGAGTCAGGTAACCGGAGTTCGCAGTTTTCAGTGCGGTATCCGCCAGACCTTTACGAGCAcc         >  1:1679151/1‑71 (MQ=255)
                                                         aCCAGACGACGAGTCAGGTAACCGGAGTTCGCAGTTTTCAGTGCGGTATCCGCCAGACCTTTACGAGCAcc         >  1:555729/1‑71 (MQ=255)
                                                         aCCAGACGACGAGTCAGGTAACCGGAGTTCGCAGTTTTCAGTGCGGTATCCGCCAGACCTTTACGAGCAcc         >  1:1608306/1‑71 (MQ=255)
                                                         aCCAGACGACGAGTCAGGTAACCGGAGTTCGCAGTTTTCAGTGCGGTATCCGCCAGACCTTTACAAGCAcc         >  1:2003846/1‑71 (MQ=255)
                                                          ccAGACGACGGGTCAGGTAACCCGAGTTCGCGGTCTTCAGCGCCGTATCCGCCAGACCTTTACGAGCGCCg        <  1:2221972/71‑1 (MQ=255)
                                                          ccAGACGACGGGTCAGGTAACCCGAGTTCGCGGTCTTCAGCGCCGTATCCGCCAGACCTTTACGAGCGCCg        <  1:2281818/71‑1 (MQ=255)
                                                          ccAGACGACGAGTCAGGTAACCGGAGTTCGCAGTTTTc                                         <  1:3292688/38‑1 (MQ=255)
                                                          ccAGACGACGAGTCAGGTAACCGGAGTTCGCAGTTTTc                                         <  1:634127/38‑1 (MQ=255)
                                                          ccAGACGACGAGTCAGGTAACCGGAGTTCGCAGTTTTCAGTGCGGTATCCGCCAGACCTTTACGAGCAcc         >  1:2607307/1‑70 (MQ=255)
                                                          ccAGACGACGAGTCAGGTAACCGGAGTTCGCAGTTTTCAGTGCGGTATCCGCCAGACCTTTACGAGCACCg        >  1:1747440/1‑71 (MQ=255)
                                                           cAGACGACGAGTCAGGTAACCGGAGTTCGCAGTTTTCAGTGCGGTATCCGCCAGACCtt                   <  1:2052106/59‑1 (MQ=255)
                                                           cAGACGACGAGTCAGGTAACCGGAGTTCGCAGTTTTCAGTGCGGTATCCGCCAGACCTTTACGAGCACCgt       >  1:1785586/1‑71 (MQ=255)
                                                            aGACGACGGGTCAGGTAACCCGAGTTCGCGGTCTTCAGCGCCGTATCCGCCAGACCTTTACGAGCGCCg        <  1:2369359/69‑1 (MQ=255)
                                                            aGACGACGAGTCAGGTAACCGGAGTTCGCAGTTTTCAGTGCGGTATCCGCCAGACCTTTACGAGCAcc         >  1:1043161/1‑68 (MQ=255)
                                                            aGACGACGAGTCAGGTAACCGGAGTTCGCAGTTTTCAGTGCGGTATCCGCCAGACCTTTACGAGCACCgtg      >  1:3337733/1‑71 (MQ=255)
                                                             gacgacGAGTCAGGTAACCGGAGTTCGCAGTTTTCAGTg                                     <  1:3069493/39‑1 (MQ=255)
                                                             gacgacGAGTCAGGTAACCGGAGTTCGCAGTTTTCAGTGCGGTATCCGCCAGACCTTTACGAGCACCgtgg     >  1:283928/1‑71 (MQ=255)
                                                                   gaGTCAGGTAACCGGAGTTCGCAGTTTTCAGTGCGGTATCCGCCAGACCTTTACGAGCACCgtgggtg  >  1:2829276/1‑68 (MQ=255)
                                                                   gaGTCAGGTAACCGGAGTTCGCAGTTTTCAGTGCGGTATCCGCCAGACCTTTACGAGCACCgtggg    >  1:390718/1‑66 (MQ=255)
                                                                   gaGTCAGGTAACCGGAGTTCGCAGTTTTCAGTGCGGTATCCGCCAGACCTTTACGAGCACCgtggg    >  1:921506/1‑66 (MQ=255)
                                                                    |                                                                  
ATACCTTCATGGGTACCACAATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGTTTTCAGTGCGGTATCCGCCAGACCTTTACGAGCACCGTGGGTG  >  minE/2254579‑2254713

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: