Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	1,234,615	+T	intergenic (‑8/+138)	dsbB ← / ← nhaB	oxidoreductase that catalyzes reoxidation of DsbA protein disulfide isomerase I/sodium:proton antiporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	1,234,612	1	.	T	83.8%	106.4 / 16.2	37	intergenic (‑5/+141)	dsbB/nhaB	oxidoreductase that catalyzes reoxidation of DsbA protein disulfide isomerase I/sodium:proton antiporter
Reads supporting (aligned to +/- strand):  ref base . (0/6);  new base T (31/0);  total (31/6)
Fisher's exact test for biased strand distribution p-value = 4.30e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01

AGCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCA‑TTTCCCTGCAATA  >  W3110S.gb/1234561‑1234625                                                     |              aGCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCA‑TTTCCCTGCAata  <  1:2086436/65‑1 (MQ=255) aGCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCA‑TTTCCCTGCAata  <  1:1117857/65‑1 (MQ=255) aGCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCA‑TTTCCCTGCAata  <  1:831697/65‑1 (MQ=255) aGCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCA‑TTTCCCTGCAata  <  1:2200865/65‑1 (MQ=255)  gCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTTCCCTGCAata  >  1:331710/1‑65 (MQ=255)  gCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTTCCCTGCAata  >  1:2086588/1‑65 (MQ=255)  gCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTTCCCTGCAata  >  1:2121439/1‑65 (MQ=255)  gCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTTCCCTGCAata  >  1:2121742/1‑65 (MQ=255)  gCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTTCCCTGCAata  >  1:2157714/1‑65 (MQ=255)  gCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTTCCCTGCAata  >  1:2205909/1‑65 (MQ=255)  gCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTTCCCTGCAata  >  1:251061/1‑65 (MQ=255)  gCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTTCCCTGCAata  >  1:2068511/1‑65 (MQ=255)  gCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTTCCCTGCAata  >  1:422776/1‑65 (MQ=255)  gCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTTCCCTGCAata  >  1:585021/1‑65 (MQ=255)  gCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTTCCCTGCAata  >  1:699239/1‑65 (MQ=255)  gCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTTCCCTGCAata  >  1:712532/1‑65 (MQ=255)  gCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTTCCCTGCAata  >  1:828868/1‑65 (MQ=255)  gCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTTCCCTGCAata  >  1:901732/1‑65 (MQ=255)  gCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTTCCCTGCAata  >  1:1475727/1‑65 (MQ=255)  gCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTTCCCTGCAata  >  1:1116033/1‑65 (MQ=255)  gCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTTCCCTGCAata  >  1:1118375/1‑65 (MQ=255)  gCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTTCCCTGCAata  >  1:1233914/1‑65 (MQ=255)  gCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTTCCCTGCAata  >  1:1307736/1‑65 (MQ=255)  gCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTTCCCTGCAata  >  1:1315542/1‑65 (MQ=255)  gCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTTCCCTGCAata  >  1:1347865/1‑65 (MQ=255)  gCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTTCCCTGCAata  >  1:138147/1‑65 (MQ=255)  gCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTTCCCTGCAata  >  1:1442647/1‑65 (MQ=255)  gCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTTCCCTGCAata  >  1:1646213/1‑65 (MQ=255)  gCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTTCCCTGCAata  >  1:1739257/1‑65 (MQ=255)  gCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTTCCCTGCAata  >  1:1797594/1‑65 (MQ=255)  gCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTTCCCTGCAata  >  1:1845959/1‑65 (MQ=255)  gCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTTCCCTGCAata  >  1:1988918/1‑65 (MQ=255)  gCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTTCCCTGCAata  >  1:20614/1‑65 (MQ=255)  gCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTTCCCTGCAat   >  1:1567452/1‑64 (MQ=255)  gCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCATTTTCCCTGCAat   >  1:591531/1‑64 (MQ=255)  gCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCA‑TTTCCCTGCAata  <  1:1100413/64‑1 (MQ=255)     aCGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCA‑TTTCCCTGCAata  <  1:2261393/61‑1 (MQ=255)                                                     |              AGCCACGCGCCCCGGCCTTGTGAACATTGGTTCAAAAATCGCAACATAATCA‑TTTCCCTGCAATA  >  W3110S.gb/1234561‑1234625

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