Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,215,023	C→T	F248F (TTC→TTT)	metA →	homoserine O‑succinyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,215,023	0	C	T	100.0%	32.5 / NA	12	F248F (TTC→TTT)	metA	homoserine O‑succinyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (8/4);  total (8/4)

GATGATTCATTCCTGGCACCGCATTCGCGCTATGCTGACTTTCCGGCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTG  >  NC_000913/4214829‑4215216                                                                                                                                                                                                   |                                                                                                                                                                                                  natgatTCATTCCTGGCACCGCATTCGCGCTATGCTGACTTTCCGGCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTTCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGcaca                                                                                                                                            >  1:151559/2‑250 (MQ=255)        cattcCTGGCACCGCATTCGCGCTATGCTGACTTTCCGGCAGCGTTGATTCGTGATTACACGGATCTGGAAATTCTGGCAGAGACGGCAGAAGGGGATGCATCTCTGTTTGCCAGTAACGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGACTTTTTTCGCGATGTGGAAg                                                                                                                                                                                      >  1:184957/1‑201 (MQ=255)                gCACCGCATTCGCGCTATGCTGACTTTCCGGCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATTTGATGCGCAAACGCTGGCGCAGGAATTTTTTCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATac                                                                                                                            >  1:252403/1‑251 (MQ=255)                                                             gATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCAGAATATGATGCGCAAACGCTGGCGCAGGACTTTTTTCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCa                                                                                >  1:39053/1‑250 (MQ=255)                                                                                                         tGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGACCTCCCGGATATGATGCGCAAACGATGGCGCAGGAATTTTTTCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATc                                                                                                                                      >  1:123456/1‑152 (MQ=255)                                                                                                                     tGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAAAGCTGGCGCAGGAATTTTTTCGCGCTGTGGAAGCCGGCCTAGACCCGGCTGTACCGCATAACTATTTCCCGCACAATGATCCg                                                                                                                                    >  1:146935/1‑142 (MQ=255)                                                                                                                          gTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATACCGAAGATGATGAGCAAACGCTGGCGCAGGACTTTTTTCGCGATGTGGCAGcc                                                                                                                                                                                    >  1:171488/1‑89 (MQ=255)                                                                                                                                      gcATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTTCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTgn                                                                                                            <  1:55492/149‑2 (MQ=255)                                                                                                                                           gCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCGGGCGCAGGAATTTTTTCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCgcgagc                                                                                                                   <  1:114097/137‑1 (MQ=255)                                                                                                                                                           nnnnnnnnnntgnnnNNCAANCGCTGGCGCAGGAAGTTTTTCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAACTTACTGTTTAc                                                                                  <  2:151559/144‑1 (MQ=255)                                                                                                                                                                            cAAACGCTGGCGCAGGAATTTTTTCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCa                                                                       >  1:147750/1‑148 (MQ=255)                                                                                                                                                                                   tGGCGCAGGAATTTTTTCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGTCTCAACTAGTACGTCTACCAGATCACGCTATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCtg  <  1:119547/210‑1 (MQ=255)                                                                                                                                                                                                   |                                                                                                                                                                                                  GATGATTCATTCCTGGCACCGCATTCGCGCTATGCTGACTTTCCGGCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTG  >  NC_000913/4214829‑4215216

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