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	97.0%	104.3 / NA	33	S68F (TCT→TTT)	kdsA	3‑deoxy‑D‑manno‑octulosonate 8‑phosphate synthase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (20/12);  minor base G (0/1);  total (20/13)
Fisher's exact test for biased strand distribution p-value = 3.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

AATCGACATTTTACTCAAGATTAAGGCGATCCTATGAAACAAAAAGTGGTTAGCATTGGCGACATCAACGTAGCAAATGACCTGCCGTTCGTACTGTTTGGCGGTATGAACGTGTTGGAATCTCGCGATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTCTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCCAGGAGTTGAAGCAGACTTTTGGCGTGAAAATTATCACCGACGTTCACGAACCAAGTCAGGCACAGCCCGTTGCTGATGTCGTGGATGTGATTCAGTTGCCGGCGTTTCTTGCTCGCCAGACTGACCTGGTTGAAGCCATGGCGAAAACCGGTGCGGTAATTAACGTCAAGAAACCACAGTTTGTCAGCCC  >  NC_000913/1268132‑1268598                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                        aatCGACATTTTACTCAAGATTAAGGCGATCCTATGAAACAAAAAGTGGTTAGCATTGGCGACATCAACGTAGCAAATGACCTGCCGTTCGTACTGTTTGGCGGTATGAACGTGTTGGAATCTCGCGATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTTTTATCGTGGACCgg                                                                                                                                                                                                                           <  2:409698/250‑1 (MQ=255)      aCATTTTACTCAAGATTAAGGCGATCCTATGAAACAAAAAGTGGTTAGCATTGGCGACATCAACGTAGCAAATGACCTGCCGTTCGTACTGTTTGGCGGTATGAACGTGTTGGAATCTCGCGATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTTTTATCGTGGACCGGGCCTgg                                                                                                                                                                                                                     <  2:52790/251‑1 (MQ=255)                cAAGATTAAGGCGATCCTATGAAACAAAAAGTGGTTAGCATTGGCGACATCAACGTAGCAAATGACCTGCCGTTCGTACTGTTTGGCGGTATGAACGTGTTGGAATCTCGCGATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTTTTATCGTGGACCGGGCCTGGAAGAAGGGAt                                                                                                                                                                                                           >  1:431417/1‑251 (MQ=255)                                      aaCAAAAAGTGGTTAGCATTGGCGACATCAACGTAGCAAATGACCTGCCGTTCGTACTGTTTGGCGGTATGAACGTGTTGGAATCTCGCGATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTTTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCCAGGAGTTGAAg                                                                                                                                                                                     <  2:221153/251‑1 (MQ=255)                                        cAAAAAGGGGTTAGCATTGGCGACATCACCGTAGCAAATGACCTGCCGTTCGTACTGTTTGGCGGTATGAACGTGTTGAAATCTCGCGCTCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTTTTATCGTGGACCGGGCCTGGAAGAAGGGATGACAATCTTCCAGGAGTTGAAGcc                                                                                                                                                                                   >  2:477248/1‑250 (MQ=255)                                         aaaaaGTGGTTAGCATTGGCGACATCAACGTAGCAAATGACCTGCCGTTCGTACTGTTTGGCGGTATGAACGTGTTGGAATCTCGCGATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTTTTATCGTGGACCGGGCCTGGAAGAAGGGCTGAAAATCTTCCAGGAGTTGAAGc                                                                                                                                                                                    >  2:283304/1‑249 (MQ=255)                                           aaaGTGGTTAGCATTGGCGACATCAACGTAGCAAATGACCTGCCGTTCGTACTGTTTGGCGGTATGACCGTGTTGGAATCTCGCGATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTCTCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTTTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCAAGGAGTTGAAGCAGAc                                                                                                                                                                                >  2:513853/1‑251 (MQ=255)                                           aaaGTGGTTAGCATTGGCGACATCAACGTAGCAAATGACCTGCCGTTCGTACTGTTTGGCGGTATGAACGTGTTGGAATCTCGCGATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTTTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCCAGGAGTTGAAGCAGAc                                                                                                                                                                                >  1:564819/1‑251 (MQ=255)                                                ggTTAGCATTGGAGACATCAACGTAGAAAATTACCTGCCGTTCGTACTTTTTGGCGGTATGAACGTGTTGGAATCTCGCGATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTGTTATCGTGGACCGGGCCTTGAAGAAGGGATGAAAATCTTCCAGGAGTTGAAGCAGACTTTTg                                                                                                                                                                           <  2:228364/251‑1 (MQ=255)                                                       aTTGGCGACATCAACGTAGCAAATGACCTGCCGTTCGTACTGTTTGGCGGTATCAACGTGTTGGAATCTAGCGATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAACGCCAACCGCTCCTCCATCCACTTTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCCAGGAGTTGAAGCAGACTTTTGGCGTGa                                                                                                                                                                     >  1:85875/1‑250 (MQ=255)                                                                 tCAACGTAGCAAATGACCTGCCGTTCGTACTGTTTGGCGGTATGAACGTGTTGGAATCTCGCGATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTTTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCCAGGAGTTGAAGCAGACTTTTGGCGTGAAAATTATCAcc                                                                                                                                                          >  1:282372/1‑251 (MQ=255)                                                                   aaCGTAGCAAATGACCTGCCGTTCGTACTGTTTGGCGGTATGAACGTGTTGGAATCTCGCGATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTTTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCCAGGAGTTGAAGCAGACTTTTGGCGTGAAAATTATCACCg                                                                                                                                                         >  1:305743/1‑250 (MQ=255)                                                                    aCGTAGCAAATGACCTGCCGTTCGTACTGGTTGGCGGTATGAACGTGTTGGAATCTCGCGATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTTTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCCAGGAGTTGAAGCAGACTTTTGGCGTGAAAATTATCACCGAc                                                                                                                                                       <  2:305743/251‑1 (MQ=255)                                                                          cAAATGTCCTGCCGTTCGTACTGTTTGGCGGTATGAACTTGTTGGAAGCTCGCGATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTTTTAACGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCCAGTAGATGAAGCAGACTTTTGGCGTGAAAATTATCACCGACGTTccc                                                                                                                                                 <  2:442591/251‑3 (MQ=255)                                                                                             aCTGTTTGGCGGTATGAACGTGTTGGAATCTCGCGATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGGATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTTTTATCGTGGACCGGGCCTGGAAGACGGGATGAAAATCTTCCAGGAGTTGAAGCAGACTTTTGGCGTGAAAATTATCACCGCCGTTCACGAACCAAGCCAGGCACAGc                                                                                                                              >  2:330081/1‑251 (MQ=255)                                                                                               tGTTTGGCGGTATGAACGTGTTGGAATCTCGCGATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTTTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCCAGGAGTTGAAGCAGACTTTTGGCGTGAAAATTATCACCGACGTTCACGAACCAAGTCAGGCACAGccc                                                                                                                            >  1:558267/1‑251 (MQ=255)                                                                                                          aTGAACGTGTTGGAATCTAGCGATCTGGAGCTGCCCATTTGCGAACACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTTTTATCGTGGACCGGACCTGGAAGAAGGGATGAACATCTTCCAGGAGTTGAAGCAGACTTTTGGCGTGAAAATTATCACCGCCGTTCACGAACCAAGTCAGGCACAGCCCGTTGCTGATGt                                                                                                                 >  2:184945/1‑251 (MQ=255)                                                                                                                            cgcgATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTTTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCCAGGAGTTGAAGCAGACTTGTGGCGTGAAAATTATCACCGACGTTCACGAACCAAGTCAGGCACAGCCCGTTGCTGATGTCGTGGATGTGATTCAg                                                                                                 <  1:283304/249‑1 (MQ=255)                                                                                                                              cgATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCACCAACTGGGTATCCCTTACGTGTTAAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCAACTTTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCCAGGAGTTGACGCAGACTTTTGGCGTGAAAATTATCACCGACGTTCACGAACCAAGTCAGGCACAGCCCGTTGCTGATGCCGTGGATGAGATTCAGTTGc                                                                                             >  1:253649/1‑251 (MQ=255)                                                                                                                                 tCTGGCGATGCGCATTTGCGCGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTTTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCCAGGAGTTGAAGCAGACTTTTGGCGTGAAAATTATCACCGACGTTCACGAACCAAGTCAGGCACAGCCCGTTGCTGATGTCGTGGATGTGCTTCAGTTGCCgg                                                                                          >  2:407026/1‑251 (MQ=255)                                                                                                                                  cTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTTTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCCAGGAGTTGAAGCAGACTTTTGGCGTGAAAATTATCACCGACGTTCACGAACCAAGTCAGGCACAGCCCGTTGCTGATGTCGTGGATGTGATTCAGTTGCCGGc                                                                                         >  1:578369/1‑251 (MQ=255)                                                                                                                                    ggCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTTTTATCGTGGACCGGGCCTTGAAGAAGGGATGAAAATCTTCCAGGAGTTGAAGCAGACTTTTGGCGTGAAAATTATCACCGACGTTCACGAACCAAGTCAGGCACAGCCCGTTGCTGATGTCGTGGATGTGATTCAGTTGCCGGc                                                                                         <  1:300997/249‑1 (MQ=255)                                                                                                                                                            cGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTTTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCCa                                                                                                                                                                                               >  1:267147/1‑123 (MQ=255)                                                                                                                                                            cGTAACTGTGACCCAGAAACTGGGTATCCCGTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTTTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCCa                                                                                                                                                                                               <  2:267147/123‑1 (MQ=255)                                                                                                                                                               aaCTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTTTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCCAGGAGTTGAAGCAGACTTTTGGCGTGAAAATTATCACCGACGTTCACGAACCAAGTCAGGCACAGCCCGTTGCTGATGTCGTGGATGTGATTCAGTTGCCGGCGTTTCTTGCTCGCCAGACTGACCTGGtt                                                             <  1:513853/250‑1 (MQ=255)                                                                                                                                                                       cccAGCAACTGGGTATCCCTTACGTGTTCAACGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTTTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCCAGGCGTTGAAGCAGACTTTTGGCGTGAAAATTATCACCGACGTTCACGAACCAAGTCAGGCACAGCCCGTTGCTGATGTCGTGGATGTGATTCAGTTGCCGGCGTTTCTTGCTCGCCAGACTGACCTGGTTGAAGCCATg                                                    >  2:351619/1‑251 (MQ=255)                                                                                                                                                                                     aTCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTTTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCAAGGAGTTGAAGCAGACTTTTGGCGTAAAAATTATCACCGACGTTCACGAACCAAGACAGAAACAGCCCGTTACTGATGTCGTGGATGTGATTCAGTTGCCGGCGTTTCTTGCTCGCCAGACTGCCCTGTTTGAAGCCATGGCGAAAACCAGTGc                                      >  1:548859/1‑251 (MQ=255)                                                                                                                                                                                            aCGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTTTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCCAGGAGTTGAAGCAGACTTTTGGCGTGAAAATTATCACCGACGTTCACGAACCAAGTCAGGCACAGCCCGTTGCTGATGTCGTGGATGTGATTCAGTTGCCGGCGTTTCTTGCTCGCCAGACTGACCTGGTTGAAGCCATGGCGAAAACCGGTGCGGTAAt                                <  2:462225/250‑1 (MQ=255)                                                                                                                                                                                              gtgtTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTTTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCCAGGAGTTGAAGCAGACTTTTGGCGTGAAAATTATCACCGACGTTCACGAACCAAGTCAGGCACAGCCCGTTGCTGATGTCGTGGATGTGATTCAGTTGCCGGCGTTTCTTGCTCGCCAGACTGCCCTGGTTGAAGCCATGGCGAAAACCGGTGCGGTAATTa                              >  1:391668/1‑250 (MQ=255)                                                                                                                                                                                                         ctctTTTGACAAAGCCAACCGCTCCTCCATCCACTTTTATCGTGGTCCGGGCCTGGAAGAAGGGATGAAAATCTTCCAGGAGTTGAACCAGACTTTTGGCGTGAAAATTATCACCGACGTTCACGAACCAAGTCAGGCCAAGCCCGTTGCTGATTTCGTGGATGTGATTCAGTTGCCGGCGGTTCTTGCTCGCCAGACTGACCTGGTTGAAGCCATGGCGAAAACTGGTGCGGTGATTAACGTCAAGAAAc                  >  2:231951/1‑251 (MQ=255)                                                                                                                                                                                                             tttGACAAAGCCAACCGCTCCTCCATCCACTTTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCCAGGAGTTGAAGCAGACTTTTGGCGTGAAAATTATCACCGACGTTCACGAACCAAGTCAGGCACAGCCCGTTGCTGATGTCGTGGATGTGATTCAGTTGCCGGCGTTTCTTGCTCGCCAGACTGACCTGGTGGAAGCCATGGCGAAAACCGGTGCGGTAATTAACGTCAAGAAACca                <  2:282372/249‑1 (MQ=255)                                                                                                                                                                                                               tGACAAAGCCAACCGCTCCTCCATCCACTTTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCCAGGAGTTGAAGCAGACTTTTGGCGTGAAAATTATCACCGACGTTCACGAACCAAGTCAGGCACAGCCCGTTGCTGATGTCGTGGATGTGATTCAGTTGCCGGCGTTTCTTGCTCGCCAGACTGACCTGGTTGAAGCCATGGCGAAAACCGGTGCGGTAATTAACGTCAAGAAACCACAGt            >  1:17169/1‑251 (MQ=255)                                                                                                                                                                                                                         aaCCGCTCCTCCATCCACTTGTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCCAGGAGTTGAAGCAGACTTTTGGCGTGAAAATTATCACCGACGTTCACGAACCAAGTCAGGCACAGCCCGTTGCTGATGTCGTGGATGTGATTCAGTTGCCGGCGTTTCTTGCTCGCCAGACTGACCTGGTTGAAGCCATGGCGAAAACCGGTGCGGTAATTAACGTCAAGAAACCACAGTTTGTCAGccc  <  2:17169/251‑1 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                        AATCGACATTTTACTCAAGATTAAGGCGATCCTATGAAACAAAAAGTGGTTAGCATTGGCGACATCAACGTAGCAAATGACCTGCCGTTCGTACTGTTTGGCGGTATGAACGTGTTGGAATCTCGCGATCTGGCGATGCGCATTTGCGAGCACTACGTAACTGTGACCCAGAAACTGGGTATCCCTTACGTGTTCAAAGCCTCTTTTGACAAAGCCAACCGCTCCTCCATCCACTCTTATCGTGGACCGGGCCTGGAAGAAGGGATGAAAATCTTCCAGGAGTTGAAGCAGACTTTTGGCGTGAAAATTATCACCGACGTTCACGAACCAAGTCAGGCACAGCCCGTTGCTGATGTCGTGGATGTGATTCAGTTGCCGGCGTTTCTTGCTCGCCAGACTGACCTGGTTGAAGCCATGGCGAAAACCGGTGCGGTAATTAACGTCAAGAAACCACAGTTTGTCAGCCC  >  NC_000913/1268132‑1268598

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