Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	1,155,387	+G	coding (511/810 nt)	pabC →	4‑amino‑4‑deoxychorismate lyase component of para‑aminobenzoate synthase multienzyme complex

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	1,155,384	1	.	G	97.5%	146.8 / ‑3.2	40	coding (508/810 nt)	pabC	4‑amino‑4‑deoxychorismate lyase component of para‑aminobenzoate synthase multienzyme complex
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (0/39);  total (0/40)
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.99e-01

GAGGCGCTGGTCCTTGACAGCGAAGGGT‑GG‑GTTACGGAATGCTGTGCGGCTAATTTGTTCTGG  >  W3110S.gb/1155357‑1155419                             |                                     gaGTCGCTGGTCCTTGACAGCGAAGGGTGGG‑GTTACGGAATGCTGTGCGGCTAATTTGTTCTgg  <  1:2430363/64‑1 (MQ=255) gaGGCGCTGGTCCTTGACAGCGAAGGGTGGG‑GTTACGGAATGCTGTGCGGCTAATTTGTTCTgg  <  1:414170/64‑1 (MQ=255) gaGGCGCTGGTCCTTGACAGCGAAGGGTGGG‑GTTACGGAATGCTGTGCGGCTAATTTGTTCTgg  <  1:1114318/64‑1 (MQ=255) gaGGCGCTGGTCCTTGACAGCGAAGGGTGGG‑GTTACGGAATGCTGTGCGGCTAATTTGTTCTgg  <  1:2266919/64‑1 (MQ=255) gaGGCGCTGGTCCTTGACAGCGAAGGGTGGG‑GTTACGGAATGCTGTGCGGCTAATTTGTTCTgg  <  1:2338210/64‑1 (MQ=255) gaGGCGCTGGTCCTTGACAGCGAAGGGTGGG‑GTTACGGAATGCTGTGCGGCTAATTTGTTCTgg  <  1:2369495/64‑1 (MQ=255) gaGGCGCTGGTCCTTGACAGCGAAGGGTGGG‑GTTACGGAATGCTGTGCGGCTAATTTGTTCTgg  <  1:2432208/64‑1 (MQ=255) gaGGCGCTGGTCCTTGACAGCGAAGGGTGGG‑GTTACGGAATGCTGTGCGGCTAATTTGTTCTgg  <  1:293944/64‑1 (MQ=255) gaGGCGCTGGTCCTTGACAGCGAAGGGTGGG‑GTTACGGAATGCTGTGCGGCTAATTTGTTCTgg  <  1:338020/64‑1 (MQ=255) gaGGCGCTGGTCCTTGACAGCGAAGGGTGGG‑GTTACGGAATGCTGTGCGGCTAATTTGTTCTgg  <  1:401485/64‑1 (MQ=255) gaGGCGCTGGTCCTTGACAGCGAAGGGTGGG‑GTTACGGAATGCTGTGCGGCTAATTTGTTCTgg  <  1:2189973/64‑1 (MQ=255) gaGGCGCTGGTCCTTGACAGCGAAGGGTGGG‑GTTACGGAATGCTGTGCGGCTAATTTGTTCTgg  <  1:427695/64‑1 (MQ=255) gaGGCGCTGGTCCTTGACAGCGAAGGGTGGG‑GTTACGGAATGCTGTGCGGCTAATTTGTTCTgg  <  1:553633/64‑1 (MQ=255) gaGGCGCTGGTCCTTGACAGCGAAGGGTGGG‑GTTACGGAATGCTGTGCGGCTAATTTGTTCTgg  <  1:607035/64‑1 (MQ=255) gaGGCGCTGGTCCTTGACAGCGAAGGGTGGG‑GTTACGGAATGCTGTGCGGCTAATTTGTTCTgg  <  1:746205/64‑1 (MQ=255) gaGGCGCTGGTCCTTGACAGCGAAGGGTGGG‑GTTACGGAATGCTGTGCGGCTAATTTGTTCTgg  <  1:821919/64‑1 (MQ=255) gaGGCGCTGGTCCTTGACAGCGAAGGGTGGG‑GTTACGGAATGCTGTGCGGCTAATTTGTTCTgg  <  1:894853/64‑1 (MQ=255) gaGGCGCTGGTCCTTGACAGCGAAGGGTGGG‑GTTACGGAATGCTGTGCGGCTAATTTGTTCTgg  <  1:907336/64‑1 (MQ=255) gaGGCGCTGGTCCTTGACAGCGAAGGGTGGG‑GTTACGGAATGCTGTGCGGCTAATTTGTTCTgg  <  1:97898/64‑1 (MQ=255) gaGGCGCTGGTCCTTGACAGCGAAGGGTGGG‑GTTACGGAATGCTGTGCGGCTAATTTGTTCTgg  <  1:1788180/64‑1 (MQ=255) gaGGCGCTGGTCCTTGACAGCGAAGGGTGGG‑GTTACGGAATGCTGTGCGGCTAATTTGTTCTgg  <  1:1226250/64‑1 (MQ=255) gaGGCGCTGGTCCTTGACAGCGAAGGGTGGG‑GTTACGGAATGCTGTGCGGCTAATTTGTTCTgg  <  1:1318109/64‑1 (MQ=255) gaGGCGCTGGTCCTTGACAGCGAAGGGTGGG‑GTTACGGAATGCTGTGCGGCTAATTTGTTCTgg  <  1:1343174/64‑1 (MQ=255) gaGGCGCTGGTCCTTGACAGCGAAGGGTGGG‑GTTACGGAATGCTGTGCGGCTAATTTGTTCTgg  <  1:1456502/64‑1 (MQ=255) gaGGCGCTGGTCCTTGACAGCGAAGGGTGGG‑GTTACGGAATGCTGTGCGGCTAATTTGTTCTgg  <  1:1486116/64‑1 (MQ=255) gaGGCGCTGGTCCTTGACAGCGAAGGGTGGG‑GTTACGGAATGCTGTGCGGCTAATTTGTTCTgg  <  1:1661458/64‑1 (MQ=255) gaGGCGCTGGTCCTTGACAGCGAAGGGTGGG‑GTTACGGAATGCTGTGCGGCTAATTTGTTCTgg  <  1:1737846/64‑1 (MQ=255) gaGGCGCTGGTCCTTGACAGCGAAGGGTGGG‑GTTACGGAATGCTGTGCGGCTAATTTGTTCTgg  <  1:1740285/64‑1 (MQ=255) gaGGCGCTGGTCCTTGACAGCGAAGGGTGGG‑GTTACGGAATGCTGTGCGGCTAATTTGTTCTgg  <  1:1954501/64‑1 (MQ=255) gaGGCGCTGGTCCTTGACAGCGAAGGGTGGG‑GTTACGGAATGCTGTGCGGCTAATTTGTTCTgg  <  1:1958513/64‑1 (MQ=255) gaGGCGCTGGTCCTTGACAGCGAAGGGTGGG‑GTTACGGAATGCTGTGCGGCTAATTTGTTCTgg  <  1:1960406/64‑1 (MQ=255) gaGGCGCTGGTCCTTGACAGCGAAGGGTGGG‑GTTACGGAATGCTGTGCGGCTAATTTGTTCTgg  <  1:2043685/64‑1 (MQ=255) gaGGCGCTGGTCCTTGACAGCGAAGGGTGGG‑GTTACGGAATGCTGTGCGGCTAATTTGTTCTgg  <  1:2059957/64‑1 (MQ=255) gaGGCGCTGGTCCTTGACAGCGAAGGGTGGG‑GTTACGGAATGCTGTGCGGCTAATTTGTTCTgg  <  1:2082262/64‑1 (MQ=255) gaGGCGCTGGTCCTTGACAGCGAAGGGTGGG‑GTTACGGAATGCTGTGCGGCTAATTTGTTCTgg  <  1:2082918/64‑1 (MQ=255) gaGGCGCTGGTCCTTGACAGCGAAGGGTGGG‑GTTACGGAATGCTGTGCGGCTAATTTGTTCTgg  <  1:2117734/64‑1 (MQ=255) gaGGCGCTGGTCCTTGACAGCGAAGGGTGGG‑GGTACGGAATGCTGTGCGGCTAATTTGTTCTgg  <  1:739241/64‑1 (MQ=255) gaGGCGCTGGTCCTTGACAGCGAAGGGT‑GGTGTTACGGAATGCTGTGCGGCTAATTTGTTCTgg  <  1:2119144/64‑1 (MQ=255) gaGGCGCTGGTCCTTGACAACGAAGGGTGGG‑GTTACGGAATGCTGTGCGGCTAATTTGTTCTgg  <  1:1996333/64‑1 (MQ=255)     cgcTGGTCCTTGACAGCGAAGTGTGGG‑GTTACGGAATGCTGTGCGGCTAATTTGTTCTgg  <  1:782591/60‑1 (MQ=255)                             |                                     GAGGCGCTGGTCCTTGACAGCGAAGGGT‑GG‑GTTACGGAATGCTGTGCGGCTAATTTGTTCTGG  >  W3110S.gb/1155357‑1155419

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