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	81.8%	128.4 / 23.8	55	K207M (AAG→ATG)	iadA	isoaspartyl dipeptidase
Reads supporting (aligned to +/- strand):  ref base T (10/0);  new base A (30/15);  total (40/15)
Fisher's exact test for biased strand distribution p-value = 4.70e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

CTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAA  >  W3110S.gb/4563534‑4563641                                                        |                                                     cTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTAc                                                <  1:105208/62‑1 (MQ=255) cTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTAc                                                <  1:4316143/62‑1 (MQ=255) cTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTAc                                                <  1:4191955/62‑1 (MQ=255) cTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTAc                                                <  1:2333384/62‑1 (MQ=255)      tGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTc                                            <  1:3038488/61‑1 (MQ=255)          cGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCa                                       >  1:2835677/1‑62 (MQ=255)          cGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCa                                       <  1:732539/62‑1 (MQ=255)          cGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCa                                       <  1:4247735/62‑1 (MQ=255)          cGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCa                                       >  1:3998207/1‑62 (MQ=255)          cGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCa                                       >  1:3329145/1‑62 (MQ=255)          cGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCa                                       >  1:2768431/1‑62 (MQ=255)          cGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCa                                       >  1:2248980/1‑62 (MQ=255)          cGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCa                                       >  1:2159017/1‑62 (MQ=255)          cGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCa                                       >  1:1525365/1‑62 (MQ=255)                         ttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTAc                                                <  1:4055307/38‑1 (MQ=255)                             agcagGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGc                     <  1:2368863/61‑1 (MQ=255)                             agcagGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGc                     <  1:3765462/61‑1 (MQ=255)                                agGTCATAAATAGGCTGTAACGCCATTTTACTGTCGcc                                         >  1:205974/1‑38 (MQ=255)                                agGTCATAAATAGGCTGTAACGCCATTTTACTGTCGcc                                         >  1:835537/1‑38 (MQ=255)                                agGTCATAAATAGGCTGTAACGCCATTTTACTGTCGcc                                         >  1:1386022/1‑38 (MQ=255)                                agGTCATAAATAGGCTGTAACGCCATTTTACTGTCGcc                                         >  1:3472815/1‑38 (MQ=255)                                agGTCATAAATAGGCTGTAACGCCATTTTACTGTCGcc                                         >  1:3395387/1‑38 (MQ=255)                                agGTCATAAATAGGCTGTAACGCCATTTTACTGTCGcc                                         >  1:2890439/1‑38 (MQ=255)                                        aaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGc                     <  1:2177401/50‑1 (MQ=255)                                        aaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGc                     <  1:718266/50‑1 (MQ=255)                                        aaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGc                     <  1:1208070/50‑1 (MQ=255)                                        aaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGc                     <  1:3857492/50‑1 (MQ=255)                                        aaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGc                     <  1:298953/50‑1 (MQ=255)                                          tAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAac                               >  1:1855047/1‑38 (MQ=255)                                          tAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAac                               >  1:1219067/1‑38 (MQ=255)                                          tAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAac                               >  1:250853/1‑38 (MQ=255)                                          tAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAac                               >  1:3774844/1‑38 (MQ=255)                                          tAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAac                               >  1:3720817/1‑38 (MQ=255)                                          tAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAac                               >  1:3715609/1‑38 (MQ=255)                                          tAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAac                               >  1:2763844/1‑38 (MQ=255)                                          tAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAac                               >  1:3406661/1‑38 (MQ=255)                                          tAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAac                               >  1:4561713/1‑38 (MQ=255)                                          tAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAac                               >  1:3386818/1‑38 (MQ=255)                                          tAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAac                               >  1:2817820/1‑38 (MQ=255)                                           aGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTAccgccg       >  1:2536957/1‑61 (MQ=255)                                           aGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTAccgccg       >  1:1130156/1‑61 (MQ=255)                                           aGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTAccgccg       >  1:2449688/1‑61 (MQ=255)                                           aGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTAccgccg       >  1:4486917/1‑61 (MQ=255)                                             gCTGTAACGCCTTTTTGCTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGa      >  1:1854455/1‑60 (MQ=255)                                             gCTGTAACGCCTTTTTGCTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAg     >  1:4561203/1‑61 (MQ=255)                                             gCTGTAACGCCTTTTTGCTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAg     >  1:935550/1‑61 (MQ=255)                                             gCTGTAACGCCTTTTTGCTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAg     >  1:4328161/1‑61 (MQ=255)                                             gCTGTAACGCCTTTTTGCTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAg     >  1:2348520/1‑61 (MQ=255)                                             gCTGTAACGCCTTTTTGCTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAg     >  1:131734/1‑61 (MQ=255)                                             gCTGTAACGCCTTTTTGCTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAg     >  1:2442086/1‑61 (MQ=255)                                             gCTGTAACGCCTTTTTGCTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAg     >  1:178707/1‑61 (MQ=255)                                             gCTGTAACGCCTTTTTGCTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAg     >  1:3292467/1‑61 (MQ=255)                                             gCTGTAACGCCTTTTTGCTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAg     >  1:1737214/1‑61 (MQ=255)                                               tGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCaa  >  1:283309/1‑62 (MQ=255)                                               tGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCaa  >  1:1011466/1‑62 (MQ=255)                                                        |                                                     CTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAA  >  W3110S.gb/4563534‑4563641

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