Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,520	A→T	P240P (CCA→CCT)	thrB →	homoserine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,520	0	A	T	100.0%	94.9 / NA	33	P240P (CCA→CCT)	thrB	homoserine kinase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (16/17);  total (16/17)

CATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCAGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCAGGCTTCCGGCAGGCGCGGCAGGCGGTCGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCAC  >  NC_000913/3277‑3760                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                 catcatAAGCCAGCAAGGGCCAGGGTTTGATGAGTGGCTGTGGGAGCAGGCGTATACGGGGATTAAAGTCTAGACGGCAGAAGAAAGGGCTATTTTAACGGCGCAGTATCGACGACAGGATTGAATTGCGCACGGGAGACATCTGGCTGGCTTCATTAACGCCAGCGATTCCAGACAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTGATCGCTGAACCATACCGAGAACGGTTACTGCCTGGCTTcc                                                                                                                                                                                                                                           >  1:165595/1‑251 (MQ=255)                 gTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCcggcaggcgcggca                                                                                                                                                                                                                              <  1:136492/248‑1 (MQ=255)                 gTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGACAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCcggcaggcgcggcaggc                                                                                                                                                                                                                           >  1:30218/1‑251 (MQ=255)                 gTGCCAGGGTTTGATGAGTGGCTGTGGGCGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGCTGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCcggcaggcgcggcacg                                                                                                                                                                                                                            >  2:136492/1‑248 (MQ=255)                                 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)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                 CATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCAGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCAGGCTTCCGGCAGGCGCGGCAGGCGGTCGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCAC  >  NC_000913/3277‑3760

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