Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	86,315	C→T	T229M (ACG→ATG)	ilvI →	acetolactate synthase/acetohydroxybutanoate synthase, catalytic subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	86,315	0	C	T	96.2%	72.1 / NA	26	T229M (ACG→ATG)	ilvI	acetolactate synthase/acetohydroxybutanoate synthase, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (17/8);  total (18/8)
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.97e-01

GTGGTCGCCCAGGACCAGTAGTCGTTGATTTACCGAAAGATATTCTTAATCCGGCGAACAAATTACCCTATGTCTGGCCGGAGTCGGTCAGTATGCGTTCTTACAATCCCACTACTACCGGACATAAAGGGCAAATTAAGCGTGCTCTGCAAACGCTGGTAGCGGCAAAAAAACCGGTTGTCTACGTAGGCGGTGGGGCAATCACGGCGGGCTGCCATCAGCAGTTGAAAGAAACGGTGGAGGCGTTGAATCTGCCCGTTGTTTGCTCATTGATGGGGCTGGGGGCGTTTCCGGCAACGCATCGTCAGGCACTGGGCATGCTGGGAATGCACGGTACCTACGAAGCCAATATGACGATGCATAACGCGGATGTGATTTTCGCCGTCGGGGTACGATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGATCCTACTTCCATTTC  >  NC_000913/86081‑86561                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                       gTGGTCGCCCAGGACCAGTAGTCGTTGATTTACCGAAAGATATTCTTAATCCGGCGAACAAATTACCCTATGTCTGGCCGGAGTCGGTCAGTATGCGTTCTTACAATCCCACTACTACCGGACATAAAGGGCAAATTAAGCGTGCTCTGCAAACGCTGGTGGCGGCAAAAAAACCGGTTGTCTACGTAGGTGGTGGGGCAATCACGGCGGGCTGCCATCAGCAGTTGAAAGAAATGGTTGAGGCGTTGAAt                                                                                                                                                                                                                                        >  1:253021/1‑251 (MQ=255)                cAGTAGTCGTTGATTTACCGAAAGATATTCTTAATCCGGCGAACAAATTACCCTATGTCTGGCCGGAGTCGGTCAGTATGCGTTCTTACAATCCCACTACTACCGGACATAAAGGGCAAATTAAGCGTGCTCTGCAAACGCTGGTGGCGGCAAAAAAACCGGTTGTCTACGTAGGTGGTGGGGCAATCACGGCGGGCTGCCATCAGCAGTTGAAAGAAATGGTGGAGGCGTTGAATCTGCCCGTTGtttcc                                                                                                                                                                                                                         <  1:365076/251‑3 (MQ=255)                      tCGTTGATTTACCGAAAGATATTCTTAATCCGGCGAACAAATTACCCTATGTCTGGCCGGAGTCGGTCAGTATGCGTTCTTACAATCCCACTACTACCGGACATAAAGGGCAAATTAAGCGTGCTCTGCAAACGCTGGTGGCGGCAAAAAAACCGGTTGTCTACGTAGGTGGGGGGGCAATCACGGCGGGGTGACATCCGCAGTTGAAAGAAACGGGGGGGGGGGTGGAACTGCCCGTTGTTTCCTCAttt                                                                                                                                                                                                                   >  2:263503/1‑250 (MQ=255)                                            tCTTAATCCGGCGAACAAATTACCCTATGTCTGGCCGGAGTCGGTCAGTATGCGTTCTTACAATCCCACTACTACCGGACATAAAGGGCAAATTAAGCGTGCTCTGCAAACGCTGGTGGCGGCAAAAAAACCGGTTGTCTACGTAGGTGGTGGGGCAATCACGGCGGGCTGCCATCAGCAGTTGAAAGAAATGGTGGAGGCGTTGAATCTGCCCGTTGTTTCCTCATTGATGGGGCTGGGGGCGTTTCCg                                                                                                                                                                                              <  2:329330/250‑1 (MQ=255)                                              ttAATCCGGCGAACAAATTACCCTATGTCTGGCAGGAGTCGGTCAGTATGCGTTCTTACAATCCCCCTAATACCGGACATAAAGGGCAAATTAAGCGTGCTCTGCAAACGCTGGTGGCGGCAAAAAAACCGGTTGTCTACGTAGGTGGTGGGGCAATAACGGCGGACTGCAATCAGCAGTTGAAAGAAATGGTGGAGGCGTTTAATATGCCCGTTGGTTCCTCATTGATGGGGCTGGGGGCGTTTCCGGc                                                                                                                                                                                            >  1:363308/1‑250 (MQ=255)                                               tAATCCGGCGAACAAATTACCCTATGTCTGGCCGGAGTCGTTCAGTATGCGTTCTTACAATCCCACTACTACCGGACATAAAGGGCAAATTAAGCGTGCTCTGCAAACGCTGGTGGCGGCAAAAAAACCGGTTGTCTACGTAGGTGGTGGGGCAATCACGGCGGGCTGCCATCAGCAGTTGAAAGAAATGGTGGAGGCGTTGAATCTGCCCGTTGTTTCCTCATTGATGGGGCTGGGGGCGTTTCCgg                                                                                                                                                                                             <  2:361254/248‑1 (MQ=255)                                                  tCCGGCGAACAAATTACCCTATGTCTGGCCGGAGTCGGTCAGTATGCGTTCTTACAATCCCACTACTACCGGACATAAAGGGCAAATTAAGCGTGCTCTGCAAACGCTGGTGGCGGCAAAAAAACCGGTTGTCTACGTAGGTGGTGGGGCAATCACGGCGGGCTGCCCTCAACAGTTGAAAGAAATGGGGGAGGGGGTGAATCTGCCCGTTGTTTCCTCCTTGATGGGGCTGGGGGCGttttcgcccacgg                                                                                                                                                                                       >  1:357163/1‑241 (MQ=255)                                                  tCCGGCGAACAAATTACCCTATGTCTGGCCGGAGTCGGTCAGTATGCGTTCTTACAATCCCACTACTACCGGACATAAAGGGCAAATTAAGCGTGCTCTGCAAACGCTGGTGGCGGCAAAAAAACCGGTTGTCTACGTAGGTGGTGGGGCAATCACGGCGGGCTGCCATCAGCAGTTGAAAGAAATGGTGGAGGGGTTGAATCTGGCCGTTGTTTCCTCATTGATGGGGCTGGGGGCGGTTCCGGCAACGc                                                                                                                                                                                       >  2:160985/1‑251 (MQ=255)                                                       cGAACAAATTAACCTATGTCTGGCCGGAGTCGGTCAGTATGCGTTCTTACAATCCCACTACTACCGGACATAAAGGGCAAATTAAGCGTGCTCTGCAAACGCTGGTGGCGGCAAAAAAACCGGTTGTCTACGTAGGTGGTGGGGCAATCACGGCGGGCTGCCATCAGCAGTTGAAAGAAATGGGGGAGGCGTTGAATCTTCCCGTTGTTTCCTCATTTATGGGGCTGGGGGCGTTTCCGGCAACGCATCgg                                                                                                                                                                                  >  2:316927/1‑250 (MQ=255)                                                                           tGGCCGGAGTCGGTCAGTATGCGTTCTTACAATCCCACTACTACCGGACATAAAGGGCAAATTAAGCGTGCTCTGCAAACGCTGGTGGCGGCAAAAAAACCGGTTGTCTACGTAGGTGGTGGGGCAATCACGGCGGGCTGCCATCAGCAGTTGAAAGAAATGGTGGAGGCGTTGAATCTGCCCGTTGTTTCCTCATTGATGGGGCTGGGGGCGTTTCCGGCAAAGCATCGTCAGGCACTGGGCATGCTggg                                                                                                                                                              >  1:372367/1‑251 (MQ=255)                                                                             gCCGGAGTCGGTCAGTATGCGTTCTTACAATCCCACTACTACCGGACATAAAGGGCAAATTAAGCGTGCTCTGCAAACGCTGGTGGCGGCAAAAAAACCGGTTGTCTACGTAGGTGGTGGGGCAATCACGGCGGGCTGCCATCAGCAGTTGAAAGAAATGGGGGAGGCGTTGAATCTGCCCGTTGTTTCCTCATTGATGGGGCTGGGGGCGTTTCCGGCAACGCATCGTCAGGCACTGGGCATGCTGGGa                                                                                                                                                             >  1:276712/1‑250 (MQ=255)                                                                                 gAGTCGGTCAGTATGCGTTCTTACAATCCCACTACTACCGGACATAAAGGGCAAATTAAGCGTGCTCTGCAAACGCTGGTGGCGGCAAAAAAACCGGTTGTCTACGTAGGTGGTGGGGCAATCACGGCGGGCTGCCATCAGCAGTTGAAAGAAATGGTGGAGGCGTTGAATCTGCCCGTTGTTTCCTCATTGATGGGGCTGGGGGCGTTTCCGGCAACGCATCGTCAGGCACTGGGCATGCTGGGAATGc                                                                                                                                                         >  1:39454/1‑250 (MQ=255)                                                                                       gTCAGTATGCGTTCTTACAATCCCACTACTACCGGACATACAGGGCAAATTAAGCGTGCTCTGCAAACGCTGGTGGCGGCAAAAAAACCGGTTGTCTACGTAGGTGGTGGGGCAATCACGGCGGGCTGCCCTCAGCAGTTGAAAGAAATTGGGGGGGGGTTGAATCTGCCCGTTGTTTTCTCTTTGATGGGGCTGGGGGGGTTTCCGGCAACGCATCGTCAGGCACTGGGCATGCTGGGAATGCAAGGTAc                                                                                                                                                  >  2:115751/1‑251 (MQ=255)                                                                                         cAGTATGCGTTCTTACAATCCCACTACTACCGGACATAAAGGGCAAATTAAGCGTGCTCTGCAAACGCTGGTGGCGGCAAAAAAACCGGTTGTCTACGTAGGTGGTGGGGCAATCACGGCGGGCTGCCATCAGCAGTTGAAAGAAATGGTGGAGGCGTTGAATCTGCCCGTTGTTTCCTCATTGATGGGGCTGGGGGCGTTTCCGGCAACGCATCGTCAGGCACTGGGCATGCTGGGAATGCACGGTAcc                                                                                                                                                 <  1:8882/250‑1 (MQ=255)                                                                                            tATGCGTTCTTACAATCCCACTACTACCGGACATAAAGGGCAAATTAAGCGTGCTCTGCAAACGCTGGTGGCGGCAAAAAAACCGGTTGTCTACGTAGGTGGTGGGGCAATCACGGCGGGCAGTCCTCAGCAGTTTAAAGAAATGGGGGAGGCGTTGGATATGCCCGTTGTTTCCTAATTGATGGGGATGGGGGCGTTTTCGGCAACGCATCGGCAGGCACTGGGCATGGTGGGAATGCACGGGACCTACg                                                                                                                                             >  2:25718/1‑251 (MQ=255)                                                                                                               actactACCGGACATAAAGGGCAAATTAAGCGTGCTCTGCAAACGCTGGTGGCGGCAAAAAAACCGGTTGTCTACGTAGGTGGTGGGGCAATCACGGCGGGCTGCCATCAGCAGTTGAAAGAAATGGTGGAGGCGTTGAATCTGCCCGTTGTTTCCTCATTGATGGGGCTGGGGGCGTTTCCTGCAACGCAGCGTCAGGCACTGGGGATGCTGGGAATGCTCGGTACCTACGAAGCCAATATGACGATGCa                                                                                                                          <  1:33696/251‑1 (MQ=255)                                                                                                                                 gggCAAATTAAGCGTGCTCTGCAAACGCTGGTGGCGGCAAAAAAACCGGTTGTCTACGTAGGTGGTGGGGCAATCACGGCGGGCCGCCATCAGCAACTGAAAGAAATGGGGGCGGCGTTGAATCTGACCGTTGTTTCCTCATTCATGGGGCCGGGGGCGTTTCCGGCAACGCATCGTCAGGCACTGGGCATGCTGGGGATGCACGGTACCTACGAAGCCCATATGACGATGCATAACGCGGATTTGAtttt                                                                                                        >  2:179753/1‑251 (MQ=255)                                                                                                                                             cGTGCTCTGCAAACGCTGGTGGCGGCAAAAAAACCGGTTGTCTCCGTAGGTGGTGGGGCAATCACGGCGGGCCGCCATCAGCAAGTGAAAGAAATGGGGGCGGCGTTGAATCTGCCCGCTGTTTTCTCATTGGTGGGGCTGGGGGCGTTTCCGGCCACGCATCGTCAGGCAATGGGCATGCTGGGAATGCACGGTACCTACGAAGCCAATATGACGATGCATAAAGCGGATGTGATTTTCGCCGTCGGGGt                                                                                            >  1:230474/1‑251 (MQ=255)                                                                                                                                               tGCTCTGCAAAAGCTGGTGGCGGCAAAACAACCGGTTGTCTACGTAGGTGGTGGGGCAATCACGGCGGGCTGCCATCAGCAGGTGAAAGAAATGGTGGAGGCGTTGAATCTGCCCGTTGTTTCCTCATTGATGGGGCTGGGGGCGTTTCCGGCAACGCATCGTCAGGCACTGGGCATGCTGGGAATGCACGGTACCTACGAAGACAATATGACGATGCATAACGCGGATGTGATTTTCGCCGTCgggggac                                                                                          >  1:83381/1‑248 (MQ=255)                                                                                                                                                     gCAAACGCTGGTGGCGGCAAAAAAACCGGTTGTCTACGTAGGTGGTGGGGCAATCACGGCGGGCTGCCATCAGCAGTTGAAAGAAATGGTGGAGGCGTTGAATCTGCCCGTTGTTTCCTCATTGATGGGGCTGGGGGCGTTTCCGGCAACGCATCGTCAGGCACTGGGCATGCTGGGAATGCACGGTACCTACGAAGCCAATATGACGATGCATAACGCGGATGTGATTTTCGCCGTCGGGGTACGAtt                                                                                      <  1:309866/249‑1 (MQ=255)                                                                                                                                                                                             ggtgggggggCAATGACGGCGGGCTGCCATCAGCAGTTGAAAGAAATGGTGGAGGCGTTGAAGCTGCCCGTTGTGTCCCCATTGATGGGGCTGGGGGCGTTTCCGGCAACGCATCGTCAGGCACTGGGCATGATGGTAATGCACGGTACCTACGAAGCCAATATGACGATGCATAACGCGGATGGGATTTTCGCCGTCGGGGTACGATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCaaa                                            <  2:304149/245‑1 (MQ=255)                                                                                                                                                                                                         aTCACGGCGGGCTGCCATCAGCAGTTGAAAGAAATGGTGGAGGCGTTGAATCTGCCCGTTGTTTCCTCATTGATGGGGCTGGGGGCGTTTCCGGCAACGCATCGTCAGGCACTGGGCATGCTGGGAATGCACGGTACCTACGAAGCCAATATGACGATGCATAACGCGGATGTGATTTTCGCCGTCGGGGTACGATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCt                                <  2:372367/251‑1 (MQ=255)                                                                                                                                                                                                            aCGGCGGGCTGCCATCAGCAGTTGAAAGAAATGGTGGAGGCGTTGAATCTGCCCGTTGTTTCCTCATTGATGGGGCTGGGGGCGTTTCCGGCAACGCATCGTCAGGCACTGGGCATGCTGGGAATGCACGGTACCTACGAAGCCAATATGACGATGCATAACGCGGATGTGATTTTCGCCGTCGGGGGACGATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGcc                             >  2:104666/1‑250 (MQ=255)                                                                                                                                                                                                                      gCCATCAGCAGTTGAAAGAAATGGTGGAGGCGTTGAATCTGCCCGTTGTTTCCTCATTGATGGGGCTGGGGGCGTTTCCGGCAACGCATCGTCAGGCACTGGGCATGCTGGGAATGCACGGTACCTACGAAGCCAATATGACGATGCATAACGCGGATGTGATTTTCGCCGTCGGGGTACGATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATAtt                   >  2:153620/1‑251 (MQ=255)                                                                                                                                                                                                                             gcagTTGAAAGAAATGGTGGAGGCGTTGAATCTGCCCGTTGTTTCCTCATTGATGGGGCTGGGGGCGTTTCCGGCAACGCATCGTCAGGCACTGGGCATGCTGGGAATGCACGGTACCTACGAAGCCAATATGACGATGCATAACGCGGATGTGATTTTCGCCGTCGGGGTACGATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGATCCt             >  1:151675/1‑250 (MQ=255)                                                                                                                                                                                                                                       gaaaTGGTGGAGGCGTTGAATCTGCCCGTTGTTTCCTCATTGATGGGGCTGGGGGCGTTTCCGGCAACGCATCGTCAGGCACTGGGCATGCTGGGAATGCACGGTACCTACGAAGCCAATATGACGATGCATAACGCGGATGTGATTTTCGTCGTCGGGGTACGATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGATCCTACTTCCATTTc  >  1:330133/1‑251 (MQ=255)                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                       GTGGTCGCCCAGGACCAGTAGTCGTTGATTTACCGAAAGATATTCTTAATCCGGCGAACAAATTACCCTATGTCTGGCCGGAGTCGGTCAGTATGCGTTCTTACAATCCCACTACTACCGGACATAAAGGGCAAATTAAGCGTGCTCTGCAAACGCTGGTAGCGGCAAAAAAACCGGTTGTCTACGTAGGCGGTGGGGCAATCACGGCGGGCTGCCATCAGCAGTTGAAAGAAACGGTGGAGGCGTTGAATCTGCCCGTTGTTTGCTCATTGATGGGGCTGGGGGCGTTTCCGGCAACGCATCGTCAGGCACTGGGCATGCTGGGAATGCACGGTACCTACGAAGCCAATATGACGATGCATAACGCGGATGTGATTTTCGCCGTCGGGGTACGATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGATCCTACTTCCATTTC  >  NC_000913/86081‑86561

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: ‑