Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	721,604	G→A	R44C (CGT→TGT)	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,604	0	G	A	100.0%	94.5 / NA	32	R44C (CGT→TGT)	kdpE	DNA‑binding transcriptional activator KdpE
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/19);  total (13/19)

CGGTGGTGGCAGAGTGGCGGCGTAATGCGACGCGCAGACGGGCCTGCAATTCGCCAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCT‑CTTCGCTGCGTGCGGAAAGCACAATC‑ACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACGGGTTGCCGCTTCCAGCAAGCCGCGTTGCAG‑CGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGTTTGCTGGGGAAGTGTAACACGAAAACAGGCACCACCTTCCGGTCGGTTGAACGCGGTAATAGTGCCCCCGTGTACATC  >  NC_000913/721364‑721846                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    cgGTGGTGGCAGAGTGGCGGCGTAATGCGACGCGCAGACGGGCCTGCAATTCGCCAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCT‑CTTCGCTGCGTGCGGAAAGCACAATC‑ACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCt                                                                                                                                                                                                                                           <  1:55665/251‑1 (MQ=255)   gtggtgGCAGAGTGGCGGCGTAATGCGACGCGCAGACGGGCCTGCAATTCGCCAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCT‑CTTCGCTGCGTGCGGAAAGCACAATC‑ACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTc                                                                                                                                                                                                                                         <  1:315274/251‑1 (MQ=255)           agagTGGCGGCGTAATGCGACGCGCAGACGGGCCTGCAATTCGCCAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCT‑CTTCGCTGCGTGCGGAAAGCACAATC‑ACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGc                                                                                                                                                                                                                                 <  1:77742/251‑1 (MQ=255)                                            cTGCAATTCGCCAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCT‑CTTCGCTGCGTGCGGAAAGCACAATC‑ACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAG‑CGTTTCGGCCTCAAAGACgcgc                                                                                                                                                                                                >  1:21757/1‑250 (MQ=255)                                                 aTTCGCCAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCT‑CTTCGCTGCGTGCGGAAAGCACAATC‑CCCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAG‑CGTTTCGGCCTCAAAGACGCGCATCCCg                                                                                                                                                                                          <  1:274510/251‑1 (MQ=255)                                                         aTGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCT‑CTTCGCTGCGTGCGGAAAGCACAATC‑ACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAG‑CGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTc                                                                                                                                                                                  >  1:132257/1‑251 (MQ=255)                                                            ccAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCT‑CTTCGCTGCGTGCGGAAAGCACAATC‑ACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAG‑CGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTcc                                                                                                                                                                                 >  1:83110/1‑249 (MQ=255)                                                                 cGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCT‑CTTCGCTGCGTGCGGAAAGCACAATC‑ACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGtt                                                                                                                                                                                                                           >  1:29179/1‑203 (MQ=255)                                                                                       gCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCT‑CTTCGCTGCGTGCGGAAAGCACAATC‑ACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAG‑CGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGc                                                                                                                                                    >  1:200300/1‑251 (MQ=255)                                                                                          ccGGCATCCAGCGCGGCGATTTTGTCGCTCT‑CTTCGCTGCGTGCGGAAAGCACAATC‑ACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAG‑CGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTg                                                                                                                                                 <  1:220944/251‑1 (MQ=255)                                                                                               aTCCAGCGCGGCGATTTTGTCGCTCT‑CTTCGCTGCGTGCGGAAAGCACAATC‑ACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAG‑CGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCAt                                                                                                                                            >  1:128650/1‑251 (MQ=255)                                                                                                   aGCGCGGCGATTTTGTCGCTCT‑CTTCGCTGCGTGCGGAAAGCACAATC‑ACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAG‑CGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTc                                                                                                                                        <  1:160391/251‑1 (MQ=255)                                                                                                   aGCGCGGCGATTTTGTCGCTCT‑CTTCGATGCGTGCGGAAAGCTCAATC‑ACCGGCACCGCGCTCCACTGGCGCATGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAG‑CGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATATCCTGTTCATCTTc                                                                                                                                        <  1:288341/251‑1 (MQ=255)                                                                                                    cgcctgccaTTTTGTCG‑TCTTCTTGGCTGCTGGCGGAAAGCACA‑TCTACCGTAACCGAGTTCCCCTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCACCGGGCAGCCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAG‑CGTGACGACCTCAAAGACGCGCAGCCCTTCGCCCTCCAGCGCCGTGCGCAGAAAGCAACGAATCGCCTGTTCATCTTca                                                                                                                                       <  2:300414/243‑1 (MQ=255)                                                                                                           gATTTTGTCGCTCT‑CTTCGCTGCGTGCGGAAAGCACAATC‑ACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAG‑CGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCa                                                                                                                                <  1:269908/251‑1 (MQ=255)                                                                                                                 gTCGCTCT‑CTTCGCTGCGTGCGGAAAGCACAATC‑ACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAG‑CGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACg                                                                                                                           <  1:352199/250‑1 (MQ=255)                                                                                                                               cTGCGTGCGGAAAGCACAATC‑ACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAG‑CGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATcc                                                                                                             >  1:333912/1‑251 (MQ=255)                                                                                                                                 gcgtgcgGAAAGCACAATC‑ACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAG‑CGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCt                                                                                                            <  1:279500/250‑1 (MQ=255)                                                                                                                                                          cACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAG‑CGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGc                                                                                   <  2:333912/251‑1 (MQ=255)                                                                                                                                                                cctcctcTGGCGCATGGCGCGGATAAACTCAATCCCCTCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCCGGCTGACAGGTTGCCGCTTCCAGCAAGCCTCGTTGCAG‑CGTTTCGGCCTCAAAGCCGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCATGGGCAGGTTg                                                                             <  2:148521/250‑1 (MQ=255)                                                                                                                                                                   ccACTGGCGCAGGTCTCGAATAAACTCTATCCCATCACCATCGGGCAGGCCGAGATCTAGAATAATCAAATCTGGCTTACAGGTTGCCTCTTCCAGCAAGCCTCGTTGTAGTCGTT‑CGGCCTCAAGTACGCGCATCCCGTCGCCCTCCATCTCCGTGCGCAGAAATCGACGAAGAGCTTGTTCATCTTCAACAATCATAACGGTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGGAGTTTGCTg                                                                          <  2:363966/251‑1 (MQ=255)                                                                                                                                                                    cACTGGCGCAGGTGGCGGATAAACTCAACCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTGCCAGCAAGCCGCGTTGCAG‑CGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAAGCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGTTTGCTgg                                                                         <  2:21757/251‑1 (MQ=255)                                                                                                                                                                        ggCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAG‑CGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGTTTGCTGGGGa                                                                      <  1:254074/250‑1 (MQ=255)                                                                                                                                                                                          aCTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAG‑CGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTGCAGGGGCAGTTTGCTGGGGAAGTGTAACACGAAAACAgg                                                   <  1:283160/251‑1 (MQ=255)                                                                                                                                                                                                           tCGGGCAGGCCGAGATCGCGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAG‑CGTTTCGGCCTCAACGACGAGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGACCGTTTGTCACATATCCTAATGAAATTCTTCACGTTCAGGGGCAGTTTGCTGGGGAAGGGTAACACGACAACAGGCACACCCTTCaggtcgg                                  >  2:155594/1‑251 (MQ=255)                                                                                                                                                                                                            cGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAG‑CGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGTTTGCTGGGGAAGTGTAACACGAAAACAGGCACCACCTTCcggtcggt                                 >  1:138686/1‑251 (MQ=255)                                                                                                                                                                                                                       agaTCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAG‑CGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGTTTGCTGGGGAAGTGTAACACGAAAACAGGCACCACCTTCCGGTCGGTTGAACGCGGTa                      <  1:67879/251‑1 (MQ=255)                                                                                                                                                                                                                        gaTCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCAAGCAAGCCGCGTTGCAG‑CGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCCGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGTTTGCTGGGGAAGTGTAACACGAAAACCGGCACCACCTTCCGGTCGGTTGAACGCGGTa                      >  1:167574/1‑250 (MQ=255)                                                                                                                                                                                                                         aTCGAGAAGAACCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAG‑CGTTGCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCATCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGTTTGCTGGGGAAGTGTAACACGAAAACAGGCACCACCTTCCGGTCGGTTGAACGCGGTAAt                    <  2:132257/251‑1 (MQ=255)                                                                                                                                                                                                                                   atCAAATATAGCTGACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAG‑CGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCACGTTCAGGGGCAGTTTGCTGGGGAAGTGTAACACGAAAACAGTCACCACCTTCCGGTCGGTTGAACGCGGTAATAGTGCCCCCg          >  1:313430/1‑251 (MQ=255)                                                                                                                                                                                                                                    tCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAG‑CGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGTTTGCTGGGGAAGTGTAACACGAAAACAGGCACCACCTTCCGGTCGGTTGAACGCGGTAATAGTGCCCCCgt         >  1:295874/1‑251 (MQ=255)                                                                                                                                                                                                                                            ggCTTACAGGTTGCCGCTTCCAGAAAGCCGCGTTGCAG‑CGTTTCGGCCTCAACGACGCGAATCCCGTCGCCCTCCAGCGCCGTGCGCAGACAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGACGTTTGCTGGGGAAGTGTAACACGAAAACATGCACCACCTTCCGGTCGGTTGACCGCGGTAATAGAGCCCCCGTGTACATc  >  2:83582/1‑250 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    CGGTGGTGGCAGAGTGGCGGCGTAATGCGACGCGCAGACGGGCCTGCAATTCGCCAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCT‑CTTCGCTGCGTGCGGAAAGCACAATC‑ACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACGGGTTGCCGCTTCCAGCAAGCCGCGTTGCAG‑CGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGTTTGCTGGGGAAGTGTAACACGAAAACAGGCACCACCTTCCGGTCGGTTGAACGCGGTAATAGTGCCCCCGTGTACATC  >  NC_000913/721364‑721846

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