Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	832,300	0	G	.	50.0%	15.4 / 18.9	12	coding (758/3744 nt)	narG	nitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (3/3);  new base . (6/0);  total (9/3)
Fisher's exact test for biased strand distribution p-value = 1.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGT  >  minE/832237‑832358                                                                |                                                           ggggCGAGCAAACCGACGTACCGGAATCTGCCGACTGGTATAACTCCAGCTACATCATCGCCTGGGGGTc                                                      >  1:41085/1‑70 (MQ=255) ggggCGAGCAAACCGACGTACCGGAATCTGCCGACTGGTATAACTCCAGCTACATCATCGCCTGGGGGTc                                                      >  1:68459/1‑70 (MQ=255)               gACGTACCGGAATCTGCCGACTGGTATAACTCCAGCTACATCATCGCCTGGGGGTCTAACGTACCGCAGAc                                       <  1:52815/71‑1 (MQ=255)               gACGTACCGGAATCTGCCGACTGGTATAACTCCAGCTACATCATCGCCTGGGGGGCTAACGTACCGCAGAc                                       <  1:114564/71‑1 (MQ=255)                               ctgACTGGTACAACTCCAGCTACATCATCGCCT‑GGGGTCAAACGTGCCGCAGa                                        >  1:4488/1‑53 (MQ=255)                                       tACAACTCCAGCTACATCATCGCCT‑GGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCActtctt               >  1:168879/1‑70 (MQ=255)                                       tACAACTCCAGCTACATCATCGCCT‑GGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCActtctt               >  1:241712/1‑70 (MQ=255)                                       tACAACTCCAGCTACATCATCGCCT‑GGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCActtctt               >  1:82368/1‑70 (MQ=255)                                              ccAGCTACATCATCGCCT‑GGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACGTCTTTACTGAAg       >  1:254677/1‑71 (MQ=255)                                                   tACATCATCGCCT‑GGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGt  >  1:272084/1‑71 (MQ=255)                                                       tcatcGCCTGGGGCTCCAACGTGCCGCAGACGCGGACCCCGGACGCGCACTTCTTCACCGagg       <  1:156111/63‑3 (MQ=255)                                                         atcGCCTGGGGCTCCAACGTGCCGCAGACGCGGACCCCGGACGCGCACTTCTTCACCGAGGTg     >  1:131092/1‑63 (MQ=255)                                                                |                                                           GGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGT  >  minE/832237‑832358

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