Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	4,563,725	+T	coding (484/1173 nt)	iadA ←	isoaspartyl dipeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	4,563,722	1	.	T	89.7%	115.0 / 7.7	39	coding (487/1173 nt)	iadA	isoaspartyl dipeptidase
Reads supporting (aligned to +/- strand):  ref base . (4/0);  new base T (35/0);  total (39/0)
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

GGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCA‑TTTCACGCCAATCACACGATCGATA  >  W3110S.gb/4563677‑4563747                                               |                          ggTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCA‑TTTCCCGCCAATcac            >  1:1414976/1‑61 (MQ=255) ggTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCA‑TTTCACg                    >  1:2864916/1‑53 (MQ=255) ggTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCA‑TTTCACGCCAATcac            >  1:2649119/1‑61 (MQ=255) ggTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCA‑TTTCACGCCAATcac            >  1:387596/1‑61 (MQ=255)           gTCCGGTGCGGCAGAACGGTGATCTGAGATGGCGCATTTTCACGCCAATGACACGATCGATa  >  1:2039761/1‑62 (MQ=255)           gTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTTCACGCCACTCACACGATCGATa  >  1:2730598/1‑62 (MQ=255)           gTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTTCACGCCAAt               >  1:2477502/1‑49 (MQ=255)           gTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTTCACGCCAATCACAcgatcgat   >  1:2678557/1‑61 (MQ=255)           gTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTTCACGCCAATCACAcgatcgat   >  1:2185582/1‑61 (MQ=255)           gTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTTCACGCCAATCACACGATCGATa  >  1:327076/1‑62 (MQ=255)           gTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTTCACGCCAATCACACGATCGATa  >  1:2679103/1‑62 (MQ=255)           gTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTTCACGCCAATCACACGATCGATa  >  1:2794504/1‑62 (MQ=255)           gTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTTCACGCCAATCACACGATCGATa  >  1:2913986/1‑62 (MQ=255)           gTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTTCACGCCAATCACACGATCGATa  >  1:2488385/1‑62 (MQ=255)           gTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTTCACGCCAATCACACGATCGATa  >  1:625874/1‑62 (MQ=255)           gTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTTCACGCCAATCACACGATCGATa  >  1:665918/1‑62 (MQ=255)           gTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTTCACGCCAATCACACGATCGATa  >  1:751412/1‑62 (MQ=255)           gTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTTCACGCCAATCACACGATCGATa  >  1:871511/1‑62 (MQ=255)           gTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTTCACGCCAATCACACGATCGATa  >  1:904405/1‑62 (MQ=255)           gTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTTCACGCCAATCACACGATCGATa  >  1:917021/1‑62 (MQ=255)           gTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTTCACGCCAATCACACGATCGATa  >  1:935040/1‑62 (MQ=255)           gTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTTCACGCCAATCACACGATCGATa  >  1:945449/1‑62 (MQ=255)           gTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTTCACGCCAATCACACGATCGATa  >  1:1706972/1‑62 (MQ=255)           gTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTTCACGCCAATCACACGATCGATa  >  1:1059459/1‑62 (MQ=255)           gTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTTCACGCCAATCACACGATCGATa  >  1:1089451/1‑62 (MQ=255)           gTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTTCACGCCAATCACACGATCGATa  >  1:1237843/1‑62 (MQ=255)           gTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTTCACGCCAATCACACGATCGATa  >  1:132868/1‑62 (MQ=255)           gTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTTCACGCCAATCACACGATCGATa  >  1:1403843/1‑62 (MQ=255)           gTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTTCACGCCAATCACACGATCGATa  >  1:1624057/1‑62 (MQ=255)           gTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTTCACGCCAATCACACGATCGATa  >  1:1638008/1‑62 (MQ=255)           gTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTTCACGCCAATCACACGATCGATa  >  1:1654115/1‑62 (MQ=255)           gTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTTCACGCCAATCACACGATCGATa  >  1:2662644/1‑62 (MQ=255)           gTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTTCACGCCAATCACACGATCGATa  >  1:2006381/1‑62 (MQ=255)           gTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTTCACGCCAATCACACGATCGATa  >  1:2015879/1‑62 (MQ=255)           gTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTTCACGCCAATCACACGATCGATa  >  1:2036838/1‑62 (MQ=255)           gTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTTCACGCCAATCACACGATCGATa  >  1:2225912/1‑62 (MQ=255)           gTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTTCACGCCAATCACACGATCGATa  >  1:2403921/1‑62 (MQ=255)           gTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTTCACGCCAATCACACGATCGATa  >  1:1049090/1‑62 (MQ=255)           gTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTTCACGCCAATCACACGATCGATa  >  1:2633517/1‑62 (MQ=255)                                               |                          GGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCA‑TTTCACGCCAATCACACGATCGATA  >  W3110S.gb/4563677‑4563747

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