Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 608,922 C→G 19.0% L444L (CTG→CTC)  pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE608,9220CG19.0% 190.9 / 36.6 84L444L (CTG→CTC) pflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base C (40/28);  new base G (8/8);  total (48/36)
Fisher's exact test for biased strand distribution p-value = 5.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATC  >  minE/608852‑608985
                                                                      |                                                               
cTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTAc                                                                 <  1:1548144/71‑1 (MQ=255)
cTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTAc                                                                 <  1:255345/71‑1 (MQ=255)
   tGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAg                                                               >  1:447005/1‑70 (MQ=255)
   tGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAg                                                               >  1:271456/1‑70 (MQ=255)
   tGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAg                                                               >  1:3899270/1‑70 (MQ=255)
   tGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAg                                                               >  1:344849/1‑70 (MQ=255)
   tGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAg                                                               >  1:3777530/1‑70 (MQ=255)
     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:508069/1‑71 (MQ=255)
            cAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTc                                                     >  1:962584/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:2792016/71‑1 (MQ=255)
                      aCCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTg                                           <  1:1038274/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:2681248/71‑1 (MQ=255)
                         aaCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAc                                        <  1:2910320/71‑1 (MQ=255)
                         aaCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAc                                        <  1:1016534/71‑1 (MQ=255)
                          aCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAcg                                       >  1:1329070/1‑71 (MQ=255)
                           ccTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAcgc                                      >  1:2043102/1‑71 (MQ=255)
                             tGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAg                                                               <  1:1021816/44‑1 (MQ=255)
                             tGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAg                                                               <  1:2522047/44‑1 (MQ=255)
                             tGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAg                                                               <  1:1945260/44‑1 (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:102482/70‑1 (MQ=255)
                                     cAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaag                             <  1:27440/70‑1 (MQ=255)
                                     cAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCAccaa                               <  1:3865116/68‑3 (MQ=255)
                                     cAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCAccaa                               <  1:1916967/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:3160660/70‑1 (MQ=255)
                                      aGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaaga                            <  1:2728383/70‑1 (MQ=255)
                                       gTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAAc                          <  1:2532648/71‑1 (MQ=255)
                                        tttttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGAACGCGCACCGAAGAACt                         >  1:2008479/1‑71 (MQ=255)
                                         ttttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTg                        >  1:3852655/1‑71 (MQ=255)
                                         ttttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTg                        >  1:3491701/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:1831769/1‑45 (MQ=16)
                                            tCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAgg                                               >  1:3681981/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:725150/69‑1 (MQ=255)
                                                   ccccgccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTa              <  1:2023686/69‑1 (MQ=255)
                                                    gccgccGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCa                                         >  1:2535825/1‑43 (MQ=255)
                                                    gccgccGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCa                                         >  1:1648793/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:1823578/71‑1 (MQ=255)
                                                               ttGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatcatc  <  1:910646/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)
                                                                      |                                                               
CTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATC  >  minE/608852‑608985

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: ‑