Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,956,510	G→A	E394K (GAA→AAA)	ilvA →	threonine deaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,956,510	0	G	A	100.0%	21.7 / NA	10	E394K (GAA→AAA)	ilvA	threonine deaminase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (7/3);  total (7/3)

CAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGACCATTCCGGAAGAAAAAGGCAGCTTCCTCAAATTCTGCCAACTGCTTGGCGGGCGTTCGGTCACCGAGTTCAACTACCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAGAAATTTTGCAGATGCTCAACGACGGCGGCTACAGCGTGGTTGATCTCTCCGACGACGAAATGGCGAAGCTACACGTGCGCTATATGGTCGGCGGACGTCCATCGCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTCCACTATCGCAGCCATGGCACCGACTACGGGCGCGT  >  NC_000913/3956313‑3956754                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     cAGAACGCTGCGAACTGGGAGAACAGCGTGAAGCGTTGTTGGCGGTGACCATTCCGGAAGAAAAAGGCAGCTTCCTCAACTTCTGCCAACTGCTTGGCGGGCGTTCGGTCACCGAGTTCAACTACCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAAAAATTTTGCAGATGCTCAACGACGGCGGCTACAGCGTGGTTGATCTCTCcgac                                                                                                                                                                                                 >  1:44541/1‑251 (MQ=255)                                                   aTTCCGGAAGAAAAAGGCAGCTTCCTCAAATTCTGCCAACTGCTTGGCGGGCGTTCGGTCACCGAGTTCAACTACCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTACGCCTGAGCCGCGCCCTCGAAGAGCGCAAAAAAATTTTGCAGATGATCAACGACGGAGGATACAGCGTGGTTGATCTCTCCGACGACGAAATGGCGAAGATACACGTGCGCTATATGGTCGGCGGACGTCCATc                                                                                                                                               >  1:84092/1‑251 (MQ=255)                                                              aaaaGGCAGCTTCCTCAAATTCTGCCAACTGCTTGGCGGGCGTTCGGTCACCGAGTTCAACTACCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAAAAATTTTGCAGATGCTCAACGACGGCGGCTACAGCGTGGTTGATCGCTCCGACGACGAAATGGCGAAGCTACACGTGCGCTATATGGGCGGCGGACGTCCATCGCATCCGTTGc                                                                                                                                    <  1:143946/251‑1 (MQ=255)                                                                     aGCTTCCTCAAATTCTGCCAACTGCTTGGCGGGCGTTCGGTCACCGAGTTCAACTACCGTCTTGCAGATGCCAAAAACGCCTGCATCTTTGTCGGTGAGCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAAAAATTTTGCAGATGCTCAACGACGGCGGCTACAGCGTGGTTGCTCTATCCGCCGACGAAATGGCGAAGCTACACGTGCGCTATATGGTCGGCGGACGTCCATCGCATCCGTTGCAGGAACg                                                                                                                             >  1:163401/1‑251 (MQ=255)                                                                                                               aCCGAGTTCAACTACCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTAGGTGTGCGCCTGCGCCGCGGCCTCGAAGAGCGCAAAAAAATTTTGCAGATGCTCAACGACGGCGGCTACAGCGTGGTTGATCTCTCCGACGACGAAATGGCGAAGCTACACGTGCGCTATATGGTCGGCGGACGTCCATCGCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCGCgctgct                                                                                   >  1:149481/1‑251 (MQ=255)                                                                                                                 cGAGTTCAACTACCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAAAAATTTTGCAGATGCTCAACGACGGCGGCTACAGCGTGGTTGATCTCTCCGACGACGAAATGGCGAAGCTACACGTGCGCTATATGGTCGGCGGACGTCCATCGCATCCGTTGCAGGAACGCCTCTTCAGCTTCGAATTCCCGGAATCACCGGGCGCGCtgctgc                                                                                 <  1:4686/251‑1 (MQ=255)                                                                                                                                                                               gCCGCGGCCTCGAAGAGCGCAAAAAAATTTTGCAGATNCTCANCNNCGNCNGCNACAGc                                                                                                                                                                                                                   <  2:80004/59‑1 (MQ=255)                                                                                                                                                                                          gAAGAGCGCAAAAAACTTTTGCAGATGCTCAACGACGGCAGCTACAGCGTAGTTGATCTCTCCAACGCCGAAATGGCGAAGCTACACGTGAGCTATATGGTCGGCGGAAGTCCATCGCATCCGTTGAAGGAACGCATCTACAGCTTCGAATTCCCGGACTCACCGGGCGCGCTGCTGCGCTTCCTCAACACACTGGGTACGTACTGGAACATTTCTTTGTTCCACTATCGCAGCAATGGCACCGACTACgg        >  1:115891/1‑251 (MQ=255)                                                                                                                                                                                                 ncAAAAAAATTTTGCAGATGCTCAACGACGGCGGCTACAGCGTGGTTGATCTCTCCGACGCCGAACTGGCGAAGCTACACGTGCGCTATATGGTCGGCGGACGTCCATCGCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTCCTGGAACATTTCTTTGTTCAACTATCGCAGCCATGGCAACGACTACGGGCGCGt  >  1:16981/2‑250 (MQ=255)                                                                                                                                                                                                 gcAAAAAAATTTTGCAGATGCTCAACGACGGCGGCTAAAGCGTGGTTGATCTCTCCGACGACGAAATGGCGAAGCTACACGTGCGCTCTATGGTCGGCGGACGTCCATCGCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCAGGCATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTGCGTACTGGAGCATTTCTTTGTTcc                                    >  1:69163/1‑216 (MQ=255)                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     CAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGACCATTCCGGAAGAAAAAGGCAGCTTCCTCAAATTCTGCCAACTGCTTGGCGGGCGTTCGGTCACCGAGTTCAACTACCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAGAAATTTTGCAGATGCTCAACGACGGCGGCTACAGCGTGGTTGATCTCTCCGACGACGAAATGGCGAAGCTACACGTGCGCTATATGGTCGGCGGACGTCCATCGCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTCCACTATCGCAGCCATGGCACCGACTACGGGCGCGT  >  NC_000913/3956313‑3956754

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