Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,604,413	A→G	53.9%	L86L (CTT→CTC)	yhiP ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,604,413	0	A	G	53.9%	2.5 / 29.1	39	L86L (CTT→CTC)	yhiP	predicted transporter
Reads supporting (aligned to +/- strand):  ref base A (9/9);  new base G (18/3);  total (27/12)
Fisher's exact test for biased strand distribution p-value = 3.46e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

CTTAAGTAGCGACATGCCGGTCATGAAGTAGCCAATCGCCAGCACAAGTGCTCCAAGAACAATGGTGCGTTTGGTCCCCAGCAGGTGGTCGCCGACATAGCCGCCAATGGAAATGAGGCCA  >  minE/2604359‑2604479                                                       |                                                                   cTTAAGTAGCGACATGCCGGTCATGAAGTAGCCAATCGCCAGCACAAGTGCTCCAAGAACAATGGTGCg                                                      >  1:3016541/1‑69 (MQ=255) cTTAAGTAGCGACATGCCGGTCATGAAGTAGCCAATCGCCAGCACAAGTGCTCCAAGAACAATGGTGCg                                                      >  1:71048/1‑69 (MQ=255) cTTAAGTAGCGACATGCCGGTCATGAAGTAGCCAATCGCCAGCACAAGTGCTCCAAGAACAATGGTGCg                                                      >  1:2951141/1‑69 (MQ=255) cTTAAGTAGCGACATGCCGGTCATGAAGTAGCCAATCGCCAGCACAAGTGCTCCAAGAACAATGGTGCg                                                      >  1:3157745/1‑69 (MQ=255) cTTAAGTAGCGACATGCCGGTCATGAAGTAGCCAATCGCCAGCACAAGTGCTCCAAGAACAATGGTGCg                                                      >  1:1665514/1‑69 (MQ=255) cTTAAGTAGCGACATGCCGGTCATGAAGTAGCCAATCGCCAGCACAAGTGCTCCAAGAACAATGGTGCg                                                      >  1:2876843/1‑69 (MQ=255) cTTAAGTAGCGACATGCCGGTCATGAAGTAGCCAATCGCCAGCACAAGTGCTCCAAGAACAATGGTGCg                                                      >  1:3230798/1‑69 (MQ=255)                          aaGTAGCCAATCGCCAGCACAAGTGCTCCAAGAACAATGGTGCGTTTGGTCCCCAGCAGGTGGTCGcc                              >  1:1720247/1‑68 (MQ=255)                          aaGTAGCCAATCGCCAGCACAAGTGCTCCAAGAACAATGGTGCGTTTGGTCCCCAGCAGGTGGTCGcc                              >  1:2522217/1‑68 (MQ=255)                           aGTAGCCAATCGCCAGCACAAGTGCTCCGAGAACAATGGTGCGTTTGGTTCCCAGCAGGTGGTCGccaa                            >  1:3209699/1‑67 (MQ=39)                           aGTAGCCAATCGCCAGCACAAGTGCTCCGAGAACAATGGTGCGTTTGGTTCCCAGCAGGTGGTCGccaa                            >  1:3054706/1‑67 (MQ=39)                           aGTAGCCAATCGCCAGCACAAGTGCTCCGAGAACAATGGTGCGTTTGGTTCCCAGCAGGTGGTCGccaa                            >  1:2997594/1‑67 (MQ=39)                           aGTAGCCAATCGCCAGCACAAGTGCTCCGAGAACAATGGTGCGTTTGGTTCCCAGCAGGTGGTCGccaa                            >  1:2853703/1‑67 (MQ=39)                           aGTAGCCAATCGCCAGCACAAGTGCTCCGAGAACAATGGTGCGTTTGGTTCCCAGCAGGTGGTCGccaa                            >  1:28181/1‑67 (MQ=39)                           aGTAGCCAATCGCCAGCACAAGTGCTCCGAGAACAATGGTGCGTTTGGTTCCCAGCAGGTGGTCGccaa                            >  1:2726536/1‑67 (MQ=39)                           aGTAGCCAATCGCCAGCACAAGTGCTCCGAGAACAATGGTGCGTTTGGTTCCCAGCAGGTGGTCGccaa                            >  1:2676177/1‑67 (MQ=39)                           aGTAGCCAATCGCCAGCACAAGTGCTCCGAGAACAATGGTGCGTTTGGTTCCCAGCAGGTGGTCGccaa                            >  1:2407054/1‑67 (MQ=39)                           aGTAGCCAATCGCCAGCACAAGTGCTCCGAGAACAATGGTGCGTTTGGTTCCCAGCAGGTGGTCGccaa                            >  1:2063789/1‑67 (MQ=39)                           aGTAGCCAATCGCCAGCACAAGTGCTCCGAGAACAATGGTGCGTTTGGTTCCCAGCAGGTGGTCGccaa                            >  1:1852276/1‑67 (MQ=39)                           aGTAGCCAATCGCCAGCACAAGTGCTCCGAGAACAATGGTGCGTTTGGTTCCCAGCAGGTGGTCGccaa                            >  1:3380640/1‑67 (MQ=39)                           aGTAGCCAATCGCCAGCACAAGTGCTCCGAGAACAATGGTGCGTTTGGTTCCCAGCAGGTGGTCGccaa                            >  1:1676848/1‑67 (MQ=39)                           aGTAGCCAATCGCCAGCACAAGTGCTCCGAGAACAATGGTGCGTTTGGTTCCCAGCAGGTGGTCGccaa                            >  1:571187/1‑67 (MQ=39)                           aGTAGCCAATCGCCAGCACAAGTGCTCCGAGAACAATGGTGCGTTTGGTTCCCAGCAGGTGGTCGccaa                            >  1:1288784/1‑67 (MQ=39)                           aGTAGCCAATCGCCAGCACAAGTGCTCCGAGAACAATGGTGCGTTTGGTTCCCAGCAGGTGGTCGccaa                            >  1:788328/1‑67 (MQ=39)                           aGTAGCCAATCGCCAGCACAAGTGCTCCGAGAACAATGGTGCGTTTGGTTCCCAGCAGGTGGTCGccaa                            >  1:842540/1‑67 (MQ=39)                           aGTAGCCAATCGCCAGCACAAGTGCTCCGAGAACAATGGTGCGTTTGGTTCCCAGCAGGTGGTCGcca                             >  1:2431032/1‑67 (MQ=39)                           aGTAGCCAATCGCCAGCACAAGTGCTCCGAGAACAATGGTGCGTTTGGTTCCCAGCAGGTGGTCGcca                             >  1:54209/1‑67 (MQ=39)                            gTAGCCAATCGCCAGCACAAGTGCTCCAAGAACAATGGTGCGTTTGGTCCCCAGCAGGTGGTCGCCGa                            <  1:846437/68‑1 (MQ=255)                            gTAGCCAATCGCCAGCACAAGTGCTCCAAGAACAATGGTGCGTTTGGTCCCCAGCAGGTGGTCGCCGa                            <  1:2919968/68‑1 (MQ=255)                                             caAGTGCTCCAAGAACAATGGTGCGTTTGGTCCCCAGCAGGTGGTCGCCGACATAGCCGCCAATGGa            <  1:2094379/67‑1 (MQ=255)                                             caAGTGCTCCAAGAACAATGGTGCGTTTGGTCCCCAGCAGGTGGTCGCCGACATAGCCGCCAATGGa            <  1:1269550/67‑1 (MQ=255)                                                gTGCTCCGAGAACAATGGTGCGTTTGGTTCCcagc                                         <  1:1715051/35‑1 (MQ=37)                                                gTGCTCCGAGAACAATGGTGCGTTTGGTTCCcagc                                         <  1:1271455/35‑1 (MQ=37)                                                gTGCTCCGAGAACAATGGTGCGTTTGGTTCCcagc                                         <  1:746818/35‑1 (MQ=37)                                                gTGCTCCAAGAACAATGGTGCGTTTGGTCCCCAGCAGGTGGTCGCCGACATAGCCGCCAATGGAAAt         <  1:221070/67‑1 (MQ=255)                                                     ccAAGAACAATGGTGCGTTTGGTCCCCAGCAGGTGGTCGCCGACATAGCCGCCAATGGAAATGAGGCCa  <  1:3371670/69‑1 (MQ=255)                                                     ccAAGAACAATGGTGCGTTTGGTCCCCAGCAGGTGGTCGCCGACATAGCCGCCAATGGAAATGAGGCCa  <  1:1428053/69‑1 (MQ=255)                                                      cAAGAACAATGGTGCGTTTGGTCCCCAGCAGGTGGTCGCCGACATAGCCGCCAATGGa            <  1:1252353/58‑1 (MQ=255)                                                      cAAGAACAATGGTGCGTTTGGTCCCCAGCAGGTGGTCGCCGACATAGCCGCCAATGGAAATGAGGCCa  <  1:1059630/68‑1 (MQ=255)                                                       aaGAACAATGGTGCGTTTGGTCCCCAGCAGGTGGTCGCCGACATAGCCGCCAATGGa            <  1:1548630/57‑1 (MQ=255)                                                       |                                                                   CTTAAGTAGCGACATGCCGGTCATGAAGTAGCCAATCGCCAGCACAAGTGCTCCAAGAACAATGGTGCGTTTGGTCCCCAGCAGGTGGTCGCCGACATAGCCGCCAATGGAAATGAGGCCA  >  minE/2604359‑2604479

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