Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,645	G→A	W282* (TGG→TAG)	thrB →	homoserine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,645	0	G	A	96.5%	95.9 / NA	29	W282* (TGG→TAG)	thrB	homoserine kinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (4/24);  minor base T (0/1);  total (4/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

GCGCACGGGCGACATCTGGCAGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCAGGCTTCCGGCAGGCGCGGCAGGCGGTCGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCA  >  NC_000913/3404‑3890                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      gcgcACGGGCGACATCTGGCAGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCAGGCTTCCGGCAGGCGCGGCAGGCGGTCGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTAGTTGGGTa                                                                                                                                                                                                                                               >  1:206871/1‑250 (MQ=255)            aCATCTGGCAGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCAGGCTTCCGGCAGGCGCGGCAGGCGGTCGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTAGTTGGGTAAGAACTACCTGc                                                                                                                                                                                                                                   <  1:109012/251‑1 (MQ=255)              aTCTGGCAGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCAGGCTTCCGGCAGGCGCGGCAGGCGGTCGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCACGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTAGTTGGGTAAGAACTACCTGCa                                                                                                                                                                                                                                  >  1:309616/1‑250 (MQ=255)                    cAGGGTTCATTCACGCCTGTTATTCCCGTCAGCCTGGGCTTGCCGCGAAGCTGAGGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCAGGCTTCCGGCAGGCGCGGCAGGCGGTCGCGGAAATCGGCGCGGTAGCGAGCGGTAGCTCCGGCTCCGGCCCGACCTTGTTCGCGCTGTGTGACAAGCCGGAAACCGCCCAGCGCTTTGCCGACGAGTTGGGTAAGAACTACCTGCAAAATCAg                                                                                                                                                                                                                           <  2:409940/251‑1 (MQ=255)                              tcaCGCCTGCTATTCCCGTCAGCCTGAGCTCGCCGCGAAGCAGACGAAAGCTGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCAGGCTTCCGGAAGGCGCGGCAGGCGGTCGCGGACATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCCGCCCGAACTTGTTCGATCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGGCTAGTTGGGTACGAACTACCTGCCAAATCAGGACGGTTTTg                                                                                                                                                                                                                 >  1:412782/1‑251 (MQ=255)                               caCGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGACCCCTACCGTGAACGGTTACTGCCAGGCTGCCGGCAGGCGCGGCAGGCGGTCGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCGCAGCGCGTTGCCGACTAGTTGGGTAAGAACTTCCTGCAAAAGCAGGAAGGTTTTGt                                                                                                                                                                                                                <  2:20399/251‑1 (MQ=255)                                   ccTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCAGGCTTCCGGCAGGCGCGGCAGGCGGTCGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTAGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTGTGTTCat                                                                                                                                                                                                            <  1:90375/251‑1 (MQ=255)                                         aTTCCCGGCAGCCTGAGCGTGCCGCGAAGCTGATGAAAGATGGTATCGCTGAACCCTACCGTGAACGGTTTCTGCCAGGCTTCCGGCAGGCGCTGCAGGCGGTCGCGGAAATCTGCGCGGTAGCGATCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCGGTTTTACAAGGCTGAAACCGCCCAGCGGGTTGCCGACTTGGTGGGTAAGAACTACCTGCGAAATCAGGAAGTTTGTGTTCATAttt                                                                                                                                                                                                        <  2:52847/249‑1 (MQ=255)                                                                           gaAAGATGTTATCGCTGATCCCTACCGTGAACGGTTACTGCCAGGCTTCCGGCAGGCGCGGCAGGCGGTCGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTAGTTGGGTAAGAACGACCTGCAAAATCAGGAAGGTTGTGTTCAGATTTGCCGGCTGGATTCGGCGGGCGCTCGAGTACTGGaaa                                                                                                                                                                    <  2:363951/251‑1 (MQ=255)                                                                           gaAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCAGGCTTCCGGCAGGCGCGGCAGGCGGTCGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTAGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGaa                                                                                                                                                                     <  1:318908/250‑1 (MQ=255)                                                                            aaaGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCAGGCTTCCGGCAGGCGCGGCAGGCGGGCGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTAGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGaaaa                                                                                                                                                                   <  2:361129/251‑1 (MQ=255)                                                                               gATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCAGGCTTCCGGCAGGCGCGGCAGGCGGTCGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTTGGTGACAAGCCGGAAACCGCCCATCGCGTTGCCGACGAGTTGGGTAAGAACTACCGGCAAAATCAGGAAGGTTTGGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTa                                                                                                                                                                <  1:345834/251‑1 (MQ=255)                                                                                attgtATCGCTGACCCCTACCGTGAACGGTTACTGCCAGGGTTCCGTCAGGCGCGGCAGGCGGTCGCGGAAATCGGCGCGGTAGCGAGCGGTTTCACCGGCTCCGGCCCGACCTTGGTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTAGTTTGGTAAGACCTACCTGGAAAATCAGGAAGGTTTTGTTCATATGTGCCGGCTGGATACGGCGGTCGCACGAGTACTGGAAAACtta                                                                                                                                                               <  2:107480/247‑3 (MQ=255)                                                                                       cGCTGAACCCTACCGTGAACGGTTACTGCCAGGCTTCCGGCAGGCGCGGCAGGCGGTCGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTAGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCCTATTTGCCGGCTGGATACGGCGGgcgccggggagctggggaaaaaaaaaaaa                                                                                                                                                         >  2:241711/1‑225 (MQ=255)                                                                                        gCTGAACCCTACCGTGAACGGTTACTGCCAGGCTTCCGGCAGGCGCGGCAGGCGGTCGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTAGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAAct                                                                                                                                                       <  1:387661/251‑1 (MQ=255)                                                                                                                                     cggcaggcgGTCGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTAGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGc                                                                                                          <  1:122307/251‑1 (MQ=255)                                                                                                                                            cgGTCGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTAGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCCGTAAcc                                                                                                   <  2:206871/251‑1 (MQ=255)                                                                                                                                                          ggCGCGTTAGCGAGCGGTATCCCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTAGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCCCAACGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCAGGGGTTGGGCaa                                                                                     <  2:298583/251‑1 (MQ=255)                                                                                                                                                                   gcgagcgGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTAGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCAGGGGTTGGGCAAAAATCAggg                                                                            <  1:241711/251‑1 (MQ=255)                                                                                                                                                                       gcggtcgctccggctccggCCCGACCTGGTTCGCGCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTAGTTGGGTAAGAACTCCCTGCAAACTCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAACCTCGACAATCTGAAAGATCACAACGATCATGTCAGCTTTGCGCAAGCCGTAACCCAGGGGTTGGGCAAAAATCAGGGGCTg                                                                        <  2:92659/244‑1 (MQ=255)                                                                                                                                                                        cggtttccccggcgccggcCCGACCTTGTTCGGGCGGTGTGACAAGCCGGAAACCGCCCAGCGCGGTGCGGCCTAGTTGGGTAAGAACTACCGGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGAGTCGGGGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGAGCACAACGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCAGGTGTTGGGCAAAAATCAGGGGCTGt                                                                       <  2:393495/243‑1 (MQ=255)                                                                                                                                                                                ccggctccggcCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTAGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCAGGGGTTGGGCAAAAATCAGGGGCTGttttt                                                                   <  1:321147/247‑1 (MQ=255)                                                                                                                                                                                        ggcCCGTCCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTAGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGTGCACGAGTACTGGAAAACTAAATGAATCTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTGGCACAAGCCGTAACCCAGGGGCTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCt                                                       <  1:150516/251‑1 (MQ=255)                                                                                                                                                                                                            gtgggACAAGCCGGAAACCGCCCAGCGTGTTGCCGACTAGTTGGGTAAGAACTACCGGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACGCTACAATCTGAAAGATCACAACGGGCAGGTCAGCTTGGCGCAAGTCGTAACCCAGGGGTTGGGCAAAAATCAGGGTCGGTTTTGTCCGCACGACCTGCCGGAATTCAGCctgac                                     <  2:290157/245‑1 (MQ=255)                                                                                                                                                                                                              gtgACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTAGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGaa                                 <  2:47461/251‑1 (MQ=255)                                                                                                                                                                                                                        ggAAACCGCCCAGCGCGTTGCCGACTAGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGTCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAAGGAAACTCTATAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGGCCTGCCGGAATTCAGCCTGACTGAAATTGATgaga                       <  1:253675/251‑1 (MQ=255)                                                                                                                                                                                                                            cccGCCCAGCGCGTTGCCGACTAGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGATATGCt                   <  1:328326/250‑1 (MQ=255)                                                                                                                                                                                                                                           gCCGACTAGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAAGTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGt    <  1:277963/251‑1 (MQ=255)                                                                                                                                                                                                                                             cGACTAGTTGGGTAAGAACTACCGGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCGACAAGCGGAAAGATCGCAACGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCATGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGGACGGCCTGCCGGAATTCAGCCTGACTGAGATTTATGAGAGGCTGAAGCTTGATTTTGTCa  <  2:255873/251‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      GCGCACGGGCGACATCTGGCAGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCAGGCTTCCGGCAGGCGCGGCAGGCGGTCGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCA  >  NC_000913/3404‑3890

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