Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,718,641	G→A	A321V (GCC→GTC)	anmK ←	anhydro‑N‑acetylmuramic acid kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,718,641	0	G	A	100.0%	60.8 / NA	19	A321V (GCC→GTC)	anmK	anhydro‑N‑acetylmuramic acid kinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (11/8);  total (11/8)

AATAGCCCCCAGTACCGTCTCCTGGCTTGCGCCAGTGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATGCCGCCAGACGCGCCATGAGTAGCGGATTACGACTACCTCCACCACATACCATCAATCGTTCGCAGCCACCGCTCAACAAAACTTGTTCAGAAATGGT  >  NC_000913/1718511‑1718774                                                                                                                                   |                                                                                                                                      aaTAGCCCCCAGTACCGTCTCCTGGCTTGCGCCAGTGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGACATCggtg                                                                                                                               <  2:505852/139‑1 (MQ=255)                 gTCTCCTGGCTTGCGCCAGTGACGGAAGGCAGATTTCCTGTTAATCCCGCCAGCGTCCGCCAAGCAAGCCAGGCGAAAGCCAATGCTTCCAAGTCATCGCCACTAATGCCAACGACATCGGTGGTGGTGACTTCTGTGc                                                                                                               >  2:415997/1‑139 (MQ=255)                    tcCTGGCTTGCGCCAGTGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGACATCGGTGGTGGTGACTTCTGTGCCTg                                                                                                            >  2:231840/1‑139 (MQ=255)                              cgcCAGTGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGACATCGGTGGTGGTGACTTCtg                                                                                                                  <  1:315178/123‑1 (MQ=255)                              cgcCAGTGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGACATCGGTGGTGGTGACTTCtg                                                                                                                  >  2:315178/1‑123 (MQ=255)                                                   ttCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGACATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATGCCGCCAGACGCGCCATGAGTAg                                                                             <  2:8428/139‑1 (MQ=255)                                                       tGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGACATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATGCCGCCAGACGCGCCATGAGTAGCGGa                                                                         >  2:675189/1‑139 (MQ=255)                                                             tCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGACATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATGCCGCCAGACGCGCCATGAGTAGCGGATTACGa                                                                   <  1:231840/139‑1 (MQ=255)                                                                        tCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGACATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATGCCGCCAGAc                                                                                          >  1:93294/1‑105 (MQ=255)                                                                        tCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGACATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATGCCGCCAGAc                                                                                          <  2:93294/105‑1 (MQ=255)                                                                           gCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGACATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATGCCGCCAGACGCGCCATGAGTAGCGGATTACGACTACCTCCACcaca                                                     <  1:675189/139‑1 (MQ=255)                                                                                          cGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGACATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATGCCGCCAGACGCGCCATGAGTAGCGGATTACGACTACCTCCACCACATACCATCAATCGTTc                                      >  2:164181/1‑139 (MQ=255)                                                                                                                  cGCCACTAATGCCAACGACATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATGCCGCCAGACGCGCCATGAGTAGCGGATTACGACTACCTCCACCACATACCATCa                                             >  1:420181/1‑108 (MQ=255)                                                                                                                  cGCCACTAATGCCAACGACATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATGCCGCCAGACGCGCCATGAGTAGCGGATTACGACTACCTCCACCACATACCATCa                                             <  2:420181/108‑1 (MQ=255)                                                                                                                  cGCCACTAATGCCAACGACATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATGCCGCCAGACGCGCCATGAGTAGCGGATTACGACTACCTCCACCACATACCATCAATCGTTCGCAGCCACCGCTCAACAAAACTTg              >  1:458984/1‑139 (MQ=255)                                                                                                                  cGCCACTAATGCCAACGACATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATGCCGCCAGACGCGCCATGAGTAGCGGATTACGACTACCTCCACCACATACCATCAATCGTTCGCAGCCACCGCTCAACAAAACTTg              >  2:49663/1‑139 (MQ=255)                                                                                                                  cGCCACTAATGCCAACGACATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATGCCGCCAGACGCGCCATGAGTAGCGGATTACGACTACCTCCACCACATACCATCAATCGTTCGCAGCCACCGCTCAACAAAACTTg              >  2:694968/1‑139 (MQ=255)                                                                                                                        tAATGCCAACGACATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATGCCGCCAGACGCGCCATGAGTAGCGGATTACGACTACCTCCACCACATACCATCAATCGTTCGCAGCCACCGCTCAACAAAACTTGTTCAGa        >  1:359561/1‑139 (MQ=255)                                                                                                                         aaTGCCAACGACATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATGCCGCCAGACGCGCCATGAGTAGCGGATTACGACTACCTCCACCACATACCATCAATCGTTCGCAGCCACCGCTCAACAAAACTTGTTCAGaa       <  2:458984/139‑1 (MQ=255)                                                                                                                              cAACGACATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATGCCGCCGGACGCGCCATGAGTGGCAGATTACGACTACCTCCACCACATACCATCAATCGTTCGCATCCACTACTCATCCAATCTTGTTCATAAATGGt  <  1:415997/139‑1 (MQ=255)                                                                                                                                   |                                                                                                                                      AATAGCCCCCAGTACCGTCTCCTGGCTTGCGCCAGTGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATGCCGCCAGACGCGCCATGAGTAGCGGATTACGACTACCTCCACCACATACCATCAATCGTTCGCAGCCACCGCTCAACAAAACTTGTTCAGAAATGGT  >  NC_000913/1718511‑1718774

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