Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,513,727	1	.	A	43.5%	4.9 / 95.8	46	coding (152/570 nt)	ubiX	3‑octaprenyl‑4‑hydroxybenzoate carboxy‑lyase
Reads supporting (aligned to +/- strand):  ref base . (26/0);  new base A (20/0);  total (46/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 = 9.67e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GTGCGTGACATCGGCTAATGCCTGCACTTCGCGCAGAG‑AAAAATCCGTTTCGAGGGATAAGGTCTGGCGCGCTGC  >  minE/1513690‑1513764                                       |                                      gtgcgtgACATCGGCTAATGCCTGCACTTCGCGCGGAG‑AAAAATCCGTTTCGAGGGATAAgg               >  1:431381/1‑62 (MQ=255) gtgcgtgACATCGGCTAATGCCTGCACTTCGCGCAGAG‑AAAAATCCGTTTCGAGGGAt                   >  1:1028224/1‑58 (MQ=255) gtgcgtgACATCGGCTAATGCCTGCACTTCGCGCAGAG‑AAAAATCCGTTTCGAGGGATAAgg               >  1:16344/1‑62 (MQ=255) gtgcgtgACATCGGCTAATGCCTGCACTTCGCGCAGAG‑AAAAATCCGTTTCGAGGGATAAgg               >  1:510004/1‑62 (MQ=255) gtgcgtgACATCGGCTAATGCCTGCACTTCGCGCAGAG‑AAAAATCCGTTTCGAGGGATAAgg               >  1:586735/1‑62 (MQ=255) gtgcgtgACATCGGCTAATGCCTGCACTTCGCGCAGAG‑AAAAATCCGTTTCGAGGGATAAgg               >  1:606996/1‑62 (MQ=255) gtgcgtgACATCGGCTAATGCCTGCACTTCGCGCAGAG‑AAAAATCCGTTTCGAGGGATAAgg               >  1:621751/1‑62 (MQ=255) gtgcgtgACATCGGCTAATGCCTGCACTTCGCGCAGAG‑AAAAATCCGTTTCGAGGGATAAgg               >  1:377333/1‑62 (MQ=255) gtgcgtgACATCGGCTAATGCCTGCACTTCGCGCAGAG‑AAAAATCCGTTTCGAGGGATAAgg               >  1:351495/1‑62 (MQ=255) gtgcgtgACATCGGCTAATGCCTGCACTTCGCGCAGAG‑AAAAATCCGTTTCGAGGGATAAgg               >  1:624973/1‑62 (MQ=255) gtgcgtgACATCGGCTAATGCCTGCACTTCGCGCAGAG‑AAAAATCCGTTTCGAGGGATAAgg               >  1:251392/1‑62 (MQ=255) gtgcgtgACATCGGCTAATGCCTGCACTTCGCGCAGAG‑AAAAATCCGTTTCGAGGGATAAgg               >  1:233143/1‑62 (MQ=255) gtgcgtgACATCGGCTAATGCCTGCACTTCGCGCAGAG‑AAAAATCCGTTTCGAGGGATAAgg               >  1:689229/1‑62 (MQ=255) gtgcgtgACATCGGCTAATGCCTGCACTTCGCGCAGAG‑AAAAATCCGTTTCGAGGGATAAgg               >  1:17231/1‑62 (MQ=255) gtgcgtgACATCGGCTAATGCCTGCACTTCGCGCAGAG‑AAAAATCCGTTTCGAGGGATAAgg               >  1:1195888/1‑62 (MQ=255) gtgcgtgACATCGGCTAATGCCTGCACTTCGCGCAGAG‑AAAAATCCGTTTCGAGGGATAAgg               >  1:1180763/1‑62 (MQ=255) gtgcgtgACATCGGCTAATGCCTGCACTTCGCGCAGAG‑AAAAATCCGTTTCGAGGGATAAgg               >  1:117049/1‑62 (MQ=255) gtgcgtgACATCGGCTAATGCCTGCACTTCGCGCAGAG‑AAAAATCCGTTTCGAGGGATAAgg               >  1:114355/1‑62 (MQ=255) gtgcgtgACATCGGCTAATGCCTGCACTTCGCGCAGAG‑AAAAATCCGTTTCGAGGGATAAgg               >  1:721738/1‑62 (MQ=255) gtgcgtgACATCGGCTAATGCCTGCACTTCGCGCAGAG‑AAAAATCCGTTTCGAGGGATAAgg               >  1:1136914/1‑62 (MQ=255) gtgcgtgACATCGGCTAATGCCTGCACTTCGCGCAGAG‑AAAAATCCGTTTCGAGGGATAAgg               >  1:72876/1‑62 (MQ=255) gtgcgtgACATCGGCTAATGCCTGCACTTCGCGCAGAG‑AAAAATCCGTTTCGAGGGATAAgg               >  1:793705/1‑62 (MQ=255) gtgcgtgACATCGGCTAATGCCTGCACTTCGCGCAGAG‑AAAAATCCGTTTCGAGGGATAAgg               >  1:1037606/1‑62 (MQ=255) gtgcgtgACATCGGCTAATGCCTGCACTTCGCGCAGAG‑AAAAATCCGTTTCGAGGGATAAgg               >  1:1031677/1‑62 (MQ=255) gtgcgtgACATCGGCTAATGCCTGCACTTCGCGCAGAG‑AAAAATCCGTTTCGAGGGATAAgg               >  1:820891/1‑62 (MQ=255) gtgcgtgACATCGGCTAATGCCAGCACTTCGCGCAGAG‑AAAAATCCGTTTCGAGGGATAAg                >  1:968902/1‑61 (MQ=255)               cTAATGCCTGCACTTCGCGCAGAGAAAAAATCCGTTTCGAGGGATAAGGTCTGTCGCGCTGc  >  1:648378/1‑62 (MQ=255)               cTAATGCCTGCACTTCGCGCAGAGAAAAAATCCGTTTCGAGGGATAAGGTCTGGCGCGCTGc  >  1:228703/1‑62 (MQ=255)               cTAATGCCTGCACTTCGCGCAGAGAAAAAATCCGTTTCGAGGGATAAGGTCTGGCGCGCTGc  >  1:1016557/1‑62 (MQ=255)               cTAATGCCTGCACTTCGCGCAGAGAAAAAATCCGTTTCGAGGGATAAGGTCTGGCGCGCTGc  >  1:839976/1‑62 (MQ=255)               cTAATGCCTGCACTTCGCGCAGAGAAAAAATCCGTTTCGAGGGATAAGGTCTGGCGCGCTGc  >  1:10272/1‑62 (MQ=255)               cTAATGCCTGCACTTCGCGCAGAGAAAAAATCCGTTTCGAGGGATAAGGTCTGGCGCGCTGc  >  1:817655/1‑62 (MQ=255)               cTAATGCCTGCACTTCGCGCAGAGAAAAAATCCGTTTCGAGGGATAAGGTCTGGCGCGCTGc  >  1:1052976/1‑62 (MQ=255)               cTAATGCCTGCACTTCGCGCAGAGAAAAAATCCGTTTCGAGGGATAAGGTCTGGCGCGCTGc  >  1:793279/1‑62 (MQ=255)               cTAATGCCTGCACTTCGCGCAGAGAAAAAATCCGTTTCGAGGGATAAGGTCTGGCGCGCTGc  >  1:740353/1‑62 (MQ=255)               cTAATGCCTGCACTTCGCGCAGAGAAAAAATCCGTTTCGAGGGATAAGGTCTGGCGCGCTGc  >  1:1116569/1‑62 (MQ=255)               cTAATGCCTGCACTTCGCGCAGAGAAAAAATCCGTTTCGAGGGATAAGGTCTGGCGCGCTGc  >  1:1140921/1‑62 (MQ=255)               cTAATGCCTGCACTTCGCGCAGAGAAAAAATCCGTTTCGAGGGATAAGGTCTGGCGCGCTGc  >  1:494126/1‑62 (MQ=255)               cTAATGCCTGCACTTCGCGCAGAGAAAAAATCCGTTTCGAGGGATAAGGTCTGGCGCGCTGc  >  1:675255/1‑62 (MQ=255)               cTAATGCCTGCACTTCGCGCAGAGAAAAAATCCGTTTCGAGGGATAAGGTCTGGCGCGCTGc  >  1:289945/1‑62 (MQ=255)               cTAATGCCTGCACTTCGCGCAGAGAAAAAATCCGTTTCGAGGGATAAGGTCTGGCGCGCTGc  >  1:390988/1‑62 (MQ=255)               cTAATGCCTGCACTTCGCGCAGAGAAAAAATCCGTTTCGAGGGATAAGGTCTGGCGCGCTGc  >  1:429764/1‑62 (MQ=255)               cTAATGCCTGCACTTCGCGCAGAGAAAAAATCCGTTTCGAGGGATAAGGTCTGGCGCGCTGc  >  1:596000/1‑62 (MQ=255)               cTAATGCCTGCACTTCGCGCAGAGAAAAAATCCGTTTCGAGGGATAAGGTCTGGCGCGCTGc  >  1:595264/1‑62 (MQ=255)               cTAATGCCTGCACTTCGCGCAGAGAAAAAATCCGTTTCGAGGGATAAGGTCTGGCGCGCTGc  >  1:482200/1‑62 (MQ=255)               cTAATGCCTGCACTTCGCGCAGAGAAAAAATCCGTTTCGAGGGATAAGGTCTGGCGCGCTGc  >  1:516860/1‑62 (MQ=255)                                       |                                      GTGCGTGACATCGGCTAATGCCTGCACTTCGCGCAGAG‑AAAAATCCGTTTCGAGGGATAAGGTCTGGCGCGCTGC  >  minE/1513690‑1513764

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