Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,506,221	G→A	R312K (AGA→AAA)	nanY →	2,7‑anhydro‑N‑acetylneuraminate hydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,506,221	0	G	A	100.0%	24.9 / NA	10	R312K (AGA→AAA)	nanY	2,7‑anhydro‑N‑acetylneuraminate hydratase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (2/8);  total (2/8)

CAGCAGGGTCACTTAGGATTGGCGGTCAGACAAAGCATTTTTTGGTCCATGAAACACAAGAAGAAGATGATGATCGTCGGAAAGGCAATATGACCTCAGAAATGGATGGCGCTATAGCATATGGTCATCCAGGAAAAAAAACACCATTATGGCTTGCCAGTTTAATAAGAAAGGAGACGTTATTCCTCCATAATATCCTCTGTGGTGCAAAACCTGAAGAAGATTATATTGACCTTCTCAATGGTGAGGCGGCCATGTCGGCGATTGCTACTGCTGATGCTGCCACTCTTTCAAGATCGCAGGACAGGAAAGTGAAAATCAGTGAGATCATTAAACATACATCAGTAATGTAAAAACTATAC  >  NC_000913/4506053‑4506414                                                                                                                                                                         |                                                                                                                                                                                                  cagcagGGTCACTTAGGATTGGCGGTCAGACAAAGCATTTTT‑GGTCCATGAAACACAAGAAGAAGATGATGATCGTCGGAAAGGCAATATGACCTCAGAAATGGATGGCGCTATAGCATATGGTCATCCAGGAAAAAAAACACCATTATGGCTTGCCAGTTTAATAAAAAAGGAGACGTTATTCCTCCATAATATCCTCTGTGGTGCAAAACCTGAAGAAGATTATATTGACCTTCTCAATGGTGAggcgn                                                                                                                <  1:63725/251‑2 (MQ=255)          cACTTAGGATTGGCGGTCAGACAAAGCATTTTTTGGTCCATGAAACACAAGAAGAAGATGATGATCGTCGGAAAGGCAATATGACCTCAGAAATGGATGGCGCTATAGCATATGGTCATCCAGGAAAAAAAACACCATTATGGCTTGCCAGTTTAATAAAAAAGGAGACGTTATTCCTCCATAATATCCTCTGTGGTGCAAAACCTGAAGAAGATTATATTGACCTGCTCAATGGTGAggcggn                                                                                                               <  1:77283/244‑2 (MQ=255)                 gATTGGCGGTCAGACAAAGCATTTTTTTGTCCATGAAACACAAGAAGAAGATGATGATCGTCGGAAAGGCAATATGACCTCAGAAATGGATGGCGCTATAGCATATGGTCATCCAGGAAAAAAAACACCATTATGGCTTGCCAGTTTAATAAAAAAGGAGACGTTATTCCTCCATAATATCCTCTGTGGTGCAAAACCTGAAGAAGATTATATTGACCTTCTCAATGGTGAGGCGGCCATGTCGGCGAttn                                                                                                 <  1:160485/251‑2 (MQ=255)                           cAGACAAAGCATTTTTTGGTCCATGAAACACAAGAAGAAGATGATGATCGTCGGAAAGGCAATCTGACCTCAGAAATGGATGGCGCTATAGCATATGGTCATACAGGAAAACAAACACCATTATGGATTGCTCGTTTAATAAAAAAGGAGACGTTATTCCTCCATAATATCCTCTGTGGTGCAAAAcc                                                                                                                                                      >  1:181892/1‑188 (MQ=255)                                   gCATGTTTTGGTCCATGAAACACAAGAAGAAGATGATGATCGTCGGAAAGGCAATATGACCTCAGAAATGGATGGCGCTATAGCATATGGTCATCCAGGAAAAAAAACACCAGTATGGCTTGCCAGTTTAATAAAAAAGGAGACGTTATTCCTCCATAATATCCTCTGTGGTGCAAAACCTGAAGAAGATTATATTGACCTtctc                                                                                                                             <  1:78458/205‑1 (MQ=255)                                             gttcATGAAACACAAGAAGAAGATGATGATCGTCGGAAANGCCATATTANNTCNGAAAGNGNNNGNNNTANAGCATATGGACATCCAGGANAANNAACACCNTTATGGCTTGCCAGTTTAATaaaaaa                                                                                                                                                                                                <  2:60552/125‑1 (MQ=255)                                             gTCCATGACACACACGAAGAAGATGATGATCGTCGGAACGGCAATATGACATCAGAAATGGATGGCGCTATAGCATATGGTCATCCCGGAAAAAAACCACCATTATGGCTTGCCCGTTTAATaaaaaa                                                                                                                                                                                                >  1:60552/1‑128 (MQ=255)                                                                           cgtcgGAAAGGCAATATGACCTCAGAAATGGATGGCGCTATAGCATATGGTCATCCAGGAAAAAAAACACCATTATGGCTTGCCAGTTTAATAAAAAAGGAGACGTGATTCCTCCATAATATCCTCTGTGGTGCAAAAc                                                                                                                                                       <  1:88624/139‑1 (MQ=255)                                                                                                                cTATAGCATATGGTCATCCAGGAAAAAAAACACCATTATGGCTTGCCAGTTTAATAAAAAAGGAGACGTTATTCCTCCATAATATCCTCTGTGGTGCAAAACCTGAAGAAGATTATATTGACCTTCTCAATGGTGAGGCGGCCATGTCGGCGATTGCTCCTGCTGATGCTGCCACTCTTTCAAGATCGCAGGACAGGACAGTGAAAATCAGTGAGATCATTAAACATACATCAGTAATGTAAAAACTATAc  >  1:128060/1‑251 (MQ=255)                                                                                                                        tatGGTCATCCAGGAAAAAAAACACCATTATGGCTTGCCAGTTTAATAAAAAAGGAGACGTTATTCCTCCATAATATCCTCTGTGGTGCAAAACCTGAAGAAGATGATATTGACCTTCTCAATGGTGTGGCGGCCATGTCGTCGATTGCTACTGCTGATGCTGCCACTCTTTCAAGATcn                                                                 <  1:53208/180‑2 (MQ=255)                                                                                                                                       aaaaaaaCACCATTATGGCTTGCCAGTTTAATAAAAAAGGAGACGTTATTCCTCCATAATATCCTCTGTGGTGCAAAACCTGAAGAAGATTATATTGACCTTCTCAATGGTGAGGCGGCCATGTCGGCGATTGCTACTGCTGATGCTGCCACTCTTTCAAGATCGCAGGACAGGAAAGTGAAAATCAGTGAGATCATTAAACATACATCAGTAATGTAAAAACtatan  <  1:30612/228‑2 (MQ=255)                                                                                                                                       aaaaaaaCACCATTATCGCTTGCCAGTTTAATAAAAAAGGAGAGGATATTCCGCCATAATATCCTCTTTGGTGCAAAACCTGAAGAAGATtata                                                                                                                                        <  1:117391/94‑1 (MQ=255)                                                                                                                                                                         |                                                                                                                                                                                                  CAGCAGGGTCACTTAGGATTGGCGGTCAGACAAAGCATTTTTTGGTCCATGAAACACAAGAAGAAGATGATGATCGTCGGAAAGGCAATATGACCTCAGAAATGGATGGCGCTATAGCATATGGTCATCCAGGAAAAAAAACACCATTATGGCTTGCCAGTTTAATAAGAAAGGAGACGTTATTCCTCCATAATATCCTCTGTGGTGCAAAACCTGAAGAAGATTATATTGACCTTCTCAATGGTGAGGCGGCCATGTCGGCGATTGCTACTGCTGATGCTGCCACTCTTTCAAGATCGCAGGACAGGAAAGTGAAAATCAGTGAGATCATTAAACATACATCAGTAATGTAAAAACTATAC  >  NC_000913/4506053‑4506414

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