Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	740,595	1	.	T	34.4%	72.3 / 92.8	61	coding (212/2190 nt)	fhuE	ferric‑rhodotorulic acid outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base . (0/40);  new base T (21/0);  total (21/40)
Fisher's exact test for biased strand distribution p-value = 8.21e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TAACAATAGTGACCGACTGAGGAATATCACGTTGAGTCATCTGCA‑TTTTGGTACCCGCAGAGGTAGACGTTACGCTGT  >  minE/740551‑740628                                              |                                  tAACAATAGTGACCGACTGAGGAATATCACGTTGAGTCATCTGCA‑TTTTGGTACCCGCag                    <  1:1028102/60‑1 (MQ=255) tAACAATAGTGACCGACTGAGGAATATCACGTTGAGTCATCTGCA‑TTTTGGTACCCGCag                    <  1:954629/60‑1 (MQ=255) tAACAATAGTGACCGACTGAGGAATATCACGTTGAGTCATCTGCA‑TTTTGGTACCCGCag                    <  1:935594/60‑1 (MQ=255) tAACAATAGTGACCGACTGAGGAATATCACGTTGAGTCATCTGCA‑TTTTGGTACCCGCag                    <  1:827444/60‑1 (MQ=255) tAACAATAGTGACCGACTGAGGAATATCACGTTGAGTCATCTGCA‑TTTTGGTACCCGCag                    <  1:624711/60‑1 (MQ=255) tAACAATAGTGACCGACTGAGGAATATCACGTTGAGTCATCTGCA‑TTTTGGTACCCGCag                    <  1:59419/60‑1 (MQ=255) tAACAATAGTGACCGACTGAGGAATATCACGTTGAGTCATCTGCA‑TTTTGGTACCCGCag                    <  1:564044/60‑1 (MQ=255) tAACAATAGTGACCGACTGAGGAATATCACGTTGAGTCATCTGCA‑TTTTGGTACCCGCag                    <  1:557341/60‑1 (MQ=255) tAACAATAGTGACCGACTGAGGAATATCACGTTGAGTCATCTGCA‑TTTTGGTACCCGCag                    <  1:530781/60‑1 (MQ=255) tAACAATAGTGACCGACTGAGGAATATCACGTTGAGTCATCTGCA‑TTTTGGTACCCGCag                    <  1:526141/60‑1 (MQ=255) tAACAATAGTGACCGACTGAGGAATATCACGTTGAGTCATCTGCA‑TTTTGGTACCCGCag                    <  1:45157/60‑1 (MQ=255) tAACAATAGTGACCGACTGAGGAATATCACGTTGAGTCATCTGCA‑TTTTGGTACCCGCag                    <  1:439821/60‑1 (MQ=255) tAACAATAGTGACCGACTGAGGAATATCACGTTGAGTCATCTGCA‑TTTTGGTACCCGCag                    <  1:403126/60‑1 (MQ=255) tAACAATAGTGACCGACTGAGGAATATCACGTTGAGTCATCTGCA‑TTTTGGTACCCGCag                    <  1:327478/60‑1 (MQ=255) tAACAATAGTGACCGACTGAGGAATATCACGTTGAGTCATCTGCA‑TTTTGGTACCCGCag                    <  1:322532/60‑1 (MQ=255) tAACAATAGTGACCGACTGAGGAATATCACGTTGAGTCATCTGCA‑TTTTGGTACCCGCag                    <  1:318458/60‑1 (MQ=255) tAACAATAGTGACCGACTGAGGAATATCACGTTGAGTCATCTGCA‑TTTTGGTACCCGCag                    <  1:188878/60‑1 (MQ=255) tAACAATAGTGACCGACTGAGGAATATCACGTTGAGTCATCTGCA‑TTTTGGTACCCGCag                    <  1:172141/60‑1 (MQ=255) tAACAATAGTGACCGAATGAGGAATATCACGTTGAGTCATCTGCA‑TTTTGGTACCCGCag                    <  1:742770/60‑1 (MQ=255)                aCTGAGGAATATCACGTTGAGTCATCTGCA‑TTTTGGTACCCGCAGAGGTAGACGTTACGCTg   <  1:516680/62‑1 (MQ=255)                aCTGAGGAATATCACGTTGAGTCATCTGCA‑TTTTGGTACCCGCAGAGGTAGACGTTACGCTg   <  1:52220/62‑1 (MQ=255)                aCTGAGGAATATCACGTTGAGTCATCTGCA‑TTTTGGTACCCGCAGAGGTAGACGTTACGCTg   <  1:200297/62‑1 (MQ=255)                aCTGAGGAATATCACGTTGAGTCATCTGCA‑TTTTGGTACCCGCAGAGGTAGACGTTACGCTg   <  1:523949/62‑1 (MQ=255)                aCTGAGGAATATCACGTTGAGTCATCTGCA‑TTTTGGTACCCGCAGAGGTAGACGTTACGCTg   <  1:1163933/62‑1 (MQ=255)                aCTGAGGAATATCACGTTGAGTCATCTGCA‑TTTTGGTACCCGCAGAGGTAGACGTTACGCTg   <  1:1138966/62‑1 (MQ=255)                aCTGAGGAATATCACGTTGAGTCATCTGCA‑TTTTGGTACCCGCAGAGGTAGACGTTACGCTg   <  1:1080178/62‑1 (MQ=255)                aCTGAGGAATATCACGTTGAGTCATCTGCA‑TTTTGGTACCCGCAGAGGTAGACGTTACGCTg   <  1:606076/62‑1 (MQ=255)                aCTGAGGAATATCACGTTGAGTCATCTGCA‑TTTTGGTACCCGCAGAGGTAGACGTTACGCTg   <  1:1074614/62‑1 (MQ=255)                aCTGAGGAATATCACGTTGAGTCATCTGCA‑TTTTGGTACCCGCAGAGGTAGACGTTACGCTg   <  1:738646/62‑1 (MQ=255)                aCTGAGGAATATCACGTTGAGTCATCTGCA‑TTTTGGTACCCGCAGAGGTAGACGTTACGCTg   <  1:773230/62‑1 (MQ=255)                aCTGAGGAATATCACGTTGAGTCATCTGCA‑TTTTGGTACCCGCAGAGGTAGACGTTACGCTg   <  1:817043/62‑1 (MQ=255)                aCTGAGGAATATCACGTTGAGTCATCTGCA‑TTTTGGTACCCGCAGAGGTAGACGTTACGCTg   <  1:1050019/62‑1 (MQ=255)                aCTGAGGAATATCACGTTGAGTCATCTGCA‑TTTTGGTACCCGCAGAGGTAGACGTTACGCTg   <  1:1008816/62‑1 (MQ=255)                aCTGAGGAATATCACGTTGAGTCATCTGCA‑TTTTGGTACCCGCAGAGGTAGACGTTACGCTg   <  1:3031/62‑1 (MQ=255)                aCTGAGGAATATCACGTTGAGTCATCTGCA‑TTTTGGTACCCGCAGAGGTAGACGTTACGCTg   <  1:1198770/62‑1 (MQ=255)                aCTGAGGAATATCACGTTGAGTCATCTGCA‑TTTTGGTACCCGCAGAGGTAGACGTTACGCTg   <  1:416274/62‑1 (MQ=255)                aCTGAGGAATATCACGTTGAGTCATCTGCA‑TTTTGGTACCCGCAGAGGTAGACGTTACGCTg   <  1:407171/62‑1 (MQ=255)                aCTGAGGAATATCACGTTGAGTCATCTGCA‑TTTTGGTACCCGCAGAGGTAGACGTTACGCTg   <  1:1178343/62‑1 (MQ=255)                 cTGAGGAATATCACGTTGAGTCATCTGCA‑TTTTGGTACCCGCAGAGGTAGACGTTACGCTg   <  1:481833/61‑1 (MQ=255)                  tGAGGAATATCACGTTGAGTCATCTGCATTTTTGGTACCCGCTGAGGTAGACGTTACGCTGt  >  1:533078/1‑62 (MQ=255)                  tGAGGAATATCACGTTGAGTCATCTGCATTTTTGGTACCCGCAGAGGTAGACGtt         >  1:423837/1‑55 (MQ=255)                  tGAGGAATATCACGTTGAGTCATCTGCATTTTTGGTACCCGCAGAGGTAGACGTTACGCTg   >  1:42829/1‑61 (MQ=255)                  tGAGGAATATCACGTTGAGTCATCTGCATTTTTGGTACCCGCAGAGGTAGACGTTACGCTg   >  1:1130860/1‑61 (MQ=255)                  tGAGGAATATCACGTTGAGTCATCTGCATTTTTGGTACCCGCAGAGGTAGACGTTACGCTg   >  1:886812/1‑61 (MQ=255)                  tGAGGAATATCACGTTGAGTCATCTGCATTTTTGGTACCCGCAGAGGTAGACGTTACGCTg   >  1:874637/1‑61 (MQ=255)                  tGAGGAATATCACGTTGAGTCATCTGCATTTTTGGTACCCGCAGAGGTAGACGTTACGCTGt  >  1:296210/1‑62 (MQ=255)                  tGAGGAATATCACGTTGAGTCATCTGCATTTTTGGTACCCGCAGAGGTAGACGTTACGCTGt  >  1:356358/1‑62 (MQ=255)                  tGAGGAATATCACGTTGAGTCATCTGCATTTTTGGTACCCGCAGAGGTAGACGTTACGCTGt  >  1:106006/1‑62 (MQ=255)                  tGAGGAATATCACGTTGAGTCATCTGCATTTTTGGTACCCGCAGAGGTAGACGTTACGCTGt  >  1:844047/1‑62 (MQ=255)                  tGAGGAATATCACGTTGAGTCATCTGCATTTTTGGTACCCGCAGAGGTAGACGTTACGCTGt  >  1:870507/1‑62 (MQ=255)                  tGAGGAATATCACGTTGAGTCATCTGCATTTTTGGTACCCGCAGAGGTAGACGTTACGCTGt  >  1:295138/1‑62 (MQ=255)                  tGAGGAATATCACGTTGAGTCATCTGCATTTTTGGTACCCGCAGAGGTAGACGTTACGCTGt  >  1:230258/1‑62 (MQ=255)                  tGAGGAATATCACGTTGAGTCATCTGCATTTTTGGTACCCGCAGAGGTAGACGTTACGCTGt  >  1:790763/1‑62 (MQ=255)                  tGAGGAATATCACGTTGAGTCATCTGCATTTTTGGTACCCGCAGAGGTAGACGTTACGCTGt  >  1:299147/1‑62 (MQ=255)                  tGAGGAATATCACGTTGAGTCATCTGCATTTTTGGTACCCGCAGAGGTAGACGTTACGCTGt  >  1:17988/1‑62 (MQ=255)                  tGAGGAATATCACGTTGAGTCATCTGCATTTTTGGTACCCGCAGAGGTAGACGTTACGCTGt  >  1:381668/1‑62 (MQ=255)                  tGAGGAATATCACGTTGAGTCATCTGCATTTTTGGTACCCGCAGAGGTAGACGTTACGCTGt  >  1:144495/1‑62 (MQ=255)                  tGAGGAATATCACGTTGAGTCATCTGCATTTTTGGTACCCGCAGAGGTAGACGTTACGCTGt  >  1:1129417/1‑62 (MQ=255)                  tGAGGAATATCACGTTGAGTCATCTGCATTTTTGGTACCCGCAGAGGTAGACGTTACGCTGt  >  1:455218/1‑62 (MQ=255)                  tGAGGAATATCACGTTGAGTCATCTGCATTTTTGGTACCCGCAGAGGTAGACGTTACGCTGt  >  1:452687/1‑62 (MQ=255)                                  gAGTCATCTGCA‑TTTTGGTCCCCGCAGAGGTAGACGTTACGCTg   <  1:683336/44‑1 (MQ=255)                                              |                                  TAACAATAGTGACCGACTGAGGAATATCACGTTGAGTCATCTGCA‑TTTTGGTACCCGCAGAGGTAGACGTTACGCTGT  >  minE/740551‑740628

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