Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,742,599	Δ1 bp	intergenic (‑67/‑134)	yiaJ ← / → yiaK	transcriptional repressor for the yiaKLMNO‑lyxK‑sgbHUE operon/2,3‑diketo‑L‑gulonate reductase, NADH‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,742,598	0	C	.	93.8%	56.2 / ‑1.8	16	intergenic (‑66/‑135)	yiaJ/yiaK	transcriptional repressor for the yiaKLMNO‑lyxK‑sgbHUE operon/2,3‑diketo‑L‑gulonate reductase, NADH‑dependent
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base . (7/8);  total (8/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.12e-01

GGCTTCCGGCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCAGATTGCCTTGT  >  NC_000913/3742463‑3742645                                                                                                                                        |                                                ggCTTCCGGCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTAT‑ACCAAt                                             >  1:815953/1‑139 (MQ=255)     tCCGGCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGAAAGATAATTTATAA‑CAATTGAt                                         >  1:315804/1‑139 (MQ=255)        ggCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAA‑CAATTGAtttt                                      >  1:201057/1‑139 (MQ=255)        ggCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGACACGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGTTAATTTATAA‑CAATTGAtttt                                      >  1:588619/1‑139 (MQ=255)                      tttCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAA‑CAATTGATTTTTATGTCTTTTGaaa                        <  2:485251/139‑1 (MQ=255)                      tttCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAA‑CAATTGATTTTTATGTCTTTTGaaa                        <  2:994237/139‑1 (MQ=255)                            gcgcCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAA‑CAATTGATTTTTATGTCTTTTGAAATtcatc                  >  2:759769/1‑139 (MQ=255)                             cgcCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAA‑CAATTGATTTTTATGTCTTTTGAAATtcatca                 <  1:635888/139‑1 (MQ=255)                             cgcCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAA‑CAATTGATTTTTATGTCTTTTGAAATtcatca                 <  2:201057/139‑1 (MQ=255)                                     cGTTCTCTTATTTTCCCATCACTTCTTTCCCCATCTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAA‑CAATTGATTTTTATGTCTTTTGAAATTCatcaatcctgtt         <  2:647118/139‑6 (MQ=255)                                       ttCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAA‑CAATTGATTTTTATGTCTTTTGAAATTCatcaatcctgttga       <  2:588619/139‑8 (MQ=255)                                            tttttttCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAA‑CAATTGATTTTTATGTCTTTTGAAATTCatcaatcctgttgacaatt  >  1:923364/1‑127 (MQ=255)                                                                                           gcgcAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAA‑CAATTGATTTTTATGTCTTTTGAAATtcatca                 >  1:369363/1‑77 (MQ=255)                                                                                           gcgcAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAA‑CAATTGATTTTTATGTCTTTTGAAATtcatca                 <  1:645774/77‑1 (MQ=255)                                                                                           gcgcAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAA‑CAATTGATTTTTATGTCTTTTGAAATtcatca                 <  2:369363/77‑1 (MQ=255)                                                                                           gcgcAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAA‑CAATTGATTTTTATGTCTTTTGAAATtcatca                 >  2:645774/1‑77 (MQ=255)                                                                                                                                        |                                                GGCTTCCGGCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCAGATTGCCTTGT  >  NC_000913/3742463‑3742645

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑