Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,316,059	C→A	Q858K (CAG→AAG)	rcsD →	RcsD phosphotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,316,059	0	C	A	100.0%	66.5 / NA	19	Q858K (CAG→AAG)	rcsD	RcsD phosphotransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (8/11);  total (8/11)

ACAATTAATTGAAGTGCAACTGGCGCAGGAAGAGGTGACAGAATCGCCTCTGGGCGGAGATGAAAATGCGCAACTCCATGCCAGCGGCTATTATGCGCTCTTTGTAGACACAGTACCGGATGATGTTAAGAGGCTGTATACTGAAGCAGCAACCAGTGACTTTGCTGCGTTAGCCCAAACGGCTCATCGTCTTAAAGGCGTATTTGCCATGCTAAATCTGGTACCCGGCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAGCAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAATTTGAAGACTCTACAGCACTGAT  >  NC_000913/2315827‑2316301                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                   aCAATTAATTGAAGTGCAACTGGCGCAGGAAGAGGTGACAGAATCGCCTCTGGGCGGAGATGAAAATGCGCAACTCCATGCCAGCGGCTATTATGCGCTCTTTGTAGACACAGTACCGGATGATGTTAAGAGGCTGTATACTGAAGCAGCAACCAGTGACTTTGCTGCGTTAGCCCAAACGGCTCATCGTCTTAAAGGCGTATTTGCCATGCTAAATCTGGTACCCGGCAAGAAGTTATGTGAAACGCTgg                                                                                                                                                                                                                                  <  2:25301/251‑1 (MQ=255)          tGAAGTGCAACTGGCGCAGGAAGAGGTGACAGAATCGCCTCTGGGCGGAGATGAAAATGCGCAACTCCATGCCAGCGGCTATTATGCGCTCTTTGTAGACACAGTACCGGATGATGTTAAGAGGCTGTATACTGAAGCAGCAACCAGTGACTTTGCTGCGTTAGCCCAAACGGCTCATCGTCTTAAAGGCGTATTTGCCATGCTAAATCTGGTACCCGGCAAGAAGTTATGTGAAACGCTGGAACATCTg                                                                                                                                                                                                                          <  1:41216/250‑1 (MQ=255)            aaGTGCAACTGGCGCAGGAAGAGGTTACATAATCGCCTCTGGGCGGAGATGAAAATGCGCAACTCCATGCCAGCGGCTATTATGCGCTCTTTTTGGACACAGTACCGGATGATGTTAAGAGGCTGTATACTGAAGCAGCAACCAGTGACTGTGCTGCGTTAGCCCAAACGGCTCATCGTCTTAAATGCGTATGTGCCAGGCTAAATCTGTTACCCGGCAATAAGTGATGGGAAACGCTCGATCATCTGAtt                                                                                                                                                                                                                       <  2:109665/251‑1 (MQ=255)                    cTGGCGCAGGAAGAGGTGACAGAATCGCCTCTGGGCGGAGATGAAAATGCGCAACTCCATGCCAGCGGCTATTATGCGCTCTTTGTAGACACAGTACCGGATGATGTTAAGAGGCTGTATACTGAAGCAGCAACCAGTGACTTTGCTGCGTTAGCCCAAACGGCTCATCGTCTTAAAGGCGTATTTGCCATGCTAAATCTGGTACCCGGCAAGAAGTTATGTGAAACGCTGGAACATCTGATTCGTgaga                                                                                                                                                                                                                >  2:34659/1‑250 (MQ=255)                                gagGTGACAGAATCGCCTCTGGGCGGAGATGAAAATGCGCAACTCCATGCCAGCGGCTATTATGCGCTCTTTGTAGACACAGTACCGGATGATGTTAAGAGGCTGTATACTGAAGCAGCAACCAGTGACTTTGCTGCGTTAGCCCAAACGGCTCATCGTCTTAAAGGCGTATTTGCCATGCTAAATCTGGTACCCGGCAAGAAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAgg                                                                                                                                                                                                   >  1:94942/1‑251 (MQ=255)                                        aGAATCGCCTCTGGGCGGAGATGCAAATGCGCAACTCCATGCCAGCGGCTATTATGCGCTCTTTGTAGACACAGTACCGGATGATGTTAAGAGGCTGTATACTGAAGAAGCAACCAGTGACTTTGCTGCGTTAGCCCAAACGGCTCATCGTCTTAAAGGCGTATTTGCCATGCTAAATCTGGTACCCGGCAAGAAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGaa                                                                                                                                                                                            >  2:138875/1‑250 (MQ=255)                                                         gagaTGAAAATGCGCAACGCCATGCCAGCGGCTATTATGCGCTCTTTGTAGACACAGTACCGGATGATGTTAAGAGGCTGTATACTGAAGCAGCAACCAGTGACTTTGCTGCGTTAGCCCAAACGGCTCATCGTCTTAAAGGCGTATTTGCCATGCTAAATCTGGTACCCGGCAAGAAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTGCCAGGAATAGAAAAAGACATCAGCGACAtt                                                                                                                                                                          <  2:184262/251‑1 (MQ=255)                                                                            ccATGCCAGCGGCTATTATGCGCTCTTTGTAGACACAGTACCGGATGATGTTAAGAGGCTGTATACTGAAGCAGCAACCAGTGACTTTGCTGCGTTAGCCCAAACGGCTCATCGTCTTAAAGGCGTATTTGCCATGCTAAATCTGGTACCCGGCAAGAAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGc                                                                                                                                                        <  2:157676/250‑1 (MQ=255)                                                                                          attatGCTCTCGTTGTAGACACAGTACCGGATGATGTTAAGAGGATGTATACTTAATCAGCAACCAGTGACTTTGCTGCGTTAGTCCAAACTGCTCATCGGCGTAAAGGCGTATTTTCCATTCTAAATCTGGTACCCGGCAAGAAGTTATGTGAATCGCTGGAACATCTTATGCTTGAGAATGATGTTCCAGGAATAGAAAAAGACATCAGTTACATTGACAGGTTTTTCAAGAGCTTGCGGTACCAAgtt                                                                                                                                         <  2:161355/251‑3 (MQ=255)                                                                                                     tttGTAGACACAGTACCGGATGATGTTAAGAGGCTGTATACTGAAGCAGCAACCAGTGACTTTGCTGCGTTAGCCCAAACGGCTCATCGTCTTAAAGGCGTATTTGCCATGCTAAATCTGGTACCCGGCAAGAAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAGCAAGGTAGCCTATTAc                                                                                                                               >  1:210108/1‑250 (MQ=255)                                                                                                              acaGTACCGGATGATGTTAAGAGGCTGTATACTGAAGCAGCAACCAGTGACTTTGCTGCGTTAGCCCAAACGGCTCATCGTCTTAAAGGCGTATTTGCCATGCTAAATCTGGTACCCGGCAAGAAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAGCAAGGTAGCCTATTACATGAACAat                                                                                                                      >  2:109493/1‑250 (MQ=255)                                                                                                                     cGGATGATGTTAAGAGGCTGTATACTGAAGCAGCAACCAGTGACTTTGCTGCGTTAGCCCAAACGGCTCATCGTCTTAAAGGCGTATTTGCCATGCTAAATCTGGTACCCGGCAAGAAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAGCAAGGTAGCCTATTACATGAACAATATGAACGt                                                                                                              <  2:115496/251‑1 (MQ=255)                                                                                                                                               gAAGCAGCAACCAGTGACTTTGCTGCGTTAGCCCAAACGGCTCATCGTCTTAAAGGCGTATTTGCCATGCTAAATCTGGTACCCGGCAAGAAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAGCAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAg                                                                                    <  1:164765/251‑1 (MQ=255)                                                                                                                                                  gcagcaACCAGTGACTTTGCTGCGTTAGCCCAAACGGCTCATCGTCTTAAAGGCGTATTTGCCATGCTAAATCTGGTACCCGGCAAGAAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAGCAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCt                                                                                 >  2:55524/1‑251 (MQ=255)                                                                                                                                                                                 aaaCGGCTCATCGTCTTAAAGGCGTATTTGCCATGCTAAATCTGGTACCCGGCAAGAAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAGCAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAAt                                                  >  2:57573/1‑251 (MQ=255)                                                                                                                                                                                                              gCCATGCTAAATCTGGTACCCGGCAAGAAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAGCAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAAttt                      <  1:2227/250‑1 (MQ=255)                                                                                                                                                                                                                 aTGCTAAATCTGGTACCCGGCAAGAAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAGCAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAATTTGAAg                  <  1:112268/251‑1 (MQ=255)                                                                                                                                                                                                                      aaaTCTGGTACCCGGCAAGAAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAGCAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCTAATTTGAAGAct               <  1:41152/249‑1 (MQ=255)                                                                                                                                                                                                                                 ccGGCAAGAAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAGCAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAATTTGAAGACTCTACAGCACTGAt  >  1:128589/1‑251 (MQ=255)                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                   ACAATTAATTGAAGTGCAACTGGCGCAGGAAGAGGTGACAGAATCGCCTCTGGGCGGAGATGAAAATGCGCAACTCCATGCCAGCGGCTATTATGCGCTCTTTGTAGACACAGTACCGGATGATGTTAAGAGGCTGTATACTGAAGCAGCAACCAGTGACTTTGCTGCGTTAGCCCAAACGGCTCATCGTCTTAAAGGCGTATTTGCCATGCTAAATCTGGTACCCGGCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAGCAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAATTTGAAGACTCTACAGCACTGAT  >  NC_000913/2315827‑2316301

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