Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,219,752	C→T	S386F (TCC→TTC)	aceK →	isocitrate dehydrogenase kinase/phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,219,752	0	C	T	96.8%	104.9 / ‑5.3	32	S386F (TCC→TTC)	aceK	isocitrate dehydrogenase kinase/phosphatase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (17/13);  minor base G (0/1);  total (18/14)
Fisher's exact test for biased strand distribution p-value = 4.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01

ACAGGGCTGTAATGAGCAGTTCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGG‑AACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGCCAGCTTGCCGCTGCTAACATTTTCCCTGGCGACATGCTGTTTAAAAACTTCGGTGTCACCC  >  NC_000913/4219510‑4219994                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    aCAGGGCTGTAATGAGCAGTTCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCgg                                                                                                                                                                                                                                                <  1:144064/248‑1 (MQ=255) aCAGGGCTGTAAGGAGCAGTTCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCa                                                                                                                                                                                                                                              <  1:462787/250‑1 (MQ=255)     gtctgaaaTGAGCAGTACATTGAAGCGCCGGGTATTCGTGGAATGGGGATGTTGGAGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACACAGGAGTTTGAAAACTTTGAGCTGGAGAAGCGGCATAGTTTCCCGGCAGTAAt                                                                                                                                                                                                                                         >  2:315274/7‑251 (MQ=255)            aTGAGCAGATCATTGAAGAGACGGGCATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACCCCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCTTATTTTACCGGCATTACTGGAAtt                                                                                                                                                                                                                                   >  2:325516/1‑250 (MQ=255)               aGCAGTTCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTg                                                                                                                                                                                                                               <  1:6055/251‑1 (MQ=255)                   gTTCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTc                                                                                                                                                                                                                           <  1:148705/251‑1 (MQ=255)                              cgcCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCg                                                                                                                                                                                                                >  1:370683/1‑251 (MQ=255)                                          gTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATATTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGATAAAATCAcc                                                                                                                                                                                                    >  2:6784/1‑251 (MQ=255)                                           tGGAATGGTGATGTTGGTGTTTACGATGCCGGGCTTTGATCGGGTATTCAAAGTCCTCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCCACTGGTGAAAGAGCACGATCGCGTGTGCCGAATGGCGGACAACCAGGAGTTTGTCAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCCTTCATGGAATTACTGCTTCTGGAAGCAGCGGAAAAAATCACCg                                                                                                                                                                                                   >  2:203810/1‑251 (MQ=255)                                                       gTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAAc                                                                                                                                                                                       <  1:325516/251‑1 (MQ=255)                                                                     tgtcGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACAACCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGAGAACAAATTGTGATTCGc                                                                                                                                                                         >  2:369294/4‑251 (MQ=255)                                                                                   cGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTg                                                                                                                                                            >  2:277099/1‑250 (MQ=255)                                                                                         ttCAAAGTCATCAAAGACAAGTTCGCGCCGCAGAACGAGATGTACGCCGCTCACGTTCGTGCGTGCTATCAAATGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACAACCAGGTGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACAGATCTCGGCGACCAAATTGTGATTTGCCATCTTTATATTGTgcggcg                                                                                                                                                     >  2:407169/1‑251 (MQ=255)                                                                                            aaaGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTAACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGACAAAATCACCGATCTCGGCGAACAAATTGGGATTCGCCATCTTTATATTGAGCGGCGGAt                                                                                                                                                  >  1:169268/1‑251 (MQ=255)                                                                                             aaGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTTCCTGCTATCAACTGGTGAAAGAGCACGAGCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCATCGGAAAAAATCACCTATCTCGGCGAACAAATTGTTATTCGCCATCTTTATATTGAGTGTCGGATg                                                                                                                                                 <  1:372235/251‑1 (MQ=255)                                                                                                 catcaacaaCAGGTTCGCGCCGCAGAAAGAGATGTCTGACGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGc                                                                                                                                             >  2:56522/9‑251 (MQ=255)                                                                                                                     cgcAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATg                                                                                                                                                                                                                                        <  1:404571/140‑1 (MQ=255)                                                                                                                     cgcAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATg                                                                                                                                                                                                                                        >  2:404571/1‑140 (MQ=255)                                                                                                                                         gCTCACGTTAGAGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGAGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGGGCCGCTCAATATCTGGCTGG‑AAAAAGTGGAAGGTCAGCAGtt                                                                                                    >  2:333039/1‑251 (MQ=255)                                                                                                                                                     gCCTGCTATCAACTGGTGAAAGAGCCCGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGCAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGTAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGG‑AACAAGTGGACGGTCAGCAGGTGCGCGACGCCAt                                                                                        >  2:333520/1‑251 (MQ=255)                                                                                                                                                      ccTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGG‑AACAAGTGGAAGGTCAGCAGTTGCGCGACGCCAtt                                                                                       <  2:457791/251‑1 (MQ=255)                                                                                                                                                        tGCTATCAACTGGTGAAAGAGCGAGCTCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGGGAAGCGGCATATTTTCCCGGCATTACTGGCATTACTGCTTCAGGAAGCAGCGGAAAAAATc                                                                                                                                                                                                       >  2:288618/1‑138 (MQ=255)                                                                                                                                                        ggCTATCAACTGGTGAAAGAGCACGAGCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTGCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATc                                                                                                                                                                                                       <  1:288618/137‑1 (MQ=255)                                                                                                                                                                   ggTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGGGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCCCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGG‑AACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATAAGGTa                                                                          <  1:343134/251‑1 (MQ=255)                                                                                                                                                                                   cgcgTGGGCCGAAAGGCGGACACCAAGGACCTTGAAAACTTTGTGCTCGAGAAGCGGCCTATTACGCCGGCGTTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCATCGCTATCGGCGCACAAATTGCGATTCGCAATCTTTATATTGAGCGGCGGATGGGGCCGAGCAATATCAGGCTGGGAAAAAG‑GGACGGTCAGCCGTTGCCCGCCGCCCTTTAAGAATAGGGTCACGCCATTCggccgct                                                          >  2:263059/1‑245 (MQ=255)                                                                                                                                                                                                               atttagAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTACTGGAAGTACTGATTCAGGAAGCAGCGGAAAAAATCACCGAGCTCGGCGAACAAATTGTGATGCGACATCTGTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGG‑AACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATGCGCCAGCTTGCCGCTGCTATCATTTTCCCTGGCGAc                              <  2:338986/246‑1 (MQ=255)                                                                                                                                                                                                                 tttGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGG‑CACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGCCAGCTTGCCGCTGCTAACATTTTCCCTGGCGACAt                            >  2:106613/1‑251 (MQ=255)                                                                                                                                                                                                                            tgtgCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGG‑AACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGCCAGCTTGCCGCTGCTAACATTTTCCCTGGCGACATGCTGTTTaaa                  >  1:285087/1‑250 (MQ=255)                                                                                                                                                                                                                             gtgCTGGAAAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGG‑AACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGCCAGCTTGCCGCTGCTAACATTTTCCCTGGCGACATGCTGTTTaaaaa                <  2:169268/251‑1 (MQ=255)                                                                                                                                                                                                                                  ggAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACATATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGG‑AACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGCCAGCTTGCCGCTGCTAACATTTTCCCTGGCGACATGCTGTTTAAAAACTTCg           <  1:99317/251‑1 (MQ=255)                                                                                                                                                                                                                                          ggCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAATCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCAGCTCAATATCTGGCTGG‑AACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGCCAGCTTGCCGCTGCTAACATATTCCCTGGCGACATGCTGTTTAAAAACTTCGGTGTCAcc   >  2:368380/1‑251 (MQ=255)                                                                                                                                                                                                                                           gCATATTTTCCCGGCATTAATGTAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAAGTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGG‑AACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAATAATACGGTAACGCTATTCGCCAGCTTGCCGCTGCTAACATTTTCCCTGGCGACATGCTGTTTAAAAACTTCGGTGTCAccc  <  1:208644/251‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    ACAGGGCTGTAATGAGCAGTTCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGG‑AACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGCCAGCTTGCCGCTGCTAACATTTTCCCTGGCGACATGCTGTTTAAAAACTTCGGTGTCACCC  >  NC_000913/4219510‑4219994

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