Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	2,812,966	T→A	I142F (ATT→TTT)	luxS ←	S‑ribosylhomocysteinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,812,966	0	T	A	80.0%	52.1 / 9.2	25	I142F (ATT→TTT)	luxS	S‑ribosylhomocysteinase
Reads supporting (aligned to +/- strand):  ref base T (0/5);  new base A (20/0);  total (20/5)
Fisher's exact test for biased strand distribution p-value = 1.88e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

GTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTC  >  W3110S.gb/2812945‑2813026                      |                                                             gTTCCAGAATGCTACGCGCAAAATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCcacag                   >  1:2843250/1‑64 (MQ=255) gTTCCAGAATGCTACGCGCAAAATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCacac                   >  1:2149050/1‑65 (MQ=255) gTTCCAGAATGCTACGCGCAAAATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCacac                   >  1:41013/1‑65 (MQ=255) gTTCCAGAATGCTACGCGCAAAATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCacac                   >  1:346213/1‑65 (MQ=255) gTTCCAGAATGCTACGCGCAAAATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCacac                   >  1:2916266/1‑65 (MQ=255) gTTCCAGAATGCTACGCGCAAAATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCacac                   >  1:2885606/1‑65 (MQ=255) gTTCCAGAATGCTACGCGCAAAATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCacac                   >  1:2517110/1‑65 (MQ=255) gTTCCAGAATGCTACGCGCAAAATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCacac                   >  1:2418850/1‑65 (MQ=255) gTTCCAGAATGCTACGCGCAAAATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCacac                   >  1:2393480/1‑65 (MQ=255) gTTCCAGAATGCTACGCGCAAAATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCacac                   >  1:2238832/1‑65 (MQ=255) gTTCCAGAATGCTACGCGCAAAATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCacac                   >  1:2045647/1‑65 (MQ=255) gTTCCAGAATGCTACGCGCAAAATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCacac                   >  1:1894258/1‑65 (MQ=255) gTTCCAGAATGCTACGCGCAAAATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCacac                   >  1:1566241/1‑65 (MQ=255) gTTCCAGAATGCTACGCGCAAAATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCacac                   >  1:443505/1‑65 (MQ=255) gTTCCAGAATGCTACGCGCAAAATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCacac                   >  1:648130/1‑65 (MQ=255) gTTCCAGAATGCTACGCGCAAAATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCacac                   >  1:1332401/1‑65 (MQ=255) gTTCCAGAATGCTACGCGCAAAATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCacac                   >  1:1223999/1‑65 (MQ=255) gTTCCAGAATGCTACGCGCAAAATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCacac                   >  1:1150137/1‑65 (MQ=255) gTTCCAGAATGCTACGCGCAAAATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCacac                   >  1:815129/1‑65 (MQ=255) gTTCCAGAATGCTACGCGCAAAATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGAGCCacac                   >  1:1089586/1‑65 (MQ=255)                  gcAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTc  <  1:1516104/65‑1 (MQ=255)                  gcAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTc  <  1:581709/65‑1 (MQ=255)                  gcAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTc  <  1:1450212/65‑1 (MQ=255)                  gcAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTc  <  1:1040777/65‑1 (MQ=255)                  gcAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTc  <  1:937205/65‑1 (MQ=255)                      |                                                             GTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTC  >  W3110S.gb/2812945‑2813026

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