Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,301	G→C	28.5%	S651S (TCC→TCG)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,301	0	G	C	28.5%	108.5 / 54.1	77	S651S (TCC→TCG)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base G (29/26);  new base C (16/6);  total (45/32)
Fisher's exact test for biased strand distribution p-value = 1.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

GTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGC  >  minE/608231‑608370                                                                       |                                                                      gTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACg                                                                       >  1:5103641/1‑71 (MQ=255)  tGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACgg                                                                      <  1:666583/71‑1 (MQ=255)  tGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACg                                                                       >  1:3162435/1‑70 (MQ=255)  tGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACg                                                                       >  1:2595823/1‑70 (MQ=255)    cGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACgg                                                                      >  1:726501/1‑69 (MQ=255)    cGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGaa                                                                    <  1:184186/71‑1 (MQ=255)         ggAACGATTGAGAAGGTGTACGAGATCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCa                                                                >  1:3755348/1‑70 (MQ=255)         ggAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCag                                                               <  1:1390292/71‑1 (MQ=255)         ggAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCag                                                               <  1:5653632/71‑1 (MQ=255)           aaCGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCa                                                                >  1:4607131/1‑68 (MQ=255)            aCGATTGAGAAGGTGTACGAGATCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCaat                                                              >  1:513774/1‑67 (MQ=255)            aCGATTGAGAAGGTGTACGAGATCCCATATTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCaat                                                              >  1:4234740/1‑67 (MQ=255)              gATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGaa                                                                    <  1:5255274/61‑1 (MQ=255)              gATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGaa                                                                    <  1:4060079/61‑1 (MQ=255)               aTAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGc                                                         <  1:1327759/71‑1 (MQ=255)                ttgagaAGGTGTACGAGATCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGc                                                         >  1:5678493/3‑70 (MQ=255)                tAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGc                                                         >  1:3348759/1‑70 (MQ=255)                tAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGc                                                         >  1:1762916/1‑70 (MQ=255)                tAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGc                                                         <  1:4135293/70‑1 (MQ=255)                 agagAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTa                                                       >  1:4378990/1‑71 (MQ=255)                 agagAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTa                                                       >  1:5573347/1‑71 (MQ=255)                 agagAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTa                                                       >  1:2969125/1‑71 (MQ=255)                    gaAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTg                                                    <  1:920718/71‑1 (MQ=255)                     aaGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGc                                                   <  1:3033685/71‑1 (MQ=255)                     aaGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGc                                                   <  1:3026919/71‑1 (MQ=255)                        gtgtAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGc                                                         <  1:15115/62‑1 (MQ=255)                        gtgtAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCAcc                                                >  1:397857/1‑71 (MQ=255)                         tgtAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCa                                                                <  1:5068583/54‑1 (MQ=255)                         tgtAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCa                                                                <  1:2288181/54‑1 (MQ=255)                          gtAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCtt                                              >  1:267150/1‑71 (MQ=255)                           tAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCaga                                                              >  1:324205/1‑54 (MQ=255)                             ggAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCagaga                                                            >  1:5311233/1‑54 (MQ=255)                             ggAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCt                                           >  1:415944/1‑71 (MQ=255)                             cgagaTCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCa                                                                <  1:813430/49‑1 (MQ=255)                              gagaTCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCtttgcg                                          >  1:4114886/1‑68 (MQ=255)                              gagaTCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCtttgcg                                          >  1:1693457/1‑68 (MQ=255)                              gagaTCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCtttgcg                                          <  1:834402/71‑4 (MQ=255)                              gagaTACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTa                                                       >  1:4826911/1‑58 (MQ=255)                              gagaTACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTa                                                       >  1:3596066/1‑58 (MQ=255)                                 aTACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGt                                                                  <  1:3959556/44‑1 (MQ=255)                                   aCCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGa                                                                     >  1:585966/1‑39 (MQ=255)                                   aCCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCAc                                     >  1:5691076/1‑71 (MQ=255)                                    ccATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCAc                                     >  1:3321072/1‑70 (MQ=255)                                    ccATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCAc                                     >  1:3762863/1‑70 (MQ=255)                                    ccATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTgg                                         <  1:1633976/66‑1 (MQ=255)                                    ccATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCAc                                     >  1:526394/1‑70 (MQ=255)                                     cATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGa                                   >  1:489970/1‑71 (MQ=255)                                     cATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGa                                   >  1:196632/1‑71 (MQ=255)                                     cATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCaga                                                              <  1:3901048/44‑1 (MQ=255)                                      aTCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGAc                                  <  1:1258203/71‑1 (MQ=255)                                       tCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGATCACGAc                                  >  1:432341/1‑70 (MQ=255)                                        cTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGAcca                                >  1:5093733/1‑70 (MQ=255)                                         tttAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCaga                                                              >  1:4230043/1‑40 (MQ=255)                                         tttAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCaga                                                              <  1:1066885/40‑1 (MQ=255)                                         tttAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCag                                                               >  1:3917764/1‑39 (MQ=255)                                          ttAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGt                               >  1:3836571/1‑70 (MQ=255)                                          ttAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTg                              >  1:1192396/1‑71 (MQ=255)                                          ttAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCCTg                              >  1:4201025/1‑71 (MQ=255)                                           tAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGAcc                                 <  1:635675/67‑1 (MQ=255)                                           tAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGt                               <  1:2513253/69‑1 (MQ=255)                                           tAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTg                              >  1:3385579/1‑70 (MQ=255)                                              cGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATc                          <  1:3595044/71‑1 (MQ=255)                                               gTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCg                         >  1:5471197/1‑71 (MQ=255)                                                tAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCg                         <  1:2509113/70‑1 (MQ=255)                                                  aGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGAc                                  >  1:4945308/1‑59 (MQ=255)                                                  aGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGAc                                  >  1:3140718/1‑59 (MQ=255)                                                  aGCAAACGGCAGTTTAGAAACGGAAGTCAGAGAGGCTACTCCACCTTTCTGGTAACGACCGTGCATCGGGt                      >  1:4693906/1‑71 (MQ=255)                                                      aaCGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGCATCGGGTTAGc                  >  1:2058665/1‑71 (MQ=255)                                                       aCGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGCATCGGGTTAGcg                 >  1:321985/1‑70 (MQ=255)                                                       aCGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGCATCGGGTTAGcg                 <  1:945225/71‑2 (MQ=255)                                                       aCGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGCATCGGGTTAGcg                 >  1:2844125/1‑70 (MQ=255)                                                       aCGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGCATCGGGTTAGcg                 >  1:3138578/1‑70 (MQ=255)                                                       aCGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGt                                        >  1:2937546/1‑48 (MQ=255)                                                       aCGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGc                  >  1:1711072/1‑70 (MQ=255)                                                         ggCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGc                             <  1:3407023/57‑1 (MQ=255)                                                         ggCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGc                             <  1:2386775/57‑1 (MQ=255)                                                          gCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGCATCGGGTTAGCGccc              <  1:1432672/71‑1 (MQ=255)                                                            aGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCgg            <  1:3073683/71‑1 (MQ=255)                                                            aGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCgg            <  1:1762056/71‑1 (MQ=255)                                                                  gCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGaa      <  1:2120626/71‑1 (MQ=255)                                                                   cAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCAc                                     >  1:1451238/1‑39 (MQ=255)                                                                      cGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGgc  <  1:5113031/71‑1 (MQ=255)                                                                       ggAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCAccc              >  1:2368787/1‑58 (MQ=255)                                                                       ggAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCAccc              >  1:1922208/1‑58 (MQ=255)                                                                       |                                                                      GTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGC  >  minE/608231‑608370

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