Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,984,916	0	A	.	33.3%	85.1 / 53.4	42	coding (656/1137 nt)	yggW	predicted oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (28/0);  new base . (14/0);  total (42/0)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GACGATGACGCGTTGTGGGATATATTCGAACAGGGGCATCAGTTATTAACCGCAGCGGGTTATCAGC  >  minE/1984897‑1984963                    |                                                gACGATGCCGCGTTGTGGG‑‑ATATTCGAACAGGGGCATCAGTTATTAACCGCAGCGGGTTATc     >  1:1097642/1‑62 (MQ=255) gACGATGACGCGTTGTGGG‑‑ATATTCGAACAGGGGCATCAg                           >  1:843672/1‑40 (MQ=37) gACGATGACGCGTTGTGGG‑‑ATATTCGAACAGGGGCATCAGTTATTAACCGCAGCGGGTTAt      >  1:3385450/1‑61 (MQ=255) gACGATGACGCGTTGTGGG‑‑ATATTCGAACAGGGGCATCAGTTATTAACCGCAGCGGGTTATc     >  1:3014688/1‑62 (MQ=255) gACGATGACGCGTTGTGGG‑‑ATATTCGAACAGGGGCATCAGTTATTAACCGCAGCGGGTTATc     >  1:1019642/1‑62 (MQ=255) gACGATGACGCGTTGTGGG‑‑ATATTCGAACAGGGGCATCAGTTATTAACCGCAGCGGGTTATc     >  1:2653551/1‑62 (MQ=255) gACGATGACGCGTTGTGGG‑‑ATATTCGAACAGGGGCATCAGTTATTAACCGCAGCGGGTTATc     >  1:2554570/1‑62 (MQ=255) gACGATGACGCGTTGTGGG‑‑ATATTCGAACAGGGGCATCAGTTATTAACCGCAGCGGGTTATc     >  1:2428798/1‑62 (MQ=255) gACGATGACGCGTTGTGGG‑‑ATATTCGAACAGGGGCATCAGTTATTAACCGCAGCGGGTTATc     >  1:2422238/1‑62 (MQ=255) gACGATGACGCGTTGTGGG‑‑ATATTCGAACAGGGGCATCAGTTATTAACCGCAGCGGGTTATc     >  1:3521363/1‑62 (MQ=255) gACGATGACGCGTTGTGGG‑‑ATATTCGAACAGGGGCATCAGTTATTAACCGCAGCGGGTTATc     >  1:3554186/1‑62 (MQ=255) gACGATGACGCGTTGTGGG‑‑ATATTCGAACAGGGGCATCAGTTATTAACCGCAGCGGGTTATc     >  1:1650694/1‑62 (MQ=255) gACGATGACGCGTTGTGGG‑‑ATATTCGAACAGGGGCATCAGTTATTAACCGCAGCGGGTTATc     >  1:56841/1‑62 (MQ=255) gACGATGACGCGTTGTGGG‑‑ATATTCGAACAGGGGCATCAGTTATTAACCGCAGCGGGTTATc     >  1:1250980/1‑62 (MQ=255)       gACGCGTTGTGGGATATATTCGAACAGGGGCATCAGTTATTACCCGCAGCGGGTTATCAGc  >  1:2436150/1‑61 (MQ=255)       gACGCGTTGTGGGATATATTCGAACAGGGGCATCAGTTATTAACCGc                >  1:1282902/1‑47 (MQ=255)       gACGCGTTGTGGGATATATTCGAACAGGGGCATCAGTTATTAACCGCa               >  1:563845/1‑48 (MQ=255)       gACGCGTTGTGGGATATATTCGAACAGGGGCATCAGTTATTAACCGCAGCGGGTTATCAg   >  1:274508/1‑60 (MQ=255)       gACGCGTTGTGGGATATATTCGAACAGGGGCATCAGTTATTAACCGCAGCGGGTTATCAGc  >  1:504834/1‑61 (MQ=255)       gACGCGTTGTGGGATATATTCGAACAGGGGCATCAGTTATTAACCGCAGCGGGTTATCAGc  >  1:345512/1‑61 (MQ=255)       gACGCGTTGTGGGATATATTCGAACAGGGGCATCAGTTATTAACCGCAGCGGGTTATCAGc  >  1:630610/1‑61 (MQ=255)       gACGCGTTGTGGGATATATTCGAACAGGGGCATCAGTTATTAACCGCAGCGGGTTATCAGc  >  1:695313/1‑61 (MQ=255)       gACGCGTTGTGGGATATATTCGAACAGGGGCATCAGTTATTAACCGCAGCGGGTTATCAGc  >  1:772433/1‑61 (MQ=255)       gACGCGTTGTGGGATATATTCGAACAGGGGCATCAGTTATTAACCGCAGCGGGTTATCAGc  >  1:3377910/1‑61 (MQ=255)       gACGCGTTGTGGGATATATTCGAACAGGGGCATCAGTTATTAACCGCAGCGGGTTATCAGc  >  1:3255518/1‑61 (MQ=255)       gACGCGTTGTGGGATATATTCGAACAGGGGCATCAGTTATTAACCGCAGCGGGTTATCAGc  >  1:3215720/1‑61 (MQ=255)       gACGCGTTGTGGGATATATTCGAACAGGGGCATCAGTTATTAACCGCAGCGGGTTATCAGc  >  1:2726825/1‑61 (MQ=255)       gACGCGTTGTGGGATATATTCGAACAGGGGCATCAGTTATTAACCGCAGCGGGTTATCAGc  >  1:3023661/1‑61 (MQ=255)       gACGCGTTGTGGGATATATTCGAACAGGGGCATCAGTTATTAACCGCAGCGGGTTATCAGc  >  1:2992583/1‑61 (MQ=255)       gACGCGTTGTGGGATATATTCGAACAGGGGCATCAGTTATTAACCGCAGCGGGTTATCAGc  >  1:283483/1‑61 (MQ=255)       gACGCGTTGTGGGATATATTCGAACAGGGGCATCAGTTATTAACCGCAGCGGGTTATCAGc  >  1:2780351/1‑61 (MQ=255)       gACGCGTTGTGGGATATATTCGAACAGGGGCATCAGTTATTAACCGCAGCGGGTTATCAGc  >  1:2577273/1‑61 (MQ=255)       gACGCGTTGTGGGATATATTCGAACAGGGGCATCAGTTATTAACCGCAGCGGGTTATCAGc  >  1:2272590/1‑61 (MQ=255)       gACGCGTTGTGGGATATATTCGAACAGGGGCATCAGTTATTAACCGCAGCGGGTTATCAGc  >  1:2206142/1‑61 (MQ=255)       gACGCGTTGTGGGATATATTCGAACAGGGGCATCAGTTATTAACCGCAGCGGGTTATCAGc  >  1:1867260/1‑61 (MQ=255)       gACGCGTTGTGGGATATATTCGAACAGGGGCATCAGTTATTAACCGCAGCGGGTTATCAGc  >  1:1835600/1‑61 (MQ=255)       gACGCGTTGTGGGATATATTCGAACAGGGGCATCAGTTATTAACCGCAGCGGGTTATCAGc  >  1:1498673/1‑61 (MQ=255)       gACGCGTTGTGGGATATATTCGAACAGGGGCATCAGTTATTAACCGCAGCGGGTTATCAGc  >  1:1423634/1‑61 (MQ=255)       gACGCGTTGTGGGATATATTCGAACAGGGGCATCAGTTATTAACCGCAGCGGGTTATCAGc  >  1:1406379/1‑61 (MQ=255)       gACGCGTTGTGGGATATATTCGAACAGGGGCATCAGTTATTAACCGCAGCGGGTTATCAGc  >  1:1229710/1‑61 (MQ=255)       gACGCGTTGTGGGATATATTCGAACAGGGGCATCAGTTATTAACCGCAGCGGGTTATCAGc  >  1:1101724/1‑61 (MQ=255)       gACGCGTTGTGGGATATATTCGAACAGGGGCATCAGTTATTAACCGCAGCGGGTTATCAGc  >  1:1073966/1‑61 (MQ=255)                    |                                                GACGATGACGCGTTGTGGGATATATTCGAACAGGGGCATCAGTTATTAACCGCAGCGGGTTATCAGC  >  minE/1984897‑1984963

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