Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,309	G→A	11.6%	L649L (CTG→TTG)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,309	0	G	A	11.6%	162.6 / 10.4	60	L649L (CTG→TTG)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base G (25/28);  new base A (2/5);  total (27/33)
Fisher's exact test for biased strand distribution p-value = 4.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

GGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGC  >  minE/608239‑608378                                                                       |                                                                      ggAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCag                                                                       <  1:694873/71‑1 (MQ=255) ggAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCag                                                                       <  1:3333443/71‑1 (MQ=255) ggAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCag                                                                       <  1:2044591/71‑1 (MQ=255)  gAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCag                                                                       <  1:3827985/70‑1 (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)        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)             aaGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTg                                                            >  1:1123783/1‑70 (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)                           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:3643701/1‑70 (MQ=255)                                                        tAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTcc                 >  1:878034/1‑70 (MQ=255)                                                        tAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTcc                 >  1:2954541/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)                                                                 aaGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGc                          <  1:751043/52‑1 (MQ=255)                                                                 aaGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGc                          <  1:1505466/52‑1 (MQ=255)                                                                  aGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGgcgcgc      >  1:779293/1‑71 (MQ=255)                                                                   gTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCg                     <  1:2161746/55‑1 (MQ=255)                                                                      agagagGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGtt                             >  1:546803/1‑44 (MQ=255)                                                                      agagagGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCa    <  1:1961746/69‑1 (MQ=255)                                                                      agagagGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGc  <  1:2220336/71‑1 (MQ=255)                                                                       |                                                                      GGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGC  >  minE/608239‑608378

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