Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	608,903	C→A	A137S (GCA→TCA)	entD ←	phosphopantetheinyltransferase component of enterobactin synthase multienzyme complex

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	608,903	0	C	A	75.0%	23.2 / 12.1	20	A137S (GCA→TCA)	entD	phosphopantetheinyltransferase component of enterobactin synthase multienzyme complex
Reads supporting (aligned to +/- strand):  ref base C (3/2);  new base A (0/15);  total (3/17)
Fisher's exact test for biased strand distribution p-value = 8.77e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTTGGCGGAAAATGCCAGTGTCAGCGCCAGAGAAAAGGCTAAACCGCAGTCTGCGAGTCGCTCGTGTTCCGCTGGTGTAATAATGTT  >  W3110S.gb/608849‑608936                                                       |                                  ctTTGGCGGAAAATGCCAGTGTCAGCGCCAGAGAAAAGGCTAAACCGCAGTCTGCGAGTCg                             >  1:1259545/1‑61 (MQ=255) ctTTGGCGGAAAATGCCAGTGTCAGCGCCAGAGAAAAGGCTAAACCGCAGTCTGCGAGTCGc                            >  1:2392823/1‑62 (MQ=255) ctTTGGCGGAAAATGCCAGTGTCAGCGCCAGAGAAAAGGCTAAACCGCAGTCTGCGAGTCGc                            >  1:2423996/1‑62 (MQ=255)               gCCAGTGTCAGCGCCAGAGAAAAGGCTAAACCGCAGTCTGCGAGTCGCTCGTGTTCCGCTgg              <  1:2067195/62‑1 (MQ=255)               gCCAGTGTCAGCGCCAGAGAAAAGGCTAAACCGCAGTCTGCGAGTCGCTCGTGTTCCGCTgg              <  1:2614748/62‑1 (MQ=255)                                         tAAACCGCAGTCTGAGAGTCGCTCGTGTTCCGCTGGTGTAATAAtgtt  <  1:847307/48‑1 (MQ=255)                                         tAAACCGCAGTCTGAGAGTCGCTCGTGTTCCGCTGGTGTAATAAtgtt  <  1:1247501/48‑1 (MQ=255)                                         tAAACCGCAGTCTGAGAGTCGCTCGTGTTCCGCTGGTGTAATAAtgtt  <  1:527361/48‑1 (MQ=255)                                         tAAACCGCAGTCTGAGAGTCGCTCGTGTTCCGCTGGTGTAATAAtgtt  <  1:329973/48‑1 (MQ=255)                                         tAAACCGCAGTCTGAGAGTCGCTCGTGTTCCGCTGGTGTAATAAtgtt  <  1:3050005/48‑1 (MQ=255)                                         tAAACCGCAGTCTGAGAGTCGCTCGTGTTCCGCTGGTGTAATAAtgtt  <  1:2763169/48‑1 (MQ=255)                                         tAAACCGCAGTCTGAGAGTCGCTCGTGTTCCGCTGGTGTAATAAtgtt  <  1:2716794/48‑1 (MQ=255)                                         tAAACCGCAGTCTGAGAGTCGCTCGTGTTCCGCTGGTGTAATAAtgtt  <  1:2713346/48‑1 (MQ=255)                                         tAAACCGCAGTCTGAGAGTCGCTCGTGTTCCGCTGGTGTAATAAtgtt  <  1:2443852/48‑1 (MQ=255)                                         tAAACCGCAGTCTGAGAGTCGCTCGTGTTCCGCTGGTGTAATAAtgtt  <  1:2439012/48‑1 (MQ=255)                                         tAAACCGCAGTCTGAGAGTCGCTCGTGTTCCGCTGGTGTAATAAtgtt  <  1:215/48‑1 (MQ=255)                                         tAAACCGCAGTCTGAGAGTCGCTCGTGTTCCGCTGGTGTAATAAtgtt  <  1:2038890/48‑1 (MQ=255)                                         tAAACCGCAGTCTGAGAGTCGCTCGTGTTCCGCTGGTGTAATAAtgtt  <  1:1777197/48‑1 (MQ=255)                                         tAAACCGCAGTCTGAGAGTCGCTCGTGTTCCGCTGGTGTAATAAtgtt  <  1:1649111/48‑1 (MQ=255)                                          aaaCCGCAGTCTGAGAGTCGCTCGTGTTCCGCTGGTGTAATAAtgtt  <  1:2698006/47‑1 (MQ=255)                                                       |                                  CTTTGGCGGAAAATGCCAGTGTCAGCGCCAGAGAAAAGGCTAAACCGCAGTCTGCGAGTCGCTCGTGTTCCGCTGGTGTAATAATGTT  >  W3110S.gb/608849‑608936

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