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%	108.4 / NA	36	K207M (AAG→ATG)	iadA	isoaspartyl dipeptidase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (24/12);  total (24/12)

GGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACC  >  W3110S.gb/4563522‑4563644                                                                    |                                                        ggTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCattt                                                      >  1:3167307/1‑71 (MQ=255) ggTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCattt                                                      >  1:155645/1‑71 (MQ=255)           agcTTGCTGATCGGCACCTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCAATTTACTGTCGcc                                            >  1:2076704/1‑71 (MQ=255)           agcTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGcc                                            >  1:2336410/1‑71 (MQ=255)           agcTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGcc                                            >  1:2197333/1‑71 (MQ=255)           agcTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGcc                                            >  1:2575821/1‑71 (MQ=255)           agcTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGcc                                            >  1:2010944/1‑71 (MQ=255)           agcTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGcc                                            >  1:1808786/1‑71 (MQ=255)           agcTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGcc                                            >  1:3232387/1‑71 (MQ=255)           agcTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGcc                                            >  1:1163449/1‑71 (MQ=255)           agcTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGcc                                            >  1:1077504/1‑71 (MQ=255)               tGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtg                                        >  1:3697538/1‑71 (MQ=255)                   gATCGGCATATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtgtg                                      >  1:2408584/1‑69 (MQ=255)                   gATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtgtg                                      >  1:3727236/1‑69 (MQ=255)                   gATCGGCACATCGCAGTTTTCCAGCAGGTAAAAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtgtg                                      >  1:3289910/1‑69 (MQ=255)                          acaTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACgg                              >  1:2330905/1‑70 (MQ=255)                          acaTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGt                             >  1:3428993/1‑71 (MQ=255)                          acaTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGt                             >  1:963581/1‑71 (MQ=255)                          acaTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGt                             >  1:1608212/1‑71 (MQ=255)                             tCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGAc                          <  1:3242779/71‑1 (MQ=255)                                       ccAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGttt                 <  1:2324799/70‑1 (MQ=255)                                          gcagGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGt                             <  1:1470643/55‑1 (MQ=255)                                          gcagGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGt                             <  1:810283/55‑1 (MQ=255)                                          gcagGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGt                             <  1:1439857/55‑1 (MQ=255)                                          gcagGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGt                             <  1:2389830/55‑1 (MQ=255)                                          gcagGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGt                             <  1:2966671/55‑1 (MQ=255)                                          gcagGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGt                             <  1:2180696/55‑1 (MQ=255)                                                cATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAg        >  1:447250/1‑70 (MQ=255)                                                 aTAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGc       >  1:3836782/1‑70 (MQ=255)                                                 aTAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGc       >  1:1132144/1‑70 (MQ=255)                                                    aaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTg                            >  1:2658015/1‑46 (MQ=255)                                                    aaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGa                           >  1:1760228/1‑47 (MQ=255)                                                    aaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGa                           >  1:2916491/1‑47 (MQ=255)                                                    aaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGa                           >  1:1387661/1‑47 (MQ=255)                                                     aTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAAcc  <  1:2917697/71‑1 (MQ=255)                                                     aTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAAcc  <  1:4090411/71‑1 (MQ=255)                                                     aTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAAcc  <  1:410968/71‑1 (MQ=255)                                                       aGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGt                             <  1:1457834/42‑1 (MQ=255)                                                                    |                                                        GGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACC  >  W3110S.gb/4563522‑4563644

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