Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,254,635 A→C 9.8% V801V (GTT→GTG rpoC ← RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,254,6350AC 9.8% 255.2 / 20.6 92V801V (GTT→GTGrpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base A (33/50);  new base C (3/6);  total (36/56)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGAGTCATCATGATACCTTCATGGGTACCACAATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGTTTTCAGTGCGGTATCCGCCAGACCTTTACGAGCAC  >  minE/2254567‑2254705
                                                                    |                                                                      
ggAGTCATCATGATACCTTCATGGGTACCACAATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAAc                                                                      <  1:95799/71‑1 (MQ=255)
ggAGTCATCATGATACCTTCATGGGTACCACAATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAAc                                                                      <  1:3140067/71‑1 (MQ=255)
ggAGTCATCATGATACCTTCATGGGTACCACAATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAAc                                                                      <  1:3059152/71‑1 (MQ=255)
 gAGTCATCATGATACCTTCATGGGTACCACAATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAAcc                                                                     <  1:2241685/71‑1 (MQ=255)
       tcatGATACCTTCATGGGTACCACAATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAgacga                                                               >  1:933175/1‑71 (MQ=255)
         atGATACCTTCATGGGTACCACAATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAgacgac                                                              >  1:2681904/1‑70 (MQ=255)
         atGATACCTTCATGGGTACCACAATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGacgacg                                                             <  1:2885991/71‑1 (MQ=255)
         atGATACCTTCATGGGTACCACAATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGacgacg                                                             <  1:617799/71‑1 (MQ=255)
         atGATACCTTCATGGGTACCACAATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGacgacg                                                             <  1:29454/71‑1 (MQ=255)
          tGATACCTTCATGGGTACCACAATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAgac                                                                 <  1:2115601/66‑1 (MQ=255)
          tGATACCTTCATGGGTACCACAATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAgac                                                                 <  1:3021307/66‑1 (MQ=255)
          tGATACCTTCATGGGTACCACAATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAgac                                                                 <  1:3220793/66‑1 (MQ=255)
            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:3122633/1‑63 (MQ=255)
                  tCATGGGTACCACAATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGAcgacga                                                            >  1:3084721/1‑63 (MQ=255)
                  tCATGGGTACCACAATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGAcgacga                                                            >  1:1817848/1‑63 (MQ=255)
                  tCATGGGTACCACAATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGt                                                     >  1:1600744/1‑70 (MQ=255)
                  tCATGGGTACCACAATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGt                                                     >  1:2967108/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)
                            cacaATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAAc                                                                      <  1:3060457/43‑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:1391619/57‑1 (MQ=255)
                                aTCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTa                                                    <  1:411812/57‑1 (MQ=255)
                                aTCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAg                                      <  1:3231330/71‑1 (MQ=255)
                                aTCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAg                                      <  1:746673/71‑1 (MQ=255)
                                 tCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGt                                     >  1:1741901/1‑71 (MQ=255)
                                 tCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGt                                     >  1:2539298/1‑71 (MQ=255)
                                    tCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTa                                                    <  1:145764/53‑1 (MQ=255)
                                    tCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGtttt                                  <  1:1208562/71‑1 (MQ=255)
                                    tCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGtttt                                  <  1:3479414/71‑1 (MQ=255)
                                      ttCGGTAACCACCAGGTCCTGCGCCACGTCAACCAgacgac                                                              >  1:2283414/1‑41 (MQ=255)
                                      ttCGGTAACCACCAGGTCCTGCGCCACGTCAACCAgacgac                                                              >  1:427744/1‑41 (MQ=255)
                                      ttCGGTAACCACCAGGTCCTGCGCCACGTCAACCAgacgac                                                              >  1:2151196/1‑41 (MQ=255)
                                      ttCGGTAACCACCAGGTCCTGCGCCACGTCAACCAgacgac                                                              >  1:454034/1‑41 (MQ=255)
                                      ttCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAAc                                                  >  1:281974/1‑53 (MQ=255)
                                        cGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAg                                                       >  1:1319297/1‑46 (MQ=255)
                                          gTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGTTTTCAGTg                             <  1:846624/70‑1 (MQ=255)
                                          gTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGTTTTCAGTg                             <  1:2128214/70‑1 (MQ=255)
                                           tcaccaccAGATCTTGCGTCACGTCCACCAGACGACGGGTc                                                         <  1:3434053/39‑1 (MQ=16)
                                           tcaccaccAGATCTTGCGTCACGTCCACCAGACGACGGGTc                                                         <  1:606162/39‑1 (MQ=16)
                                           tcaccaccAGATCTTGCGTCACGTCCACCAGACGACGGGTc                                                         <  1:196598/39‑1 (MQ=16)
                                           tcaccaccAGATCTTGCGTCACGTCCACCAGACGACGGGTc                                                         <  1:3329271/39‑1 (MQ=16)
                                           tcaccaccAGATCTTGCGTCACGTCCACCAGACGACGGGTc                                                         <  1:691904/39‑1 (MQ=16)
                                           tcaccaccAGATCTTGCGTCACGTCCACCAGACGACGGGTc                                                         <  1:801822/39‑1 (MQ=16)
                                           tAACCACCAGGTCCTGCGCCACGTCAACCAGAcgacga                                                            <  1:428827/38‑1 (MQ=255)
                                           tAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCgg                                               <  1:2749328/51‑1 (MQ=255)
                                           tAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCgg                                               <  1:1200365/51‑1 (MQ=255)
                                           tAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCgg                                               <  1:2208723/51‑1 (MQ=255)
                                           tAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCgg                                               <  1:1947422/51‑1 (MQ=255)
                                           tAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCgg                                               <  1:2330950/51‑1 (MQ=255)
                                           tAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCgg                                               <  1:2366973/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:2577199/1‑41 (MQ=19)
                                             accaccAGATCTTGCGTCACGTCCACCAGACGACGGGTCAg                                                       >  1:2836454/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:2045586/1‑68 (MQ=255)
                                                           cgcCACGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGTTTTCAGTGCGGTATCCGCCAGAc              >  1:1407179/1‑68 (MQ=255)
                                                           cgcCACGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGTTTTCAGTGCGGTATCCGCCAGAc              >  1:2030528/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:3264060/45‑1 (MQ=255)
                                                             ccACGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGttt                                   <  1:3006562/45‑1 (MQ=255)
                                                              cACGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGTTTTCAGTGCGGTATCCGCCAGACCTTTAc        <  1:3055801/71‑1 (MQ=255)
                                                               aCGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGTTTTc                                 >  1:943406/1‑45 (MQ=255)
                                                               aCGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGTTTTc                                 >  1:1068827/1‑45 (MQ=255)
                                                               aCGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGTTTTc                                 >  1:1354788/1‑45 (MQ=255)
                                                               aCGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGTTTTCa                                >  1:279901/1‑46 (MQ=255)
                                                               aCGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGTTTTCAGTGCGGTATCCGCCAGACCTTTa         >  1:218428/1‑69 (MQ=255)
                                                                cGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGTTTTCAGTGCGGTATCCGCCAGACCttt          <  1:2718938/67‑1 (MQ=255)
                                                                cGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGTTTTCAGTGCGGTATCCGCCAGACCttt          <  1:267509/67‑1 (MQ=255)
                                                                cGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGTTTTCAGTGCGGTATCCGCCAGACCttt          <  1:2648786/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)
                                                                    |                                                                      
GGAGTCATCATGATACCTTCATGGGTACCACAATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGTTTTCAGTGCGGTATCCGCCAGACCTTTACGAGCAC  >  minE/2254567‑2254705

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: