Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	406,920	C→T	I172I (ATC→ATT)	aroL →	shikimate kinase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	406,920	0	C	T	100.0%	146.0 / NA	42	I172I (ATC→ATT)	aroL	shikimate kinase 2
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (19/23);  total (19/23)

CTGAGCGAAGAAGTTCAGGAAGTGCTGGAAGAACGCGATGCGCTATATCGCGAAGTTGCGCATATTATCATCGACGCAACAAACGAACCCAGCCAGGTGATTTCTGAAATTCGCAGCGCCCTGGCACAGACGATCAATTGTTGATTTTCGAGCGCCTATACTTAACGTTCATCCCGTGAAATAAGGAAGAACGATGCCAACGAAACCGCCTTATCCTCGTGAAGCATATATAGTGACGATTGAAAAAGGAAAGCCAGGACAGACGGTAA  >  NC_000913/406786‑407054                                                                                                                                       |                                                                                                                                       cTGAGCGAAGAAGTTCAGGAAGTGCTGGAAGAACGCGATGCGCTATATCGCGAAGTTGCGCATATTATCATCGACGCAACAAACGAACCCAGCCAGGTGATTTCTGAAATTCGCAGCGCCCTGGCACAGACGATTAAttg                                                                                                                                   >  1:156603/1‑140 (MQ=255)             gTTCAGGAAGTGCTGGAAGAACGCGATGCGCTATATCGCGAAGTTGCGCATATTATCATCGACGCAACAAACGAACCCAGCCAGGTGATTTCTGAAATTCGCAGCGCCCTGGCACAGACGATTAATTGTTGATTTTCGAg                                                                                                                       <  2:419812/140‑1 (MQ=255)             gTTCAGGAAGTGCTGGAAGAACGCGATGCGCTATATCGCGAAGTTGCGCATATTATCATCGACGCAACAAACGAACCCAGCCAGGTGATTTCTGAAATTCGCAGCGCCCTGGCACAGACGATTAATTGTTGATTTTCGAg                                                                                                                       <  1:467891/140‑1 (MQ=255)                cAGGAAGTGCTGGAAGAACGCGATGCGCTATATCGCGAAGTTGCGCATATTATCATCGACGCAACAAACGAACCCAGCCAGGTGATTTCTGAAATTCGCAGCGCCCTGGCACAGACGATTAATTGTTGATTTTCGAgcgc                                                                                                                    >  2:14109/1‑140 (MQ=255)                cAGGAAGTGCTGGAAGAACGCGATGCGCTATATCGCGAAGTTGCGCATATTATCATCGACGCAACAAACGAACCCAGCCAGGTGATTTCTGAAATTCGCAGCGCCCTGGCACAGACGATTAATTGTTGATTTTCGAgcgc                                                                                                                    <  2:327868/140‑1 (MQ=255)                cAGGAAGTGCTGGAAGAACGCGATGCGCTATATCGCGAAGTTGCGCATATTATCATCGACGCAACAAACGAACCCAGCCAGGTGATTTCTGAAATTCGCAGCGCCCTGGCACAGACGATTAATTGTTGATTTTCGAgcgc                                                                                                                    <  1:545376/140‑1 (MQ=255)                 aGGAAGTGCTGGAAGAACGCGATGCGCTATATCGCGAAGTTGCGCATATTATCATCGACGCAACAAACGAACCCAGCCAGGTGATTTCTGAAATTCGCAGCGCCCTGGCACAGACGATTAATTGTTGATTTTCGAGCGcc                                                                                                                   <  1:528290/140‑1 (MQ=255)                   gAAGTGCTGGAAGAACGCGATGCGCTATATCGCGAAGTTGCGCATATTATCATCGACGCAACAAACGAACCCAGCCAGGTGATTTCTGAAATTCGCAGCGCCCTGGCACAGACGATTAATTGTTGATTTTCGAGCGCCta                                                                                                                 <  1:95664/140‑1 (MQ=255)                       tGCTGGAAGAACGCGATGCGCTATATCGCGAAGTTGCGCATATTATCATCGACGCAACAAACGAACCCAGCCAGGTGATTTCTGAAATTCGCAGCGCCCTGGCACAGACGATTAATTGTTGATTTTCGAGCGCCTATACt                                                                                                             >  2:268952/1‑140 (MQ=255)                          tGGAAGAACGCGATGCGCTATATCGCGAAGTTGCGCATATTATCATCGACGCAACAAACGAACCCAGCCAGGTGATTTCTGAAATTCGCAGCGCCCTGGCACAGACGATTAATTGTTGATTTTCGAGCGCCTATACTTaa                                                                                                          <  1:268952/140‑1 (MQ=255)                                 aCGCGATGCGCTATATCGCGAAGTTGCGCATATTATCATCGACGCAACAAACGAACCCAGCCAGGTGATTTCTGAAATTCGCAGCGCCCTGGCACAGACGATTAATTGTTGATTTTCGAGCGCCtat                                                                                                                <  1:343157/127‑1 (MQ=255)                                 aCGCGATGCGCTATATCGCGAAGTTGCGCATATTATCATCGACGCAACAAACGAACCCAGCCAGGTGATTTCTGAAATTCGCAGCGCCCTGGCACAGACGATTAATTGTTGATTTTCGAGCGCCtat                                                                                                                >  2:343157/1‑127 (MQ=255)                                   gcgATGCGCTATATCGCGAAGTTGCGCATATTATCATCGACGCAACAAACGAACCCAGCCAGGTGATTTCTGAAATTCGCAGCGCCCTGGCACAGACGATTAATTGTTGATTTTCGAGCGCCTATACTTAACGTTCATcc                                                                                                 >  1:500956/1‑140 (MQ=255)                                     gATGCGCTATATCGCGAAGTTGCGCATATTATCATCGACGCAACAAACGAACCCAGCCAGGTGATTTCTGAAATTCGCAGCGCCCTGGCACAGACGATTAATTGTTGATTTTCGAGCGCCTATACTTAACGTTCATCCCg                                                                                               <  1:550451/140‑1 (MQ=255)                                      aTGCGCTATATCGCGAAGTTGCGCATATTATCATCGACGCAACAAACGAACCCAGCCAGGTGATTTCTGAAATTCGCAGCGCCCTGGCACAGACGATTAATTGTTGATTTTCGAGCGCCTATACTTAACGTTCATCCCGt                                                                                              >  2:375394/1‑140 (MQ=255)                                            tataTCGCGAAGTTGCGCATATTATCATCGACGCAACAAACGAACCCAGCCAGGTGATTTCTGAAATTCGCAGCGCCCTGGCACAGACGATTAATTGTTGATTTTCGAGCGCCTATACTTAACGTTCATCCCGTGAAATa                                                                                        >  1:758956/1‑140 (MQ=255)                                                 cgcgAAGTTGCGCATATTATCATCGACGCAACAAACGAACCCAGCCAGGTGATTTCTGAAATTCGCAGCGCCCTGGCACAGACGATTAATTGTTGATTTTCGAGCGCCTATACTTAACGTTCATCCCGTGAAATAAGgaa                                                                                   >  2:755572/1‑140 (MQ=255)                                                    gAAGTTGCGCATATTATCATCGACGCAACAAACGAACCCAGCCAGGTGATTTCTGAAATTCGCAGCGCCCTGGCACAGACGATTAATTGTTGATTTTCGAGCGCCTATACTTAACGTTCATCCCGTGAAATAAGgaagaa                                                                                <  1:14109/140‑1 (MQ=255)                                                        ttGCGCATATTATCATCGACGCAACAAACGAACCCAGTCAGGTGATTTCTGAAATTCGCAGCGCCCTGGCACAGACGATTAATTGTTGATTTTCGAGCGCCTATACTTAACGTTCATCCCGTGAAATAAGGAAGAACGAt                                                                            >  1:353844/1‑140 (MQ=255)                                                          gcgcATATTATCATCGACGCAACAAACGAACCCAGCCAGGTGATTTCTGAAATTCGCAGCGCCCTGGCACAGACGATTAATTGTTGATTTTCGAGCGCCTATACTTAACGTTCATCCCGTGAAATAAGGAAGAACGATGc                                                                          >  2:57218/1‑140 (MQ=255)                                                          gcgcATATTATCATCGACGCAACAAACGAACCCAGCCAGGTGATTTCTGAAATTCGCAGCGCCCTGGCACAGACGATTAATTGTTGATTTTCGAGCGCCTATACTTAACGTTCATCCCGTGAAATAAGGAAGAACGATGc                                                                          >  2:467076/1‑140 (MQ=255)                                                                     catcGACGCAACAAACGAACCCAGCCAGGTGATTTCTGAAATTCGCAGCGCCCTGGCACAGACGATTAATTGTTGATTTTCGAGCGCCTATACTTAACGTTCATCCCGTGAAATAAGGAAGAACGATGCCAACGaaa                                                                  <  1:548550/137‑1 (MQ=255)                                                                     catcGACGCAACAAACGAACCCAGCCAGGTGATTTCTGAAATTCGCAGCGCCCTGGCACAGACGATTAATTGTTGATTTTCGAGCGCCTATACTTAACGTTCATCCCGTGAAATAAGGAAGAACGATGCCAACGaaa                                                                  >  2:548550/1‑137 (MQ=255)                                                                          aCGCAACAAACGAACCCAGCCAGGTGATTTCTGAAATTCGCAGCGCCCTGGCACAGACGATTAATTGTTGATTTTCGAGCGCCTATACTTAACGTTCATCCCGTGAAATAAGGAAGAACGATGCCAACGAAACCGCCTTa                                                          >  2:15775/1‑140 (MQ=255)                                                                          aCGCAACAAACGAACCCAGCCAGGTGATTTCTGAAATTCGCAGCGCCCTGGCACAGACGATTAATTGTTGATTTTCGAGCGCCTATACTTAACGTTCATCCCGTGAAATAAGGAAGAACGATGCCAACGAAACCGCCTTa                                                          >  2:365274/1‑140 (MQ=255)                                                                              aacaaACGAACCCAGCCAGGTGATTTCTGAAATTCGCAGCGCCCTGGCACAGACGATTAATTGTTGATTTTCGAGCGCCTATACTTAACGTTCATCCCGTGAAATAAGGAAGAACGATGCCAACGAAACCGCCTTATCCt                                                      >  1:491406/1‑140 (MQ=255)                                                                               acaaACGAACCCAGCCAGGTGATTTCTGAAATTCGCAGCGCCCTGGCACAGACGATTAATTGTTGATTTTCGAGCGCCTATACTTAACGTTCATCCCGTGAAATAAGGAAGAACGATGCCAACGAAACCGCCTTATCCTc                                                     <  2:500956/140‑1 (MQ=255)                                                                                 aaaCGAACCCAGCCAGGTGATTTCTGAAATTCGCAGCGCCCTGGCACAGACGATTAATTGTTGATTTTCGAGCGCCTATACTTAACGTTCATCCCGTGAAATAAGGAAGAACGATGCCAACGAAACCGCCTTATCCTCGt                                                   <  1:28402/140‑1 (MQ=255)                                                                                 aaaCGAACCCAGCCAGGTGATTTCTGAAATTCGCAGCGCCCTGGCACAGACGATTAATTGTTGATTTTCGAGCGCCTATACTTAACGTTCATCCCGTGAAATAAGGAAGAACGATGCCAACGAAACCGCCTTATCCTCGt                                                   <  2:381333/140‑1 (MQ=255)                                                                                 aaaCGAACCCAGCCAGGTGATTTCTGAAATTCGCAGCGCCCTGGCACAGACGATTAATTGTTGATTTTCGAGCGCCTATACTTAACGTTCATCCCGTGAAATAAGGAAGAACGATGCCAACGAAACCGCCTTATCCTCGt                                                   <  1:15775/140‑1 (MQ=255)                                                                                   aCGAACCCAGTCAGGTGATTTCTGAAATTCGCAGCGCCCTGGCACAGACGATTAATTGTTGATTTTCGAGCGCCTATACTTAACGTTCATCCCGTGAAATAAGGAAGAACGATGCCAACGAAACCGCCTTATCCTCGTGa                                                 <  2:353844/140‑1 (MQ=255)                                                                                   aCGAACCCAGCCAGGTGATTTCTGAAATTCGCAGCGCCCTGGCACAGACGATTAATTGTTGATTTTCGAGCGCCTATACTTAACGTTCATCCCGTGAAATAAGGAAGAACGATGCCAACGAAACCGCCTTATCCTCGTGa                                                 <  2:591312/140‑1 (MQ=255)                                                                                   aCGAACCCAGCCAGGTGATTTCTGAAATTCGCAGCGCCCTGGCACAGACGATTAATTGTTGATTTTCGAGCGCCTATACTTAACGTTCATCCCGTGAAATAAGGAAGAACGATGCCAACGAAACCGCCTTATCCTCGTGa                                                 <  1:474718/140‑1 (MQ=255)                                                                                   aCGAACCCAGCCAGGTGATTTCTGAAATTCGCAGCGCCCTGGCACAGACGATTAATTGTTGATTTTCGAGCGCCTATACTTAACGTTCATCCCGTGAAATAAGGAAGAACGATGCCAACGAAACCGCCTTATCCTCGTGa                                                 <  2:491406/140‑1 (MQ=255)                                                                                      aaCCCAGCCAGGTGATTTCTGAAATTCGCAGCGCCCTGGCACAGACGATTAATTGTTGATTTTCGAGCGCCTATACTTAACGTTCATCCCGTGAAATAAGGAAGAACGATGCCAACGAAACCGCCTTATCCTCGTGAAGc                                              <  2:775341/140‑1 (MQ=255)                                                                                      aaCCCAGCCAGGTGATTTCTGAAATTCGCAGCGCCCTGGCACAGACGATTAATTGTTGATTTTCGAGCGCCTATACTTAACGTTCATCCCGTGAAATAAGGAAGAACGATGCCAACGAAACCGCCTTATCCTCGTGAAGc                                              <  1:755572/140‑1 (MQ=255)                                                                                      aaCCCAGCCAGGTGATTTCTGAAATTCGCAGCGCCCTGGCACAGACGATTAATTGTTGATTTTCGAGCGCCTATACTTAACGTTCATCCCGTGAAATAAGGAAGAACGATGCCAACGAAACCGCCTTATCCTCGTGAAGc                                              <  1:57218/140‑1 (MQ=255)                                                                                                ggTGATTTCTGAAATTCGCAGCGCCCTGGCACAGACGATTAATTGTTGATTTTCGAGCGCCTATACTTAACGTTCATCCCGTGAAATAAGGAAGAACGATGCCAACGAAACCGCCTTATCCTCGTGaa                                                >  2:359435/1‑128 (MQ=255)                                                                                                ggTGATTTCTGAAATTCGCAGCGCCCTGGCACAGACGATTAATTGTTGATTTTCGAGCGCCTATACTTAACGTTCATCCCGTGAAATAAGGAAGAACGATGCCAACGAAACCGCCTTATCCTCGTGaa                                                <  1:359435/128‑1 (MQ=255)                                                                                                                  cAGCGCCCTGGCACAGACGATTAATTGTTGATTTTCTAGCGCCTATACTTAACGTTCATCCCGTGAAATAAGGAAGAACGATGCCAACGAAACCGCCTTATCCTCGTGAAGCATATATAGTGACGATTGAAAAAGGAAAg                  >  1:69912/1‑140 (MQ=255)                                                                                                                       cccTGGCACAGACGATTAATTGTTGATTTTCGAGCGCCTATACTTAACGTTCATCCCGTGAAATAAGGAAGAACGATGCCAACGAAACCGCCTTATCCTCGTGAAGCATATATAGTGACGATTGAAAAAGGAAAGCCAgg             >  2:208970/1‑140 (MQ=255)                                                                                                                                  aCGATTAATTGTTGATTTTCGAGCGCCTATACTTAACGTTCATCCCGTGAAATAAGGAAGAACGATGCCAACGAAACCGCCTTATCCTCGTGAAGCATATATAGTGACGATTGAAAAAGGAAAGCCAGGACAGACGGTaa  >  2:353080/1‑140 (MQ=255)                                                                                                                                       |                                                                                                                                       CTGAGCGAAGAAGTTCAGGAAGTGCTGGAAGAACGCGATGCGCTATATCGCGAAGTTGCGCATATTATCATCGACGCAACAAACGAACCCAGCCAGGTGATTTCTGAAATTCGCAGCGCCCTGGCACAGACGATCAATTGTTGATTTTCGAGCGCCTATACTTAACGTTCATCCCGTGAAATAAGGAAGAACGATGCCAACGAAACCGCCTTATCCTCGTGAAGCATATATAGTGACGATTGAAAAAGGAAAGCCAGGACAGACGGTAA  >  NC_000913/406786‑407054

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑