Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	721,372	G→A	A121V (GCC→GTC)	kdpE ←	DNA‑binding transcriptional activator KdpE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	721,372	0	G	A	100.0%	114.9 / NA	36	A121V (GCC→GTC)	kdpE	DNA‑binding transcriptional activator KdpE
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (17/19);  total (17/19)

GAGCAGCACCGCCAGCAGGCGGAACTCAATTGGTGTGAGATGCACCTCTTCCTCACCCCGGTGAATCACGCGGGCGGCTAAATCGACGGTAACATCGGAAAATTTTACCAGCGGATCGGGCGCGGTGGTGGCAGAGTGGCGGCGTAATGCGACGCGCAGACGGGCCTGCAATTCGCCAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGG  >  NC_000913/721242‑721498                                                                                                                                   |                                                                                                                               gAGCAGCACCGCCAGCAGGCGGAACTCAATTGGTGTGAGATGCACCTCTTCCTCACCCCGGTGAATCACGCGGGCGGCTAAATCGACGGTAACATCGGAAAATTTTACCAGCGGATCGGGCGCGGTGGTGACa                                                                                                                              <  1:394116/133‑1 (MQ=255)                            aaTTGGTGTGAGATGCACCTCTTCCTCACCCCGGTGAATCACGCGGGCGGCTAAATCGACGGTAACATCGGAAAATTTTACCAGCGGATCGGGCGCGGTGGTGACAGAGTGGCGGCGTAATGCGACGCGCAGa                                                                                                   <  1:432242/133‑1 (MQ=255)                                 gtgtGAGATGCACCTCTTCCTCACCCCGGTGAATCACGCGGGCGGCTAAATCGACGGTAACATCGGAAAATTTTACCAGCGGATCGGGCGCGGTGGTGACAGCGTGGCGGGGTATTGGGACGCGGAGACGGGc                                                                                              >  1:165954/1‑133 (MQ=255)                                 gtgtGAGATGCACCTCTTCCTCACCCCGGTGAATCACGCGGGCGGCTAAATCGACGGTAACATCGGAAAATTTTACCAGCGGATCGGGCGCGGTGGTGACAGAGTGGCGGCGTAATGCGACGCGCAGACGGGc                                                                                              >  2:116714/1‑133 (MQ=255)                                 gtgtGAGATGCACCTCTTCCTCACCCCGGTGAATCACGCGGGCGGCTAAATCGACGGTAACATCGGAAAATTTTACCAGCGGATCGGGCGCGGTGGTGACAGAGTGGCGGCGT‑ATGCGACGCGCAGACGGGcc                                                                                             >  2:381697/1‑133 (MQ=255)                                      agaTGCACCTCTTCCTCACCCCGGTGAATCACGCGGGCGGCTAAATCGACGGTAACATCGGAAAATTTTACCAGCGGATCGGGCGCGGTGGTGACAGAGTGGCGGCGTAATGCGACGCGCAGACGGGCCTGCa                                                                                         >  2:442307/1‑133 (MQ=255)                                      agaTGCACCTCTTCCTCACCCCGGTGAATCACGCGGGCGGCTAAATCGACGGTAACATCGGAAAATTTTACCAGCGGATCGGGCGCGGTGGTGACAGAGTGGCGGCGTAATGCGACGCGCAGACGGGCCTGCa                                                                                         >  1:373644/1‑133 (MQ=255)                                      agaTGCACCTCTTCCTCACCCCGGTGAATCACGCGGGCGGCTAAATCGACGGTAACATCGGAAAATTTTACCAGCGGATCGGGCGCGGTGGTGACAGAGTGGCGCCGTGATGCGACGTGCAGTCGGGCCTGCa                                                                                         >  2:281588/1‑133 (MQ=255)                                       gaTGCACCCCTTCCTGACTCCGGTGAATCAAGCGGGCGGCTAAATCGACGGTAACATCGGAAAATTTTACCAGCGGATCGGGCGCGGTGGTGACAGAGTGGCGGCGTAATGCGACGCGCAGACGGGCCTGCaa                                                                                        <  1:247913/133‑1 (MQ=255)                                            agctctTCCACACCCAGGTGAATGCCGCGGGCAGCTAAATCCACGGTAACATCGGAAAATTTTACCAGCGGATCGGGCGCGGTGGTGACAGAGTGGCGGCGTAATGCGACGCGCAGACGGGCCTGCAATTCGc                                                                                   <  2:298000/131‑1 (MQ=255)                                            aCCTCTTCCTCACCCCGGTGAATCACGCGGGCGGCTAAATCGACGGTAACATCGGAAAATTTTACCAGCGGATCGGGCGCGGTGGTGACAGCGTGGCGGCATAAGGCGACGCGCAGAAGGGCCTGCAATTCGc                                                                                   >  2:247913/1‑133 (MQ=255)                                            aCCTCTTCCTCACCCCGGTGAATCACGCGGGCGGCTAAATCGACGGTAACATCGGAAAATTTTACCAGCGGATCGGGCGCGGTGGTGACAGAGTGGCGGCGTAATGCGACGCGCAGACGGGCCTGCAAATCGc                                                                                   >  1:619301/1‑133 (MQ=255)                                              ctctTCCTCTCCCCGGTGAATCACGCGTGCGGCTAAATCGACGGTAACATCGGAAAATTTTACCAGCGGATCGGGCGCGGTGGTGACAGAGTGGCGGCGTAATGCGACGCGCAGACGGGCCTGCAATTCGCCa                                                                                 <  1:395890/133‑1 (MQ=255)                                                ctTCCTCACCCCGGTGAATCACGCGGGCGGCTAAATCGACGGTAACATCGGAAAATTTTACCAGCGGATCGGGCGCGGTGGTGACAGAGTGGCGGCGTAATGCGACGCGCAGACGGGCCTGCAATTCGCCAAt                                                                               >  1:87515/1‑133 (MQ=255)                                                     tCACCCCGGTGAATCACGCGGGCGGCTAAATCGACGGTAACATCGGAAAATTTTACCAGCGGATCGGGCGCGGTGGTGACAGAGTGGCGGCGTAATGCGACGCGCAGACGGGCCTGCAATTCGCCAATGCCaa                                                                          <  1:452324/133‑1 (MQ=255)                                                      cACCCCGGTGAATCACGCGGGCGGCTAAATCGACGGTAACATCGGAAAATTTTACCAGCGGATCGGGCGCGGTGGTGACAGAGTGGCGGCGTAATGCGACGCGCAGACGGGCCTGCAATTCGc                                                                                   <  2:583348/123‑1 (MQ=255)                                                      cACCCCGGTGAATCACGCGGGCGGCTAAATCGACGGTAACATCGGAAAATTTTACCAGCGGATCGGGCGCGGTGGTGACAGAGTGGCGGCGTAATGCGACGCGCAGACGGGCCTGCAATTCGc                                                                                   >  1:583348/1‑123 (MQ=255)                                                        ccccGGTGAATCACGCGGGCGGCTAAATCGACGGTAACATCGGAAAATTTTACCAGCGGATCGGGCGCGGTGGTGACAGAGTGGCGGCGTAATGCGACGCGCTGACGGGCGTGCCATACGCCAAGGCCAAACg                                                                       >  1:256599/1‑133 (MQ=255)                                                             gTGAATCATGCGGGCTGCTAAACCGACGGTCACATCGGAAAATTTTACCAGCGGATCGGGCGCGGTGGTGACAGAGTGGCGGCGTAATGCGACGCGCAGACGGGCCTGCAATTCGCCAATGCCAAACGGCTTa                                                                  <  1:184409/133‑1 (MQ=255)                                                             gTGAATCACGCGGGCGGCTAAATCGACGGTAACATCGGAAAATTTTACCAGCGGATCGGGCGCGGTGGTGACAGAGTGGCGGCGTAATGCGGCGCGCAGAGGGGACTGGAATTCGCCTATGCAAAACGGCTTa                                                                  >  2:176459/1‑133 (MQ=255)                                                             gTGAATCACGCGGGCGGCTAAATCGACGGTAACATCGGAAAATTTTACCAGCGGATCGGGCGCGGTGGTGACAGAGTGGCGGCGTAATGCGACGCGCAGACGGGCCTGCAATTCGCCAATGCCAAACGGCTTa                                                                  <  1:146278/133‑1 (MQ=255)                                                             gTGAATCACGCGGGCGGCTAAATCGACGGTAACATCGGAAAATTTTACCAGCGGATCGGGCGCGGTGGTGACAGAGTGGCGGCGTAATGCGACGCGCAGACGGGCCTGCAATTCGCCAATGCCAAACGGCTTa                                                                  <  1:442307/133‑1 (MQ=255)                                                              tGAATCACGCGGGCGGCTAAATCGACGGTAACATCGGAAAATTTTACCAGCGGATCGGGCGCGGTGGTGACAGAGTGGCGGCGTAATGCGACGCGCAGACGGGCCTGCAATTCGCCAATGCCAAACGGCTTAc                                                                 <  1:480827/133‑1 (MQ=255)                                                               gAATCACGCGGGCGGCTAAATCGACGGTAACATCGGAAAATTTTACCAGCGGATCGGGCGCGGTGGTGACAGAGTGGCGGCGTAATGCGACGCGCAGACGGGCCTGCAATTCGCCAATGCCAAACGGCTTACt                                                                <  2:160151/133‑1 (MQ=255)                                                                   cACGCGGGCGGCTAAATCGACGGTAACATCGGAAAATTTTACCAGCGGATCGGGCGCGGTGGTGACAGAGTGGCGGCGTAATGCGACGCGCAGACGGGCCTGCAATTCGCCAATGCCAAACGGCTTACTCAGa                                                            >  2:349448/1‑133 (MQ=255)                                                                            ggcTAAATCGACGGTAACATCGGAAAATTTTACCAGCGGAGCGGGCGCGGTGGTGACAGAGTGGCGGCGT‑ATGCGACGCGCAGACGGGCCTGCAATTCGCCAATGCTAAACGGCTTACTCAGATAATCATCCg                                                  <  1:381697/133‑1 (MQ=255)                                                                             gcTAAATCGACGGTAACATCGGAAAATTTTACCAGCGGATCGGGCGCGGTGGTGACAGAGTGGCGGCGTAATGCGACGCGCAGACGGGCCTGCAATTCGCCAATGCCAAACGGCTTACTCAGATAATCATCCg                                                  <  1:116714/133‑1 (MQ=255)                                                                                      aCGGTAACATCGGAAAATTTTACCAGCGGATCGGGCGCGGTGGTGACAGAGTGGCGGCGTAATGCGACGCGCAGACGGGCCTGCAATTCGCCAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCAt                                         >  2:56824/1‑133 (MQ=255)                                                                                             cATCGGAAAATTTTACCAGCGGATCGGGCGCGGTGGTGACAGAGTGGCGGCGTAATGCGACGCGCAGACGGGCCTGCAATTCGCCAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAgcgc                                  >  1:27833/1‑133 (MQ=255)                                                                                                  gAAAATTTTACCAGCGGATCGGGCGCGGTGGTGACAGAGTGGCGGCGTAATGCGACGCGCAGACGGGCCTGCAATTCGCCAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGa                             >  2:353670/1‑133 (MQ=255)                                                                                                       ttttACCAGCGGATCGGGCGCGGTGGTGACAGAGTGGCGGCGTAATGCGACGCGCAGACGGGCCTGCAATTCGCCAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTg                        <  1:353670/133‑1 (MQ=255)                                                                                                         ttACCAGCGGCTCGGACGCGGTGGCGACTGAGCGGCGGCGTCATGCGACGCGCAGACGGGCCTGCAATTCGCCAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTc                      <  2:165954/133‑1 (MQ=255)                                                                                                         ttACCAGCGGATCGGGCGCGGTGGTGACAGAGTGGCGGCGTAATGCGACGCGCAGACGGGCCTGCAATTCGCCAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTc                      <  1:149753/133‑1 (MQ=255)                                                                                                            ccAGCGGATCGGGCGTGGAGGTTACAGAGTCGCGGCGTAATGCGACGCGCAGACGGGCCTGCAATTCGCCAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGct                   <  1:176459/133‑1 (MQ=255)                                                                                                                         cgcgGTGGTGACAGAGTGGCGGCGTAATGCGACGCGCAGACGGGCCTGCAATTCGCCAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCtgcgt      <  1:349448/133‑1 (MQ=255)                                                                                                                             gtggtgACAGAGTGGCGGCGTAATGCGACGCGCAGACGGGCCTGCAATTCGCCAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCgg  >  2:16357/1‑133 (MQ=255)                                                                                                                                   |                                                                                                                               GAGCAGCACCGCCAGCAGGCGGAACTCAATTGGTGTGAGATGCACCTCTTCCTCACCCCGGTGAATCACGCGGGCGGCTAAATCGACGGTAACATCGGAAAATTTTACCAGCGGATCGGGCGCGGTGGTGGCAGAGTGGCGGCGTAATGCGACGCGCAGACGGGCCTGCAATTCGCCAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGG  >  NC_000913/721242‑721498

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