Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,301	G→C	12.3%	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	12.3%	191.5 / 15.5	73	S651S (TCC→TCG)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base G (33/31);  new base C (3/6);  total (36/37)
Fisher's exact test for biased strand distribution p-value = 4.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCG  >  minE/608231‑608371                                                                       |                                                                       gTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACg                                                                        <  1:3339879/71‑1 (MQ=255)  tGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACgg                                                                       <  1:269442/71‑1 (MQ=255)  tGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACgg                                                                       <  1:1722830/71‑1 (MQ=255)    cGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACgg                                                                       >  1:3867587/1‑69 (MQ=255)    cGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACgg                                                                       >  1:3834838/1‑69 (MQ=255)     gTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAg                                                                    >  1:369184/1‑71 (MQ=255)       tCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGt                                                                   >  1:361587/1‑70 (MQ=255)        cGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGt                                                                   >  1:3140270/1‑69 (MQ=255)         ggAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCag                                                                <  1:694873/71‑1 (MQ=255)         ggAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCag                                                                <  1:2044591/71‑1 (MQ=255)         ggAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCag                                                                <  1:3333443/71‑1 (MQ=255)          gAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCag                                                                <  1:3827985/70‑1 (MQ=255)           aaCGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCa                                                                 >  1:3352267/1‑68 (MQ=255)           aaCGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCa                                                                 >  1:1412528/1‑68 (MQ=255)           aaCGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCa                                                                 >  1:1982606/1‑68 (MQ=255)            aCGATTGAGAAGGTGTACGAGATCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCaat                                                               >  1:2992609/1‑67 (MQ=255)            aCGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCagaga                                                             <  1:324297/71‑1 (MQ=255)            aCGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCaga                                                               >  1:3302705/1‑69 (MQ=255)            aCGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCaga                                                               >  1:3283272/1‑69 (MQ=255)               attgagaAGGTGTACGAGATCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGc                                                          <  1:2742794/68‑1 (MQ=255)               attgagaAGGTGTACGAGATCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGc                                                          <  1:267180/68‑1 (MQ=255)               aTAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGc                                                          <  1:2445083/71‑1 (MQ=255)                tAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACgg                                                                       >  1:1607299/1‑57 (MQ=255)                tAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGc                                                          >  1:428323/1‑70 (MQ=255)                 agagAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCagag                                                              >  1:2221826/1‑65 (MQ=255)                   agaAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGc                                                          <  1:2212037/67‑1 (MQ=255)                   agaAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACt                                                      >  1:2131970/1‑71 (MQ=255)                     aaGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGaa                                                                     >  1:523529/1‑54 (MQ=255)                     aaGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTg                                                     >  1:1123783/1‑70 (MQ=255)                         tgtACGAGATCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGac                                                                     <  1:2837786/50‑2 (MQ=255)                          gtAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCtt                                               >  1:953217/1‑71 (MQ=255)                           tacgagaTCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCttt                                              <  1:1753908/68‑1 (MQ=255)                             ggAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTc                                             <  1:2401660/70‑1 (MQ=255)                                 aTACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTc                                        <  1:1654394/71‑1 (MQ=255)                                   aCCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGt                                                                   <  1:2029132/42‑1 (MQ=255)                                   aCCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTc                                        >  1:1331851/1‑69 (MQ=255)                                     cATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGa                                    >  1:720036/1‑71 (MQ=255)                                      aTCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGc                                                          <  1:3942664/48‑1 (MQ=255)                                      aTCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACATTTCTGGTCACGAc                                   >  1:3260619/1‑71 (MQ=255)                                       tCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGAc                                   >  1:1051499/1‑70 (MQ=255)                                         tttAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGt                                >  1:1980820/1‑71 (MQ=255)                                           tAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGTACCTTTCTGGTCACGACCGTg                               >  1:87529/1‑70 (MQ=255)                                           tAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGAcc                                  <  1:2037412/67‑1 (MQ=255)                                              cGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATc                           <  1:3250035/71‑1 (MQ=255)                                              cGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTAACGACCGTGc                              <  1:913851/68‑1 (MQ=255)                                              cGTAAGAAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATc                           <  1:3164314/71‑1 (MQ=255)                                                tAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTc                                             <  1:3559758/51‑1 (MQ=255)                                                tAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCgg                         >  1:1019022/1‑71 (MQ=255)                                                tAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTAATg                                                     <  1:355408/43‑1 (MQ=255)                                                  aGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGc                                                    >  1:2251034/1‑42 (MQ=255)                                                  aGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCATCTTTCTGGTCACGACCGTGCATCGGGt                       <  1:3696793/71‑1 (MQ=255)                                                  aGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGt                       >  1:850934/1‑71 (MQ=255)                                                  aGCAAACGGCAGTTTAGCAACGGAAGGCAGCGCGGCTACTGCACCTTTCTGGTCACGACCGTGCATCggg                        <  1:3344432/70‑1 (MQ=255)                                                   gtgaaCGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGCATCGGGtt                      <  1:2517230/68‑1 (MQ=255)                                                   gCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGt                                         <  1:1104843/52‑1 (MQ=255)                                                   gCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCgg                         >  1:1086037/1‑68 (MQ=255)                                                   gCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGtt                      <  1:3646045/71‑1 (MQ=255)                                                      aaCGGCAGTTTAGCAACGGAAGTCAGAGAGGCTAc                                                       >  1:1135848/1‑35 (MQ=255)                                                      aaCGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCa                                                   >  1:1794823/1‑39 (MQ=255)                                                      aaCGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATc                           <  1:1799456/63‑1 (MQ=255)                                                       aCGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCa                  <  1:552843/71‑1 (MQ=255)                                                        cGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGCATCGGGTTAGcgc                 >  1:2072388/1‑69 (MQ=255)                                                         ggCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCAcc                <  1:1770690/71‑1 (MQ=255)                                                            aGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCgg             <  1:1867519/71‑1 (MQ=255)                                                            aGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCgg             <  1:539325/71‑1 (MQ=255)                                                             gTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCCTGCATCGGGTTAGCGCCCggc            <  1:3263540/71‑2 (MQ=255)                                                              tttAGCATCGGAAGTCTGAGAGGCTACTGCACCTTTc                                             <  1:2397951/37‑1 (MQ=255)                                                              tttAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGAcc                                  <  1:2470539/48‑1 (MQ=255)                                                              tttAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTc           >  1:3031956/1‑71 (MQ=255)                                                                tAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTcc          >  1:2954541/1‑70 (MQ=255)                                                                tAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTcc          >  1:3643701/1‑70 (MQ=255)                                                                tAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTcc          >  1:878034/1‑70 (MQ=255)                                                                tAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTc           >  1:2297856/1‑69 (MQ=255)                                                                 aGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCAcc                >  1:323299/1‑63 (MQ=255)                                                                  gCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGa        <  1:1527459/70‑1 (MQ=255)                                                                    aaCGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACg     <  1:2364158/71‑1 (MQ=255)                                                                     aCGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACgg    >  1:3186160/1‑71 (MQ=255)                                                                     aCGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACg     <  1:3366397/70‑1 (MQ=255)                                                                      cGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGgc   >  1:1922/1‑71 (MQ=255)                                                                       ggAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGgcg  >  1:3165444/1‑71 (MQ=255)                                                                       |                                                                       GTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCG  >  minE/608231‑608371

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