Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	4,563,589	T→A	100%	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%	95.0 / NA	34	K207M (AAG→ATG)	iadA	isoaspartyl dipeptidase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (22/12);  total (22/12)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGA  >  W3110S.gb/4563522‑4563656                                                                    |                                                                    ggTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCattt                                                                  <  1:1344892/71‑1 (MQ=255) ggTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCattt                                                                  >  1:1649156/1‑71 (MQ=255) ggTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCattt                                                                  <  1:1715893/71‑1 (MQ=255)     ggCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACt                                                              <  1:1981586/71‑1 (MQ=255)       cagcagCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGt                                                            >  1:3791442/1‑71 (MQ=255)           agcTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGcc                                                        <  1:3705913/71‑1 (MQ=255)           agcTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGcc                                                        <  1:1668492/71‑1 (MQ=255)             cTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCa                                                      >  1:2149559/1‑71 (MQ=255)             cTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCa                                                      <  1:1931746/71‑1 (MQ=255)              ttGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAt                                                     >  1:400812/1‑71 (MQ=255)                gCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtgt                                                   >  1:31007/1‑71 (MQ=255)                 cTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtgtg                                                  <  1:2057646/71‑1 (MQ=255)                    aTCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGTCTGTAACGCCATTTTACTGTCGCCCATGTgg                                                 <  1:3169175/69‑1 (MQ=255)                    aTCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTgg                                                 <  1:138724/69‑1 (MQ=255)                    aTCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTgg                                                 <  1:253022/69‑1 (MQ=255)                          acaTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGt                                         <  1:395740/71‑1 (MQ=255)                             tCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCatttt                                                                 <  1:4056632/44‑1 (MQ=255)                              cGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGaa                                               >  1:2240864/1‑61 (MQ=255)                              cGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGt                                         >  1:1293025/1‑67 (MQ=255)                              cGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGt                                         >  1:962004/1‑67 (MQ=255)                              cGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGa                                       >  1:2668953/1‑69 (MQ=255)                              cGCAGTTATCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGt                                         >  1:1243967/1‑67 (MQ=255)                               gCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTAc                                                               >  1:148563/1‑44 (MQ=255)                               gCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGc                                    >  1:606482/1‑71 (MQ=255)                                                cATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGa                                       >  1:915574/1‑51 (MQ=255)                                                cATAAATAGGCTGTAACGCCATTTTACTGTCG‑‑CATGTGGaa                                               >  1:3430920/1‑41 (MQ=25)                                                   aaaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCaaa                >  1:3588778/1‑71 (MQ=255)                                                   aaaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCaaa                >  1:496748/1‑71 (MQ=255)                                                    aaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCaaa                >  1:535760/1‑70 (MQ=255)                                                    aaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAAc               >  1:853277/1‑71 (MQ=255)                                                      tAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCg             >  1:1084495/1‑71 (MQ=255)                                                      tAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCg             >  1:1978889/1‑71 (MQ=255)                                                       aGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAAc               >  1:2699141/1‑68 (MQ=255)                                                        ggCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGcc           <  1:3757258/71‑1 (MQ=255)                                                        ggCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGc            <  1:1946632/70‑1 (MQ=255)                                                               aCGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCgg    >  1:193691/1‑71 (MQ=255)                                                               aCGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCgg    >  1:742004/1‑71 (MQ=255)                                                                 gCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGa  <  1:1972632/71‑1 (MQ=255)                                                                    |                                                                    GGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGA  >  W3110S.gb/4563522‑4563656

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