Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,131,375	G→A	R514H (CGC→CAC)	metL →	fused aspartate kinase/homoserine dehydrogenase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,131,375	0	G	A	100.0%	31.7 / NA	13	R514H (CGC→CAC)	metL	fused aspartate kinase/homoserine dehydrogenase 2
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (5/8);  total (5/8)

CACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGTATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGTTTGTGCTGGCAGGTGTGGTGGACAGCCGCCGCAGCCTGTTGAGCTATGACGGGCTGGACGCCAGCCGCGCGTTAGCCTTCTTCAACGATGAAGCGGTTGAGCAGGATGAAGAGTCGTTGTTCCTGTGGATGCGCGCCCATCCGTATGATGATTTAGTGGTGCTGGACGTTACCGCCAGCCAGCAGCTTGCTGATCAGTATCTTGATTTCGCCAGCCACGGTTTCCACGTTATCAGCGCCAACAAACTGGCGGGAGCCAGCGACAGC  >  NC_000913/4131145‑4131574                                                                                                                                                                                                                                       |                                                                                                                                                                                                        cACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGTATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGTTTGTGCTGGCAGGTGTGGTGGACAGCCGCCGCAGCCTGTTGAGCTATGACGGGCTGGACGCCAGCCACGCGTTAGCCTTCTTCAAcn                                                                                                                                                                                     <  1:32187/251‑2 (MQ=255)            ccgaAAGCCTGATTCAGGGGCTTCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGTATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGTTTGTGCTGGCAGGTGTGGTGGACAGCCGCCGCAGCCTGTTGAGCTATGACGGGCTGGACGCCAGCCACGCGTTAGCCGTCTTCAACGAGGAAGCGGtt                                                                                                                                                                          <  1:76256/251‑1 (MQ=255)                   ccTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAACCGCATCGGCCTGGTATTGTTAGGTACGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCAAGCTTTCGGCACGTACCGGCTTTGAGTTTGTGCTGGCAGGTGTGGTGGACAGACGCCGCAGCCTGTTGAGCTATGACGGGCTGGACGCCAGCCACGCGTTAGCCTTCTTCAACGATGAAGCGGTTGAGCAgg                                                                                                                                                                   >  1:46103/1‑251 (MQ=255)                       aTTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCGGGTATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGTTTGTGCTGGCAGGTGTGGTGGACAGCCGCCGCAGCCTGTTGAGCTATGACGGGCTGGACGCCAGCCACGCGTTAGCCTTCTTCAACGATGAAGCGGTTGAGCAGGATGa                                                                                                                                                               <  1:16151/251‑1 (MQ=255)                                          ccGTCTTCCGCGCAGAAAAACGCATCGGCCTGGTATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGTTTGTGCTGGCAGGTGTGGTGGACAGCCGCCGCAGCCTGTTGAGCTATGACGGGCTGGACGCCAGCCACGCGTTAGCCTTCTTCAACGATGAAGCGGTTGAGCAGGATGAAGAGGCGTTGTTCCtgtg                                                                                                                                             <  1:95756/250‑1 (MQ=255)                                                 ccGCGCAGAAAAACGCATCGGCCTGGTATTGTTCGGTAAGGGCAATATCGGTTCCCGGTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGTTTGTGCTGGCAGGTGTGGTGGACAGCCGCCGCAGCCTGTTGAGCTATGACGGGCTGGACGCCAGCCACGCGTTAGCCTTCTTCAACGATGAAGCGGTTGAGCAGGATGAAGAGTCGGTGTTCCTGTGGATgcgcg                                                                                                                                     <  1:283908/251‑1 (MQ=255)                                                                                          gCAATATCGGTTCCCGTTGGATGGAACTGTTCGCCCGTGAGCAGAGCAGGCTTTGGGCACGTACAGGCTTTGAGATTGTGCTGGCAGGAGTGGTGGAAAGCCGCCGCAGGCTGTTGCGCTATGACGGGCTGGAAGCAAGCCACGCGTTAGCCTTCTTCAACGATGAAGCGGTTGAGCAGGATGAAGAGTCGTTGTTCCTGTGGATGCGCGCCCATCCGGATGATGATTTAGTGGTGCTGGACGTTACCgcc                                                                                            >  1:340335/1‑251 (MQ=255)                                                                                           cAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCCGAGCACGCTTTCGGCACGTCCCGGCTTTGAGTTTGTGCTGGCAGGTGTGGTGGCCAGCCGCCGCAGCCTGTTGAGCTATGACGGGCTGGACGCCAGCCACGCGTTAGCCTTCTTCAACGATGAAGCGGTTGAGCAGGATGAAGAGTCGTTGTTCCTGTGGATGCGCGCCCATCCGTATGATGATTTAGTGGTGCTGGACGTTACCgcc                                                                                            >  1:348835/1‑250 (MQ=255)                                                                                             ataccgggtCCCGTTGGCTTGCCCTGTTCGCCCGTGAGCATAACTCGCTTTCTGCACGTCCCGTCTTTGAGTTTGTGCTGGCAGGGGTGGTGAGCAGCCGCCGCAGCCTGTTGAGCTATGACGGGCTGGACGCCATCCACGCGTTAGCCTTCTTCAACGATGAAGCGTTTGAGCAGGATTAAGAGTCGTTGTTCCTGTGGATGCGCGCCCACCCGTATGATGATTTAGTGGTGCTGGACGTTACCgccag                                                                                          <  1:256063/242‑1 (MQ=255)                                                                                                                                                        aCCGGCTTTGAGGTCGTGCCGGCAGGTGGGGTGGACAGCCCCCGCAGATTGTTGAGGTATTAGGGGCCGGACGCGAGCCACGCGTTTGCCTTCTTCAACGATGAAGCGGTTGAGCAGGATGAagag                                                                                                                                                           <  1:142005/126‑1 (MQ=255)                                                                                                                                                                                   gtggtgGACAGCCGCCGCCGCCTGTTGAGCTATGACGGGCTGGACGCCAGCCACGCGTTAGCCTTCTTCAACGATGAAGCGGTTGAGCAGGATGAAGAGTCGTTGTTCCTGTGGATGCGCGCCCATCCGTATGATGATTTAGTGGTGCTGGACGTTACCGCCAGCCAGCAGCTTGCTGATCAGTATCTTGATTTCGCCAGCCACGGTTTCCACGTTATCAGCGCCAACAAACTGGCGGGAGCCAGCGAc     >  1:98337/1‑249 (MQ=255)                                                                                                                                                                                     ggtggACAGCCGCCGCAGCCTGTTGAGCTATGACGGGCTGGACGCCAGCCACGCGTTAGCCTTCTTCAACGATGAAGCGGTTGAGCAGGATGAAGAGTCGTTGTTCCTGTGGATGCGCGCCCATCCGTATGATGATTTAGTGGTGCTGGCCGTTACCGCCAGCCAGCAGCTTGCTGATCAGTATCTTGATTTCGCCAGCCACGGTTTCCACGTTATCAGCGCCAACAAACTGGCGGGAGCCAGCGACAGc  >  1:234553/1‑250 (MQ=255)                                                                                                                                                                                                                            ggACGCCAGCCACGCGTTAGCCTTCTTCAACGATGAAGCGGTTGAGCAGGATGAAGAGTCGTTGTTCCTGTGGATGCGCGCCCATCCGTATGATGATTTAGTGGTGCTGGACGTTACCGCCAGCCAGCAGCTTGCTGATCAGTAGCTTGATTTCGCCAGCCACGGTTTcn                                           <  1:101942/170‑2 (MQ=255)                                                                                                                                                                                                                                       |                                                                                                                                                                                                        CACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGTATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGTTTGTGCTGGCAGGTGTGGTGGACAGCCGCCGCAGCCTGTTGAGCTATGACGGGCTGGACGCCAGCCGCGCGTTAGCCTTCTTCAACGATGAAGCGGTTGAGCAGGATGAAGAGTCGTTGTTCCTGTGGATGCGCGCCCATCCGTATGATGATTTAGTGGTGCTGGACGTTACCGCCAGCCAGCAGCTTGCTGATCAGTATCTTGATTTCGCCAGCCACGGTTTCCACGTTATCAGCGCCAACAAACTGGCGGGAGCCAGCGACAGC  >  NC_000913/4131145‑4131574

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