Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,052,885	G→A	intergenic (+871/‑758)	yeeJ → / → shiA	inverse autotransporter adhesin/shikimate:H(+) symporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,052,885	0	G	A	100.0%	109.0 / NA	32	intergenic (+871/‑758)	yeeJ/shiA	inverse autotransporter adhesin/shikimate:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (14/18);  total (14/18)

TTCTGCAATCGATAATGAAAATGTACGTTCCTGCGAGCTACCTTCTAAATTGAACGCAATGTAAGATCCTAACGATCGCATTTCCTCGCGCACCTCATCGAGTACATCCTCACTTAGTGGCAATTCATATATTGGCCTGACTGCTGGAAAACCCGCCTCACGCATCATAAATGCCCATGTCATAGGTACGGGAGCCCGGAGTTTCTGATCCATACTGGACGCGTTCTTGCACAAAGGGGAGAAGCAATTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTACGCATCCCCTCAATGCAAAGATCCGGCGTACCGTATTGCTGTGTTATGTTCTTTGCTAAATCTTTTATTTCTTTTAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAA  >  NC_000913/2052645‑2053127                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                   ttCTGCAATCGATAATTAAAATGTACGTTCCTGCGAGCTACCTTCTAAAGTGAACGCAATGTAAGATCCTAACGCTCGCATTTCCTCGCGCTCCTCAGCGAGTCCAGCCTCACTGAGTGGCAATTCAAATATTGGCCTCACTGCTGGAAAACCCGCCTCACGCATCATAAATGCCCATGTCATAGGTACGGGAGCCCGGAGTTTCTTTTCCATACTGGACGCGTTCTTGCACAAAGGGGAAAATCAATTc                                                                                                                                                                                                                                           <  1:714725/250‑1 (MQ=255)     gCAATCGATAATGAAAATGTACGTTCCTGCGAGCTACCTTCTAAATTGAACGCAATGTAAGATCCTAACGATCGCATTTCCTCGCGCACCTCATCGAGTACATCCTCACTTAGTGGCAATTCATATATTGGCCTGACTGCTGGAAAACCCGCCTCACGCATCATAAATGCCCATGTCATAGGTACGGGAGCCCGGAGTTTCTGATCCATACTGGACGCGTTCTTGCACAAAGGGGAAAAGCAATTCAt                                                                                                                                                                                                                                         <  1:606417/248‑1 (MQ=255)                                cgctatcTACCTTCAAAATTGATCGCAATTTAAGATCCTAACAATCTCATTTCCCCGCGCTCCTCATCGATTTCATCCTCACTTATTGGCAATGCATATATGTGCCGGACTGCTGAAAAACCCGCCTCACGCATCATAAATGCCCATGTCATAGGTACGGGAGCCCGGAGTTTCTGTTCCATTCTGGACGCGTACTTGCACAAAGGGGAAAAGCAATTCATGGTTATACCAACAACCTGAAAATCCGtttt                                                                                                                                                                                                           <  2:440280/245‑1 (MQ=255)                                  cGAGCTACCTTCTAAATTGAACTCAATGTAAGATCCTAACGACCGCATTTCCTCGCGCACCTCATCGAGTACATCCTCACTTAGTGGCAATTCAGATATTGGCCTGACTGCTGGAAAACCCGCCTCACGCATCATAAATGCCCATGTCATAGGGACGGGAGCCCGGAGTTTCTGATCCATACTGGACGCGTTCTTGCACAAAGGGGAAAAGCAATTCATGGTTATGCCAACAACCTGAAAATTCGTTTTTg                                                                                                                                                                                                         <  1:592644/251‑1 (MQ=255)                                                 aTTGAACGCAATGTAAGATCCTAACGATCGCATTTCCTCGCGCACCTCATCGAGTACATCCTCACTTAGTGGCAATTCATATATTGGCCTGACTGCTGCAAAACCCGCCTCACGCATCATAAATGCCCATGTCATAGGTACGGGAGCCCGGAGTTTCTGATCCATACTGGCCGCGTTCTTGCCCAAAGGGGAAAAGCAATTCATGGTTATACCAACAACCTGAAACTTCGTTTTTGCTTTCAactgactg                                                                                                                                                                                           >  1:744362/1‑250 (MQ=255)                                                                    ccTAACGATCGCATTTCCTCGCGCACCTCATCGAGTACCTCCTCACTTAGTGGCAATTCATATATTGGCCTGACTGCTGGAAAACCCGCCTCACGCATCATAAATGCCCATGTCATAGGTACGGGAGCCCGGAGTTTCTGATCCATACTGGACGCGTTCTTGCACACAGGGGAAAAGCAATTCATGGTTATACCAACAACCTGAAACTTCGATTTTGCTTTCAACTGAATGATAAATAACATCGTTTTc                                                                                                                                                                         >  2:446204/1‑249 (MQ=255)                                                                             cGCATTTCCTCGCGCACCTCACCGAGTCCACCCTCACTTAGTGCCAATTCATATATTGCCCTCACTGCTGAAAAACCCGCCTCACGCATCATAAATGCCCATGTCATAGGTACGGGAGCCCGGAGTTTCTGATCCATACTGGACGCGTTCTTGCACAAAGGGGAAAAGCAATTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCttttt                                                                                                                                                              <  1:425715/251‑1 (MQ=255)                                                                                aTTTCCTCGCGCACCTCATCGAGTACATCCTCACTTAGTGGCAATTCATATATTGGCCTGACTGCTGGAAAACCCGCCTCACGCATCATAAATGCCCATGTCATAGGTACGGGAGCCCGGAGTTTCTGATCCATACTGGACGCGTTCTTGCACAAAGGGGAAAAGCAATTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCCGGTTATTTTTACg                                                                                                                                                           >  2:328856/1‑251 (MQ=255)                                                                                         cgcACCTCATCGATTACATCCTCACTTAGTGGCAATTCATATATTGGCCTGACTGCTGGAAAACCCGCCTCACGCATCATAAATGCCCATGTCATAGGTACGGGAGCCCGGAGTTTCTGATCCATACTGGACGCGTTCTTGCACAAAGGGGAAAAGCAATTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTACGCATCCCCTc                                                                                                                                                  <  1:739484/251‑1 (MQ=255)                                                                                            cccTCATCGAGTACATCCCCACTTAGAGGCAATTCATATATTGGCCTGACTGCTGGAAAAACCGCCTCACGCATCATAAATGCCCATGTCATAGGAACGGGCGCCCGGAGTTTCTGATCCATACTGAACGCGTTCTTGCACAAAGGGGAAAAGCAATTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATACCATCGTTTTCAGGTTCTTTTTACGCATCCCCTCAAt                                                                                                                                               >  2:566898/2‑251 (MQ=255)                                                                                                     aGTACAGCCTCACTTAGTGGCAATTCATATATTGGCCTGACTGCTGGAAAACCCGCCTCACGCATCATAAATGCCCATGTCATAGGTACGGGAGCCCGGAGTTTCTGATCCATACTGGACGCGTTCTTGCACAAAGGGGAAAAGCAATTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTACGCATCCCCTCAATGCaaa                                                                                                                                          <  2:743702/247‑1 (MQ=255)                                                                                                            ccTCACTTAGTGGCAATTCATATATTGGCCTGACTGCTGGAAAACCCGCCTCACGCATCATAAATGCCCATGTCATAGGTACGGGGGCCCGGAGTTTCTGACCCATACTGGACGCGTTCTTGCACAAAGGGGAAAAGCAATTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTACGTATCCCCTCATTGCAAAGATCCGGctt                                                                                                                                >  2:554462/1‑248 (MQ=255)                                                                                                               cACTTAGTGGCAATTCATATATTGGCCTGACTGCTGGAAAACCCGCCTCACGCATCATAAATGCCCATGTCATAGGTACGGGAGCCCGGAGTTTCTGATCCATACTGGACGCGTTCTTGCACAAAGGGGAAAAGCAATTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTACGCATCCCCTCAATGCAAAGATCCGGCGTAcc                                                                                                                             <  1:107441/250‑1 (MQ=255)                                                                                                                        gCAATTCATATATGGGCCTTACTGCTGGAAAACCCGCCTCACGCATCATAAATGCCCATGTCATAGGTACGGGAGCCCGGAGTTTCTGATCCATACTGGACGCGTTCTTGCACAAAGGGGAAAAGCAATTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTACGCATCCCCTCAATGCAAAGATCCGGCGTACCGTATTGCtgt                                                                                                                   <  1:670088/251‑1 (MQ=255)                                                                                                                                 aaattggccatagttcgggAAAACCCGCCTCACGCATCATAAATGCCCATGTCATAGGTACGGGAGCCCTGATATTCTGATCCATACTGGACGCGTTCTTGCAAAATGGGGAAAAGCAATTCATGGTTATACCAACAATCAGAAAAGTCGTGTTTGCTTTCAACTTACTGATAAATAACATCGTTTTCAGGTGTTTTTTACGCATCCCCTCAATGCAAAGATCCGGCGGACCGTATTGCTGTGTTAGGTTc                                                                                                          <  2:591453/234‑1 (MQ=255)                                                                                                                                               gctgTAAACCCCGCCTCACGCAGCATAAATGCCCATGTCATAGGTACGGGAGCCCGGAGTTTCTGATCCATACTGGACGCGTTCTTGCACAAAGGGGAAAAGCAATTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTACGCATCCCCTCAATGCAAAGATCCGGCGTACCGTATTGCTGTGTTATGTTCTTTGCTAAATCt                                                                                              <  1:446204/249‑1 (MQ=255)                                                                                                                                                ctgGAAAACCCGCCTCACGCATCATAAATGCCCATGTCATAGGTACGGGAGCCCGGAGTTTCTGATCCATACTGGACGCGTTCTTGCACAAAGGGGAAAAGCAATTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTACGCATCCCCTCAATGCAAAGATCCGGCGTACCGTATTGCTGTGTTATGTTCTTTGCTAAATCtttt                                                                                           >  1:493635/1‑251 (MQ=255)                                                                                                                                                   gAAAACCCGCCTCACGCATCATAAATGCCCATGCCATAGGTACGGGAGCCCGGCGTTTCTGATCCATACTGGACGCGTTCTTGCACAAAGGGGAAAAGCAATTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCTTTTTCAGGTTCTTTTTACGAATCCCCTCACTGCAAAGATCCGTATTACCGTATTGCTGTGTTATGTTATTTGCTAAATCTTTTAtt                                                                                        >  2:417296/1‑251 (MQ=255)                                                                                                                                                     aaaCCCGCCTCACGCATCATAAATGCCCATGTCATAGGTACGGGAGCCCGGAGTTTCTGATCCATACTGGACGCGTTCTTGCACAAAGGGGAAAAGCAATTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTACGCATCCCCTCAATGCAAAGATCCGGCGTACCGTATTGCTGTGTTATGTTCTTTGCTAAATCTTTTATTTc                                                                                      <  2:223943/251‑1 (MQ=255)                                                                                                                                                      aaCCCGCCTCACGCATCATAAATGCCCATGTCATAGGTACGGGAGCCCGGAGTTTCTGATCCATACTGGACGCGTTCTTGCACAAAGGGGAAAAGCAATTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTACGCATCCCCTCAATGCAAAGATCCGGCGTACCGTATTGCTGTGTTATGTTCTTTGCTAAATCTTTTATTTCt                                                                                     >  1:210099/1‑251 (MQ=255)                                                                                                                                                      aaCCCGCCTCACGCATCATAAAGGCCCATGTCATAGGTACGGGAGCCCGGAGTTTCTGATCCATACTGGACGCGTTCTTGCACAAAGGGGAAAAGCAATTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTACGCATCCCCTCAATGCAAAGATCCGGCGTACCGTATTGCTGTGTTATGTTCTTTGCTAAATCTTTTATTTCt                                                                                     <  1:453010/251‑1 (MQ=255)                                                                                                                                                              tCACGCATCATAAATGCCCATGTCATAGGTACGGGAGCCCGGAGTTTCTGATCCATACTGGACGCGTTCTTGCACAAAGGGGAAAAGCAATTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTACGCATCCCCTCAATGCAAAGATCCGGCGTACCGTATTGCTGTGTTATGTTCTTTGCTAAATCTTTTATTTCTTTTAAt                                                                               >  1:254204/1‑249 (MQ=255)                                                                                                                                                              tCACGCATCATAAATGCCCATGTCATAGGTACGGGAGCCCGGAGGTTCTGATCCATACTGGACGCGTTCTTGCACAAAGGGGAAAAGCAATTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTACGCATCCCCTCAATGCAAAGATCCGGCGTACCGTATTGCTGTGTTATGTTCTTTGCTAAATCTTTTATTTCTTTTAAt                                                                               <  2:254204/249‑1 (MQ=255)                                                                                                                                                                      catAAATGCCCATGTCATAGGTACGGGAGCCCGGAGTTTCTGATCCATACTGGACGCGTTCTTGCACAAAGGGGAAAAGCAATTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTACGCATCCCCTCAATGCAAAGATCCGGCGTACCGTATTGCTgtgt                                                                                                                 <  2:173313/207‑1 (MQ=255)                                                                                                                                                                      catAAATGCCCATGTCATAGGTACGGGAGCCCGGAGTTTCTGATCCATACTGGACGCGTTCTTGCACAAAGGGGAAAAGCAATTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTACGCATCCCCACAATGCAAAGATCCGGCGTACCGTATTGCTgtgt                                                                                                                 >  1:173313/1‑207 (MQ=255)                                                                                                                                                                                 aTGTCATAGGTACGGGAGCCCGGAGTTTCTGATCCATACTTGACGCGTTCTTGCACAAAGGGGAAAAGCAATTCATGGTTATTCCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCGTTTTACGCATCCCCTCAATGCAAAGATCCGGCGTACCGTATTGCTGTGTTATGTTCTTTGCTAAATCTTTTATGTCTTTTCATGTTGCGTGATCCTGCATAGTc                                                          <  2:193466/251‑1 (MQ=255)                                                                                                                                                                                         ggTACGGGAGCCCGGAGTTTCTGATCCATACTGGACGCGTTCTTGCACAAAGGGGAAAAGCAATTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTACGCATCCCCTCAATGCAAAGATCCGGCGTACCGTATTGCTGTGTTATGTTCTTTGCTAAATCTTTTATTTCTTTTAATGTTGCGTGATCCTGCATAGTCATTGTGAc                                                  <  1:328856/251‑1 (MQ=255)                                                                                                                                                                                                           ttCTGATCCATACTGGACGCGTTCTTGCACAAAGGGGAAAAGCAATTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTACGCATCCCCTCAATGCAAAGATCCGGCGTACCGTATTGCTGTGTTATGTTCTTTGCTAAATCTTTTATTTCTTTTAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCt                                 <  2:716222/250‑1 (MQ=255)                                                                                                                                                                                                                aTCCATACTGGACGCGTTCTTGCACAAAGGGGAAAAGCAATTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTACGCATCCCCTCAATGCAAAGATCAGGCGAACCGTATTGCTGTGTTATGTTCTTTGCTAAATCTTTTATTTCTTTTAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTCGTCTGTTCa                            >  1:561986/1‑250 (MQ=255)                                                                                                                                                                                                                             gcgTTCTTGCACAAAGGGGAAAAGCAATTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTACGCATCCCCTCAATGCAAGGATCCGGCGTACCGTATTGCTGTGTTATGTTCTTTGCTAAATCTTTTATTTCTTTTAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCAtt              >  1:149527/1‑251 (MQ=255)                                                                                                                                                                                                                               gTTCTTGCACAAAGGGGAAAAGCAATTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTACGCATCCCCTCAATGCAAAGATCCGGCGTACCGTATTGCTGTGTTATGTTCTTTGCTAAATCTTTTATTTCTTTTAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTa             >  1:257671/1‑250 (MQ=255)                                                                                                                                                                                                                                          aaGGGGAAAAGCAATTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTACGCATCCCCTCAATGCAAAGATCCGGCGTACCGTATTGCTGTGTTATGTTCTTTGCTAAATCTTTTATTTCTTTTAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTaa  >  1:322834/1‑250 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                   TTCTGCAATCGATAATGAAAATGTACGTTCCTGCGAGCTACCTTCTAAATTGAACGCAATGTAAGATCCTAACGATCGCATTTCCTCGCGCACCTCATCGAGTACATCCTCACTTAGTGGCAATTCATATATTGGCCTGACTGCTGGAAAACCCGCCTCACGCATCATAAATGCCCATGTCATAGGTACGGGAGCCCGGAGTTTCTGATCCATACTGGACGCGTTCTTGCACAAAGGGGAGAAGCAATTCATGGTTATACCAACAACCTGAAAATTCGTTTTTGCTTTCAACTGACTGATAAATAACATCGTTTTCAGGTTCTTTTTACGCATCCCCTCAATGCAAAGATCCGGCGTACCGTATTGCTGTGTTATGTTCTTTGCTAAATCTTTTATTTCTTTTAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAA  >  NC_000913/2052645‑2053127

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