Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,124	C→T	12.3%	V710V (GTG→GTA)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,124	0	C	T	12.3%	219.8 / 17.6	81	V710V (GTG→GTA)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base C (36/35);  new base T (2/8);  total (38/43)
Fisher's exact test for biased strand distribution p-value = 9.42e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.43e-01

ACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCA  >  minE/608056‑608193                                                                     |                                                                      acacGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCAc                                                                     >  1:1149639/1‑71 (MQ=255)  cacGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCAc                                                                     <  1:251882/70‑1 (MQ=255)     ggATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGttg                                                                 <  1:1812077/71‑2 (MQ=255)     ggATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGtt                                                                  >  1:2612834/1‑70 (MQ=255)           tCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTc                                                           <  1:928150/71‑1 (MQ=255)            cAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGtt                                                            >  1:254837/1‑69 (MQ=255)             aGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAg                                                         <  1:3241491/71‑1 (MQ=255)             aGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAg                                                         >  1:1010417/1‑71 (MQ=255)             aGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAg                                                         >  1:3155460/1‑71 (MQ=255)              gCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtcatc                                                                       >  1:2916952/1‑56 (MQ=255)              gCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCa                                                          <  1:1912623/69‑1 (MQ=255)               cTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGt                                                       <  1:2538351/71‑1 (MQ=255)               cTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGt                                                       >  1:40621/1‑71 (MQ=255)               cTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGt                                                       <  1:3583132/71‑1 (MQ=255)                 taggATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTgt                                                     >  1:1934025/3‑70 (MQ=255)                 gCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGc                                                     >  1:3690568/1‑71 (MQ=255)                    gATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGa                                                  <  1:1929460/71‑1 (MQ=255)                    gATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGa                                                  >  1:616967/1‑71 (MQ=12)                       atTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGt                                                                   <  1:1328372/51‑1 (MQ=255)                       atTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGacca                                               >  1:564393/1‑71 (MQ=255)                        tttttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGaccac                                              <  1:2149654/71‑1 (MQ=255)                         ttttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGaccacc                                             >  1:855116/1‑71 (MQ=255)                         ttttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGaccacc                                             >  1:534188/1‑71 (MQ=255)                         ttttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGaccacc                                             >  1:3251582/1‑71 (MQ=255)                          tttCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGaccgcc                                             <  1:3278536/70‑1 (MQ=255)                              cGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCga                                        <  1:1278827/71‑1 (MQ=255)                              cGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCga                                        <  1:868636/71‑1 (MQ=255)                               gggTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgat                                       <  1:1627687/71‑1 (MQ=255)                                 gTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatg                                      >  1:1553947/1‑70 (MQ=255)                                  ttttCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAg                                                         <  1:237535/50‑1 (MQ=255)                                  ttttCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAg                                                         <  1:1953432/50‑1 (MQ=255)                                  ttttCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatgga                                    >  1:1659397/1‑71 (MQ=255)                                    ttCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatggat                                   >  1:1977977/1‑70 (MQ=255)                                     tCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGc                                 >  1:3689004/1‑71 (MQ=255)                                      ccATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTGAGGTGCTGACCACCTTCGATGGATGc                                 <  1:3419300/70‑1 (MQ=255)                                      ccATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGc                                 <  1:1213134/70‑1 (MQ=255)                                       cATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAgg                                                        <  1:904983/46‑1 (MQ=255)                                       cATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCtt                               >  1:3171065/1‑71 (MQ=255)                                        aTCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCtt                                           >  1:1568133/1‑58 (MQ=255)                                        aTCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCtt                                           >  1:1577473/1‑58 (MQ=255)                                        aTCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTc                              >  1:350886/1‑71 (MQ=255)                                         tCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCg                             >  1:136118/1‑71 (MQ=255)                                          cgcgTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAg                                                         >  1:1356558/1‑42 (MQ=255)                                           gcgTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTg                                                      >  1:3904503/1‑44 (MQ=255)                                           gcgTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTg                                                      >  1:2517574/1‑44 (MQ=255)                                           gcgTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTg                                                      >  1:1872647/1‑44 (MQ=255)                                             gTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGa                                                  >  1:3597212/1‑46 (MQ=255)                                             gTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCgtggt                         >  1:2158204/1‑71 (MQ=255)                                              tCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCga                                        <  1:2830291/55‑1 (MQ=255)                                                gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCt                                                    >  1:3094586/1‑41 (MQ=255)                                                gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTc                              <  1:37952/63‑1 (MQ=255)                                                gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGtggtgg                       <  1:1361846/70‑1 (MQ=255)                                                gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGtggtgg                       >  1:1551415/1‑70 (MQ=255)                                                gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGtggtgg                       <  1:1775066/70‑1 (MQ=255)                                                gAGCAGCAATTCACGGTTCATCACGTTAACGTTCAg                                                         >  1:1227772/1‑36 (MQ=255)                                                cagcagcATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGtggtgg                       <  1:3562883/69‑1 (MQ=255)                                                cagcagcATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGa                      <  1:1099635/70‑1 (MQ=255)                                                cagcagcATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGa                      <  1:1307701/70‑1 (MQ=255)                                                 agcagcATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGaccgcc                                             <  1:1502184/47‑1 (MQ=255)                                                 agcagcATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGaa                     <  1:2845961/71‑1 (MQ=12)                                                 agcagcATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGa                                                  >  1:2861022/1‑42 (MQ=255)                                                  gcagcaTTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGaaa                    <  1:2209198/71‑2 (MQ=255)                                                   cagcaTTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGt                   <  1:410139/71‑1 (MQ=255)                                                     gcaTTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTaa                 >  1:686513/1‑71 (MQ=255)                                                     gcaTTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTaa                 >  1:690240/1‑71 (MQ=255)                                                       aTTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAcc               <  1:3187811/71‑1 (MQ=255)                                                        tttCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCgtggt                         <  1:3859320/60‑1 (MQ=255)                                                        tttCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAcca              >  1:653763/1‑71 (MQ=255)                                                         ttCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccat             >  1:1077683/1‑71 (MQ=255)                                                          tCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatggat                                   <  1:2052940/48‑1 (MQ=255)                                                           cACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCgtggt                         >  1:923323/1‑57 (MQ=255)                                                            aCGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatcc           <  1:597418/70‑1 (MQ=255)                                                            aCGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatcc           <  1:1859894/70‑1 (MQ=255)                                                             cGGTTCATCACGTTAACGTTCAGGTGCTTACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatccat         <  1:3467065/71‑1 (MQ=255)                                                             cGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatccat         <  1:103763/71‑1 (MQ=255)                                                             cGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatccat         <  1:1082883/71‑1 (MQ=255)                                                               gTTCATCACGTTAACGTTCAGGTTCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCa       <  1:1865490/71‑1 (MQ=255)                                                               gTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatccat         >  1:2550594/1‑69 (MQ=255)                                                               gTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCa       <  1:3007658/71‑1 (MQ=255)                                                               gTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCa       <  1:743175/71‑1 (MQ=255)                                                               gTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCa       <  1:1093349/71‑1 (MQ=255)                                                                ttCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCa       <  1:400357/70‑1 (MQ=255)                                                                  catcaCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAc    >  1:2114380/1‑71 (MQ=255)                                                                  catcaCGTTAACGATCAGGTGCTGACCACCTACgatgga                                    >  1:2430576/1‑39 (MQ=37)                                                                   atcaCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAc    <  1:1054824/70‑1 (MQ=255)                                                                    tcaCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAcca  <  1:3408674/71‑1 (MQ=255)                                                                     caCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAg                    >  1:3902437/1‑52 (MQ=255)                                                                     |                                                                      ACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCA  >  minE/608056‑608193

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