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%	143.8 / NA	50	K207M (AAG→ATG)	iadA	isoaspartyl dipeptidase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (24/26);  total (24/26)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGG  >  W3110S.gb/4563522‑4563654                                                                    |                                                                  ggTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCattt                                                                <  1:2237072/71‑1 (MQ=255)  gTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCatttt                                                               <  1:3092359/71‑1 (MQ=255)            gcTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGccc                                                     <  1:3665425/71‑1 (MQ=255)             cTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCa                                                    >  1:3248203/1‑71 (MQ=255)             cTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCa                                                    <  1:1991296/71‑1 (MQ=255)              ttGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAt                                                   >  1:2966524/1‑71 (MQ=255)               tGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCg                                                        >  1:360091/1‑65 (MQ=255)                 cTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtgtg                                                >  1:590740/1‑71 (MQ=255)                    aTCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTgg                                               <  1:3673928/69‑1 (MQ=255)                    aTCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTgg                                               <  1:2295739/69‑1 (MQ=255)                    aTCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGaa                                             <  1:3681597/71‑1 (MQ=255)                    aTCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACCCCATTTTACTGTCGCCCATGTgg                                               <  1:1868945/69‑1 (MQ=255)                     tCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTa                                                              >  1:544251/1‑53 (MQ=255)                      cGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtg                                                  >  1:1526131/1‑64 (MQ=255)                      cGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAaca                                           >  1:1401070/1‑71 (MQ=255)                       ggCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAaca                                           <  1:1623598/70‑1 (MQ=255)                         cacaTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGaa                                             >  1:1686776/1‑66 (MQ=255)                         cacaTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACgg                                        <  1:2281627/71‑1 (MQ=255)                          acaTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGt                                       <  1:2096476/71‑1 (MQ=255)                          acaTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGt                                       <  1:3318892/71‑1 (MQ=255)                                  gTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTgg                                               <  1:605382/55‑1 (MQ=255)                                   ttttCCAGCAGGTCATAAATAGGCTGTAACGCCatttt                                                               <  1:1964384/38‑1 (MQ=255)                                     ttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGtt                            >  1:2784697/1‑71 (MQ=255)                                     ttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGtt                            >  1:878999/1‑71 (MQ=255)                                     ttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGtt                            >  1:1938486/1‑71 (MQ=255)                                      tCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGtt                            <  1:1356715/70‑1 (MQ=255)                                      tCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGtt                            <  1:2958757/70‑1 (MQ=255)                                       ccAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTa                          >  1:540814/1‑71 (MQ=255)                                       ccAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTa                          >  1:1620821/1‑71 (MQ=255)                                       ccAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTa                          >  1:2285869/1‑71 (MQ=255)                                            agGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACAAGGTGACGCCAGGTTTAccgcc                     <  1:2963345/71‑1 (MQ=255)                                             ggTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTgg                                               <  1:1535395/44‑1 (MQ=255)                                               tCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtgt                                                 >  1:961366/1‑40 (MQ=255)                                               tCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtgt                                                 >  1:1917171/1‑40 (MQ=255)                                               tCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTAccgcc                     <  1:1616821/68‑1 (MQ=255)                                               tCATAAATAGGCTGTAAAGCCATTTTACTGTCGCCCAtgt                                                 >  1:3362074/1‑40 (MQ=38)                                                 aTAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGa                   >  1:2144928/1‑68 (MQ=255)                                                  tAAATAGGCTGTCACGCCATTTTACTGTCGCCCAt                                                   <  1:1952522/35‑1 (MQ=37)                                                  tAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTAccgcc                     <  1:768616/65‑1 (MQ=255)                                                  tAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCaa               <  1:1364416/71‑1 (MQ=255)                                                     aTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAAcc            <  1:21577/71‑1 (MQ=255)                                                      tAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGt                                       <  1:2087364/43‑1 (MQ=255)                                                      tAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCg           >  1:816557/1‑71 (MQ=255)                                                       aGGCTGTAAGGCCATTTTACTGTCGCCCATGTGGAACACGGTGa                                     >  1:3381317/1‑44 (MQ=255)                                                        ggCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGcc         <  1:2053508/71‑1 (MQ=255)                                                        ggCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGcc         <  1:161214/71‑1 (MQ=255)                                                         gCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGc          <  1:3433877/69‑1 (MQ=255)                                                             tAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAAcgc    <  1:2891425/71‑1 (MQ=255)                                                               aCGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGc                                  >  1:425807/1‑39 (MQ=255)                                                               aCGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAg                               >  1:1615780/1‑42 (MQ=255)                                                               aCGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCgg  >  1:980066/1‑71 (MQ=255)                                                               aCGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCgg  <  1:3607419/71‑1 (MQ=255)                                                               aCGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCgg  >  1:2763347/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: ‑