Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 608,927 T→G 20.0% M443L (ATG→CTG)  pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE608,9270TG20.0% 132.8 / 50.5 80M443L (ATG→CTG) pflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base T (38/26);  new base G (8/8);  total (46/34)
Fisher's exact test for biased strand distribution p-value = 5.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

ATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACG  >  minE/608857‑608994
                                                                      |                                                                   
aTCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCAt                                                                     <  1:3374726/71‑1 (MQ=255)
 tCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATg                                                                    >  1:3559816/1‑71 (MQ=255)
  cGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATgg                                                                   >  1:3071503/1‑71 (MQ=255)
   ggTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGt                                                                  <  1:1785360/71‑1 (MQ=255)
     ttCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATgg                                                                   <  1:1368281/68‑1 (MQ=255)
       cAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTc                                                              >  1:962584/1‑71 (MQ=255)
       cAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTc                                                              >  1:508069/1‑71 (MQ=255)
           cTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATg                                                                    >  1:1254041/1‑61 (MQ=255)
            ttCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAg                                                         >  1:3006853/1‑71 (MQ=255)
             tCGGCCCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAgg                                                        >  1:1510627/1‑71 (MQ=255)
             tCGGACCAACCTGCATTTTCAGTTTTTCGTCAACTCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAgg                                                        >  1:2552676/1‑71 (MQ=255)
                 aCCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTg                                                    <  1:1038274/71‑1 (MQ=255)
                 aCCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTg                                                    <  1:2792016/71‑1 (MQ=255)
                 aCCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTg                                                    <  1:2491041/71‑1 (MQ=255)
                 aCCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTg                                                    <  1:1971130/71‑1 (MQ=255)
                  ccgaCCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAgg                                                        >  1:3826218/4‑66 (MQ=255)
                   cgaCCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGtta                                                     >  1:1234575/3‑67 (MQ=255)
                    aaCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAc                                                 <  1:3040962/71‑1 (MQ=255)
                    aaCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAc                                                 <  1:2910320/71‑1 (MQ=255)
                    aaCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAc                                                 <  1:1016534/71‑1 (MQ=255)
                    aaCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAc                                                 <  1:2681248/71‑1 (MQ=255)
                     aCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAcg                                                >  1:1329070/1‑71 (MQ=255)
                      ccTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAcgc                                               >  1:2043102/1‑71 (MQ=255)
                        tGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAcgcg                                              >  1:24167/1‑70 (MQ=255)
                          cATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACgcgc                                             <  1:1895898/69‑1 (MQ=255)
                              ttCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCg                                         >  1:2381497/1‑69 (MQ=255)
                              ttCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCg                                         >  1:3654931/1‑69 (MQ=255)
                              ttCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCAcca                                         >  1:1494808/1‑68 (MQ=255)
                                cAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaag                                      <  1:27440/70‑1 (MQ=255)
                                cAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaag                                      <  1:102482/70‑1 (MQ=255)
                                cAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCAccaa                                        <  1:1916967/68‑3 (MQ=255)
                                cAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCAccaa                                        <  1:3865116/68‑3 (MQ=255)
                                cAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCAcca                                         <  1:3618961/67‑2 (MQ=255)
                                cAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCaaag                                      <  1:1401060/70‑1 (MQ=255)
                                cAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCaaag                                      <  1:1703505/70‑1 (MQ=255)
                                 aGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaagaa                                    >  1:595649/1‑71 (MQ=255)
                                 aGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaaga                                     <  1:2728383/70‑1 (MQ=255)
                                 aGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaaga                                     <  1:3160660/70‑1 (MQ=255)
                                  gTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAAc                                   <  1:2532648/71‑1 (MQ=255)
                                   tttttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGAACGCGCACCGAAGAACt                                  >  1:2008479/1‑71 (MQ=255)
                                    ttttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTg                                 >  1:3491701/1‑71 (MQ=255)
                                    ttttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTg                                 >  1:3852655/1‑71 (MQ=255)
                                     tttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAcg                                                >  1:1668574/1‑55 (MQ=255)
                                     tttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTg                                 >  1:3339379/1‑70 (MQ=255)
                                      ttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaag                                      >  1:229808/1‑64 (MQ=255)
                                       tCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAgg                                                        >  1:3681981/1‑45 (MQ=16)
                                       tCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAgg                                                        >  1:1831769/1‑45 (MQ=16)
                                       tCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCa                               <  1:2971132/70‑1 (MQ=255)
                                        cGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGATCTGCAt                              <  1:1270/70‑1 (MQ=255)
                                         gTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCAt                              <  1:52892/69‑1 (MQ=255)
                                          tCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTg                           <  1:2163221/71‑1 (MQ=255)
                                            aaCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGtt                         >  1:1856739/1‑71 (MQ=255)
                                             aCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGttt                        <  1:3815239/71‑1 (MQ=255)
                                              ccccgccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTa                       <  1:2023686/69‑1 (MQ=255)
                                              ccccgccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTa                       <  1:725150/69‑1 (MQ=255)
                                               gccgccGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCa                                                  >  1:1648793/1‑43 (MQ=255)
                                               gccgccGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCa                                                  >  1:2535825/1‑43 (MQ=255)
                                               gccgccGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTAc                      >  1:2708355/1‑71 (MQ=255)
                                                ccgccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCaaa                                       >  1:2585727/1‑53 (MQ=255)
                                                  gccgTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCaa                   >  1:128065/1‑71 (MQ=255)
                                                    cgTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAAc                  >  1:953683/1‑70 (MQ=255)
                                                    cgTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAAc                  >  1:3075564/1‑70 (MQ=255)
                                                    cgTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACg                 <  1:677044/71‑1 (MQ=255)
                                                       tGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCa                                                  >  1:2935016/1‑35 (MQ=255)
                                                       tGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGat               >  1:657558/1‑70 (MQ=255)
                                                        gATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAc                                                 >  1:2819634/1‑35 (MQ=255)
                                                        gATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatca             >  1:587614/1‑71 (MQ=255)
                                                        gATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatca             >  1:1852293/1‑71 (MQ=255)
                                                        gATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatc              <  1:4005748/70‑1 (MQ=255)
                                                        gATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTCGTTTACCAACGatca             >  1:685199/1‑71 (MQ=255)
                                                          ttGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatcatc           <  1:910646/71‑1 (MQ=255)
                                                          ttGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatcatc           <  1:1823578/71‑1 (MQ=255)
                                                            gCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCaaagaa                                    >  1:3238564/1‑44 (MQ=12)
                                                            gCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCa                                            >  1:2117140/1‑36 (MQ=255)
                                                                aCAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatcat            >  1:111947/1‑64 (MQ=255)
                                                                aCAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAAAGatcat            >  1:1450022/1‑64 (MQ=255)
                                                                   gCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCt      >  1:1288974/1‑67 (MQ=255)
                                                                   gCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACg  <  1:43335/71‑1 (MQ=255)
                                                                   gCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACg  <  1:2648335/71‑1 (MQ=255)
                                                                   gCATGGTTTTCGCCAGGTTTGCACCCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACg  >  1:3004311/1‑71 (MQ=255)
                                                                   gCATGGGTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCAttt                            >  1:538428/1‑45 (MQ=255)
                                                                      tGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCa                    <  1:1360880/50‑1 (MQ=255)
                                                                      tGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCgg         <  1:4123172/61‑1 (MQ=255)
                                                                      tGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCgg         <  1:1040688/61‑1 (MQ=255)
                                                                      tGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCgg         <  1:2636176/61‑1 (MQ=255)
                                                                      tGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCGTACg  >  1:3760010/1‑68 (MQ=255)
                                                                      |                                                                   
ATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACG  >  minE/608857‑608994

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