Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,424	A→T	A208A (GCA→GCT)	thrB →	homoserine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,424	0	A	T	96.9%	91.4 / NA	32	A208A (GCA→GCT)	thrB	homoserine kinase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (15/16);  minor base G (1/0);  total (16/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGGCGAGCTGGAAGGCCGTATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTCGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCAGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCAGGCTTCCGGCAGGCGCGGCAGGCGGTCGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGC  >  NC_000913/3180‑3665                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  tGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGATTTCTTGTTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCGGGGCTATTTTACCGGAGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGTCATCTGGCTGGCttc                                                                                                                                                                                                                                             >  2:188885/1‑251 (MQ=255)     cgagctggaacttcgactcaccGGCAGCATTCCTTACGACAACGTGGTACCGTGTTTTCTTGGTGGTATGCACTTGATGATCGAAGAAAACGCCATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGAGTACCCGGGGATTAAAGTCTCGACGGCCGAAGCCATGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCGGGCttcattc                                                                                                                                                                                                                                         >  2:224754/21‑251 (MQ=255)     cGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCttcattc                                                                                                                                                                                                                                         >  2:39537/1‑251 (MQ=255)                                        aCGACAACGTGGCACCGTGTTTGCTTGGTGGTATGCAGTTGATGAGCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGTCATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCcgcg                                                                                                                                                                                                      <  2:186588/251‑1 (MQ=255)                                           ccAACGTGGTAACGTGTTTTCTCGGTGGTAAGCAGCTGATGAGCGAAGAAAACGACATCACCAGCCAGTAAGTGTCAGGGTTTCATGAGTGTCTGTGGGTGCTGGCGTATCCGGGTATTAAAGTCTCGACGGCAGAAGATAGGGCTATGTTACCGGCGCAGTATCGTCTCCAGGATTGCATTGCGCACGGGCGACATCTGGATGTTTTCATTCACGCCTGCTATTCCCTTCAGCCTGAGCTTGCCGCGAAg                                                                                                                                                                                                   <  1:175843/250‑1 (MQ=255)                                             aaCGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCt                                                                                                                                                                                                 <  1:110781/251‑1 (MQ=255)                                                     aCCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCCTCAGCCAGCAAGTGACAGGGTTTGATGAGTGGATGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGAGCAGAATCGACGCCAGGATTGCCTTGCGCACGGGCGACATCTGACTGGCTTCATTCACGCCTGCTATACGCGTCAGCCTGAGCATGCCGCGAAGATGATGaa                                                                                                                                                                                           >  1:113189/1‑249 (MQ=255)                                                                 tggtggtATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCg                                                                                                                                                                                <  2:67934/247‑1 (MQ=255)                                                                            aGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGt                                                                                                                                                                  <  1:39537/251‑1 (MQ=255)                                                                             gTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGATTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTg                                                                                                                                                                 >  1:171742/1‑251 (MQ=255)                                                                                 atgatCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACg                                                                                                                                                             <  1:13705/251‑1 (MQ=255)                                                                                   gatCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAGGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGt                                                                                                                                                           >  1:131961/1‑251 (MQ=255)                                                                                      cGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGTAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCGTTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCGTGCCGCGAAGCTGATGAAAGATGTTATTGCGGAACCCTTGCGTGAACGGTTa                                                                                                                                                         <  2:22579/250‑1 (MQ=255)                                                                                          gaaAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTAGTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGcc                                                                                                                                                    <  2:231763/251‑1 (MQ=255)                                                                                             aaCGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGcctgg                                                                                                                                                 >  1:110820/1‑248 (MQ=255)                                                                                                  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)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  TGGGCGAGCTGGAAGGCCGTATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTCGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCAGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCAGGCTTCCGGCAGGCGCGGCAGGCGGTCGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGC  >  NC_000913/3180‑3665

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