Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,547	C→T	V249V (GTC→GTT)	thrB →	homoserine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,547	0	C	T	100.0%	93.3 / NA	33	V249V (GTC→GTT)	thrB	homoserine kinase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (13/20);  total (13/20)

AGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCAGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCAGGCTTCCGGCAGGCGCGGCAGGCGGTCGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCC  >  NC_000913/3309‑3787                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                 ggTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGgc                                                                                                                                                                                                                                      <  2:17294/250‑1 (MQ=255)               cTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGGAATCGGCGCGGGCGCGAGCgg                                                                                                                                                                                                                        >  1:1192/1‑251 (MQ=255)                            ggcttAAAGTCTCGGCGGCAGAGGCCAGGGCTAGTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACTCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATTTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTGCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGATCGGTATctccggctcc                                                                                                                                                                                                           <  2:77899/248‑1 (MQ=255)                              aTTAAAGTCTCGGCGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGCGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGACATCGGCGCGGTAGCGCGCGGTATctacggctccgg                                                                                                                                                                                                         >  2:82147/1‑251 (MQ=255)                                 aaaGTCTCGACGGCAGAAGCCCGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGGAATCGGCGCGGGAGCGGGCGGTGTctccggccaccgccc                                                                                                                                                                                                      >  2:64144/1‑243 (MQ=255)                                                            aTTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGGAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGcc                                                                                                                                                                           >  1:76589/1‑251 (MQ=255)                                                                tACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCCCCGGCCCGACCTTGTTCGCTCTGTGTGGCAAGCCGGat                                                                                                                                                                       >  2:127454/1‑250 (MQ=255)                                                                  ccGGCGCACTCTCGTCGCCAGGATTGCACTGCGCGCGGGCGACACCTGGCTGGCTGCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAATCTGATGAAAGATGGTATCGCTGAACCCTACCGTGAACGGTTACTTCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAAGCGGCTCGGTAGCGAGCGGTATCTCCGGCTCCTGCCCGTCCTTGTTCGCTCTGTGTGACAAGCCGGatac                                                                                                                                                                     <  1:156653/251‑4 (MQ=255)                                                                              cgccgcCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGt                                                                                                                                                         <  1:250570/251‑1 (MQ=255)                                                                                       gATTGCAGTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCGGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCCGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTAGCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCACAGCGCGTTGCCGACtt                                                                                                                                                <  2:147270/251‑2 (MQ=255)                                                                                            cattgcGCACGGGCGCCATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGGTgg                                                                                                                                           >  2:262801/1‑251 (MQ=255)                                                                                                         cGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCt                                                                                                                              <  1:130786/251‑1 (MQ=255)                                                                                                           aCATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGTCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCGCCGGCCCGACCTTGTTCGCTCTGTGGGACAAGCCGTAGACCGCCCAGCGCGTGGCCGACTGGTTGGGTAAGAACTTCCTGc                                                                                                                            <  2:94956/251‑1 (MQ=255)                                                                                                                  gctggCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGCAATCGGCGCGGGAGCGAGCGGTATCTCCGGCTCCGGCCCGCCCCTGTTACCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGGCGACTGGTTGGGTAAGAACTACCTGCCAAATc                                                                                                                      >  2:80635/4‑250 (MQ=255)                                                                                                                          cattcaCGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCCCCGGCCCGACCTTGGTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTTCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAgggt                                                                                                             >  1:82946/1‑249 (MQ=255)                                                                                                                          cattcaCGCCTGCTATTCCCGGCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGATCCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGtt                                                                                                             <  1:151885/251‑1 (MQ=255)                                                                                                                            ttcaCGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGCTGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGCTGCGGAAATCGGCGCGGGAACGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGCGTGACAAGCCGGGTACCGCCCCGCGCGGTGCCGACTGGGTGGGTAAGAAATCCCTGCAAAATCAAGAAGGtttt                                                                                                           >  2:165048/1‑251 (MQ=255)                                                                                                                                 gCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCa                                                                                                      <  1:209602/251‑1 (MQ=255)                                                                                                                                        aTTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATAttt                                                                                                 <  2:33582/249‑1 (MQ=255)                                                                                                                                           cccGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCg                                                                                             <  1:182071/250‑1 (MQ=255)                                                                                                                                           cccGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTAGCTCTGTGTGACAAGCCGGATAACGCCCAGCGCGTTGCCGACTGGTTGGGGAAGAACTACCTGCCAAATTAGGAAGGGTTTGTTCATATTTGCCgg                                                                                            >  1:121885/1‑251 (MQ=255)                                                                                                                                                             cagcgAAGCCGATGAAGTATGTTATCGCTGAACCCGGCCTTGAACGGGTACTGCCTGGCTGCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGTCGCGGTAGCGTGCGGTATCTCCGGCCCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCACAGCGCGTTGCCGACTGGTTGGGTAAGAACTCCCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGgcgcc                                                                          <  2:219285/249‑2 (MQ=255)                                                                                                                                                                       gatgaAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGGAGCGAGCGGGATCTCCGGCTCCGGCCCGACCTTGTTCCCTCTGTGTGACAAGCCGGGTAACGCCCCGCGCGTTGCCGACCGGTTGGGTAAGAACTCCCTGCAAAATTCGGGAGGTTTTTTTCATATTTGCCGGCTGGATACggcgggggggagcgtgcttt                                                                >  2:112347/1‑237 (MQ=255)                                                                                                                                                                          gaAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGGTTTGTTCATATTTGCCGGCTGGATACggcgggcggcccagtaacggga                                                              >  1:210731/1‑236 (MQ=255)                                                                                                                                                                                     tCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGGTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGGGCAccggggcatggaaaaaaaaataa                                                   >  1:255985/1‑228 (MQ=255)                                                                                                                                                                                      cGCTGAACCCTACCTTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAAGCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAAc                                                 <  2:210731/251‑1 (MQ=255)                                                                                                                                                                                                 aCCGTGAACGGTTACTGCCTGGCTGCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTg                                      <  1:188885/251‑1 (MQ=255)                                                                                                                                                                                                      gagcGGTTACTGGCTGGCTGCCGGCAGGCGCTGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGGTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCGGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGa                                 <  1:224754/248‑1 (MQ=255)                                                                                                                                                                                                          ggTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATcac                             <  1:112347/251‑1 (MQ=255)                                                                                                                                                                                                                 gcctggCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGTAAGGTTTTGTTCATATTTGCCGGCTGTATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGc                      <  2:121885/247‑1 (MQ=255)                                                                                                                                                                                                                 gcctggCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAg                       <  1:262801/246‑1 (MQ=255)                                                                                                                                                                                                                               aggcgcggcaggcgGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTgcc        <  2:1192/251‑2 (MQ=255)                                                                                                                                                                                                                                     ggcaggcgGTTGCGGAAACCGGCGCGGTAGCTAGCGGTATCTCCGGTGCCTGCCCGACCTTGTTCGCTCTGTGTGACAAGCCAGATACCGCAAAGCGCGTTGCCGACTGGTTGGTTAAGAACTACCTGCAAAATCTGGAAGCTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCTCGAGTACTGGAAAACTAAATGAAACTCTACAATCTTAAAGATCACAATGAGCAGGTCAGCTTTGCGTAAGcc  <  1:267657/251‑1 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                 AGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCAGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCAGGCTTCCGGCAGGCGCGGCAGGCGGTCGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCC  >  NC_000913/3309‑3787

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑