Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	4,257,559	C→T	intergenic (+79/+76)	ubiA → / ← plsB	p‑hydroxybenzoate octaprenyltransferase/glycerol‑3‑phosphate O‑acyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	4,257,559	0	C	T	80.0%	10.5 / ‑0.3	10	intergenic (+79/+76)	ubiA/plsB	p‑hydroxybenzoate octaprenyltransferase/glycerol‑3‑phosphate O‑acyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base T (0/8);  total (0/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGGCCTACACAACACATCGCAATTTATTGAATTTGCAGATTATGGAAGGCCGGATAAGGCGTTTTCGCCGCA  >  W3110S.gb/4257544‑4257616                |                                                          cAGGCCTACACAACATATCGCAATTTATTGAATTTGCAGATCATGGTAGGCCGAATAAgg               <  1:130026/60‑1 (MQ=25) cAGGCCTACACAACATATCGCAATTTATTGAATTTGCAGATCATGGTAGGCCGAATAAgg               <  1:1386136/60‑1 (MQ=25) cAGGCCTACACAACATATCGCAATTTATTGAATTTGCAGATCATGGTAGGCCGAATAAgg               <  1:1388211/60‑1 (MQ=25) cAGGCCTACACAACATATCGCAATTTATTGAATTTGCAGATCATGGTAGGCCGAATAAgg               <  1:359298/60‑1 (MQ=25) cAGGCCTACACAACATATCGCAATTTATTGAATTTGCAGATCATGGTAGGCCGAATAAgg               <  1:1485751/60‑1 (MQ=25) cAGGCCTACACAACATATCGCAATTTATTGAATTTGCAGATCATGGTAGGCCGAATAAgg               <  1:3000732/60‑1 (MQ=25) cAGGCCTACACAACATATCGCAATTTATTGAATTTGCAGATCATGGTAGGCCGAATAAgg               <  1:2943669/60‑1 (MQ=25) cAGGCCTACACAACATATCGCAATTTATTGAATTTGCAGATCATGGTAGGCCGAATAAgg               <  1:2205746/60‑1 (MQ=25) cAGGCCTACACAACACATCGCAATTTATTGAATTTGCAGATTATGGAAGGCCGGATAAgg               <  1:1094031/60‑1 (MQ=255) cAGGCCTACACAACACATCGCAATTTATTGAATTTGCAGATTATGGAAGGCCGGATAAgg               <  1:2334475/60‑1 (MQ=255)          acaacaccTCGCAATTTATTGAATTTGCAGATTATGGAAGGCCGGATAAgg               <  1:772115/51‑1 (MQ=255)            aacacaTCGCAATTTATTGAATTTGCAGATTATGGAAGGCCGGATAAGGCGTTTTCGCCGCa  <  1:2246985/62‑1 (MQ=255)            aacacaTCGCAATTTATTGAATTTGCAGATTATGGAAGGCCGGATAAGGCGTTTTCGCCGCa  <  1:1979163/62‑1 (MQ=255)            aacacaTCGCAATTTATTGAATTTGCAGATTATGGAAGGCCGGATAAGGCGTTTTCGCCGCa  <  1:2564043/62‑1 (MQ=255)            aacacaTCGCAATTTATTGAATTTGCAGATTATGGAAGGCCGGATAAGGCGTTTTCGCCGCa  <  1:614574/62‑1 (MQ=255)               acaTCGCAATTTATTGAATTTGCAGATTATGGAAGGCCGGATAAGGCGTTTTcgccgc   <  1:2375233/58‑1 (MQ=255)               acaTCGCAATTTATTGAATTTGCAGATTATGGAAGGCCGGATAAGGCGTTTTcgccgc   <  1:1914413/58‑1 (MQ=255)               acaTCGCAATTTATTGAATTTGCAGATTATGGAAGGCCGGATAAGGCGTTTTcgccgc   <  1:341015/58‑1 (MQ=255)               acaTCGCAATTTATTGAATTTGCAGATTATGGAAGGCCGGATAAGGCGTTTTcgccgc   <  1:1485/58‑1 (MQ=255)               acaTCGCAATTTATTGAATTTGCAGATTATGGAAGGCCGGATAAGGCGTTTTcgccgc   <  1:477422/58‑1 (MQ=255)               acaTCGCAATTTATTGAATTTGCAGATTATGGAAGGCCGGATAAGGCGTTTTcgccgc   <  1:1333085/58‑1 (MQ=255)                caTCGCAATTTATTGAATTTGCAGATTATGGAAGGCCGGATAAGGCGTTTTcgccgc   <  1:1918434/57‑1 (MQ=255)                |                                                          CAGGCCTACACAACACATCGCAATTTATTGAATTTGCAGATTATGGAAGGCCGGATAAGGCGTTTTCGCCGCA  >  W3110S.gb/4257544‑4257616

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