Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	4,563,589	T→A	K207M (AAG→ATG)	iadA ←	isoaspartyl dipeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	4,563,589	0	T	A	100.0%	159.3 / NA	48	K207M (AAG→ATG)	iadA	isoaspartyl dipeptidase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (7/41);  total (7/41)

TCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCG  >  W3110S.gb/4563542‑4563633                                                |                                             tCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGccc                                >  1:1282455/1‑62 (MQ=255) tCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGccc                                >  1:1949985/1‑62 (MQ=255) tCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGccc                                >  1:29750/1‑62 (MQ=255) tCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGccc                                >  1:1241389/1‑62 (MQ=255) tCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGccc                                >  1:429122/1‑62 (MQ=255) tCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGccc                                >  1:730926/1‑62 (MQ=255) tCGGCACATAGCAGTTTTCCAGCAGGACATAAATAGGCTGTAACGCCATTTTACTGTCGccc                                >  1:1073030/1‑62 (MQ=255)                    cagcagGTCAT‑AATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGc             <  1:2097723/61‑1 (MQ=255)                     agcagGTCATAAATAGGCTGTGACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGc             <  1:1860269/61‑1 (MQ=255)                     agcagGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGc             <  1:2665166/61‑1 (MQ=255)                     agcagGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGc             <  1:78430/61‑1 (MQ=255)                     agcagGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGc             <  1:1999295/61‑1 (MQ=255)                     agcagGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGc             <  1:2324813/61‑1 (MQ=255)                     agcagGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGc             <  1:2487021/61‑1 (MQ=255)                     agcagGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGc             <  1:2615075/61‑1 (MQ=255)                     agcagGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGc             <  1:2646442/61‑1 (MQ=255)                     agcagGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGc             <  1:2910064/61‑1 (MQ=255)                     agcagGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGc             <  1:877987/61‑1 (MQ=255)                     agcagGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGc             <  1:768520/61‑1 (MQ=255)                     agcagGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGc             <  1:742818/61‑1 (MQ=255)                     agcagGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGc             <  1:59363/61‑1 (MQ=255)                     agcagGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGc             <  1:626963/61‑1 (MQ=255)                     agcagGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGc             <  1:1891629/61‑1 (MQ=255)                     agcagGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGc             <  1:1943169/61‑1 (MQ=255)                     agcagGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGc             <  1:1011598/61‑1 (MQ=255)                     agcagGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGc             <  1:1622175/61‑1 (MQ=255)                     agcagGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGc             <  1:1107714/61‑1 (MQ=255)                     agcagGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGc             <  1:1565341/61‑1 (MQ=255)                     agcagGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGc             <  1:1386547/61‑1 (MQ=255)                     agcagGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGc             <  1:1110615/61‑1 (MQ=255)                     agcagGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGc             <  1:1193756/61‑1 (MQ=255)                      gcagGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGc             <  1:1268335/60‑1 (MQ=255)                               aaaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTAccg  <  1:1008227/62‑1 (MQ=255)                               aaaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTAccg  <  1:1073926/62‑1 (MQ=255)                               aaaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTAccg  <  1:1078641/62‑1 (MQ=255)                               aaaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTAccg  <  1:662358/62‑1 (MQ=255)                               aaaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTAccg  <  1:576038/62‑1 (MQ=255)                               aaaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTAccg  <  1:44486/62‑1 (MQ=255)                               aaaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTAccg  <  1:1305325/62‑1 (MQ=255)                               aaaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTAccg  <  1:2804454/62‑1 (MQ=255)                               aaaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTAccg  <  1:1601476/62‑1 (MQ=255)                               aaaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTAccg  <  1:2435480/62‑1 (MQ=255)                               aaaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTAccg  <  1:2305178/62‑1 (MQ=255)                               aaaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTAccg  <  1:2281184/62‑1 (MQ=255)                               aaaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTAccg  <  1:1741031/62‑1 (MQ=255)                                      cTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGc             <  1:2897819/44‑1 (MQ=255)                                        gTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTAccg  <  1:1339191/53‑1 (MQ=255)                                           aCGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGc             <  1:1115254/39‑1 (MQ=255)                                                |                                             TCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCG  >  W3110S.gb/4563542‑4563633

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