Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,956,719	C→T	F463F (TTC→TTT)	ilvA →	threonine deaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,956,719	0	C	T	100.0%	80.9 / NA	26	F463F (TTC→TTT)	ilvA	threonine deaminase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (9/17);  total (9/17)

GCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAGAAATTTTGCAGATGCTCAACGACGGCGGCTACAGCGTGGTTGATCTCTCCGACGACGAAATGGCGAAGCTACACGTGCGCTATATGGTCGGCGGACGTCCATCGCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTCCACTATCGCAGCCATGGCACCGACTACGGGCGCGTACTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTG‑GTTCGGCAAAATCTT  >  NC_000913/3956479‑3956950                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                         gcgcCTGAGCCGCGGCCTCGAAGAGCGCAAAGAAATTTTGCAGATGCTCAACGACGGCGGCTACAGCGTGGTTGATCTCTCCGACGACGAAATGGCGAAGCTACACGTGCGCTATATGGTCGGCGGACGTCCATCGCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTTCACTATCGCa                                                                                                                                                                                                                                <  1:136100/251‑1 (MQ=255)             cgGCCTCGAAGAGCGCAAAGAAATTTTGCAGATGCTCAACGACGGCGGCTACAGCGTGGTTGATCTCTCCGACGACGAAAGGGCGAAGCGACACGTGCGCTATATGGTCGGCGGACGGCCAGCGCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCGTCAACACGCTGGGTACGTACTGGACCATTTCTGTGTTTCACTATCGCAGCCATGGCACCg                                                                                                                                                                                                                    <  2:133245/251‑1 (MQ=255)                              aaGAAATTTTGCAGATGCTCAACGACGGCGGCTACAGCGTGGTTGATCTCTCCGACGACGAAATGGCGAAGCTACACGTGCGCTATATGGTCGGCGGACGTCCATCGCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTTCACTATCGCAGCCATGGCACCGACTACGGGCGCGTattg                                                                                                                                                                                                   <  1:84114/251‑4 (MQ=255)                                 aaaTTTTGCAGATGCTCAACGACGGCGGCTACAGCGTGGTTGATCTCTCCGACGACGAAATGGCGAAGCTACACGTGCGCTATATGGTCGGCGGACGTCCATCGCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTTCACTATCGCAGCCATGGCACCGACTACGGGCGCGTATTggcg                                                                                                                                                                                                <  1:221043/251‑1 (MQ=255)                                 aaaTTTTGCAGATGCTCAACGACGGCGGCTACAGCGTGGTTGATCTCTCCGACGACGAAATGGCGAAGCTACACGTGCGCTATATGGTCGGCGGACGTCCATCGCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTTCACTATCGCAGCCATGGCACCGACTACGGGCGCGTATTggcg                                                                                                                                                                                                <  2:120868/251‑1 (MQ=255)                                     gttGCTGATGCTCAACGACGGCGGCTACAGCGTGGTTGATCTCTCCGACGACGAAATGGCGAAACTACACGCGCGCTATATGGTCGGCGGACGTCCATCGCATCCGTTGCAGGAACGCCTCTACATCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTTCACTATCGCAGCCATGGCACCGACTACGGGCGCGTATTGgcggcg                                                                                                                                                                                             <  2:164418/249‑1 (MQ=255)                                         cAGATGCTCAAGGACGGCGTCTACAGCGTGGTTGAGCTCTCCGACGACGAAATGGCGAAGCTACACGTGCGCTATATGGTCGGCGGAGGTCCATCGCATCCGTTGCAGGACCGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTTCACTATCGCATCCATGGCACCGACTTCGGTCGCGTATTGGGGGCGGTCGa                                                                                                                                                                                        <  2:135725/251‑1 (MQ=255)                                                   acgacgGCGGCTACAGCGTGGTTGATCTCTCCGACGACGAAATGGCGAAGCTACACGTGCGCTATATGGTCGGCGGACGTCCATCGCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTTCACTATCGCAGCCATGGCACCGACTACGGGCGCGTATTGGCGGCGTTCGAACTTGGCGAc                                                                                                                                                                              <  1:36190/251‑1 (MQ=255)                                                                                  cgacgacgaAATGGCGAAGCTACACGTGCGCTATATGGTCGGCGGACGTCCATCGCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTTCACTATCGCAGCCATGGCACCGACTACGGGCGCGTATTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATg                                                                                                                                               <  1:56003/251‑1 (MQ=255)                                                                                          aaaTGGCGAAGCTACACGTGCGCTATATGGTCGGCGGACGTCCATCGCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTTCACTATCGCAGCCATGGCACCGACTACGGGCGCGTATTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGc                                                                                                                                       <  2:120751/251‑1 (MQ=255)                                                                                              ggCGAAGCTACACGTGCGCTATATGGTCGGCGGACGTCCATCGCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTTCACTATCGCAGCCATGGCACCGACTACGGGCGCGTATTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTAc                                                                                                                                    <  1:44300/250‑1 (MQ=255)                                                                                                       gcacGTGCGCTATATGGTCGCGGGACGTCCATCGCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTTCACTATCGCAGCCATGGCACCGACTACGGGCGCGTATTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCac                                                                                                                           <  1:15144/249‑1 (MQ=255)                                                                                                          cGTGCGCTATATGGTCGGCGGACGTCCATCGCATCCGTAGCAGGAACGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTTCACTATCGCAGCCATGGCACCGACTACGGGCGCGTATTGGCGGCGTTCGAACTTGGTGACCATGCACCGGGTTTCGAACCCCCGTTCAATGAGCTGGGCAACGATTGCCacgacg                                                                                                                       >  2:169568/1‑251 (MQ=255)                                                                                                            ggcgcTATATGGTCGGCGGACGTCCATCGCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTTCACTATCGCAGCCATGGCACCGACTACGGGCGCGTATTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGaa                                                                                                                     <  2:41346/250‑1 (MQ=255)                                                                                                                         cggcggACGTCCATCGCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCCGGGATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTTCACTATCGCAGCCATGGCACCGACTACGGGCGCGTATTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAAACCgg                                                                                                        >  2:241449/1‑251 (MQ=255)                                                                                                                               aCGTCCATCGCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTTCACTATCGCAGCCATGGCACCGACTACGGGCGCGTATTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCa                                                                                                  <  2:75571/251‑1 (MQ=255)                                                                                                                                     aTCGCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTTCACTATCGCAGCCATGGCACCGACTACGGGCGCGTATTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCt                                                                                            >  2:280271/1‑251 (MQ=255)                                                                                                                                      tCGCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCCGGCATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACCTTTCTTTGTTTCACTATCGCAGCCATGGCCCCGACTACGGGCGCGTATTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGACTGGGCGGGGCTCCGATTGCCATGACGAAACCAATAACCCGGCGTTCAGGttatt                                                                                           >  2:15328/1‑248 (MQ=255)                                                                                                                                                                      cGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTTCACTATCGCAGCCATGGCACCGACTACGGGCGCGTATTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAgcg                                                           <  1:47820/251‑1 (MQ=255)                                                                                                                                                                      cGAATTCCCGGAATCACAGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTTCACTATCGCAGCCATGGCACCGACTACGGGCGCGTATTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGAAGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGACTggtca                                                             >  2:72296/1‑245 (MQ=255)                                                                                                                                                                                         ggCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTTCACTATCGCAGCCATGGCACCGACTACGGGCGCGTATTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTAc                                        <  2:23624/251‑1 (MQ=255)                                                                                                                                                                                                              aaCACGCTGGGTACGTACTGGAACATTTCTTTGTTTCACTATCGCAGCCATGGCACCGACTACGGGCGCGTATTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTtgt                   >  1:46791/1‑251 (MQ=255)                                                                                                                                                                                                                           cgtacTGGAACATTTCTTTGTTTCACTATCGCAGCCATGGCACCGACTACGGGCGCGTATTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGGTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTG‑GTTCGGCAAAAt     >  1:246851/1‑251 (MQ=255)                                                                                                                                                                                                                             tacTGGAACATTTCTTTGTTTCACTATCGCAGCCATGGCACCGACTACGGGCGCGTATTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGGTAGGGAAAAAAGCCTTGTAGAGCTTCGCTTATTAGGCCTACCCGCGCGCCAACGTCATTT‑TGTGTTCGGCaaaaact   >  2:44632/1‑248 (MQ=255)                                                                                                                                                                                                                              ccTGGAACATTTCTTTGTTTCACTATCGCAGCCATGGCACCGACTACGGGCGCGTATTGGCGGCGTTCGAACTTGGCGACCATGACCCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTG‑GTTCGGCAGAATCt   <  2:157903/249‑1 (MQ=255)                                                                                                                                                                                                                              acTGGAACATTTCTTTGTTTCACTATCGCAGCCATGGCACCGACTACGGGCGCGTATTGGGGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGGTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGAGGGTTAGGGAAAAAAGCCTGATAGAGCTATGCTTAACAGGACTACCCGCGCGACAACGGCATTTGTG‑GTTCGGCAAAATCtt  >  1:118968/1‑251 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                         GCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAGAAATTTTGCAGATGCTCAACGACGGCGGCTACAGCGTGGTTGATCTCTCCGACGACGAAATGGCGAAGCTACACGTGCGCTATATGGTCGGCGGACGTCCATCGCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTCCACTATCGCAGCCATGGCACCGACTACGGGCGCGTACTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTG‑GTTCGGCAAAATCTT  >  NC_000913/3956479‑3956950

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