Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,309,231	0	G	.	59.7%	‑0.0 / 151.1	67	coding (666/960 nt)	ftsN	essential cell division protein
Reads supporting (aligned to +/- strand):  ref base G (27/0);  new base . (0/40);  total (27/40)
Fisher's exact test for biased strand distribution p-value = 2.44e-19
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CAGCCGTACCAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAAGCCGCAGCAAGCTGCGCC  >  minE/2309178‑2309248                                                      |                  cagcCGTACCAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAA‑CCgcagc            <  1:2674358/60‑1 (MQ=255) cagcCGTACCAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAA‑CCgcagc            <  1:1936900/60‑1 (MQ=255) cagcCGTACCAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAA‑CCgcagc            <  1:2176143/60‑1 (MQ=255) cagcCGTACCAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAA‑CCgcagc            <  1:2184689/60‑1 (MQ=255) cagcCGTACCAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAA‑CCgcagc            <  1:2202375/60‑1 (MQ=255) cagcCGTACCAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAA‑CCgcagc            <  1:45033/60‑1 (MQ=255) cagcCGTACCAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAA‑CCgcagc            <  1:2242341/60‑1 (MQ=255) cagcCGTACCAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAA‑CCgcagc            <  1:2452574/60‑1 (MQ=255) cagcCGTACCAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAA‑CCgcagc            <  1:253634/60‑1 (MQ=255) cagcCGTACCAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAA‑CCgcagc            <  1:2603690/60‑1 (MQ=255) cagcCGTACCAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAA‑CCgcagc            <  1:107226/60‑1 (MQ=255) cagcCGTACCAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAA‑CCgcagc            <  1:469786/60‑1 (MQ=255) cagcCGTACCAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAA‑CCgcagc            <  1:2687244/60‑1 (MQ=255) cagcCGTACCAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAA‑CCgcagc            <  1:3628640/60‑1 (MQ=255) cagcCGTACCAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAA‑CCgcagc            <  1:2782977/60‑1 (MQ=255) cagcCGTACCAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAA‑CCgcagc            <  1:281250/60‑1 (MQ=255) cagcCGTACCAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAA‑CCgcagc            <  1:2851714/60‑1 (MQ=255) cagcCGTACCAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAA‑CCgcagc            <  1:2888481/60‑1 (MQ=255) cagcCGTACCAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAA‑CCgcagc            <  1:2981151/60‑1 (MQ=255) cagcCGTACCAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAA‑CCgcagc            <  1:3337121/60‑1 (MQ=255) cagcCGTACCAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAA‑CCgcagc            <  1:3436984/60‑1 (MQ=255) cagcCGTACCAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAA‑CCgcagc            <  1:851154/60‑1 (MQ=255) cagcCGTACCAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAA‑CCgcagc            <  1:1869279/60‑1 (MQ=255) cagcCGTACCAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAA‑CCgcagc            <  1:1154695/60‑1 (MQ=255) cagcCGTACCAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAA‑CCgcagc            <  1:1857195/60‑1 (MQ=255) cagcCGTACCAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAA‑CCgcagc            <  1:523791/60‑1 (MQ=255) cagcCGTACCAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAA‑CCgcagc            <  1:1754932/60‑1 (MQ=255) cagcCGTACCAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAA‑CCgcagc            <  1:166385/60‑1 (MQ=255) cagcCGTACCAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAA‑CCgcagc            <  1:539151/60‑1 (MQ=255) cagcCGTACCAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAA‑CCgcagc            <  1:789610/60‑1 (MQ=255) cagcCGTACCAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAA‑CCgcagc            <  1:1550478/60‑1 (MQ=255) cagcCGTACCAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAA‑CCgcagc            <  1:120383/60‑1 (MQ=255) cagcCGTACCAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAA‑CCgcagc            <  1:596428/60‑1 (MQ=255) cagcCGTACCAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAA‑CCgcagc            <  1:625028/60‑1 (MQ=255)  agcCGTACCAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAA‑CCgcagc            <  1:3499141/59‑1 (MQ=255)    ccGTACCAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAA‑CCgcagc            <  1:1268440/57‑1 (MQ=255)         ccAGGAGCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAA‑CCgcagc            <  1:488210/52‑1 (MQ=255)          cAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAAGCCgc               >  1:1617832/1‑49 (MQ=255)          cAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAAGCCgc               >  1:1140878/1‑49 (MQ=255)          cAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAAGCCGCAGCAAGCTgcg    >  1:1900583/1‑60 (MQ=255)          cAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAAGCCGCAGCAAGCTGCGcc  >  1:3570331/1‑62 (MQ=255)          cAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAAGCCGCAGCAAGCTGCGcc  >  1:586021/1‑62 (MQ=255)          cAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAAGCCGCAGCAAGCTGCGcc  >  1:441505/1‑62 (MQ=255)          cAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAAGCCGCAGCAAGCTGCGcc  >  1:2639591/1‑62 (MQ=255)          cAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAAGCCGCAGCAAGCTGCGcc  >  1:3487363/1‑62 (MQ=255)          cAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAAGCCGCAGCAAGCTGCGcc  >  1:3445830/1‑62 (MQ=255)          cAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAAGCCGCAGCAAGCTGCGcc  >  1:854667/1‑62 (MQ=255)          cAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAAGCCGCAGCAAGCTGCGcc  >  1:866780/1‑62 (MQ=255)          cAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAAGCCGCAGCAAGCTGCGcc  >  1:3230292/1‑62 (MQ=255)          cAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAAGCCGCAGCAAGCTGCGcc  >  1:114431/1‑62 (MQ=255)          cAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAAGCCGCAGCAAGCTGCGcc  >  1:1172333/1‑62 (MQ=255)          cAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAAGCCGCAGCAAGCTGCGcc  >  1:1209068/1‑62 (MQ=255)          cAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAAGCCGCAGCAAGCTGCGcc  >  1:1413110/1‑62 (MQ=255)          cAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAAGCCGCAGCAAGCTGCGcc  >  1:1540562/1‑62 (MQ=255)          cAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAAGCCGCAGCAAGCTGCGcc  >  1:161387/1‑62 (MQ=255)          cAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAAGCCGCAGCAAGCTGCGcc  >  1:2019494/1‑62 (MQ=255)          cAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAAGCCGCAGCAAGCTGCGcc  >  1:2183604/1‑62 (MQ=255)          cAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAAGCCGCAGCAAGCTGCGcc  >  1:220639/1‑62 (MQ=255)          cAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAAGCCGCAGCAAGCTGCGcc  >  1:2542700/1‑62 (MQ=255)          cAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAAGCCGCAGCAAGCTGCGcc  >  1:2704620/1‑62 (MQ=255)          cAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAAGCCGCAGCAAGCTGCGcc  >  1:2725542/1‑62 (MQ=255)          cAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAAGCCGCAGCAAGCTGCGcc  >  1:2815578/1‑62 (MQ=255)          cAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAAGCCGCAGCAAGCTGCGcc  >  1:3057828/1‑62 (MQ=255)          cAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAAGCCGCAGCAAGCTGCGcc  >  1:3394850/1‑62 (MQ=255)               tCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAA‑CCgcagc            <  1:1124177/46‑1 (MQ=255)                 tgctgcAAACTCCTGCGCACACGACTGCGCAGTCGAA‑CCgcagc            <  1:1824780/44‑1 (MQ=39)                   ctgcAAACTCCTGCGCACACGACTGCGCAGTCGAA‑CCgcagc            <  1:1111732/42‑1 (MQ=39)                                                      |                  CAGCCGTACCAGGATCTGCTGCAAACTCCTGCGCACACGACTGCGCAGTCGAAGCCGCAGCAAGCTGCGCC  >  minE/2309178‑2309248

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