Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,622	A→T	E274D (GAA→GAT)	thrB →	homoserine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,622	0	A	T	83.9%	70.1 / 2.7	31	E274D (GAA→GAT)	thrB	homoserine kinase
Reads supporting (aligned to +/- strand):  ref base A (3/2);  new base T (7/19);  total (10/21)
Fisher's exact test for biased strand distribution p-value = 2.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.03e-01

CGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCAGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTG‑CCGCGAAGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCAGGCTTCCGGCAGGCGCGGCAGGCGGTCGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCCGT‑AACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAG  >  NC_000913/3378‑3868                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                        cgcAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTG‑CCGCGAAGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGAGGGCGGGATCTCCGGATCCGGCCCGACCTTGGTCGCTCTGGGTGTCAAGCCGGGAACCGcc                                                                                                                                                                                                                                                   >  1:265378/1‑251 (MQ=255)   cAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTG‑CCGCGAAGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGGATCTCCGGCCCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATAACGcccc                                                                                                                                                                                                                                                 >  1:338857/1‑250 (MQ=255)       aTCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTG‑CCGCGAAGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTTCGGAAATCGGCGCGGGAGCGAGCGGTATCTCCGGCTCCGGCCCGGCCTTGTTCGCTCTGTGTGACAAGCCGGAATACGCCCAgcgc                                                                                                                                                                                                                                             >  1:196253/1‑251 (MQ=255)       aTCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTG‑CCGCGAAGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGGAATAGGGGCGGGAACGAGAGGAATCTCCGGCTCCGGCCCGACCATGTTCGCTCTGTGTCACAAGACGGATACCGCCCAgcgc                                                                                                                                                                                                                                             >  1:261877/1‑251 (MQ=255)          gccgccAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTG‑CCGCGAAGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGtt                                                                                                                                                                                                                                          <  1:178858/251‑1 (MQ=255)           cctccAGGATTGCATTGCGCTCGGGCGACAGTTGGCTGGCTTCATTCACGCCGGCGATTCCCGTCAGCCTGAGCTTG‑CCGGGAAGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCGGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCTGTAGCGAGCGGTATCTCCGGCTCCGTCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTg                                                                                                                                                                                                                                         <  2:368578/248‑1 (MQ=255)                   attgcattgcGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTG‑CCGCGAAGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCCGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGGTACCGCCCAGCGCGTTGCCGGCTgg                                                                                                                                                                                                                                 >  1:122527/1‑251 (MQ=255)                      gcatggcGCACGGGTGACCTCTGGCGGGCTTCAGTCACGCCTGCTTTTCCCGGCCGCCTGAG‑TTGTCCGCGAAGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCGCCGGCCCGACCTGGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGtt                                                                                                                                                                                                                               <  2:338857/250‑1 (MQ=255)                                              ctggCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTG‑CCGCGAAGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGGAGCCAGCGGGATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATc                                                                                                                                                                                                        >  1:52712/3‑249 (MQ=255)                                                 gCTTCATGCACGCCTGATATTCTCGTCAGCCTGTGCTTG‑CCGAGTAGCTGATGAAAGA‑GTGTATCGTTGAACCCTCCCGGGAACGGTTAATGCCGGGCTTCCGGCAGGCGCTGCAGGCGGTGGCGGAAAACGGCGCGGGAGTGAGCGGTATCTCCGTCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCGCATCGCGTTGCCGACTGGTTGGGTAAGAACTTCCTGCAAATTCAGGaa                                                                                                                                                                                                   <  2:251038/251‑1 (MQ=255)                                                         caCGCCTGCTATTCCCGTCAGCCTGAGCTTG‑CCGCGAAGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGTGGTTGCGGAAATCGGCGCGGTAGTGAGCGGTATCGCCTGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGt                                                                                                                                                                                           <  1:48111/251‑1 (MQ=255)                                                           cGCCTGCTATTCCCGTCAGCCTGAGCTTG‑CCGCGAAGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGGAATCGGCGCGGTAGCGAGCGGGATCTCCGGCTCCGGCCCGACCTTGGTCGCTCTGTGTGACAAGCCGGGTACCGCCCAGCGCGATGCCGACTGGGTGGGGAAGAACTACCTGCACAAACAGGAAGGTTTTGttt                                                                                                                                                                                         >  1:137976/1‑250 (MQ=255)                                                                     tCCCGTCAGCCTGAGCTTG‑CCGCGAAGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGCCTGGGTGGGGAAGAACTACCTGCAAAATCAGGAAGGGTTTGTTCATATTTTCCg                                                                                                                                                                               >  1:78528/1‑251 (MQ=255)                                                                           cAGCCTGAGCTTG‑CCGCGAAGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTgg                                                                                                                                                                          <  1:110133/250‑1 (MQ=255)                                                                                 gtgcgtG‑CCGCGAAGCTGATGAAGGATGT‑TATCGCTGATCCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATAcg                                                                                                                                                                     <  2:78528/244‑1 (MQ=255)                                                                                           gcgAAGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCAGGCTTCCGGCAGGCGCGGCAGGCGGTCGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACggcggggggccg                                                                                                                                                          >  1:167627/1‑245 (MQ=255)                                                                                                   gatgaAAGATGG‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCAGGCTTCCGGCAGGCGCGGCAGGCGGTCGCTGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTgg                                                                                                                                                  <  2:167627/251‑1 (MQ=255)                                                                                                      gaAAGATGT‑TATCGCTGAACACTACCGTGAACGGGTACTGCCTGGCTTCCGGCATGCGCGGCCGGCGGTTGCGGAAATCGGCTCGGTAGCGCTCGGTATCTCCGGCTCCGGCCCTACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACGGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCGGGATACGGCGGGCGCACGAGTACTGGaaa                                                                                                                                               <  2:9166/251‑1 (MQ=255)                                                                                                       aaaGTTGT‑TATCGCTGAACCCTACCGTGAACGGGTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAAGCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGaaaa                                                                                                                                              <  1:12004/251‑1 (MQ=255)                                                                                                                       gAACCCTACCGTGAACGTTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCTGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTTCTGGAAAACTAAATGAAACTCTa                                                                                                                               <  1:330107/251‑1 (MQ=255)                                                                                                                          cccTACCGTGAACTGTTACTTCCCGGCTTCCGGCATGCGCGGCATGCGGTTGCGTAAAGCGGCGCGGTAGCGAGACGTATCTCCGGCTCCGGCCCGACCTCGTTCGCTCTGTGTGACAAGCCGGATTCCGCCCAGCGCGTCGCCGACTTGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACaa                                                                                                                            <  2:154166/251‑1 (MQ=255)                                                                                                                            cTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCaaaa                                                                                                                                                                                                          <  1:150559/171‑1 (MQ=255)                                                                                                                                              gcctggCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAg                                                                                                         <  1:18553/246‑1 (MQ=255)                                                                                                                                                  ggttcccGGCAGGCGCGGCAGGCGGTTGCTGAAATCGGCGCGGTGGCGAGCGGTATCTCCGGCGCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGTCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAAGCAGGAAGGTTTTGGTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGTAAAACTAAATGAAACTCTACAATCTGAATGATCACAATGAGCAGGt                                                                                                    <  2:196253/246‑1 (MQ=255)                                                                                                                                                            cggcgcggcagtagGTTGCTGAAATCGGCGCGGGAGCGATCGGTATCGCCGTCTCCGGCCCGACCTTGGTCGCTCTGTGTGCCAAGCCGGATACCGCCCATTGCGTTTCCGACTGGTTGGGTAAGAACTACCTTCAACATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTAATGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTgcg                                                                                          <  2:302839/250‑1 (MQ=255)                                                                                                                                                                          gttgcgGCAATCGGCCCGGTAGCGAGGGTTACCCCCGGCTCCGGCCCGACATGGTTTGTTCTGTGTGACAATCCGGATACCGCTCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACGAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGT‑AACCCa                                                                           <  2:4234/248‑1 (MQ=255)                                                                                                                                                                                     cGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCCGT‑AACCCAGGGGTTGGGCa                                                                <  1:208635/251‑1 (MQ=255)                                                                                                                                                                                     cGGCGCGGTAGCGAGCGGTAGCTCCGGCTCCTGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGTCCAGCGCGTTGCCGACTGGTTGGGTAATAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTCCAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGT‑AACCCAGGGGTTGGGCa                                                                <  1:48066/251‑1 (MQ=255)                                                                                                                                                                                                            ccggctccggcCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGT‑AACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTc                                           <  1:68650/249‑1 (MQ=255)                                                                                                                                                                                                            ccggctccggcACGGCCTTGTTCTCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGT‑AACCCAGGGGTTGGGCAAAAATCAGGGGCTGttttt                                             <  1:217421/247‑1 (MQ=255)                                                                                                                                                                                                                                                   cGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGGACGAGGGCTGGGAAACTAAAAGAAACTCTACCAACTTAAAGAGTACAAAAAACAGGTGCGCTTTTCGCAAGACGTCAACCCAGGGGTGGGGGCACAA‑CACGGGGTGTTTTTTCCGCCCGCCCTGCCGGAAATCAGcccgtcttgaatttttggg  >  1:65061/1‑234 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                        CGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCAGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTG‑CCGCGAAGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCAGGCTTCCGGCAGGCGCGGCAGGCGGTCGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCCGT‑AACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAG  >  NC_000913/3378‑3868

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