Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,171,233	C→T	S250S (AGC→AGT)	qseC →	sensor histidine kinase QseC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,171,233	0	C	T	100.0%	41.0 / NA	15	S250S (AGC→AGT)	qseC	sensor histidine kinase QseC
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (9/6);  total (9/6)

ATCATCATGATGGTACTACTGGGTCGTGAACTCGCGCCGCTGAACAAACTGGCGCTGGCACTACGTATGCGTGACCCTGACTCGGAAAAACCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGGCCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGATGATGGTTCGTGAACGACGCTTTACCTCCGACGCAGCTCACGAACTTCGTAGCCCGTTAACGGCGCTGAAAGTGCAAACCGAAGTTGCGCAGCTCTCTGACGATGATCCGCAGGCGCGGAAAAAAGCACTGCTCCAATTACATTCCGGGATCGATCGCGCTACTCGTCTGGTTGATCAACTGCTCACGCTATCGCGGCTGGACTCACTGGATAACCTTCAGGACGTCGCGGAGATCCCGCTTGAAGATCTCCTGCAATCGTCGGTGATGGATATTTACCAC  >  NC_000913/3171009‑3171461                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                     atcatcATGATGGTACTACTGGGTCGTGAACTCGCGCCGCTGAACAAACTGGCGCTGGCACTACGTATGCGTGACCCTGACTCGGAAAAACCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGGCCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGATGATGGTTCGTGAACGACGCTTTACCTCCGACGCAGCTCACGAACTTCGTAGTCCGTTAACGGCGCTGAAAGTGCAAAc                                                                                                                                                                                                            <  1:115002/251‑1 (MQ=255)                                       gcTGAACAAACTGGCGCTGGCACTACGTATGCGTGACCCTGACTCGGAAAAACCACTAACCGCGACTGGCGTACCCAGCGAAGTGCGGCCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGATGATGGTTCGTGAACGACGCTTTACCTCCGACGCAGCTCACGAACTTCGTAGTCCGTTAACGGCGCTGAAAGTGCAAACCGAAGTTGCGCCGCTCTCTGACGATGATCCGCAGgcgc                                                                                                                                                                      >  1:147979/1‑251 (MQ=255)                                          gAACAACCTGGCGCTGGCACTACGTCTGCGTGCCCCATACTCGGAAAAACCACTAATCGCGACTGGCGTCCCCACCGAAGTGCGGCCACTGGTTTAGTCGCTAAATCACCTGTTCGCCCGCACACATGCTATGATGGTTCGTGAACGACGCTTTACCTCCGCCGCAGCTCACGAGCTTCGTAGTCCGTTAACCGCGCTGAAAGTGCACTCCGAATTTGCGCAGCTCTCTGACGATGATCCGCAGgcgcn                                                                                                                                                                     <  1:61466/249‑2 (MQ=255)                                                                                   cGGAAAAACCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGGCCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGATGATGGTTCGTGAACGACGCTTTACCTCCGACGCAGCTCACGAACTTCGTAGTCCGTTAACGGCGCTGAAAGTGCAAACCGAAGTTGCGCAGCTCTCTGACGATGATCCGCAGGCGCGGAAAAAAGCACTGCTCCAATTACATTCCGGGATCGATCgcgc                                                                                                                           <  1:182646/250‑1 (MQ=255)                                                                                          aCCACTAAACGAGACTGGCGTACCAAGAGAAGTGCGGCCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGATGATGGTTCGTGAACGACGCTTTACCTCCGACGCAGCTCACGAACTTCGTAGTCCGTTAACGGCGCTGAAAGTGCAAACCGAAGTTGCGCAGCTCTCTGAc                                                                                                                                                                                      >  1:165414/1‑184 (MQ=255)                                                                                                    gcgACTGGCGTACCCAGCGAAGTGCGGCCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGATGATGGTTCGTGAACGACGCTTTACCTCCGACGCAGCTCACGAACTTCGTAGTCCGTTAACGGCGCTGAAAGTGCAAACCGAAGTTGCGCAGCTCTCTGACGATGATCCGCAGGCGCGGAAAAAAGCACTGCTCCAATTACATTCCGGGATCGATCGCGCTACTCGTCTGGTTGATc                                                                                                          >  1:47623/1‑250 (MQ=255)                                                                                                             ntACCCAGCGAAGTGCGGCCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCCCACATGCGATGATGGTTCGTGAACGAAGCTTTACCTCCGACGCAGCTCACGAACTTCGTAGTCCGTTAACGGAGCTGAAAGTGCAAACCGAAGTTGCGCAGCTCTCTGACGATGATCCGCAGGCGCGGAAAAAAGCACTGCTCCAATTACATTCCGGGATCGCTCGCGCTACTCGTCTGGTTGATCAACTGCTCAc                                                                                                >  1:130021/2‑251 (MQ=255)                                                                                                               aCCCAGCGAAGTGCGGCCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGATGATGGTTCGTGAACGACGCTTTACCTCCGACGCAGCTCACGAACTTCGTAGTCCGTTaan                                                                                                                                                                                                                              <  1:54080/123‑2 (MQ=255)                                                                                                                              nccACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGATGATGGATCGTGAACGACGCTTTACCTCCGACGCAGCTCACGAACTTCGTAGTCCGTTAACGGCGCTGAAAGTGCAAACCGAAGTTGCGCAGCTCTCTGACGATGATCCGCAGGCGCGGAAAAAAGCACTGCTCCAATTACATTCCGGGATCGATCGCGCTACTCGTCTGGTTGATCAACTGCTCACGCTATCGCGGCTGGACt                                                                               >  1:128641/2‑251 (MQ=255)                                                                                                                                                              nccGCACACATGCGATGATGGTTCGTGAACGACGCTTTACCTCCGACGCAGCTCACGAACTTCGTAGTCCGTTAACGGCGCTGAAAGTGCAAACCGAAGTTGCGCAGCTCTCTGACGATGATCCGCAGGCGCGGAAAAAAGCACTGCTCCAATTACATTCCGGGATCGATCGCGCTACTCGTCTGGTTGATCAACTGCTCACGCTATCGCGGCTGGACTCACTGGATAACCTTCAGGACGTCGCgg                                                    >  1:42204/2‑246 (MQ=255)                                                                                                                                                                      caTGCGATGATGGTTCGTGAACGACGNTTTACCNNCGNCGCAGCTCACGAACTTCGTNGTNNNTNNNCNGCGCTNAANNTGCAAACCGNAGTTGCGCAGCTCTCTGACGATGATCCGCAGGCGCGGAAACAAGCACTGCTCCNATTACATTCNNGGATCGATCGCGCTACTCGTCNGGTNNANNANCTGNTc                                                                                                  >  2:46459/1‑192 (MQ=255)                                                                                                                                                                      caTGCGATGATGGGTCGTGAACGACGCTTTACCTCCGACGCAGCTCACGAACTTCGTAGTCCGTTAACGGCGCTGAAAGTGCAAACCGAAGTTGCGCAGCTCTCTGACGATGATCCGCAGGCGCGGAAAAAAGCACTGCTCCAATTACATTCCGGGAGCGATCGCGCTACTCGTCTGGTTGATCAACTGCTCACGCTAGCGCGGCTgg                                                                                  <  1:46459/208‑1 (MQ=255)                                                                                                                                                                                                     ccTCCGACGCAGCTCACGAACTTCGTAGTCCGTTAACGGCGCTGAAAGTGCAAACCGAAGTTGCGCAGCTCTCTGACGATGATCCGCAGGCGCGGAAAAAAGCACTGCTCCAATTACATTCCGGGATCGATCGCGCTACTCGTCTGGTTGATCAACTGCTCACGCTATcgcg                                                                                       <  1:58648/172‑1 (MQ=255)                                                                                                                                                                                                        ccGACGCAGCTCACGAACTTCGTAGTCCGTTAACGGCGCTGAAAGTGCAAACCGAAGTTGCGCAGCTCTCTGACGATGATCCGCAGGCGCGGAAAAAAGCACTGCTCCAATTACATTCCGGGATCGATCGCGCTACTCGTCTGGTTGATCAACTGCTCACGCTATCGCGGCTGGACTCACTGGATAACCTTCAGGACGTCGCGGAGATCCCGCTTGAAGATCTCCTGCAATCGTCGGTGATGGATATTTAc     >  1:48623/1‑251 (MQ=255)                                                                                                                                                                                                            cGCAGCTCACGAACTTCGTCGTCCGTTAACGGCGCTGAAAGTGCAAACCGAAGTTGCGCAGCTCTCTGACGATGATCCGCAGGCGCGGAAAAAAGCACTGCTCCAATTACATTCCGGGATCGATCGCGCTACTCGTCTGGTTGATCAACTGCTCACGCTATCGCGGCTGGACTCACTGGATAACCTTCAGGACGTCGCGGAGATCCCGCTTGAAGATCTCCTGCAATCGTCGGTGATGGATATTTACcac  >  1:167449/1‑250 (MQ=255)                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                     ATCATCATGATGGTACTACTGGGTCGTGAACTCGCGCCGCTGAACAAACTGGCGCTGGCACTACGTATGCGTGACCCTGACTCGGAAAAACCACTAAACGCGACTGGCGTACCCAGCGAAGTGCGGCCACTGGTTGAGTCGCTAAATCAACTGTTCGCCCGCACACATGCGATGATGGTTCGTGAACGACGCTTTACCTCCGACGCAGCTCACGAACTTCGTAGCCCGTTAACGGCGCTGAAAGTGCAAACCGAAGTTGCGCAGCTCTCTGACGATGATCCGCAGGCGCGGAAAAAAGCACTGCTCCAATTACATTCCGGGATCGATCGCGCTACTCGTCTGGTTGATCAACTGCTCACGCTATCGCGGCTGGACTCACTGGATAACCTTCAGGACGTCGCGGAGATCCCGCTTGAAGATCTCCTGCAATCGTCGGTGATGGATATTTACCAC  >  NC_000913/3171009‑3171461

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