Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,604,434	C→T	55.3%	G79G (GGG→GGA)	yhiP ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,604,434	0	C	T	55.3%	3.3 / 27.4	38	G79G (GGG→GGA)	yhiP	predicted transporter
Reads supporting (aligned to +/- strand):  ref base C (2/15);  new base T (18/3);  total (20/18)
Fisher's exact test for biased strand distribution p-value = 8.56e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

AAGTAGCCAATCGCCAGCACAAGTGCTCCAAGAACAATGGTGCGTTTGGTCCCCAGCAGGTGGTCGCCGACATAGCCGCCAATGGAAATGAGGCCATAG  >  minE/2604384‑2604482                                                   |                                                 aaGTAGCCAATCGCCAGCACAAGTGCTCCAAGAACAATGGTGCGTTTGGTCCCCAGCAGGTGGTCGcc                                 >  1:2522217/1‑68 (MQ=255) aaGTAGCCAATCGCCAGCACAAGTGCTCCAAGAACAATGGTGCGTTTGGTCCCCAGCAGGTGGTCGcc                                 >  1:1720247/1‑68 (MQ=255)  aGTAGCCAATCGCCAGCACAAGTGCTCCGAGAACAATGGTGCGTTTGGTTCCCAGCAGGTGGTCGccaa                               >  1:842540/1‑67 (MQ=39)  aGTAGCCAATCGCCAGCACAAGTGCTCCGAGAACAATGGTGCGTTTGGTTCCCAGCAGGTGGTCGccaa                               >  1:2676177/1‑67 (MQ=39)  aGTAGCCAATCGCCAGCACAAGTGCTCCGAGAACAATGGTGCGTTTGGTTCCCAGCAGGTGGTCGccaa                               >  1:2407054/1‑67 (MQ=39)  aGTAGCCAATCGCCAGCACAAGTGCTCCGAGAACAATGGTGCGTTTGGTTCCCAGCAGGTGGTCGccaa                               >  1:2726536/1‑67 (MQ=39)  aGTAGCCAATCGCCAGCACAAGTGCTCCGAGAACAATGGTGCGTTTGGTTCCCAGCAGGTGGTCGccaa                               >  1:28181/1‑67 (MQ=39)  aGTAGCCAATCGCCAGCACAAGTGCTCCGAGAACAATGGTGCGTTTGGTTCCCAGCAGGTGGTCGccaa                               >  1:2063789/1‑67 (MQ=39)  aGTAGCCAATCGCCAGCACAAGTGCTCCGAGAACAATGGTGCGTTTGGTTCCCAGCAGGTGGTCGccaa                               >  1:2853703/1‑67 (MQ=39)  aGTAGCCAATCGCCAGCACAAGTGCTCCGAGAACAATGGTGCGTTTGGTTCCCAGCAGGTGGTCGccaa                               >  1:1852276/1‑67 (MQ=39)  aGTAGCCAATCGCCAGCACAAGTGCTCCGAGAACAATGGTGCGTTTGGTTCCCAGCAGGTGGTCGccaa                               >  1:2997594/1‑67 (MQ=39)  aGTAGCCAATCGCCAGCACAAGTGCTCCGAGAACAATGGTGCGTTTGGTTCCCAGCAGGTGGTCGccaa                               >  1:1676848/1‑67 (MQ=39)  aGTAGCCAATCGCCAGCACAAGTGCTCCGAGAACAATGGTGCGTTTGGTTCCCAGCAGGTGGTCGccaa                               >  1:3054706/1‑67 (MQ=39)  aGTAGCCAATCGCCAGCACAAGTGCTCCGAGAACAATGGTGCGTTTGGTTCCCAGCAGGTGGTCGccaa                               >  1:3209699/1‑67 (MQ=39)  aGTAGCCAATCGCCAGCACAAGTGCTCCGAGAACAATGGTGCGTTTGGTTCCCAGCAGGTGGTCGccaa                               >  1:3380640/1‑67 (MQ=39)  aGTAGCCAATCGCCAGCACAAGTGCTCCGAGAACAATGGTGCGTTTGGTTCCCAGCAGGTGGTCGccaa                               >  1:1288784/1‑67 (MQ=39)  aGTAGCCAATCGCCAGCACAAGTGCTCCGAGAACAATGGTGCGTTTGGTTCCCAGCAGGTGGTCGccaa                               >  1:571187/1‑67 (MQ=39)  aGTAGCCAATCGCCAGCACAAGTGCTCCGAGAACAATGGTGCGTTTGGTTCCCAGCAGGTGGTCGccaa                               >  1:788328/1‑67 (MQ=39)  aGTAGCCAATCGCCAGCACAAGTGCTCCGAGAACAATGGTGCGTTTGGTTCCCAGCAGGTGGTCGcca                                >  1:54209/1‑67 (MQ=39)  aGTAGCCAATCGCCAGCACAAGTGCTCCGAGAACAATGGTGCGTTTGGTTCCCAGCAGGTGGTCGcca                                >  1:2431032/1‑67 (MQ=39)   gTAGCCAATCGCCAGCACAAGTGCTCCAAGAACAATGGTGCGTTTGGTCCCCAGCAGGTGGTCGCCGa                               <  1:2919968/68‑1 (MQ=255)   gTAGCCAATCGCCAGCACAAGTGCTCCAAGAACAATGGTGCGTTTGGTCCCCAGCAGGTGGTCGCCGa                               <  1:846437/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)                               aGAACAATGGTGCGTTTGGTCCCCAGCAGGTGGTCGCCGACATAGCCGCCAATGGAAATGAGGCCATAg  <  1:1880374/69‑1 (MQ=255)                               aGAACAATGGTGCGTTTGGTCCCCAGCAGGTGGTCGCCGACATAGCCGCCAATGGAAATGAGGCCATAg  <  1:1284370/69‑1 (MQ=255)                               aGAACAATGGTGCGTTTGGTCCCCAGCAGGTGGTCGCCGACATAGCCGCCAATGGAAATGAGGCCATAg  <  1:743391/69‑1 (MQ=255)                                     aTGGTGCGTTTGGTCCCCAGCAGGTGGTCGCCGACATAGCCGCCAATGGAAAt            <  1:1596495/53‑1 (MQ=255)                                      tGGTGCGTTTGGTCCCCAGCAGGTGGTCGCCGACATAGCCGCCAATGGa               <  1:1719481/49‑1 (MQ=255)                                                   |                                                 AAGTAGCCAATCGCCAGCACAAGTGCTCCAAGAACAATGGTGCGTTTGGTCCCCAGCAGGTGGTCGCCGACATAGCCGCCAATGGAAATGAGGCCATAG  >  minE/2604384‑2604482

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