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

GGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGG  >  W3110S.gb/4563522‑4563654                                                                    |                                                                  ggTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCattt                                                                <  1:4345284/71‑1 (MQ=255)  gTCGGGAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCatttt                                                               <  1:2878056/71‑1 (MQ=255)  gTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCatttt                                                               <  1:4322746/71‑1 (MQ=255)  gTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCatttt                                                               <  1:1775741/71‑1 (MQ=255)     ggCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTAc                                                             <  1:813938/70‑1 (MQ=255)           agcTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGcc                                                      >  1:4499338/1‑71 (MQ=255)           agcTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTAATGTCGcc                                                      >  1:1292923/1‑71 (MQ=255)            gcTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGccc                                                     <  1:3832673/71‑1 (MQ=255)             cTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCa                                                    <  1:2221846/71‑1 (MQ=255)             cTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCa                                                    <  1:1380950/71‑1 (MQ=255)             cTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCa                                                    <  1:3570225/71‑1 (MQ=255)             cTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCa                                                    <  1:3221084/71‑1 (MQ=255)                    aTCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGaa                                             <  1:3207992/71‑1 (MQ=255)                      cGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAaca                                           >  1:2533027/1‑71 (MQ=255)                       ggCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAacac                                          <  1:455566/71‑1 (MQ=255)                       ggCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAacac                                          <  1:4406259/71‑1 (MQ=255)                          acaTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGt                                       <  1:4239806/71‑1 (MQ=255)                          acaTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGt                                       <  1:3105797/71‑1 (MQ=255)                          acaTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGt                                       <  1:3416696/71‑1 (MQ=255)                             tCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGt                                       <  1:241757/68‑1 (MQ=255)                               gCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCa                                                    >  1:408829/1‑53 (MQ=255)                               gCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGc                                  >  1:4075192/1‑71 (MQ=255)                                cAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGcc                                 <  1:1508954/71‑1 (MQ=255)                                cAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGcc                                 <  1:336817/71‑1 (MQ=255)                                   ttttCCAGCAGGTCATAAATAGGCTGTAACGCCatttt                                                               <  1:3399253/38‑1 (MQ=255)                                   ttttCCAGCAGGTCATAAATAGGCTGTAACGCCatttt                                                               <  1:2637619/38‑1 (MQ=255)                                     ttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGtt                            >  1:861319/1‑71 (MQ=255)                                        cagcagGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTa                          <  1:220136/70‑1 (MQ=255)                                        cagcagGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTa                          <  1:1881263/70‑1 (MQ=255)                                              gTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGt                                       <  1:4549417/51‑1 (MQ=255)                                               tCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAg                  >  1:1737678/1‑71 (MQ=255)                                                cATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAg                  >  1:706735/1‑70 (MQ=255)                                                cATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAg                  >  1:4586615/1‑70 (MQ=255)                                                cATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAg                  >  1:1660890/1‑70 (MQ=255)                                                  tAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAac                                            <  1:3307724/42‑1 (MQ=255)                                                  tAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAac                                            <  1:469968/42‑1 (MQ=255)                                                  tAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCaa               <  1:3924572/71‑1 (MQ=255)                                                  tAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCaa               <  1:1608994/71‑1 (MQ=255)                                                    aaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTgg                                               >  1:2908502/1‑37 (MQ=255)                                                    aaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAAc             >  1:1003586/1‑71 (MQ=255)                                                      tAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCg           >  1:3748153/1‑71 (MQ=255)                                                       aGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGc          >  1:470998/1‑71 (MQ=255)                                                          cTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGt                                       <  1:4468144/39‑1 (MQ=255)                                                          cTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGt                                       <  1:1933783/39‑1 (MQ=255)                                                          cTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGt                                       <  1:1005145/39‑1 (MQ=255)                                                             tAACGCCATTTTACTGTCGCCCAAGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAAcgc    >  1:1308914/1‑71 (MQ=255)                                                              aaCGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAAcgcg   >  1:1844774/1‑71 (MQ=255)                                                               aCGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTAc                         >  1:2215424/1‑48 (MQ=255)                                                               aCGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCgg  >  1:1927004/1‑71 (MQ=255)                                                               aCGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCgg  >  1:540615/1‑71 (MQ=255)                                                               aCGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCgg  >  1:2739848/1‑71 (MQ=255)                                                                    |                                                                  GGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGG  >  W3110S.gb/4563522‑4563654

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