Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,134,927	T→C	29.1%	N289N (AAT→AAC)	pfkB →	6‑phosphofructokinase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,134,927	0	T	C	29.1%	93.6 / 43.1	96	N289N (AAT→AAC)	pfkB	6‑phosphofructokinase II
Reads supporting (aligned to +/- strand):  ref base T (35/33);  new base C (12/16);  total (47/49)
Fisher's exact test for biased strand distribution p-value = 5.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

AAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTA  >  minE/1134864‑1134989                                                                |                                                               aaaTGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAggg                                                           >  1:1629943/1‑69 (MQ=255) aaaTGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAggg                                                           >  1:148575/1‑69 (MQ=255)              gAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGGCACACTCAATCAGGGAACACGTCtg                                                 <  1:309365/66‑1 (MQ=255)              gAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCtg                                                 <  1:3171414/66‑1 (MQ=255)              gAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCtg                                                 <  1:1088780/66‑1 (MQ=255)              gAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCtg                                                 <  1:652642/66‑1 (MQ=255)              gAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCtg                                                 <  1:534394/66‑1 (MQ=255)              gAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCtg                                                 <  1:516652/66‑1 (MQ=255)              gAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCtg                                                 <  1:2087265/66‑1 (MQ=255)              gAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCtg                                                 <  1:3016281/66‑1 (MQ=255)              gAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCtg                                                 <  1:3259856/66‑1 (MQ=255)              gAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCtg                                                 <  1:1546377/66‑1 (MQ=255)              gAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCtg                                                 <  1:2431450/66‑1 (MQ=255)              gAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGCAACACGTCtg                                                 <  1:2469919/66‑1 (MQ=255)               aaGTGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCtg                                                 <  1:3378984/65‑1 (MQ=255)                   gaTGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCtg                                                 <  1:15178/61‑1 (MQ=255)                    aTGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCAt                                        <  1:2472731/69‑1 (MQ=255)                    aTGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCAt                                        <  1:1638400/69‑1 (MQ=255)                    aTGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCAt                                        <  1:1925629/69‑1 (MQ=255)                    aTGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCAt                                        <  1:2476560/69‑1 (MQ=255)                    aTGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCAt                                        <  1:1574311/69‑1 (MQ=255)                    aTGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCAt                                        <  1:3266520/69‑1 (MQ=255)                    aTGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCAt                                        <  1:1417829/69‑1 (MQ=255)                    aTGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCAt                                        <  1:1011620/69‑1 (MQ=255)                    aTGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCAt                                        <  1:405485/69‑1 (MQ=255)                    aTGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCAt                                        <  1:45504/69‑1 (MQ=255)                     tGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCAt                                        <  1:394494/68‑1 (MQ=255)                     tGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCAt                                        <  1:2688317/68‑1 (MQ=255)                     tGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCAt                                        <  1:2062571/68‑1 (MQ=255)                                         gCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTAc                   <  1:2105545/69‑1 (MQ=255)                                         gCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTAc                   <  1:1528046/69‑1 (MQ=255)                                         gCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTAc                   <  1:1331029/69‑1 (MQ=255)                                               aGTGCAGCCACACTCAATCAGGGAACCCGTCTGTGCTCCCATGACGATACGCAAAAAAt                       >  1:2401773/1‑59 (MQ=255)                                               aGTGCAGCCACACTCAATCAGGGAACACGTCTGTTCTCCCATGACGATACGCAAAAAAtt                      >  1:2274439/1‑60 (MQ=255)                                               aGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATg                                       >  1:2818462/1‑43 (MQ=255)                                               aGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAAtt                      >  1:2256606/1‑60 (MQ=255)                                               aGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAAtt                      >  1:3152480/1‑60 (MQ=255)                                               aGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAAtt                      >  1:3105823/1‑60 (MQ=255)                                               aGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAAtt                      >  1:2406000/1‑60 (MQ=255)                                               aGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAAtt                      >  1:975656/1‑60 (MQ=255)                                               aGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAAtt                      >  1:1176754/1‑60 (MQ=255)                                               aGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAAtt                      >  1:1555209/1‑60 (MQ=255)                                               aGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAAtt                      >  1:2071074/1‑60 (MQ=255)                                               aGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAAtt                      >  1:1562478/1‑60 (MQ=255)                                               aGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAAtt                      >  1:1580307/1‑60 (MQ=255)                                               aGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAAtt                      >  1:1629012/1‑60 (MQ=255)                                               aGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTAc             <  1:2864480/69‑1 (MQ=255)                                               aGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTAc             <  1:2455424/69‑1 (MQ=255)                                               aGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTAc             <  1:2459310/69‑1 (MQ=255)                                               aGTGCAGACACACTCAATCAGGAAACACGTCTGTGCTCCCATGACGATACGCAAAAAAtt                      >  1:1857426/1‑60 (MQ=255)                                                   cAGCCACACTCAATCAGGGCACACGTCTGTGCTCCCATGACGATACGCa                             >  1:1312604/1‑49 (MQ=255)                                                   cAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCttt         >  1:1630557/1‑69 (MQ=255)                                                   cAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCttt         >  1:2961855/1‑69 (MQ=255)                                                   cAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCttt         >  1:2988986/1‑69 (MQ=255)                                                   cAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCttt         >  1:3033584/1‑69 (MQ=255)                                                   cAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCttt         >  1:1705146/1‑69 (MQ=255)                                                   cAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCttt         >  1:2205899/1‑69 (MQ=255)                                                   cAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCttt         >  1:1133640/1‑69 (MQ=255)                                                   cAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCttt         >  1:2267035/1‑69 (MQ=255)                                                   cAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCttt         >  1:2237425/1‑69 (MQ=255)                                                   cAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCttt         >  1:2569588/1‑69 (MQ=255)                                                       cacaCTCAACCAGGGAACACGCCTGTGCTCCCATGACGATACGCAAAAAATTTTCGCGTACCttt         >  1:2918282/1‑65 (MQ=255)                                                       cacaCTCAACCAGGGAACACGCCTGTGCTCCCATGACGATACGCAAAAAATTTGCGCGTACCttt         >  1:2922207/1‑65 (MQ=255)                                                       cacaCTCAACCAGGGAACACGCCTGTGCTCCCATGACGATACGCAAAAAATTTACGCGTACCttt         >  1:1314452/1‑65 (MQ=255)                                                       cacaCTCAACCAGGGAACACGCCTGTGCTCCCATGACGATACGCAAAAAATTTACGCGTACCttt         >  1:3183944/1‑65 (MQ=255)                                                       cacaCTCAACCAGGGAACACGCCTGTGCTCCCATGACGATACGCAAAAAATTTACGCGTACCttt         >  1:3134046/1‑65 (MQ=255)                                                       cacaCTCAACCAGGGAACACGCCTGTGCTCCCATGACGATACGCAAAAAATTTACGCGTACCttt         >  1:1943008/1‑65 (MQ=255)                                                       cacaCTCAACCAGGGAACACGCCTGTGCTCCCATGACGATACGCAAAAAATTTACGCGTACCttt         >  1:1535980/1‑65 (MQ=255)                                                       cacaCTCAACCAGGGAACACGCCTGTGCTCCCATGACGATACGCAAAAAATTTACGCGTACCttt         >  1:1435569/1‑65 (MQ=255)                                                       cacaCTCAACCAGGGAACACGCCTGTGCTCCCATGACGATACGCAAAAAATTTACGCGTACCttt         >  1:2624846/1‑65 (MQ=255)                                                       cacaCTCAACCAGGGAACACGCCTGTGCTCCCATGACGATACGCAAAAAATTTACGCGTACCttt         >  1:724271/1‑65 (MQ=255)                                                       cacaCTCAACCAGGGAACACGCCTGTGCTCCCATGACGATACGCAAAAAATTTACGCGTACCttt         >  1:2803302/1‑65 (MQ=255)                                                       cacaCTCAACCAGGGAACACGCCTGTGCTCCCATGACGATACGCAAAAAATTTAAGCGTACCttt         >  1:343078/1‑65 (MQ=255)                                                           cTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTa  >  1:1024998/1‑68 (MQ=255)                                                           cTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTa  >  1:408335/1‑68 (MQ=255)                                                           cTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTa  >  1:718639/1‑68 (MQ=255)                                                           cTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTa  >  1:1788760/1‑68 (MQ=255)                                                           cTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTa  >  1:2817065/1‑68 (MQ=255)                                                           cTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTa  >  1:1783615/1‑68 (MQ=255)                                                           cTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTAACCTTTCCCGCTa  >  1:2805931/1‑68 (MQ=255)                                                             cAACCAGGGAACACGCCTGTGCTCCCTTGACGATACGCAAAAAATTTACGCGTACCTTTc        <  1:3004209/60‑1 (MQ=255)                                                             cAACCAGGGAACACGCCTGTGCTCCCATGACGATACGCAAAAAATTTACGCGTACCTTTc        <  1:2923810/60‑1 (MQ=255)                                                             cAACCAGGGAACACGCCTGTGCTCCCATGACGATACGCAAAAAATTTACGCGTACCTTTc        <  1:850778/60‑1 (MQ=255)                                                             cAACCAGGGAACACGCCTGTGCTCCCATGACGATACGCAAAAAATTTACGCGTACCTTTc        <  1:2140151/60‑1 (MQ=255)                                                             cAACCAGGGAACACGCCTGTGCTCCCATGACGATACGCAAAAAATTTACGCGTACCTTTc        <  1:2648244/60‑1 (MQ=255)                                                             cAACCAGGGAACACGCCTGTGCTCCCATGACGATACGCAAAAAATTTACGCGTACCTTTc        <  1:1260874/60‑1 (MQ=255)                                                             cAACCAGGGAACACGCCTGTGCTCCCATGACGATACGCAAAAAATTTACGCGTACCTTTc        <  1:2077283/60‑1 (MQ=255)                                                             cAACCAGGGAACACGCCTGTGCTCCCATGACGATACGCAAAAAATTTACGCGTACCTTTc        <  1:2889600/60‑1 (MQ=255)                                                             cAACCAGGGAACACGCCTGTGCTCCCATGACGATACGCAAAAAATTTACGCGTACCTTTc        <  1:1997102/60‑1 (MQ=255)                                                             cAACCAGGGAACACGCCTGTGCTCCCATGACGATACGCAAAAAATTTACGCGTACCTTTc        <  1:1735143/60‑1 (MQ=255)                                                             cAACCAGGGAACACGCCTGTGCTCCCATGACGATACGCAAAAAATTTACGCGTACCTTTc        <  1:3425308/60‑1 (MQ=255)                                                             cAACCAGGGAACACGCCTGTGCTCCCATGACGATACGCAAAAAATTTACGCGTACCTTTc        <  1:1910738/60‑1 (MQ=255)                                                             cAACCAGGGAACACGCCTGTGCTCCCATGACGATACGCAAAAAATTTACGCGTACCTTTc        <  1:1938212/60‑1 (MQ=255)                                                             cAACCAGGGAACACGCCTGTGCTCCCATGACGATACGCAAAAAATTTACGCGTACCTTTc        <  1:2282932/60‑1 (MQ=255)                                                             cAACCAGGGAACACGCCTGTGCTCCCATGACGATACGCAAAAAATTTACGCGTACCTTTc        <  1:3211976/60‑1 (MQ=255)                                                             cAACCAGGGAACACGCCTGTGCTCCCATCACGATACGCAAAAAATTTACGCGTACCTTTc        <  1:1974878/60‑1 (MQ=255)                                                                |                                                               AAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTA  >  minE/1134864‑1134989

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