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

GGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGC  >  W3110S.gb/4563522‑4563652                                                                    |                                                                ggTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCattt                                                              <  1:2490577/71‑1 (MQ=255) ggTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCattt                                                              <  1:1572656/71‑1 (MQ=255)        agcagcTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTc                                                       >  1:630425/1‑71 (MQ=255)           agcTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGcc                                                    >  1:246753/1‑71 (MQ=255)                gCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtg                                                <  1:3160959/70‑1 (MQ=255)                gCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtg                                                <  1:1922332/70‑1 (MQ=255)                gCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAt                                                 <  1:2672949/69‑1 (MQ=255)                   gATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTgg                                             >  1:2901671/1‑70 (MQ=255)                    aTCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTa                                                            <  1:2848447/54‑1 (MQ=255)                       ggCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAaca                                         <  1:903063/70‑1 (MQ=255)                          acaTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAt                                                 >  1:2222173/1‑59 (MQ=255)                           caTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGt                                     >  1:1472496/1‑70 (MQ=255)                           caTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTg                                    <  1:3104079/71‑1 (MQ=255)                            aTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCattt                                                              <  1:1587147/44‑1 (MQ=255)                            aTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGa                                   >  1:657211/1‑71 (MQ=255)                             tCGCAGTTTTCCAGCATGTCATAAATAGGCTGTAACGCCatttt                                                             <  1:2753992/44‑1 (MQ=38)                             tCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCg                                                      <  1:1227439/51‑1 (MQ=255)                               gCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTa                                                            >  1:1228040/1‑43 (MQ=255)                               gCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACt                                                          >  1:390232/1‑45 (MQ=255)                                 aGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCa                              <  1:2014612/71‑1 (MQ=255)                                 aGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCa                              <  1:2762726/71‑1 (MQ=255)                                   ttttCCAGCAGGTCATAAATAGGCTGTAACGCCatttt                                                             <  1:2086026/38‑1 (MQ=255)                                   ttttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAacac                                        <  1:2808673/59‑1 (MQ=255)                                    tttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAac                                          >  1:2325003/1‑56 (MQ=255)                                     ttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTAc                                                           >  1:75042/1‑38 (MQ=255)                                     ttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGt                                     >  1:524797/1‑60 (MQ=255)                                     ttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGtt                          >  1:2227126/1‑71 (MQ=255)                                      tCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGtt                          <  1:555436/70‑1 (MQ=255)                                       ccAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGttt                         <  1:2292794/70‑1 (MQ=255)                                       ccAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTa                        >  1:516292/1‑71 (MQ=255)                                          gcagGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTAccg                     <  1:2321562/71‑1 (MQ=255)                                              gTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGa                 <  1:2511927/71‑1 (MQ=255)                                               tCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAaca                                         <  1:2174008/46‑1 (MQ=255)                                               tCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAaca                                         <  1:1261517/46‑1 (MQ=255)                                               tCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAg                <  1:939674/71‑1 (MQ=255)                                                cATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtg                                                >  1:1275843/1‑38 (MQ=255)                                                cATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAg                >  1:3099674/1‑70 (MQ=255)                                                 aTAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGcc                               <  1:2234274/54‑1 (MQ=255)                                                 aTAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCa              <  1:2191511/71‑1 (MQ=255)                                                  tAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCaa             <  1:1569225/71‑1 (MQ=255)                                                  tAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCaa             <  1:1399458/71‑1 (MQ=255)                                                    aaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGaa                                           >  1:1855675/1‑39 (MQ=255)                                                       aGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAg                             <  1:1867573/50‑1 (MQ=255)                                                        ggCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGcc       <  1:1028696/71‑1 (MQ=255)                                                        ggCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGcc       <  1:1027106/71‑1 (MQ=255)                                                         gCTGTAACGCCATTTTACTGGCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGcc       <  1:3109151/70‑1 (MQ=255)                                                          cTGTAACGCCATTTTACTGTCGCCCATGTGGAaca                                         >  1:548594/1‑35 (MQ=255)                                                             tAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCaaa            <  1:2168160/61‑1 (MQ=255)                                                             tAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAAcgc  >  1:2456969/1‑71 (MQ=255)                                                              aaCGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAAc           >  1:1837629/1‑61 (MQ=255)                                                               aCGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAAcc          >  1:2335200/1‑61 (MQ=255)                                                                    |                                                                GGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGC  >  W3110S.gb/4563522‑4563652

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