Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,619,153	0	G	C	40.0%	3.0 / 10.3	10	intergenic (+164/‑39)	nikR/rhsB	transcriptional repressor, Ni‑binding/Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base G (1/5);  new base C (2/2);  total (3/7)
Fisher's exact test for biased strand distribution p-value = 5.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TCTAAATATATTCCTAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAA  >  NC_000913/3619088‑3619170                                                                  |                  tCTAAATATATTCCTAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGaa                 <  1:617698/68‑1 (MQ=255) tCTAAATATATTCCTAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGaa                 <  2:639296/68‑1 (MQ=255) tCTAAATATATTCCTAAAAATCTTCAAATCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGaa                 <  1:27151/68‑1 (MQ=255)    aaaTATATTCCTAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGAt              <  2:631916/68‑1 (MQ=255)      atatatTCCTAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATAc            >  1:14691/1‑68 (MQ=255)                aaaaaTCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTaaa  <  1:535977/68‑1 (MQ=255)                                           aGTTTTTCTTCGCTTTTTCGTATCAAGATACTGTCATTaaa  <  1:68133/41‑1 (MQ=255)                                           aGTTTTTCTTCGCTTTTTCGTATCAAGATACTGTCATTaaa  <  1:69035/41‑1 (MQ=255)                                           aGTTTTTCTTCGCTTTTTCGTATCAAGATACTGTCATTaaa  >  2:68133/1‑41 (MQ=255)                                           aGTTTTTCTTCGCTTTTTCGTATCAAGATACTGTCATTaaa  >  2:69035/1‑41 (MQ=255)                                                                  |                  TCTAAATATATTCCTAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAA  >  NC_000913/3619088‑3619170

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