Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,369,674	C→T	E247K (GAA→AAA)	yjeH ←	L‑methionine/branched chain amino acid exporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,369,674	0	C	T	100.0%	32.3 / NA	13	E247K (GAA→AAA)	yjeH	L‑methionine/branched chain amino acid exporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (10/3);  total (10/3)

GCACCGCATTGAGGGCATTATTCGGGATATGGCGAGAAGAGAGGCGTGCCAGGTAGTGGT‑CAGGATTATGTTGCGCCTGCGACCAG‑ACCAGGCGGGCGAAGCTCTGTATATAAATGTTGAGACTGGCAAAGCAGGCCAGATAGCCAATCACGCAGGCAATCCATAACGCTCCTACACCGAACAACTGCACTACAATTTTTGGAAGCGATGCTGCCGCCGCCATTTTTTCACCATAGGCGTCGAAGTGTAAGACGACTACCGTACAGCCCCAGTAGACTAATCCTGCCAGCAGCAGACCAATCATCAAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATC  >  NC_000913/4369446‑4369837                                                                                                                                                                                                                                       |                                                                                                                                                                    gCACCGCATTGAGGGCATTATTCGGGATATGGCGAGACGAGAGGCGTGCCAGGTAGTGGTACAG‑CTTATGTTGCGCCTGCGACCCGCACC‑GGCGGGCGAAGCCCTGTATTTAAATGTTGAGACTGGCACAACAGGCCACATAGCCCATCAAGCCGGCAATCCATAACGCTCCTAACCCCAACAACTGCACTACAATTTTTTGAACCGATCCTTCCGACGCCCTTTTTTTACCCTAGGCGTCGacgattaa                                                                                                                                                >  1:125539/1‑243 (MQ=255)                    attCGGGATATGGCGAGAAGAGAGGCGTGCCAGGTAGTGGT‑CAGGATTATGTTGCGCCTGCGACCAG‑ACCAGGCGGGCGAAGCTCTGTATATAAATGTTGAGACTGGCAAAGCAGGCCAGATAGCCAATCACGCAGGCAATCCATAACGCTCCTACACCGAACAACTGCACTACAATTTTTGGAAGGGATGCTGCCGCCGCCATTTTTGTACCATAGGCGTCGAAGTGTAAGACGACTACCGTACAGccc                                                                                                                             <  1:93103/250‑1 (MQ=255)                                                 ccAGGTAGTGGT‑CAGGATTATGTTGCGCCTGCGACCAG‑ACCAGGCGGGCGAAGCTCTGTATATAAATGTTGGGACTGGCAAAGCAGGCCAGATAGCCAATCACGAAGGCAATCCATAACGCTCCTCCCCCGAACAACTGCCCTACAATTTTTGGAAGCGATGCTGCCGCCGCCATTTTTTTACCCTAGGCCTCGCAGTGTa                                                                                                                                                 >  1:78240/1‑201 (MQ=255)                                                                                                    gAAGCTCTGTATATAAATGTTGAGACTGGCAAAGCAGGACAGATAGCCACTCACGAAGGCAATCCATAACGCTCCTACACCACACAACTGCACTAAAATTTTTGGACGCGATGCTGCCGCCGCAATTTTTTTACCATAGGCGTCGAAGTGTAAGACGACTACCGTACAGCCCCAGTAGACTAATCCTGCcagcagcag                                                                                                   >  1:186792/1‑198 (MQ=255)                                                                                                                                         gCCAGATAGCCAATCACGCAGGCAATCCATAACGCTCCTACACCGAACAACTGCACTACAATTTTTGGAAGCGATGCTGCCGCCGCCATTTTTTTACCATAGGCGTCGAAGTGTAAGACGACTACCGTACAGCCCCAGTAGACTAATCCTGCCAGCAGCAGACCAATCATCAAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACAcc          >  1:93755/1‑250 (MQ=255)                                                                                                                                          ncAGATAGCCAATCACGCAGGCAATCCATAACGCTACTACAGCGAACAACTGCACTACAATTTTTGGAAGCGCTGCTGCCGCCGCCATTTTTTTACCATAGGCGTCGAAGTGTAAGACGACTACCGTACAGCCCCAGTAGACTCATCCTGCCAGCAGCAGACCAATCATCAAAGCACg                                                                                 >  1:25539/2‑178 (MQ=255)                                                                                                                                                               cAATCCATAACGCTCCTACACCGAACAACTGCACTACAATTTTTGGAAGCGATGCTGCCGCCGCAATTTTTTTACCATAGGCGTCGAAGTGTAGgacgac                                                                                                                                          >  1:239853/1‑100 (MQ=255)                                                                                                                                                                                                           ggAAGCGATGCTGCCGCCGCCATTTTTTTACCATAGGCGTCGAAGTGTAAGCCGACTCCAGTACAACCCCAGTAGACTTATCATGCCAGAAGCCGACAAATCATCAAAGCCCGAGGAAACTCACGctct                                                                 >  1:222986/1‑129 (MQ=255)                                                                                                                                                                                                             aaGCGATGCTGCCGCCGCCATTTTTTTACCATAGGCGTCGAAGTGTAAGACGACTACCGTACAGCCCCAGTAGACTAAt                                                                                                                 >  1:22308/1‑79 (MQ=255)                                                                                                                                                                                                             aaGCGATGCTGCCGCCGCCATTTTTTTACCATAGGCGTCGAAGTGTAAGACGACTACCGTACAGCCCCAGTAGACTAAt                                                                                                                 <  2:22308/79‑1 (MQ=255)                                                                                                                                                                                                              aGCGATGCTGCCGCCGCCATTTTTTTACCATAGGCGTCGAAGTGTAAGTCGACTACCGTACAGCCCCAGTAGAc                                                                                                                     >  2:16636/1‑74 (MQ=255)                                                                                                                                                                                                              aGCGATGCTGCCGCCGCCATTTTTTTACCATAGGCGTCGAAGTGTAAGACGACTACCGTACAGCCCCAGTAGAc                                                                                                                     <  1:16636/74‑1 (MQ=255)                                                                                                                                                                                                                          gccgccATTTTTTTACCATAGGCGTCGAAGTGTAAGACGACTACCGTACAGCCCCAGTAGACTAATCCTGCCAGCAGCAGACCAATCATCAAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGCCAAAACACCAGAACATc  >  1:62349/1‑177 (MQ=255)                                                                                                                                                                                                                                       |                                                                                                                                                                    GCACCGCATTGAGGGCATTATTCGGGATATGGCGAGAAGAGAGGCGTGCCAGGTAGTGGT‑CAGGATTATGTTGCGCCTGCGACCAG‑ACCAGGCGGGCGAAGCTCTGTATATAAATGTTGAGACTGGCAAAGCAGGCCAGATAGCCAATCACGCAGGCAATCCATAACGCTCCTACACCGAACAACTGCACTACAATTTTTGGAAGCGATGCTGCCGCCGCCATTTTTTCACCATAGGCGTCGAAGTGTAAGACGACTACCGTACAGCCCCAGTAGACTAATCCTGCCAGCAGCAGACCAATCATCAAAGCACGAGGAAAATCACGCTCTGGATTTTTAAATTCCGAGGCGAGATGGGCAAATGCCTCCAGACCGACAAAACACCAGAACATC  >  NC_000913/4369446‑4369837

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