Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 608,411 T→G 15.0% K615Q (AAG→CAG)  pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE608,4110TG15.0% 164.3 / 28.7 80K615Q (AAG→CAG) pflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base T (33/35);  new base G (7/5);  total (40/40)
Fisher's exact test for biased strand distribution p-value = 7.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.86e-01

GCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTG  >  minE/608342‑608477
                                                                     |                                                                  
gCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTa                                                                   >  1:4342200/1‑71 (MQ=255)
     gggTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATa                                                              <  1:1332168/71‑1 (MQ=255)
          aGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACaa                                                         >  1:4628740/1‑71 (MQ=255)
            cACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACg                                                       >  1:5220220/1‑71 (MQ=255)
            cACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACg                                                       <  1:3273572/71‑1 (MQ=255)
             aCCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGt                                                      <  1:1479695/71‑1 (MQ=255)
              cccGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGt                                                      >  1:2916497/1‑70 (MQ=255)
                cGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAg                                                   <  1:1196942/71‑1 (MQ=255)
                cGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAg                                                   <  1:4386550/71‑1 (MQ=255)
                    ccGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTAcc                                                                 >  1:321984/1‑53 (MQ=255)
                     cGAACGGTGTTCCGGCGCGACGACCGTCCGGCGTATTACCGGTTTTCTGGCCGTACACCACGTTAGAAGt                                               >  1:3633108/1‑70 (MQ=255)
                     cGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTg                                              <  1:899285/71‑1 (MQ=255)
                       aaCGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAg                                                   <  1:2862864/64‑1 (MQ=255)
                       aaCGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAg                                                   <  1:5010659/64‑1 (MQ=255)
                       aaCGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGAt                                            >  1:3351370/1‑71 (MQ=255)
                         cGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATgg                                          >  1:3361145/1‑71 (MQ=255)
                           gcgcgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGt                                         >  1:1669343/1‑70 (MQ=255)
                           gcgcgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGt                                         >  1:2058523/1‑70 (MQ=255)
                           gcgcgcCAGCACGACGACCGTCTGGGGTGTTAACCGTTTTCTTACCATACACAACGTTAGAAGTGATGGt                                         >  1:2981590/1‑70 (MQ=255)
                             gcgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATaca                                                            <  1:2340756/49‑1 (MQ=255)
                             gcgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATaca                                                            <  1:3983295/49‑1 (MQ=255)
                             gcgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGt                                                      <  1:210484/55‑1 (MQ=255)
                             gcgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGt                                                      <  1:4115035/55‑1 (MQ=255)
                              cgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAg                                      >  1:2490241/1‑70 (MQ=255)
                               gcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGaa                                    >  1:781163/1‑71 (MQ=255)
                                   gCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATacac                                                           <  1:1390336/44‑1 (MQ=255)
                                   gCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTg                                              <  1:392881/57‑1 (MQ=255)
                                   gCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGa                                >  1:1605881/1‑71 (MQ=255)
                                   gCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGa                                >  1:1599268/1‑71 (MQ=255)
                                   gCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGa                                >  1:2419006/1‑71 (MQ=255)
                                    cACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATacac                                                           <  1:3516012/43‑1 (MQ=255)
                                    cACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCa                                       >  1:2961006/1‑63 (MQ=255)
                                    cACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGAc                               <  1:5056468/71‑1 (MQ=255)
                                      cgacgaCCGTCCGGCGTATTACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTg                             <  1:985849/71‑1 (MQ=255)
                                      cgacgaCCGTCCGGCGTATTACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTg                             <  1:1131484/71‑1 (MQ=255)
                                       gacgacCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTATaa                                                 >  1:964193/1‑50 (MQ=255)
                                       gacgacCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACt                              >  1:779657/1‑69 (MQ=255)
                                        acgacCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAg                                                   <  1:5523749/47‑1 (MQ=255)
                                         cgacCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGt                          >  1:1411661/1‑71 (MQ=255)
                                          gacCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTc                         <  1:2463970/71‑1 (MQ=255)
                                           acCGGCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCg                        <  1:998088/71‑1 (MQ=255)
                                            ccGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTa                                                    >  1:2208333/1‑42 (MQ=255)
                                            ccGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCgg                       <  1:1327464/71‑1 (MQ=255)
                                            ccGTCCGGCGTATTACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTGGGTAgg                       >  1:3422674/1‑71 (MQ=255)
                                             cGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCggg                      <  1:3416414/71‑1 (MQ=255)
                                             cGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCggg                      <  1:858281/71‑1 (MQ=255)
                                             cGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCggg                      <  1:3959433/71‑1 (MQ=255)
                                              gTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACa                                  >  1:1857874/1‑58 (MQ=255)
                                              gTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACa                                  >  1:5010318/1‑58 (MQ=255)
                                                 tGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGaa                                    <  1:3540171/53‑1 (MQ=255)
                                                 tGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGaa                                    <  1:3858501/53‑1 (MQ=255)
                                                  ggggTGTTACCCGTTTTCTTACCATACACAACGTTAg                                                   <  1:1368028/37‑1 (MQ=255)
                                                        ttACCGGTTTTCTGGCCGTACACCACGTTAGAAGt                                               >  1:5545369/1‑35 (MQ=12)
                                                        ttACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGa                                             <  1:18036/37‑1 (MQ=12)
                                                        ttACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTGGGta                         >  1:2900815/1‑56 (MQ=255)
                                                        ttACCCGTTTTCTTACCATACACAACGTTAGAAGt                                               <  1:3801679/35‑1 (MQ=255)
                                                        ttACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGaa                                    <  1:1265460/46‑1 (MQ=255)
                                                        ttACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACg           >  1:4274323/1‑71 (MQ=255)
                                                        ttACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGCAt                    <  1:4533783/62‑1 (MQ=255)
                                                         tACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTGGGTAgg                       <  1:4226166/58‑1 (MQ=255)
                                                         tACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTGGGTAGGGACGGCGt               <  1:3955736/66‑1 (MQ=255)
                                                          aCCGGTTTTCTGGCCGTACACCACGTTAGAAGTGAta                                           >  1:4473590/1‑36 (MQ=12)
                                                          aCCCGTTTTCTTACCATACACAACGTTAGAAGTGa                                             >  1:3691654/1‑35 (MQ=255)
                                                          aCCCGTTTTCTTACCATACACAACGTTAGAAGTGa                                             >  1:1242288/1‑35 (MQ=255)
                                                          aCCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACgg          >  1:115970/1‑70 (MQ=255)
                                                              gTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACgg          >  1:3634530/1‑66 (MQ=255)
                                                              gTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACgg          >  1:1062633/1‑66 (MQ=255)
                                                              gTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGgt     >  1:4710746/1‑71 (MQ=255)
                                                              gTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGgt     >  1:2541054/1‑71 (MQ=255)
                                                              gTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGgt     >  1:2879597/1‑71 (MQ=255)
                                                              gTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTGGGTAGGGACGGCGTTGCGATAGgt     >  1:2102100/1‑71 (MQ=255)
                                                               ttttCTTACCATACACAACGTTAGAAGTGATGGTCAGAACa                                  <  1:5667625/41‑1 (MQ=255)
                                                               ttttCTTACCATACACAACGTTAGAAGTGATGGTCAGAACa                                  <  1:38625/41‑1 (MQ=255)
                                                               ttttCTTACCATACACAACGTTAGAAGTGATGGTCAGAACa                                  <  1:2378405/41‑1 (MQ=255)
                                                               ttttCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCAc            <  1:4317281/63‑1 (MQ=255)
                                                                tttCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACt                              >  1:4245046/1‑44 (MQ=12)
                                                                 ttCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGgtgt   <  1:4831545/70‑1 (MQ=255)
                                                                 ttCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGtgtg  >  1:1910693/1‑71 (MQ=255)
                                                                 ttCTTACCATACACAACGTTAGAAGAGATGGTCAGAACAGACTGAGTCGGGATAg                  >  1:3934166/1‑55 (MQ=255)
                                                                  tCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGa                                <  1:2128709/40‑1 (MQ=255)
                                                                     |                                                                  
GCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTG  >  minE/608342‑608477

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: