Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,268,367	C→T	S68F (TCT→TTT)	kdsA →	3‑deoxy‑D‑manno‑octulosonate 8‑phosphate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,268,367	0	C	T	100.0%	109.8 / NA	34	S68F (TCT→TTT)	kdsA	3‑deoxy‑D‑manno‑octulosonate 8‑phosphate synthase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (17/17);  total (17/17)

AATCGACATTTTACTCAAGATTAAGGCGATCCTATGAAACAAAAAGTGGTTAGCATTGGCGACATCAACGTAGCAAATGACCTGCCGTTCGTACTGTTTGGCGGTATGAACGTGTTGGAATCTCGCGATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTCTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCCAGGAGTTGAAGCAGACTTTTGGCGTGAAAATTATCACCGACGTT‑CACGAACCAAGTCAGGCACAGCCCGTTGCTGATGTCGTGGATGTGATTCAGTTGCCGGCGTTTCTTGCTCGCCAGACTGACCTGGTTGAAGCCATGGCGAAAACCGGTGCGGTAATTAACGTCAAGAAACCACAGTTTGTCAGCCCGGGACA  >  NC_000913/1268132‑1268604                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                               aatCGACATTTTACTCAAGATTAAGGCGATCCTATGAAACAAAAAGTGGTTAGCATTGGCGACATCAACGTAGCAAATGACCTGCCGTTCGTACTGTTTGGCGGTATGAACGTGTTGGAATCTCGCGATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTTTTATCGTGGACCggg                                                                                                                                                                                                                                 <  1:332540/251‑1 (MQ=255)                            gATCCTATGACACAAAAAGTGGTTAGCATTGGCGACATCAACGTAGCAACTGACCTGCCGTTCGTACTGTTTGGCGGTATGAACGTGTTGGAATCTCGCGATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTTTTATCGTGGACCGGGCTTGGAAGAAGGGCTGAAAATCTTcc                                                                                                                                                                                                       >  2:237735/1‑250 (MQ=255)                                  aTGAAACAAAATGTGGTTAGCATTGTCGACATCAACGTAGCAAATGACCTGCCGTTCGTACTGTTTGGCGGTATGAACGTGTTGGAATCTCGCGATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTGCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTTTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCCAGGAGtt                                                                                                                                                                                                <  1:286714/251‑1 (MQ=255)                                          aaaaGTGGTTAGCATTGGCGACATCAACGTAGCAAATGACCTGCCGTTCGTACTGTTTGGCGGTATGAACGTGTTGGAATCTCGCGATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTTTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCCAGGAGTTGAAGc                                                                                                                                                                                           >  1:356772/1‑248 (MQ=255)                                            aaGTGGTTAGCATTGGCGACATCAACGTAGCAAATGACCTGCCGTTCGTACTGTTTGGCGGTATGAACGTGTTGGAATCTCGCGATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTTTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCCAGGAGTTGAAGCAGACt                                                                                                                                                                                      >  1:12809/1‑251 (MQ=255)                                                           gCGACATCAACGTAGCAAATGACCTGCCGTTCGTACTGTTTGGCGGTATGAACGTGTTGGAATCTCGCGATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTTTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCCAGGAGTTGAAGCAGACTTTTGGCGTGAAAAtt                                                                                                                                                                       <  2:2847/251‑1 (MQ=255)                                                             gACATCAACGTAGCAAATGACCTGCCGTTCGTACTGTTTGGCGGTATGAACGTGTTGGAATCTCGCGATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTTTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCCAGGAGTTGAAGCAGACTTTTGGCGTGAAAATTAt                                                                                                                                                                     >  2:176839/1‑251 (MQ=255)                                                                  cAACGTAGCAAATGACCTGCCGTTCGTACTGTTTGGCGGTATGAACGTGTTGGAATCTCGCGATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTTTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCCAGGAGTTGAAGCAGACTTTTGGAGTGAACATTATCACCg                                                                                                                                                                >  2:553210/1‑251 (MQ=255)                                                                        aGCAAATGACCTGCCGTTCGTACTGTTTGGCGGTATGAACGTGTTGGAATCTCGCGATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTTTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCCAGGAGTTGAAGCAGACTTTTGGCGTGAAAATTATCACCGACGTT‑c                                                                                                                                                         <  1:494384/251‑1 (MQ=255)                                                                           aaaTGACCTGCCGTTCGTACTGTTTGGCGGTATGAACGTGTTGGAATCTCGCGATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTTTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCCAGGAGTTGAAGCAGACTTTTGACGTGAAAATTATCACCGACGTT‑CACg                                                                                                                                                      >  1:10397/1‑251 (MQ=255)                                                                                      cGTTCGTACTGTTTGGCGGTATGAACGTGTTGGAATCTCGCGATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGATCCTCCATCCACTTTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCCAGGAGTTGAAGCAGACTTTTGGCGTGAAAATTATCACCGACGTT‑CACGAACCAAGTc                                                                                                                                             >  1:358467/1‑249 (MQ=255)                                                                                      cGTTCGTACTGCTTGGCGGTATGCACGTGTTGGACTCTCGCGATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTTTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCCAGGAGTTGAAGCAGACTTTTGGCGTGAAAATTATCACCGACGTT‑CACGAACCAAGTc                                                                                                                                             >  2:386632/1‑249 (MQ=255)                                                                                           gTACTGTTTGGCGGTATGAACGTGTTGGAATCTCGCGATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTTTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCCAGGAGTTGAAGCAGACTTTTGGCGTGAAAATTATCACCGACTTT‑CACGAACCAAGTCAGGcaca                                                                                                                                      <  2:532107/251‑1 (MQ=255)                                                                                                ttttGGCGGTATGAACGTGTTGGAATCTCGCGATCTGGCGATGCGCAATTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTCCGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTTTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAAGCTTCCAGGAGTTGAAGCAGACTTTTGGCGTGAAAATTATCACCGACGTT‑CACGAACCAAGTCAGGCACAGccc                                                                                                                                  <  2:85363/249‑1 (MQ=255)                                                                                                                 tgtTGGAATCTCGCGATCTGGCGATGCGCATTTGCGAGAACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTTTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCCAGGAGTTGAAGCAGACTTTTGGCGTGAAAATTATCACCGACGTT‑CACGAACAAAGTCAGGCACAGCCCGTTGCTGATGTCGTGGAt                                                                                                                >  2:396939/1‑251 (MQ=255)                                                                                                                  gtTGGAATCTCGCGATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTTTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCCAGGAGTTGAAGCAGACTTTTGGCGTGAAAATTATCACCGACGTT‑CACGAACCAAGTCAGGCACAGCCCGTTGCTGATGTCGTGGAt                                                                                                                <  1:126913/250‑1 (MQ=255)                                                                                                                            cgcgATCTTGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTTTTATCGTGGACCGGGCCTGGAAGAAGGGCTGACAATCTTCCAGGAGTTGAAGCAGACTTTTGGCGTGAAAATTATCACCGACGTT‑CACGAACCAAGTCAGGCACAGCCCGTTGCTGCTGTCGTGGATGTGATTCAGtt                                                                                                     >  1:334458/1‑251 (MQ=255)                                                                                                                            cgcgATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTTTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCCAGGAGTTGAAGCAGACTTTTGGCGTGAAAATTATCACCGACGTT‑CACGAACCAAGTCAGGCACAGCCCGTTGCTGATGTCGTGGATGTGATTCAGtt                                                                                                     <  2:35104/251‑1 (MQ=255)                                                                                                                                    ggCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTTTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCCAGGAGTTGAAGCAGACTTTTGGCGTGAAAATTATCACCGACGTT‑CACGAACCAAGTCAGGCACAGCCCGTTGCTGATGTCGTGGATGTGATTCAGTTGCCGGCGt                                                                                             <  2:69415/251‑1 (MQ=255)                                                                                                                                      cGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTCTCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCCACCGCTCCTCCATCCACTTTTATCGTGGCCCGGGCCTGGAAGAAGGGATGAAAATCTTCCAGGAGTTGAAGCAGACTTTTGGCGTGAAAATTATCACCGACGTT‑CACGAACCAAGTCAGGCACAGCCCGTTGCTGATGTCGTGGATGTGATTCAGTTGCCGGCGttt                                                                                           >  1:366373/1‑251 (MQ=255)                                                                                                                                         tGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGCATCCCTTAAGTGTTCAAAGACTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTTTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCCAGGAGTTGAAGCAGACTTTTGGCGTGAAAATTATCACCGACGTT‑CACGAACCAAGTCAGGCACAGCCCGTTGCTGCTGTCGTGGATGTGATTCAGTTGCCGGCGTTTCtt                                                                                        >  2:302392/1‑251 (MQ=255)                                                                                                                                                        aCTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCCACCACTCCTCCATCCACTTTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCCAGGAGTTGAAGCAGACTTTTGGCGTGAAACTTATCACCGACGTT‑AACGAACCAAGTCAGGCACAGCCCGTTGCTGATGTCGTGGATGTGATTCAGTTGCCGACGTTTCTTGCTCGCCAGACTgtc                                                                         >  1:535857/1‑249 (MQ=255)                                                                                                                                                                   gtgacccaccatccgccgaTCCCTGACGTGTCCATAGCCTTGTTTGTCAAAGCCTACCGTTCCCCCATCCACTTTTATCGTCGACCGGGTCTGGAAGAATGGATGAAAATCTTCCAGGAGTTGAAGCAGATTTTTGGCTGGATAATTATCCCTGTCGTTTCA‑GAATCACGTCATGCTTCGCCCGTTTCTGATGTCGTGGATGTGATGCAGGTGCGGGTGTTTCTTCCTCGCCCGTCTCTCCTGGTTGaa                                                                <  2:359214/231‑1 (MQ=255)                                                                                                                                                                     gACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTTTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCCAGGAGTTGAAGCAGACTTTTGGCGTGAAAATTATCACCGACGTT‑CACGAACCAAGTCAGGCACAGCCCGTTGCTGATGTCGTGGATGTGATTCAGTTGCCGGCGTTTCTTGCTCGCCAGACTGACCTGGTTGAAGcc                                                             >  1:492580/1‑250 (MQ=255)                                                                                                                                                                                 gggTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTTTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCCAGGCGTTGAAGCAGACTTTTGGCGTGAAAATTATCACCGACGTT‑CACGAACCAAGTCAGGCACAGCCCGTTGCTGATGTCGTGGATGTGATTCAGTTGCCGGCGTTTCTTGCTCGCCAGCCTGACCTGGTTGAAGCCATGGCGAAAACCg                                                >  2:187491/1‑251 (MQ=255)                                                                                                                                                                                            aCGTGTTCAAAGCCTCGTTTGACAAAGCCAACCGCTCCTCCATCCACTTTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCCAGGAGTTGAAGCAGACTTTTGGCGTGAAAATTATCACCGACGTT‑CACGAACCAAGTCAGGCACAGCCCGTTGCTGATGTCGTGGATGTGATTCAGTTGCCGGCGTTTCTTGCTCGCCAGACTGACCTGGTTGAAGCCATGGCGAAAACCGGTGCGGTAAt                                      <  1:547158/250‑1 (MQ=255)                                                                                                                                                                                                gtTCAAAGCCTTTTTTGACAAAGCCAACCGCTCCTCCATCCACTTTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCCAGGAGTTGAAGCAGACTTTTGGCGTGAAAATTATCACCGACGTT‑CACGAACCAAGTCAGGCACAGCCCGTTGCTGATGTCGTGGATGTGATTCAGTTGCCGGCGTTTCTTGCTCGCCAGACTGACCTGGTTGAAGCCATGGCGAAAACCGGTGCGGTAATTAAc                                  <  2:489621/250‑1 (MQ=255)                                                                                                                                                                                                     aaGCCTCTTGTGACAAAGCCAACCGCTCCTCCATCCACTTTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCCAGGAGTTGAAGCAGACTTTTGGCGTGAAAATTATCTCCGACGTT‑CACGAACCAAGTCAGGCACAGCCCGTTGCTGATGTCGTGGATGTGATTCAGTTGCCGGCGTTTCTTGCTCGCCAGACTGACCTGGTTGAAGCCATGGCGAAAACCGGTGCGGTAATTAACGTCAAg                            <  1:237735/251‑1 (MQ=255)                                                                                                                                                                                                           ctTTTGACAAAGCCAACCGCTCCTCCATCCACTTTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCCAGGAGTTGAAGCAGACTTTTGGCGTGAAAATTATCACCGACGTT‑CACGAACCAAGTCAGGCACAGCCCGTTGCTGATGTCGTGGATGTGATGCAGTTGCCGGCGTTTCTGGCTCGCCAGACTGACCTGGTTGAGGCCATGGCGAAAACCGGTGTGTTAAGTAACGGGAAGAAAc                        <  2:358467/249‑1 (MQ=255)                                                                                                                                                                                                               tGACAAAGCCAACCGCTCCTCCATCCACTTTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCCAGGAGTTGAAGCAGACTTTTGGCGTGAAAATTATCACCGACGTT‑CACGAACCAAGTCAGGCACAGCCCGTTGCTGATGTCGTGGATGTGATTCAGTTGCCGGCGTTTCTTGCTCGCCAGACTGACCTGGTTGAAGCCATGGCGAAAACCGGTGCGGTAATTAACGTCAAGAAACCACAGt                  >  2:179166/1‑251 (MQ=255)                                                                                                                                                                                                                gACAAAGCCAACCGCTCCTCCATCCACTTTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCCAGGAGTTGAAGCAGACGTTTGGCGTGAAAATTATCACCGACGTT‑CACGAACCAAGTCAGGCACAGCCCGTTGCTGATGTCGTGGATGTGATTCAGTTGCCGGCGTTTCGTGCTCGCCAGACTGACCTGGTTGAAGCCATGGCGAAAACCGGTGCGGTAATTAACGTCAAGAAACCACAGt                  <  1:386632/250‑1 (MQ=255)                                                                                                                                                                                                                    aaGCCAACCGCTCCTCCATCCACTTTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCCAGGAGTTGAAGCAGACTTTTGGCGTGAAAATTATCACCGACGTT‑CACGAACCAAGCCAGGCACAGCCCGTTGCTGATGTCGTGGATGTGATTCAGTTGCCGGCGTTTCTTGCTCGCCAGACTGACCTGGTTGAAGCCATGGCGAAAACCGGTGCGGTAATTAACGTCAAGAAACCACAGTTTGTc             >  2:204624/1‑251 (MQ=255)                                                                                                                                                                                                                          aCCGCTCCTCCATCCACTTTGATCGTGGACCGGGCCGGGAAGAAGGGATGAAAATCTTCCAGGAGTTGAAGCAGACTTTTGGCGTGAAAATTATCACCGACGTT‑CACGAACCAAGTCAGGCACAGCCCGTTGCTGATGTCGTGGATGTGATTCAGTTGCCGGCGTTTCTTGCTCGCCAGACTGACCTGGGTGAAGCCATGGCGAAAACCGGTGCGGTAATTAACGTCAAGAAACAACAGTTTGTCAGCCCg       <  2:432690/251‑1 (MQ=255)                                                                                                                                                                                                                                cctccATCCACTTTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCCAGGAGTTGAAGCAGACTTTTGGCGTGAAAATTATCACCGACGTT‑CACGAACCAAGTCAGGCACAGCCCGTTGCTGATGTCGTGGATGTGATTCAGTTGCCGGCGTTTCTTGCTCGCCAGACTGACCTGGTTGAAGCCATGGCGAAAACCGGTGCGGTAATTAACGTCAAGAAACCACAGTTTGTCAGCCCGGGACa  <  1:290631/250‑1 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                               AATCGACATTTTACTCAAGATTAAGGCGATCCTATGAAACAAAAAGTGGTTAGCATTGGCGACATCAACGTAGCAAATGACCTGCCGTTCGTACTGTTTGGCGGTATGAACGTGTTGGAATCTCGCGATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTCTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCCAGGAGTTGAAGCAGACTTTTGGCGTGAAAATTATCACCGACGTT‑CACGAACCAAGTCAGGCACAGCCCGTTGCTGATGTCGTGGATGTGATTCAGTTGCCGGCGTTTCTTGCTCGCCAGACTGACCTGGTTGAAGCCATGGCGAAAACCGGTGCGGTAATTAACGTCAAGAAACCACAGTTTGTCAGCCCGGGACA  >  NC_000913/1268132‑1268604

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