Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,054,066	G→A	E142K (GAA→AAA)	shiA →	shikimate:H(+) symporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,054,066	0	G	A	100.0%	212.8 / NA	60	E142K (GAA→AAA)	shiA	shikimate:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (25/35);  total (25/35)

ATGCTGACCGTCTGGATGATGGGCATCGCGACAGCCTTGATTGGTATTCTTCCTTCATTCTCGACCATTGGGTGGTGGGCACCTATTTTGCTGGTGACACTGCGTGCCATTCAGGGATTTGCAGTCGGCGGCGAATGGGGAGGCGCGGCGTTGCTTTCCGTTGAAAGTGCACCGAAAAATAAAAAAGCCTTTTACAGTAGCGGTGTACAAGTTGGCTACGGTGTAGGTTTACTGCTTTCAACCGGACTGGTTTCATTGATCAGT  >  NC_000913/2053934‑2054197                                                                                                                                     |                                                                                                                                    aTGCTGACCGTCTGGATGATGGGCATCGCGACAGCCTTGATTGGTATTCTTCCTTCATTCTCGACCATTGGGTGGTGGGCACCTATTTTGCTGGTGACACTGCGTGCCATTCAGGGATTTGCAGTCGGCGGCAAATgggg                                                                                                                              <  1:738266/140‑1 (MQ=255) aTGCTGACCGTCTGGATGATGGGCATCGCGACAGCCTTGATTGGTATTCTTCCTTCATTCTCGACCATTGGGTGGTGGGCACCTATTTTGCTGGTGACACTGCGTGCCATTCAGGGATTTGCAGTCGGCGGCAAATgggg                                                                                                                              <  2:917639/140‑1 (MQ=255)      gACCGTCTGGATGATGGGCATCGCGACAGCCTTGATTGGTATTCTTCCTTCATTCTCGACCATTGGGTGGTGGGCACCTATTTTGCTGGTGACACTGCGTGCCATTCAGGGATTTGCAGTCGGCGGCAAATGGGGAGgcg                                                                                                                         <  2:320464/140‑1 (MQ=255)      gACCGTCTGGATGATGGGCATCGCGACAGCCTTGATTGGTATTCTTCCTTCATTCTCGACCATTGGGTGGTGGGCACCTATTTTGCTGGTGACACTGCGTGCCATTCAGGGATTTGCAGTCGGCGGCAAATGGGGAGgcg                                                                                                                         <  2:648901/140‑1 (MQ=255)        ccGTCTGGATGATGGGCATCGCGACAGCCTTGATTGGTATTCTTCCTTCATTCTCGACCATTGGGTGGTGGGCACCTATTTTGCTGGTGACACTGCGTGCCATTCAGGGATTTGCAGTCGGCGGCAAATGGGGAGGcgcg                                                                                                                       >  1:815423/1‑140 (MQ=255)                         aTCGCGACAGCCTTGATTGGTATTCTTCCTTCATTCTCGACCATTGGGTGGTGGGCACCTATTTTGCTGGTGACACTGCGTGCCATTCAGGGATTTGCAGTCGGCGGCAAATGGGGAGGCGCGGCGTTGCTTTCCGTTGa                                                                                                      >  1:827054/1‑140 (MQ=255)                         aTCGCGACAGCCTTGATTGGTATTCTTCCTTCATTCTCGACCATTGGGTGGTGGGCACCTATTTTGCTGGTGACACTGCGTGCCATTCAGGGATTTGCAGTCGGCGGCAAATGGGGAGGCGCGGCGTTGCTTTCCGTTGa                                                                                                      >  1:823005/1‑140 (MQ=255)                         aTCGCGACAGCCTTGATTGGTATTCTTCCTTCATTCTCGACCATTGGGTGGTGGGCACCTATTTTGCTGGTGACACTGCGTGCCATTCAGGGATTTGCAGTCGGCGGCAAATGGGGAGGCGCGGCGTTGCTTTCCGTTGa                                                                                                      <  1:13596/140‑1 (MQ=255)                                   ccTTGATTGGTATTCTTCCTTCATTCTCGACCATTGGGTGGTGGGCACCTATTTTGCTGGTGACACTGCGTGCCATTCAGGGATTTGCAGTCGGCGGCAAATGGGGAGGCGCGGCGTTGCTTTCCGTTGaaa                                                                                                    <  2:439163/132‑1 (MQ=255)                                   ccTTGATTGGTATTCTTCCTTCATTCTCGACCATTGGGTGGTGGGCACCTATTTTGCTGGTGACACTGCGTGCCATTCAGGGATTTGCAGTCGGCGGCAAATGGGGAGGCGCGGCGTTGCTTTCCGTTGaaa                                                                                                    >  1:439163/1‑132 (MQ=255)                                              aTTCTTCCTTCATTCTCGACCATTGGGTGGTGGGCACCTATTTTGCTGGTGACACTGCGTGCCATTCAGGGATTTGCAGTCGGCGGCAAATGGGGAGGCGCGGCGTTGCTTTCCGTTGAAAGTGCACCGAAAAATaaaaa                                                                                 <  2:233892/140‑1 (MQ=255)                                              aTTCTTCCTTCATTCTCGACCATTGGGTGGTGGGCACCTATTTTGCTGGTGACACTGCGTGCCATTCAGGGATTTGCAGTCGGCGGCAAATGGGGAGGCGCGGCGTTGCTTTCCGTTGAAAGTGCACCGAAAAATaaaaa                                                                                 >  1:303417/1‑140 (MQ=255)                                                  ttccttcATTCTCGACCATTGGGTGGTGGGCACCTATTTTGCTGGTGACACTGCGTGCCATTCAGGGATTTGCAGTCGGCGGCAAATGGGGAGGCGCGGCGTTGCTTTCCGTTGAAAGTGCACCGAAAAATAAAAAAGcc                                                                             >  2:273602/1‑140 (MQ=255)                                                  ttccttcATTCTCGACCATTGGGTGGTGGGCACCTATTTTGCTGGTGACACTGCGTGCCATTCAGGGATTTGCAGTCGGCGGCAAATGGGGAGGCGCGGCGTTGCTTTCCGTTGAAAGTGCACCGAAAAATAAAAAAGcc                                                                             <  2:815423/140‑1 (MQ=255)                                                   tccttcATTCTCGACCATTGGGTGGTGGGCACCTATTTTGCTGGTGACACTGCGTGCCATTCAGGGATTTGCAGTCGGCGGCAAATGGGGAGGCGCGGCGTTGCTTTCCGTTGAAAGTGCACCGAAAAATAAAAAAGCCt                                                                            <  2:303417/140‑1 (MQ=255)                                                       tcATTCTCGACCATTGGGTGGTGGGCACCTATTTTGCTGGTGACACTGCGTGCCATTCAGGGATTTGCAGTCGGCGGCAAATGGGGAGGCGCGGCGTTGCTTTCCGTTGAAAGTGCACCGAAAAATAAAAAAGCCTTTTa                                                                        >  1:229240/1‑140 (MQ=255)                                                          ttCTCGACCATTGGGTGGTGGGCACCTATTTTGCTGGTGACACTGCGTGCCATTCAGGGATTTGCAGTCGGCGGCAAATGGGGAGGCGCGGCGTTGCTTTCCGTTGAAAGTGCACCGAAAAATAAAAAAGCCTTTTACAg                                                                     >  2:697059/1‑140 (MQ=255)                                                          ttCTCGACCATTGGGTGGTGGGCACCTATTTTGCTGGTGACACTGCGTGCCATTCAGGGATTTGCAGTCGGCGGCAAATGGGGAGGCGCGGCGTTGCTTTCCGTTGAAAGTGCACCGAAAAATAAAAAAGCCTTTTACAg                                                                     <  2:108659/140‑1 (MQ=255)                                                          ttCTCGACCATTGGGTGGTGGGCACCTATTTTGCTGGTGACACTGCGTGCCATTCAGGGATTTGCAGTCGGCGGCAAATGGGGAGGCGCGGCGTTGCTTTCCGTTGAAAGTGCACCGAAAAATAAAAAAGCCTTTTACAg                                                                     <  1:292818/140‑1 (MQ=255)                                                           tctcGACCATTGGGTGGTGGGCACCTATTTTGCTGGTGACACTGCGTGCCATTCAGGGATTTGCAGTCGGCGGCAAATGGGGAGGCGCGGCGTTGCTTTCCGTTGAAAGTGCACCGAAAAATaaaa                                                                                  <  1:672191/126‑1 (MQ=255)                                                           tctcGACCATTGGGTGGTGGGCACCTATTTTGCTGGTGACACTGCGTGCCATTCAGGGATTTGCAGTCGGCGGCAAATGGGGAGGCGCGGCGTTGCTTTCCGTTGAAAGTGCACCGAAAAATaaaa                                                                                  >  2:672191/1‑126 (MQ=255)                                                            ctcGACCATTGGGTGGTGGGCACCTATTTTGCTGGTGACACTGCGTGCCATTCAGGGATTTGCAGTCGGCGGCAAATGGGGAGGCGCGGCGTTGCTTTCCGTTGAAAGTGCACCGAAAAATAAAAAAGCCTTTTACAGTa                                                                   >  2:234241/1‑140 (MQ=255)                                                             tcGACCATTGGGTGGTGGGCACCTATTTTGCTGGTGACACTGCGTGCCATTCAGGGATTTGCAGTCGGCGGCAAATGGGGAGGCGCGGCGTTGCTTTCCGTTGAAAGTGCACCGAAAAATAAAAAAGCCTTTTACAGTAg                                                                  <  1:697059/140‑1 (MQ=255)                                                             tcGACCATTGGGTGGTGGGCACCTATTTTGCTGGTGACACTGCGTGCCATTCAGGGATTTGCAGTCGGCGGCAAATGGGGAGGCGCGGCGTTGCTTTCCGTTGAAAGTGCACCGAAAAATAAAAAAGCCTTTTACAGTAg                                                                  <  1:234241/140‑1 (MQ=255)                                                             tcGACCATTGGGTGGTGGGCACCTATTTTGCTGGTGACACTGCGTGCCATTCAGGGATTTGCAGTCGGCGGCAAATGGGGAGGCGCGGCGTTGCTTTCCGTTGAAAGTGCACCGAAAAATAAAAAAGCCTTTTACAGTAg                                                                  <  2:827054/140‑1 (MQ=255)                                                               gACCATTGGGTGGTGGGCACCTATTTTGCTGGTGACACTGCGTGCCATTCAGGGATTTGCAGTCGGCGGCAAATGGGGAGGCGCGGCGTTGCTTTCCGTTGAAAGTGCACCGAAAAATAAAAAAGCCTTTTACAGTAGCg                                                                <  1:121604/140‑1 (MQ=255)                                                               gACCATTGGGTGGTGGGCACCTATTTTGCTGGTGACACTGCGTGCCATTCAGGGATTTGCAGTCGGCGGCAAATGGGGAGGCGCGGCGTTGCTTTCCGTTGAAAGTGCACCGAAAAATAAAAAAGCCTTTTACAGTAGCg                                                                <  2:681804/140‑1 (MQ=255)                                                                    ttGGGTGGTGGGCACCTATTTTGCTGGTGACACTGCGTGCCATTCAGGGATTTGCAGTCGGCGGCAAATGGGGAGGCGCGGCGTTGCTTTCCGTTGAAAGTGCACCGAAAAATAAAAAAGCCTTTTACAGTAGCGGTGTa                                                           <  2:703974/140‑1 (MQ=255)                                                                     tGGGTGGTGGGCACCTATTTTGCTGGTGACACTGCGTGCCATTCAGGGATTTGCAGTCGGCGGCAAATGGGGAGGCGCGGCGTTGCTTTCCGTTGAAAGTGCACCGAAAAATAAAAAAGCCTTTTACAGTAGCGGTGTAc                                                          <  2:229240/140‑1 (MQ=255)                                                                       ggtggtGGGCACCTATTTTGCTGGTGACACTGCGTGCCATTCAGGGATTTGCAGTCGGCGGCAAATGGGGAGGCGCGGCGTTGCTTTCCGTTGAAAGTGCACCGAAAAATAAAAAAGCCTTTTACAGTAGCGgtg                                                             <  1:960168/135‑1 (MQ=255)                                                                       ggtggtGGGCACCTATTTTGCTGGTGACACTGCGTGCCATTCAGGGATTTGCAGTCGGCGGCAAATGGGGAGGCGCGGCGTTGCTTTCCGTTGAAAGTGCACCGAAAAATAAAAAAGCCTTTTACAGTAGCGgtg                                                             >  2:960168/1‑135 (MQ=255)                                                                           gtggGCACCTATTTTGCTGGTGACACTGCGTGCCATTCAGGGATTTGCAGTCGGCGGCAAATGGGGAGGCGCGGCGTTGCTTTCCGTTGAAAGTGCACCGAAAAATAAAAAAGCCTTTTACAGTAGCGgtg                                                             <  2:601573/131‑1 (MQ=255)                                                                           gtggGCACCTATTTTGCTGGTGACACTGCGTGCCATTCAGGGATTTGCAGTCGGCGGCAAATGGGGAGGCGCGGCGTTGCTTTCCGTTGAAAGTGCACCGAAAAATAAAAAAGCCTTTTACAGTAGCGgtg                                                             >  1:601573/1‑131 (MQ=255)                                                                           gtggGCACCTATTTTGCTGGTGACACTGCGTGCCATTCAGGGATTTGCAGTCGGCGGCAAATGGGGAGGCGCGGCGTTGCTTTCCGTTGAAAGTGCACCGAAAAATAAAAAAGCCTTTTACAGTAGCGGTGTACAAGTTg                                                    >  2:673182/1‑140 (MQ=255)                                                                           gtggGCACCTATTTTGCTGGTGACACTGCGTGCCATTCAGGGATTTGCAGTCGGCGGCAAATGGGGAGGCGCGGCGTTGCTTTCCGTTGAAAGTGCACCGAAAAATAAAAAAGCCTTTTACAGTAGCGGTGTACAAGTTg                                                    >  2:25878/1‑140 (MQ=255)                                                                             gggCACCTATTTTGCTGGTGACACTGCGTGCCATTCAGGGATTTGCAGTCGGCGGCAAATGGGGAGGCGCGGCGTTGCTTTCCGTTGAAAGTGCACCGAAAAATAAAAAAGCCTTTTACAGTAGCGGTGTACAAGTTGGc                                                  <  1:273602/140‑1 (MQ=255)                                                                                   cTATTTTGCTGGTGACACTGCGTGCCATTCAGGGATTTGCAGTCGGCGGCAAATGGGGAGGCGCGGCGTTGCTTTCCGTTGAAAGTGCACCGAAAAATAAAAAAGCCTTTTACAGTAGCGGTGTACaa                                                        >  1:972235/1‑128 (MQ=255)                                                                                   cTATTTTGCTGGTGACACTGCGTGCCATTCAGGGATTTGCAGTCGGCGGCAAATGGGGAGGCGCGGCGTTGCTTTCCGTTGAAAGTGCACCGAAAAATAAAAAAGCCTTTTACAGTAGCGGTGTACaa                                                        <  2:972235/128‑1 (MQ=255)                                                                                   cTATTTTGCTGGTGACACTGCGTGCCATTCAGGGATTTGCAGTCGGCGGCAAATGGGGAGGCGCGGCGTTGCTTTCCGTTGAAAGTGCACCGAAAAATAAAAAAGCCTTTTACAGTAGCGGTGTACAAGTTGGc                                                  <  1:405235/134‑1 (MQ=255)                                                                                   cTATTTTGCTGGTGACACTGCGTGCCATTCAGGGATTTGCAGTCGGCGGCAAATGGGGAGGCGCGGCGTTGCTTTCCGTTGAAAGTGCACCGAAAAATAAAAAAGCCTTTTACAGTAGCGGTGTACAAGTTGGc                                                  >  2:405235/1‑134 (MQ=255)                                                                                      ttttGCTGGTGACACTGCGTGCCATTCAGGGATTTGCAGTCGGCGGCAAATGGGGAGGCGCGGCGTTGCTTTCCGTTGAAAGTGCACCGAAAAATAAAAAAGCCTTTTACAGTAGCGGTGTACAAGTTGGCTACGGTGTa                                         >  1:616518/1‑140 (MQ=255)                                                                                        ttGCTGGTGACACTGCGTGCCATTCAGGGATTTGCAGTCGGCGGCAAATGGGGAGGCGCGGCGTTGCTTTCCGTTGAAAGTGCACCGAAAAATAAAAAAGCCTTTTACAGTAGCGGTGTACAAGTTGGCTACGGTGTAgg                                       >  2:124923/1‑140 (MQ=255)                                                                                            tGGTGACACTGCGTGCCATTCAGGGATTTGCAGTCGGCGGCAAATGGGGAGGCGCGGCGTTGCTTTCCGTTGAAAGTGCACCGAAAAATAAAAAAGCCTTTTACAGTAGCGGTGTACAAGTTGGCTACGGTGTAGGTTTa                                   <  1:124923/140‑1 (MQ=255)                                                                                             ggTGACACTGCGTGCCATTCAGGGATTTGCAGTCGGCGGCAAATGGGGAGGCGCGGCGTTGCTTTCCGTTGAAAGTGCACCGAAAAATAAAAAAGCCTTTTACAGTAGCGGTGTACAAGTTGGCTACGGTGTAGGTTTAc                                  >  2:826370/1‑140 (MQ=255)                                                                                              gTGACACTGCGTGCCATTCAGGGATTTGCAGTCGGCGGCAAATGGGGAGGCGCGGCGTTGCTTTCCGTTGAAAGTGCACCGAAAAATAAAAAAGCCTTTTACAGTAGCGGTGTACAAGTTGGCTACGGTGTAGGTTTACt                                 >  2:205852/1‑140 (MQ=255)                                                                                                 acacTGCGTGCCATTCAGGGATTTGCAGTCGGCGGCAAATGGGGAGGCGCGGCGTTGCTTTCCGTTGAAAGTGCACCGAAAAATAAAAAAGCCTTTTACAGTAGCGGTGTACAAGTTGGCTACGGTGTAGGTTTACTGCt                              <  1:455698/140‑1 (MQ=255)                                                                                                        gTGCCATTCAGGGATTTGCAGTCGGCGGCAAATGGGGAGGCGCGGCGTTGCTTTCCGTTGAAAGTGCACCGAAAAATAAAAAAGCCTTTTACAGTAGCGGTGTACAAGTTGGCTACGGTGTAGGTTTACTGCTTTCAAcc                       >  2:268996/1‑140 (MQ=255)                                                                                                         tGCCATTCAGGGATTTGCAGTCGGCGGCAAATGGGGAGGCGCGGCGTTGCTTTCCGTTGAAAGTGCACCGAAAAATAAAAAAGCCTTTTACAGTAGCGGTGTACAAGTTGGCTACGGTGTAGGTTTACTGCTTTCAACCg                      >  1:222594/1‑140 (MQ=255)                                                                                                          gCCATTCAGGGATTTGCAGTCGGCGGCAAATGGGGAGGCGCGGCGTTGCTTTCCGTTGAAAGTGCACCGAAAAATAAAAAAGCCTTTTACAGTAGCGGTGTACAAGTTGGCTACGGTGTAGGTTTACTGCTTTCAACCgg                     <  2:616518/140‑1 (MQ=255)                                                                                                               tCAGGGATTTGCAGTCGGCGGCAAATGGGGAGGCGCGGCGTTGCTTTCCGTTGAAAGTGCACCGAAAAATAAAAAAGCCTTTTACAGTAGCGGTGTACAAGTTGGCTACGGTGTAGGTTTACTGCTTTCAACCGGACTgg                >  1:648847/1‑140 (MQ=255)                                                                                                                    gATTTGCAGTCGGCGGCAAATGGGGAGGCGCGGCGTTGCTTTCCGTTGAAAGTGCACCGAAAAATAAAAAAGCCTTTTACAGTAGCGGTGTACAAGTTGGCTACGGTGTAGGTTTACTGCTTTCAACCGGACTGGTTTCa           <  2:823005/140‑1 (MQ=255)                                                                                                                    gATTTGCAGTCGGCGGCAAATGGGGAGGCGCGGCGTTGCTTTCCGTTGAAAGTGCACCGAAAAATAAAAAAGCCTTTTACAGTAGCGGTGTACAAGTTGGCTACGGTGTAGGTTTACTGCTTTCAACCGGACTGGTTTCa           <  2:222594/140‑1 (MQ=255)                                                                                                                       ttGCAGTCGGCGGCAAATGGGGAGGCGCGGCGTTGCTTTCCGTTGAAAGTGCACCGAAAAATAAAAAAGCCTTTTACAGTAGCGGTGTACAAGTTGGCTACGGTGTAGGTTTACTGCTTTCAACCGGACTGGTTTCATTg        <  1:826370/140‑1 (MQ=255)                                                                                                                         gCAGTCGGCGGCAAATGGGGAGGCGCGGCGTTGCTTTCCGTTGAAAGTGCACCGAAAAATAAAAAAGCCTTTTACAGTAGCGGTGTACAAGTTGGCTACGGTGTAGGTTTACTGCTTTCAACCGGACTGGTTTCATTGAt      <  1:673182/140‑1 (MQ=255)                                                                                                                         gCAGTCGGCGGCAAATGGGGAGGCGCGGCGTTGCTTTCCGTTGAAAGTGCACCGAAAAATAAAAAAGCCTTTTACAGTAGCGGTGTACAAGTTGGCTACGGTGTAGGTTTACTGCTTTCAACCGGACTGGTTTCATTGAt      <  1:268996/140‑1 (MQ=255)                                                                                                                         gCAGTCGGCGGCAAATGGGGAGGCGCGGCGTTGCTTTCCGTTGAAAGTGCACCGAAAAATAAAAAAGCCTTTTACAGTAGCGGTGTACAAGTTGGCTACGGTGTAGGTTTACTGCTTTCAACCGGACTGGTTTCATTGAt      <  1:25878/140‑1 (MQ=255)                                                                                                                         gCAGTCGGCGGCAAATGGGGAGGCGCGGCGTTGCTTTCCGTTGAAAGTGCACCGAAAAATAAAAAAGCCTTTTACAGTAGCGGTGTACAAGTTGGCTACGGTGTAGGTTTACTGCTTTCAACCGGACTGGTTTCATTGAt      <  1:205852/140‑1 (MQ=255)                                                                                                                             tCGGCGGCAAATGGGGAGGCGCGGCGTTGCTTTCCGTTGAAAGTGCACCGAAAAATAAAAAAGCCTTTTACAGTAGCGGTGTACAAGTTGGCTACGGTGTAGGTTTACTGCTTTCAACCGGACTGGTTTCATTGATCAGt  >  2:240405/1‑140 (MQ=255)                                                                                                                              cggcggCAAATGGGGAGGCGCGGCGTTGCTTTCCGTTGAAAGTGCACCGAAAAATAAAAAAGCCTTTTACAGTAGCGGTGTACAAGTTGGCTACGGTGTAGGTTTACTGCTTTCAACCGGACTGGTTTCATTGAt      >  1:904726/1‑135 (MQ=255)                                                                                                                              cggcggCAAATGGGGAGGCGCGGCGTTGCTTTCCGTTGAAAGTGCACCGAAAAATAAAAAAGCCTTTTACAGTAGCGGTGTACAAGTTGGCTACGGTGTAGGTTTACTGCTTTCAACCGGACTGGTTTCATTGAt      <  2:904726/135‑1 (MQ=255)                                                                                                                                     |                                                                                                                                    ATGCTGACCGTCTGGATGATGGGCATCGCGACAGCCTTGATTGGTATTCTTCCTTCATTCTCGACCATTGGGTGGTGGGCACCTATTTTGCTGGTGACACTGCGTGCCATTCAGGGATTTGCAGTCGGCGGCGAATGGGGAGGCGCGGCGTTGCTTTCCGTTGAAAGTGCACCGAAAAATAAAAAAGCCTTTTACAGTAGCGGTGTACAAGTTGGCTACGGTGTAGGTTTACTGCTTTCAACCGGACTGGTTTCATTGATCAGT  >  NC_000913/2053934‑2054197

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