Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	812,424	A→G	40.2%	F198F (TTT→TTC)	ychM ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	812,424	0	A	G	40.2%	27.8 / 49.4	72	F198F (TTT→TTC)	ychM	predicted transporter
Reads supporting (aligned to +/- strand):  ref base A (18/25);  new base G (15/14);  total (33/39)
Fisher's exact test for biased strand distribution p-value = 4.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACCGCGCAACCAGCCAGCAAGGCCGGAAGGTGACCGGGTAAACGAATGCCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGGCAGCATCACCCACATTAATGGTCGGCAGCGC  >  minE/812361‑812489                                                                |                                                                  cACCGCGCAACCAGCCAGCAAGGCCGGAAGGTGACCGGGTAAACGAATGCCCAGACGCGGCCaaaaa                                                                >  1:2932870/1‑67 (MQ=255) cACCGCGCAACCAGCCAGCAAGGCCGGAAGGTGACCGGGTAAACGAATGCCCAGACGCGGCCaaaaa                                                                >  1:2097163/1‑67 (MQ=255) cACCGCGCAACCAGCCAGCAAGGCCGGAAGGTGACCGGGTAAACGAATGCCCAGACGCGGCCaaaaa                                                                >  1:2477242/1‑67 (MQ=255) cACCGCGCAACCAGCCAGCAAGGCCGGAAGGTGACCGGGTAAACGAATGCCCAGACGCGGCCaaaaa                                                                >  1:2512457/1‑67 (MQ=255) cACCGCGCAACCAGCCAGCAAGGCCGGAAGGTGACCGGGTAAACGAATGCCCAGACGCGGCCaaaaa                                                                >  1:277559/1‑67 (MQ=255) cACCGCGCAACCAGCCAGCAAGGCCGGAAGGTGACCGGGTAAACGAATGCCCAGACGCGGCCaaaaa                                                                >  1:2787880/1‑67 (MQ=255) cACCGCGCAACCAGCCAGCAAGGCCGGAAGGTGACCGGGTAAACGAATGCCCAGACGCGGCCaaaaa                                                                >  1:1716677/1‑67 (MQ=255) cACCGCGCAACCAGCCAGCAAGGCCGGAAGGTGACCGGGTAAACGAATGCCCAGACGCGGCCaaaaa                                                                >  1:2993258/1‑67 (MQ=255) cACCGCGCAACCAGCCAGCAAGGCCGGAAGGTGACCGGGTAAACGAATGCCCAGACGCGGCCaaaaa                                                                >  1:3121556/1‑67 (MQ=255)  aCCGCGCAACCAGCCAGCAAGGCCGGAAGGTGACCGGGTAAACGAATGCCCAGACGCGGCCaaaaa                                                                >  1:2829490/1‑66 (MQ=255)  aCCGCGCAACCAGCCAGCAAGGCCGGAAGGTGACCGGGTAAACGAATGCCCAGACGCGGCCaaaaa                                                                >  1:919948/1‑66 (MQ=255)  aCCGCGCAACCAGCCAGCAAGGCCGGAAGGTGACCGGGTAAACGAATGCCCAGACGCGGCCaaaaa                                                                >  1:2222019/1‑66 (MQ=255)  aCCGCGCAACCAGCCAGCAAGGCCGGAAGGTGACCGGGTAAACGAATGCCCAGACGCGGCCaaaaa                                                                >  1:2474276/1‑66 (MQ=255)  aCCGCGCAACCAGCCAGCAAGGCCGGAAGGTGACCGGGTAAACGAATGCCCAGACGCGGCCaaaaa                                                                >  1:2864015/1‑66 (MQ=255)  aCCGCGCAACCAGCCAGCAAGGCCGGAAGGTGACCGGGTAAACGAATGCCCAGACGCGGCCaaaaa                                                                >  1:3300146/1‑66 (MQ=255)  aCCGCGCAACCAGCCAGCAAGGCCGGAAGGTGACCGGGTAAACGAATGCCCAGACGCGGCCaaaaa                                                                >  1:2977257/1‑66 (MQ=255)  aCCGCGCAACCAGCCAGCAAGGCCGGAAGGTGACCGGGTAAACGAATGCCCAGACGCGGCCaaaaa                                                                >  1:3165318/1‑66 (MQ=255)  aCCGCGCAACCAGCCAGCAAGGCCGGAAGGTGACCGGGTAAACGAATGCCCAGACGCGGCCaaaaa                                                                >  1:3104852/1‑66 (MQ=255)  aCCGCGCAACCAGCCAGCAAGGCCGGAAGGTGACCGGGTAAACGAATGCCCAGACGCGGCCaaaa                                                                 >  1:1778127/1‑65 (MQ=255)   ccGCGCAACCAGCCAGCAAGGCCGGAAGGTGGCCGGGTAAACGTATGCCCAGACGCGGCCAGAAAAcc                                                             >  1:287222/1‑67 (MQ=255)   ccGCGCAACCAGCCAGCAAGGCCGGAAGGTGGCCGGGTAAACGTATGCCCAGACGCGGCCAGAAAAcc                                                             >  1:1184732/1‑67 (MQ=255)   ccGCGCAACCAGCCAGCAAGGCCGGAAGGTGGCCGGGTAAACGTATGCCCAGACGCGGCCAGAAAAcc                                                             >  1:33237/1‑67 (MQ=255)   ccGCGCAACCAGCCAGCAAGGCCGGAAGGTGGCCGGGTAAACGTATGCCCAGACGCGGCCAGAAAAcc                                                             >  1:595848/1‑67 (MQ=255)   ccGCGCAACCAGCCAGCAAGGCCGGAAGGTGGCCGGGTAAACGTATGCCCAGACGCGGCCAGAAAAcc                                                             >  1:274244/1‑67 (MQ=255)   ccGCGCAACCAGCCAGCAAGGCCGGAAGGTGGCCGGGTAAACGTATGCCCAGACGCGGCCAGAAAAcc                                                             >  1:2548429/1‑67 (MQ=255)   ccGCGCAACCAGCCAGCAAGGCCGGAAGGTGGCCGGGTAAACGTATGCCCAGACGCGGCCAGAAAAcc                                                             >  1:611455/1‑67 (MQ=255)   ccGCGCAACCAGCCAGCAAGGCCGGAAGGTGGCCGGGTAAACGTATGCCCAGACGCGGCCAGAAAAcc                                                             >  1:746624/1‑67 (MQ=255)   ccGCGCAACCAGCCAGCAAGGCCGGAAGGTGGCCGGGTAAACGTATGCCCAGACGCGGCCAGAAAAcc                                                             >  1:1746619/1‑67 (MQ=255)   ccGCGCAACCAGCCAGCAAGGCCGGAAGGTGGCCGGGTAAACGTATGCCCAGACGCGGCCAGAAAAcc                                                             >  1:1631323/1‑67 (MQ=255)   ccGCGCAACCAGCCAGCAAGGCCGGAAGGTGGCCGGGTAAACGTATGCCCAGACGCGGCCAGAAAAcc                                                             >  1:1668242/1‑67 (MQ=255)   ccGCGCAACCAGCCAGCAAGGCCGGAAGGTGGCCGGGTAAACGTATGCCCAGACGCGGCCAGAAAAcc                                                             >  1:189763/1‑67 (MQ=255)   ccGCGCAACCAGCCAGCAAGGCCGGAAGGTGGCCGGGTAAACGTATGCCCAGACGCGGCCAGAAAAcc                                                             >  1:1310515/1‑67 (MQ=255)   ccGCGCAACCAGCCAGCAAGGCCGGAAGGTGGCCGGGTAAACGTATGCCCAGACGCGGCCAGAAAAcc                                                             >  1:1248065/1‑67 (MQ=255)   ccGCGCAACCAGCCAGCAAGGCCGGAAGGTGGCCGGGTAAACGTATGCCCAGACGCGGCCAGAAAAc                                                              >  1:2946269/1‑67 (MQ=255)              gccagcAAGGCCGGAAGGTGACCGGGTAAACGAATGCCCAGACGCGGCCAAAAAACAAGAATACCTAGc                                                 >  1:3399329/1‑69 (MQ=255)              gccagcAAGGCCGGAAGGTGACCGGGTAAACGAATGCCCAGACGCGGCCAAAAAACAAGAATACCTAGc                                                 >  1:3418701/1‑69 (MQ=255)              gccagcAAGGCCGGAAGGTGACCGGGTAAACGAATGCCCAGACGCGGCCAAAAAACAAGAATACCTAGc                                                 >  1:1842242/1‑69 (MQ=255)              gccagcAAGGCCGGAAGGTGACCGGGTAAACGAATGCCCAGACGCGGCCAAAAAACAAGAATACCTAGc                                                 >  1:3098513/1‑69 (MQ=255)                    aaGGCCGGAAGGTGACCGGGTAAACGAATGCCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTcaca                                           >  1:3003048/1‑69 (MQ=255)                    aaGGCCGGAAGGTGACCGGGTAAACGAATGCCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTcaca                                           >  1:301105/1‑69 (MQ=255)                    aaGGCCGGAAGGTGACCGGGTAAACGAATGCCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTcaca                                           >  1:2466360/1‑69 (MQ=255)                    aaGGCCGGAAGGTGACCGGGTAAACGAATGCCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTcaca                                           >  1:2889886/1‑69 (MQ=255)                    aaGGCCGGAAGGTGACCGGGTAAACGAATGCCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTcaca                                           >  1:1865850/1‑69 (MQ=255)                    aaGGCCGGAAGGTGACCGGGTAAACGAATGCCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTca                                             >  1:2491363/1‑67 (MQ=255)                        ccGGAAGGTGACCGGGTAAACGAATGCCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTcac                                            >  1:2579757/1‑64 (MQ=255)                        ccGGAAGGTGACCGGGTAAACGAATGCCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTcac                                            >  1:3123829/1‑64 (MQ=255)                               gTGACCGGGTAAACGAATGCCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGgc                                <  1:331585/69‑1 (MQ=255)                               gTGACCGGGTAAACGAATGCCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGgc                                <  1:132444/69‑1 (MQ=255)                               gTGACCGGGTAAACGAATGCCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGgc                                <  1:318787/69‑1 (MQ=255)                               gTGACCGGGTAAACGAATGCCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGgc                                <  1:3422954/69‑1 (MQ=255)                               gTGACCGGGTAAACGAATGCCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGgc                                <  1:1086250/69‑1 (MQ=255)                               gTGACCGGGTAAACGAATGCCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGgc                                <  1:39534/69‑1 (MQ=255)                               gTGACCGGGTAAACGAATGCCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGgc                                <  1:2567122/69‑1 (MQ=255)                               gTGACCGGGTAAACGAATGCCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGgc                                <  1:2007176/69‑1 (MQ=255)                               gTGACCGGGTAAACGAATGCCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGgc                                <  1:1903701/69‑1 (MQ=255)                               gTGACCGGGTAAACGAATGCCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGgc                                <  1:871850/69‑1 (MQ=255)                               gTGACCGGGTAAACGAATGCCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGgc                                <  1:2472342/69‑1 (MQ=255)                                tGACCGGGTAAACGAATGCCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGgc                                >  1:571854/1‑68 (MQ=255)                                tGACCGGGTAAACGAATGCCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGgc                                >  1:2045735/1‑68 (MQ=255)                                tGACCGGGTAAACGAATGCCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGgc                                >  1:886558/1‑68 (MQ=255)                                tGACCGGGTAAACGAATGCCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGgc                                >  1:459599/1‑68 (MQ=255)                                tGACCGGGTAAACGAATGCCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGgc                                >  1:3373120/1‑68 (MQ=255)                                tGACCGGGTAAACGAATGCCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGgc                                >  1:277264/1‑68 (MQ=255)                                   ccGGGTAAACGAATGCCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGgc                                <  1:3017946/65‑1 (MQ=255)                                        tAAACGAATGCCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGGCAGCATCAcc                       <  1:1228475/69‑1 (MQ=255)                                        tAAACGAATGCCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGGCAGCATCAcc                       <  1:20479/69‑1 (MQ=255)                                        tAAACGAATGCCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGGCAGCATCAcc                       <  1:378520/69‑1 (MQ=255)                                        tAAACGAATGCCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGGCAGCATCAcc                       <  1:3220492/69‑1 (MQ=255)                                        tAAACGAATGCCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGGCAGCATCAcc                       <  1:1990092/69‑1 (MQ=255)                                        tAAACGAATGCCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGGCAGCATCAcc                       <  1:2826569/69‑1 (MQ=255)                                        tAAACGAATGCCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGGCAGCATCAcc                       <  1:1896559/69‑1 (MQ=255)                                              aaTGCCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGGCAGCATCACCCACAtt                 <  1:2052701/69‑1 (MQ=255)                                              aaTGCCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGGCAGCATCACCCACAtt                 <  1:1667210/69‑1 (MQ=255)                                              aaTGCCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGGCAGCATCACCCACAtt                 <  1:2084187/69‑1 (MQ=255)                                              aaTGCCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGGCAGCATCACCCACAtt                 <  1:624749/69‑1 (MQ=255)                                               aTGCCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGGCAGCATCACCCACAtt                 <  1:81351/68‑1 (MQ=255)                                               aTGCCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGGCAGCATCACCCACAtt                 <  1:1919904/68‑1 (MQ=255)                                                             ccAGAAAACCAGAATGCCAAGCGTCACAATGCCAATGGCTGCATCACCCACATTAATGGTCGGCAgcgc  <  1:997494/69‑1 (MQ=255)                                                             ccAGAAAACCAGAATGCCAAGCGTCACAATGCCAATGGCTGCATCACCCACATTAATGGTCGGCAgcgc  <  1:910872/69‑1 (MQ=255)                                                             ccAGAAAACCAGAATGCCAAGCGTCACAATGCCAATGGCTGCATCACCCACATTAATGGTCGGCAgcgc  <  1:2273510/69‑1 (MQ=255)                                                             ccAGAAAACCAGAATGCCAAGCGTCACAATGCCAATGGCTGCATCACCCACATTAATGGTCGGCAgcgc  <  1:2474093/69‑1 (MQ=255)                                                             ccAGAAAACCAGAATGCCAAGCGTCACAATGCCAATGGCTGCATCACCCACATTAATGGTCGGCAgcgc  <  1:1280675/69‑1 (MQ=255)                                                             ccAGAAAACCAGAATGCCAAGCGTCACAATGCCAATGGCTGCATCACCCACATTAATGGTCGGCAgcgc  <  1:2712542/69‑1 (MQ=255)                                                             ccAGAAAACCAGAATGCCAAGCGTCACAATGCCAATGGCTGCATCACCCACATTAATGGTCGGCAgcgc  <  1:1884568/69‑1 (MQ=255)                                                             ccAGAAAACCAGAATGCCAAGCGTCACAATGCCAATGGCTGCATCACCCACATTAATGGTCGGCAgcgc  <  1:1875244/69‑1 (MQ=255)                                                             ccAGAAAACCAGAATGCCAAGCGTCACAATGCCAATGGCTGCATCACCCACATTAATGGTCGGCAgcgc  <  1:1867518/69‑1 (MQ=255)                                                             ccAGAAAACCAGAATGCCAAGCGTCACAATGCCAATGGCTGCATCACCCACATTAATGGTCGGCAgcgc  <  1:1725242/69‑1 (MQ=255)                                                             ccAGAAAACCAGAATGCCAAGCGTCACAATGCCAATGGCTGCATCACCCACATTAATGGTCGGCAgcgc  <  1:3341644/69‑1 (MQ=255)                                                             ccAGAAAACCAGAATGCCAAGCGTCACAATGCCAATGGCTGCATCACCCACATTAATGGTCGGCAgcgc  <  1:1534562/69‑1 (MQ=255)                                                             ccAGAAAACCAGAATGCCAAGCGTCACAATGCCAATGGCTGCATCACCCACATTAATGGTCGGCAgcgc  <  1:1342549/69‑1 (MQ=255)                                                             ccAGAAAACCAGAATGCCAAGCGTCACAATGCCAATGGCTGCATCACCCACATTAATGGTCGGCAgcgc  <  1:3222822/69‑1 (MQ=255)                                                                |                                                                  CACCGCGCAACCAGCCAGCAAGGCCGGAAGGTGACCGGGTAAACGAATGCCCAGACGCGGCCAAAAAACAAGAATACCTAGCGTCACAATGCCAATGGCAGCATCACCCACATTAATGGTCGGCAGCGC  >  minE/812361‑812489

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑