Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	721,886	(C)5→4	coding (2529/2685 nt)	kdpD ←	sensory histidine kinase KdpD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	721,882	0	C	.	100.0%	116.3 / NA	28	coding (2533/2685 nt)	kdpD	sensory histidine kinase KdpD
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (14/14);  total (14/14)

CTTCCGGTCGGTTGAACGCGGTAATAGTGCCCCCGTGTACATCCACTATCGCCCGACAAATTGCCAGTCCAAGCCCTACCCCCGGTACTGCCGACTCTTTATTCCCGCGAGCAAACTTATCAAATATCGTCTGCTCCTGGCCTGGCGGAAGA  >  NC_000913/721804‑721955                                                                               |                                                                          cTTCCGGTCGGTTGAACGCGGTAATAGTGCCCCCGTGTACATCCACTATCGCCCGACAAATTGCCAGTCCAAGCCCTAccc                                                                         >  1:154312/1‑81 (MQ=255)    ccGGTCGGTTGAACGCGGTAATAGTGCCCCCGTGTACATCCACTATCGCCCGACAAATTGCCAGTCCAAGCCCTAcccc                                                                        <  1:154282/79‑1 (MQ=255)        tcggtTGAACGCGGTAATAGTGCCCCCGTGTACATCCACTATCGCCCGACAAATTGCCAGTCCAAGCCCTA‑CCCCGGTACt                                                                 >  1:154315/1‑81 (MQ=255)        tcggtTGAACGCGGTAATAGTGCCCCCGTGTACATCCACTATCGCCCGACAAATTGCCAGTCCAAGCCCTA‑CCCCGGTACt                                                                 >  2:154326/1‑81 (MQ=255)           gtTGAACGCGGTAATAGTGCCCCCGTGTACATCCACTATCGCCCGACAAATTGCCAGTCCAAGCCCTA‑CCCCGGTACTGcc                                                              <  1:154283/81‑1 (MQ=255)              gAACGCGGTAATAGTGCCCCCGTGTACATCCACTATCGCCCGACAAATTGCCAGTCCAAGCCCTA‑CCCCGGTACTGCCGAc                                                           >  1:154319/1‑81 (MQ=255)                  gcgGTAATAGTGCCCCCGTGTACATCCACTATCGCCCGACAAATTGCCAGTCCAAGCCCTA‑CCCCGGTACTGCCGACTCtt                                                       >  1:154334/1‑81 (MQ=255)                       aaTAGTGCCCCCGTGTACATCCACTATCGCCCGACAAATTGCCAGTCCAAGCCCTA‑CCCCGGTACTGCCGACTCTTTATTc                                                  >  2:154320/1‑81 (MQ=255)                        aTAGTGCCCCCGTGTACATCCACTATCGCCCGACAAATTGCCAGTCCAAGCCCTA‑CCCCGGTACTGCCGACTCTTTATTc                                                  <  1:154284/80‑1 (MQ=255)                         tAGTGCCCCCGTGTACATCCACTATCGCCCGACAAATTGCCAGTCCAAGCCCTA‑CCCCGGTACTGCCGACTCTTTATTccc                                                <  1:154286/81‑1 (MQ=255)                         tAGTGCCCCCGTGTACATCCACTATCGCCCGACAAATTGCCAGTCCAAGCCCTA‑CCCCGGTACTGCCGACTCTTTATTccc                                                <  1:154285/81‑1 (MQ=255)                          aGTGCCCCCGTGTACATCCACTATCGCCCGACAAATTGCCAGTCCAAGCCCTA‑CCCCGGTACTGCCGACTCTTTATTCCcg                                               >  1:154321/1‑81 (MQ=255)                            tGCCCCCGTGTACATCCACTATCGCCCGACAAATTGCCAGTCCAAGCCCTA‑CCCCGGTACTGCCGACTCTTTATTCCcgcg                                             <  1:154287/81‑1 (MQ=255)                               ccccGTGTACATCCACTATCGCCCGACAAATTGCCAGTCCAAGCCCTA‑CCCCGGTACTGCCGACTCTTTATTCCCGCGAg                                           >  1:154324/1‑80 (MQ=255)                                cccGTGTACATCCACTATCGCCCGACAAATTGCCAGTCCAAGCCCTA‑CCCCGGTACTGCCGACTCTTTATTCCCGCGAGCa                                         <  1:154288/81‑1 (MQ=255)                                    tgtACATCCACTATGGCCCGACAAATTGCCAGTCCAAGCCCTA‑CCCCGGTACTGCCGACTCTTTATTCCCGCGAGCAAACt                                     <  1:154289/81‑1 (MQ=255)                                        cATCCACTATCGCCCGACAAATTGCCAGTCCAAGCCCTA‑CCCCGGTACTGCCGACTCTTTATTCCCGCGAGCAAACTTAt                                  <  1:154290/80‑1 (MQ=255)                                         aTCCACTATCGCCCGACAAATTGCCAGTCCAAGCCCTA‑CCCCGGTACTGCCGACTCTTTATTCCCGCGAGCAAACTTATc                                 >  1:154323/1‑80 (MQ=255)                                            cACTATCGCCCGACAAATTGCCAGTCCAAGCCCTA‑CCCCGGTACTGCCGACTCTTTATTCCCGCGAGCAAACTTATCAAat                             <  1:154291/81‑1 (MQ=255)                                             aCTATCGCCCGGCAAATTGCCAGTCCAAGCCCTA‑CCCCGGTACTGCCGACTCTTTATTCCCGCGAGCAAACTTATCAAata                            >  2:154360/1‑81 (MQ=255)                                               tATCGCCCGACAAATTGCCAGTCCAAGCCCTA‑CCCTGGTACTGCCGACTCTTTATTCCCGCGAGCAAACTTATCAAATATc                          >  2:154325/1‑81 (MQ=255)                                                  cGCCCGACAAATTGCCAGTCCAAGCCCTA‑CCCCGGTACTGCCGACTCTTTATTCCCGCGAGCAAACTTATCAAATATCGTc                       >  2:154333/1‑81 (MQ=255)                                                   gCCCGACAAATTGCCAGTCCAAGCCCTA‑CCCCGGTACTGCCGACTCTTTATTCCCGCGAGCAAACTTATCAAATATCGTc                       <  2:154292/80‑1 (MQ=255)                                                           aaTTGCCAGTCCAAGCCCTA‑CCCCGGTACTGCCGACTCTTTATTCCCGCGAGCAAACTTATCAAATATCGTCTGCTcctgg              <  2:154293/81‑1 (MQ=255)                                                                ccAGTCCAAGCCCTA‑CCCCGGTACTGCCGACTCTTTATTCCCGCGAGCAAACTTATCAAATATCGTCTGCTcctggcctg          <  2:154294/80‑1 (MQ=255)                                                                   gTCCAAGCCCTA‑CCCCGGTACTGCCGACTCTTTATTCCCGCGAGCAAACTTATCAAATATCGTCTGCTCCTGGCCTGGCgg      <  2:154295/81‑1 (MQ=255)                                                                    tCCAAGCCCTA‑CCCCGGTACTGCCGACTCTTTATTCCCGCGAGCAAACTTATCAAATATCGTCTGCTCCTGGCCTGGCGGa     <  1:154296/81‑1 (MQ=255)                                                                     ccAAGCCCTA‑CCCCGGTACTGCCGACTCTTTATTCCCGCGAGCAAACTTATCAAATATCGTCTGCTCCTGGCCTGGCgg      >  1:154347/1‑79 (MQ=255)                                                                      cAAGCCCTA‑CCCCGGTACTGCCGACTCTTTATTCCCGCGAGCAAACTTATCAAATATCGTCTGCTCCTGGCCTGGCgg      >  1:154328/1‑78 (MQ=255)                                                                       aaGCCCTA‑CCCCGGTACTGCCGACTCTTTATTCCCGCGAGCAAACTTATCAAATATCGTCTGCTCCTGGCCTGGCGGAAGa  >  1:154335/1‑81 (MQ=255)                                                                               |                                                                          CTTCCGGTCGGTTGAACGCGGTAATAGTGCCCCCGTGTACATCCACTATCGCCCGACAAATTGCCAGTCCAAGCCCTACCCCCGGTACTGCCGACTCTTTATTCCCGCGAGCAAACTTATCAAATATCGTCTGCTCCTGGCCTGGCGGAAGA  >  NC_000913/721804‑721955

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