Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,715,661	C→T	D77N (GAT→AAT)	pdxY ←	pyridoxal kinase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,715,661	0	C	T	97.0%	106.3 / NA	33	D77N (GAT→AAT)	pdxY	pyridoxal kinase 2
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (19/13);  minor base G (0/1);  total (19/14)
Fisher's exact test for biased strand distribution p-value = 4.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCAGCTCAACCAGATTTGGCGCAATGATATCGCTGGCAGGCAAACCGTGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATGACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTCACCCTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACA‑ATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCT  >  NC_000913/1715438‑1715896                                                                                                                                                                                                                                |                                                                                                                                                                                                                                             tCCAGCGCAACCAGATTGGGCGCAATGATAGCGCTGGCAGTCAAACCGTGCCGCACATGAAACTCTGCGTCACCCGGTGCAAAGATACAGCCGTTTTCCGGAGGACCCAGTACCGGATCGCAAAAAGAGTGCTCCTGCGGTGCCGCGGCTTTCACCTGACGGACGATACCGAGGATAGGTCCACCCTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATggcgg                                                                                                                                                                                                                    <  1:296378/250‑1 (MQ=255)  ccAGCTCAACCAGATTTGGCGCAATGATATCGCTGGCAGGCAAACCGTGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATGACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTCACCCTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATggcggc                                                                                                                                                                                                                   >  2:494867/1‑250 (MQ=255)         aaCCAGATTTGGCGCAATGATATCGCTGGCAGGCAAACCGTGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATGACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTCACCCTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCtt                                                                                                                                                                                                           >  1:402632/1‑251 (MQ=255)                              ctCGCTGGCAGGCAAACCGTGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATCCAGCCTTTTTCCGGATGACCCATTACCGGATCGCACCAATATTTCGACTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTCACCCTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCTTTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGC‑CACATTTCTGTTTAATGGc                                                                                                                                                                                     >  2:209839/2‑251 (MQ=255)                                           aacccgggCCGCACATGAAACGCGGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATGACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTCACCCTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACA‑ATTTCGGTTAAATGGCTGGGCGGCATCAc                                                                                                                                                                        <  2:143456/244‑1 (MQ=255)                                                        cccgatcttctGCGACCCCCGGTTCAACGATACAGCCTTTTTCCGGATGACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATCCCGAGGATATGTTCCCCCTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACA‑ATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCAt                                                                                                                                                           <  2:72564/243‑1 (MQ=255)                                                                tctGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATGACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTCACCCTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACA‑ATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCAt                                                                                                                                                           >  1:72564/1‑243 (MQ=255)                                                                       cccccGGTGCAACAATACATCCTTTTCCCGGATGACCCATTACCGGACCGAAAAAATATTTCGCCAGCGGATTCGCTGCTTTCCCCGGCCGTACGATTCCGAGGATTTGTTCCCCCTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACA‑ATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATGTGCCGTATGGGgtg                                                                                                                                            <  2:517708/249‑1 (MQ=255)                                                                                           ccTTTTTCAGGATAACCCATTACCGGATCACAACCACATTTCCCCTGCGCATTCGCGGATTTCACCTGACCGACGATACCGAGGATCTGTTCACCCTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTACTTTATCAATGGCGGCAATTCCTTGCACA‑ATTTCGGTTAAATGGCTGGGCGGCATCTCGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTCGAAAATTGAACggg                                                                                                                        >  2:560765/1‑250 (MQ=255)                                                                                                  ccggaggccccatgtccgaaTCGCAACAATATGCCGCCTGCGAGTTCGCGGCTTTCCCTTGACGGTCGATCCGGAGTATATGTTCACCCTGCTCCGCCGATCCCAGATAGCCACTTAATTCGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACA‑ATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAg                                                                                                                  <  2:463076/231‑1 (MQ=255)                                                                                                           ccATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTCACCCTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACA‑ATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACg                                                                                                        >  1:370889/1‑251 (MQ=255)                                                                                                               tACCGGATCGCAAAACTCTTTAGCCTGCGTATTCGCGGCTTTCACATGACGGACGATACCGAGGATATGTTCACCCTGCTCAGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCATTGCACA‑ATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGACGTATTGGGCGTGATTAGAAAATTGAACGGGGATCAGCGGCCAGTCGTTcg                                                                                                    >  2:351168/1‑251 (MQ=255)                                                                                                               tACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCCCCTGACGGACGATACCGAGGATATGTTCACCCTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACA‑ATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTcg                                                                                                    >  2:73194/1‑251 (MQ=255)                                                                                                                  cGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTCACCCTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACA‑ATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCgcgc                                                                                                 >  1:379601/1‑251 (MQ=255)                                                                                                                                atTTCGCCTGCGGATTCGTGGCTTGCACCTGACGGCCGATACCGGGGATATGTTCCCCCGGCTCCGCCGATCCCAGATAGCCACGTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAAGGCCTTGCACA‑ATGTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATGAGAAAATGGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCATGCGGCGCAt                                                                                    <  2:523503/250‑1 (MQ=255)                                                                                                                                atTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTCACCCTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACA‑ATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATc                                                                                   >  1:238860/1‑251 (MQ=255)                                                                                                                                           ggctccgcgGCTTTCACCTGACGGACGATACCGAGGATATGTTCCCCCTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACA‑ATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGc                                                                        <  2:279263/246‑1 (MQ=255)                                                                                                                                                          aCCTGACGGACGATACCGAGGCTATGTTCATCCTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACA‑ATTTCGGTTAAATGGCTGGGCGGCATCACACAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGACGGCGCATCGGAAACTCTGCCGCAATGTTACTCGc                                                         >  1:454529/1‑251 (MQ=255)                                                                                                                                                                      aTACCGAGGATATGTTCAACCTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACA‑ATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCCTCGGAAACTCTGCCGCACTGTTACCCGCATGACCATaaa                                              >  1:112327/1‑250 (MQ=255)                                                                                                                                                                        aCCGAGTATATGTTCCCCCTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACA‑ATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGtt                                                                                                                     <  2:392413/177‑1 (MQ=255)                                                                                                                                                                        aCCGAGGATATGTTCACCCTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACA‑ATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGtt                                                                                                                     >  1:392413/1‑177 (MQ=255)                                                                                                                                                                             ggATATGTTCCCCCTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACA‑ATTTCGGTTAAATGGCGGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTga                                      <  2:553496/251‑1 (MQ=255)                                                                                                                                                                               atagtgccACCCTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATGACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACA‑ATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTgaga                                    <  2:171866/244‑1 (MQ=255)                                                                                                                                                                                     tCACCCGGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACA‑ATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGa                               <  2:379601/250‑1 (MQ=255)                                                                                                                                                                                        cccggccccgccgATCCCAGATAGCCACTTATTACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACA‑ATTTCGGTTAAATGGCTGGGCTGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCAGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGAt                              <  2:464207/241‑1 (MQ=255)                                                                                                                                                                                             cTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACA‑ATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGCCTGGATAGCgaga                       >  2:430432/1‑250 (MQ=255)                                                                                                                                                                                                agccgATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACA‑ATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTATTACCCGCATGACCATAAACAACGGGAGACTGGATAACGAGAatat                   >  2:422297/2‑251 (MQ=255)                                                                                                                                                                                                  ccgATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACA‑ATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATAttt                 >  1:95611/1‑251 (MQ=255)                                                                                                                                                                                                       cccAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACA‑ATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTtcat            <  1:327209/251‑1 (MQ=255)                                                                                                                                                                                                              aGCCACTTAATACGGCATTACAGGTGTGTAATTTAGCAATGGCGGCAATGCCTTGCACA‑ATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTc     <  2:310999/251‑1 (MQ=255)                                                                                                                                                                                                               gCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACA‑ATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTcc    >  1:58396/1‑251 (MQ=255)                                                                                                                                                                                                               gCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACA‑ATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGGGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTcc    >  1:31852/1‑251 (MQ=255)                                                                                                                                                                                                                 cACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACA‑ATTTCGGTTCAATGGCTGGGCGGCATCCCGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCCGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCt  >  2:557085/1‑251 (MQ=255)                                                                                                                                                                                                                                |                                                                                                                                                                                                                                             TCCAGCTCAACCAGATTTGGCGCAATGATATCGCTGGCAGGCAAACCGTGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATGACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTCACCCTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACA‑ATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCT  >  NC_000913/1715438‑1715896

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: ‑