Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	351,119	C→T	P370S (CCT→TCT)	prpC →	2‑methylcitrate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	351,119	0	C	T	93.8%	105.2 / NA	32	P370S (CCT→TCT)	prpC	2‑methylcitrate synthase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (15/15);  minor base G (0/2);  total (15/17)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

CGCAGGAAGGCGGCTCGCTGAAGATGTACAACATCGCCGATCGCCTGGAAACGGTGATGTGGGAGAGCAAAAAGATGTTCCCCAATCTCGACTGGTTCTCTGCTGTTTCCTACAACATGATGGGCGTTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAACGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGATAAGCGCCAGTAAACCTCTACGAATAACAATAAGGAAACGTACCCAATGTCAGCTCAAATCAACAACATCCGCCCGGAATTTGATCGTGAAATCGTTGATATCGTCGATTACGTCATGAACTACGAAATCAGCTCTAAAGTGGCCTACGACACCGCACATTACTGCCTGCTC  >  NC_000913/350895‑351340                                                                                                                                                                                                                                 |                                                                                                                                                                                                                              cGCAGGAAGGCGGCTCGCTGAAGATGTACAACATCGCCGATCGCCTGGAAACGGTGATGTGGGAGAGCAAAAAGATGTTCCCCAATCTCGACTGGTTCTCTGCTGTTTCCTACAACATGATGGGCGTTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAACGTCAGGACAACAAAATTATCCGTTCTTCCGCCAATTATGTTGGACCGGa                                                                                                                                                                                                      <  1:121092/250‑1 (MQ=255)   cAGGAAGGCGGCTCGCTGAAGATGTACAACATCGCCGATCGCCTGGAAACGGTGATGTGGGAGAGCAAAAAGATGTTCCCCAATCTCGACTGGTTCTCTGCTGTTTCCTACAACATGATGGGCGTTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAACGTCAGGACAACAAAATTATCCGTTCTTCCGCCAATTATGTTGGACCGGa                                                                                                                                                                                                      >  1:393901/1‑248 (MQ=255)       aaGGCGGCTCGCTGAAGATGTACAACATCGCCGATCGCCTGGAAACGGTGATGTGGGAGAGCAAAAAGATGTTCCCCAATCTCGCCTGGTTCTCTGCTGTTTCCTACAACATGATGGGCGTTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAACGTCAGGACAACAAAATTATCCGTTCTTCCGCCACTTATGTTGGACCGGAAGACCGc                                                                                                                                                                                               >  1:327242/1‑251 (MQ=255)         ggcggcTCGCTGAAGATGTACAACATCGCCGATCGCCTGGAAACGGTGATGTGGGAGAGCAAAAAGATGTTCCCCAATCTCGACTGGTTCTCTGCTGTTTCCTACAACATGATGGGCGTTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAACGTCAGGACAACAAAATTATCCGTTCTTCCGCCAATTATGTTGGACCGGAAGACCGccc                                                                                                                                                                                             <  2:121444/251‑1 (MQ=255)           cggcTCGCTGAAGATGTACAACATCGCCGATCGCCTGGAAACGGTGATGTGGGAGAGCAAAAAGATGTTCCCCAATCTCGACTGGTTCTCTGCTGTTTCCTACAACATGATGGGCGTTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAACGTCAGGACAACAAAATTATCCGTTCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGt                                                                                                                                                                                           >  1:71798/1‑251 (MQ=255)                cGCTGAAGATGTACAACATCGCCGATCGCCTGGAAACGGTGATGTGGGAGAGCAAAAAGATGTTCCCCAATCTCGACTGGTTCTCTGCTGTTTCCTACAACATGATGGGCGTTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAACGTCAGGACAACAAAATTATCCGTTCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTg                                                                                                                                                                                        <  2:16027/249‑1 (MQ=255)                agttGAATATGTACAACATCGCCGATCGCCGGGATACTGTTATGTGGGAGAGCAAAATGATGTTTCCCACTCTCGACTGGTTTTCTGCTGTTTTCTTTGTTAAGAATACTGCTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAACGTCAGGACAACAAAATTATCCGTTCTTCCGCCAATTATGTGGGACCGGAAGACCGCCCGTTTGt                                                                                                                                                                                       <  2:411680/247‑1 (MQ=255)                    gAAGATGTACAACAGCGCCGATCGCCTGGAAACGGTGATGTGGGAGAGCAAAAAGATGTTCCCCAATCTCGACTGGTTCTCTGCTGTTTCCTACAACATGATGGGCGTTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAACGTCAGGACAACAAAATTATCCGTTCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTcgc                                                                                                                                                                                    <  2:373855/249‑1 (MQ=255)                     aaGATGTACAACATCGCCGATCGCCTGGAAACGGTGATGTGGGAGAGCAAAAAGATGTTCCACAATCTCGACTGGTTCTCTGCTGTTTCCTACAACATGATGGGCGTTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAACGTCAGGACAACAAAATTATCCGTTCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCt                                                                                                                                                                                 >  1:193349/1‑251 (MQ=255)                        aTGTACAACATCGCCGATCGCCTGGAAACGGTGATGTGGGAGAGCAAAAAGATGTTCCCCAATCTCGACTGGTTCTCTGCTGTTTCCTACAACATGATGGGCGTTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAACGTCAGGACAACAAAATTATCCGTTCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGa                                                                                                                                                                              <  1:61043/251‑1 (MQ=255)                         tGTACAACATCGCCGATCGCCTGGAAACGGTGATGTGGGAGAGCAAAAAGATGTTCCCCAATCTCGACTGGTTCTCTGCTGTTTCCTACAACATGATGGGCGTTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAACGTCAGGACAACAAAATTATCCGTTCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGAt                                                                                                                                                                             >  1:368024/1‑251 (MQ=255)                            acaacaTCGCCGATCGCCTGGAAACGGTGATGTGGGAGAGCAAAAAGATGTTCCCCAATCTCGACTGGTTCTCTGCTGTTTCCTACAACATGATGGGCGTTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAACGTCAGGACAACAAAATTATCCGTGCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGATaa                                                                                                                                                                           <  1:393015/250‑1 (MQ=255)                                 aTCGCCGATCGCCTGGAAACGGTGATGTGGGAGAGCAAAAAGATGTTCCCCAATCTCGACTGGTTCTCTGCTGTTTCCTACAACATGATGGGCGTTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAACGTCAGGACAACAAAATTATCCGTTCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGATAAGCGCCa                                                                                                                                                                     <  1:74120/251‑1 (MQ=255)                                   cGCCGATCGCCTGGAAACGGTGATGTGGGAGAGCAAAAAGATGTTCCCCAATCTCGACTGGTTCTCTGCTGTTTCCTACAACATGATGGGCGTTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAACGTCAGGACAACAAAATTATCCGTTCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGATTGTCTCGCTGGATAAGCGCCAGt                                                                                                                                                                   >  2:6518/1‑251 (MQ=255)                                               ggAAACGGTGATGTGGGAGAGCAAAAAGATGTTCCCCAATCTCGACTGGTTCTCTGCTGTTTCCTACAACATGATGGGCGTTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAACGTCAGGACAACAAAATTATCCGTTCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGATAAGCGCCAGTAAACCTCTAc                                                                                                                                                         <  2:270466/249‑1 (MQ=255)                                                    cGGTGATGTGGGAGAGCAAAAAGATGTTCCCAAATCTCGACTGGTTCTCTGCTGTTTCCTACAACATGATGGGCGTTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAACGTCAGGACAACAAAATTATCCGTTCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGATAAGCGCCAGTAAACCTCTACGAATAAc                                                                                                                                                  >  1:324111/1‑251 (MQ=255)                                                             gggAGAGCAAAAAGATGTGCCCCAATATCGACTGGTTCTCTGCTGTTTCCTACAACATGATGGGCGGTCCCACCGAGATGTTCACACCGCTGTGTGTTATCGCCCGCGTCACCGTCTGGGCGTCGCACATTATCGAACAACGTCAGGACAACAAAAGTATCCGTGCTTCCGACAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCACTGGATAAGCGCCAGTAAACCTCTACGAATAACAATAAGGaa                                                                                                                                         <  2:306632/251‑1 (MQ=255)                                                                                         cGACTGGTTCTCTGCTGTTTCCTACAACATGATGGGCGTTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAACGTCAGGACAACAAAATTATCCGTTCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGATAAGCGCCAGTAAACCTCTACGAATAACAATAAGGAAACGTACCCAATGTCAGCTCAAATc                                                                                                                 <  1:322952/247‑1 (MQ=255)                                                                                          gACTGGTTCTCTGCTGTTTCCTACAACATGATGGGCGTTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAACGTCAGGACAACAAAATTATCCGTTCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGATAAGCGCCAGTAAACCTCTACGAATAACAATAAGGAAACGTACCCAATGTCAGCTCAAATcaaca                                                                                                             >  1:77928/1‑250 (MQ=255)                                                                                                         gTTTCCTACAACATGATGGGCGTTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAACGTCAGGACAACAAAATTATCCGTTCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGATAAGCGCCAGTAAACCTCTACGAATAACAATAAGGAAACGTACCCAATGTCAGCTCAAATCAACAACATCCGCCCGGaa                                                                                               <  2:258609/249‑1 (MQ=255)                                                                                                            tCCTACAACATGATGGGCGTTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAACGTCAGGACAACAAAATTATCCGTTCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGATAAGCGCCAGTAAACCTCTACGAATAACAATAAGGAAACGTACCCAATGTCAGCTCAAATCAACAACATCCGCCCGGAATTTg                                                                                           <  2:368024/250‑1 (MQ=255)                                                                                                                 caacaTGATGGGCGTTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAACGTCAGGACAACAAAATTATCCGTTCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGATAAGCGCCAGTAAACCTCTACGAATAACAATAAGGAAACGTACCCAATGTCAGCTCAAATCAACAACATCCGCCCGGAATTTGATCGTg                                                                                     >  1:69065/1‑251 (MQ=255)                                                                                                                     agggtgGGCGCTCCCACCGGGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCGGGGCGGCGCACATTATCGAACAACGTCAGGACAACAAAATTATCCGTTCTTCCGCCAATTATGTTGGACCGGGAGACCGCCCGTTTGTCGCGCTGGATAAGCGCCAGTAAACCTCTACGAATAACAATAAGGAAACTTACCCAATGTCAGCTCAAAGCAACAACATCCGCCCTGAATTTGATCGGGAAAt                                                                                 <  2:56351/247‑1 (MQ=255)                                                                                                                          gggCGTTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAACGTCAGGACAACAAAATTATCCGTTCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGATAAGCGCCAGTAAACCTCTACGAATAACAATAAGGAAACGTACCCAATGTCAGCTCAAATCAACAACATCCGCCCGGAATTTGATCGTGAAATCGTTg                                                                            >  1:42848/1‑251 (MQ=255)                                                                                                                                 cccACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAACGTCAGGACAACAAAATTATCCGTTCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGATAAGCGCCAGTAAACCTCTACGAATAACAATAAGGAAACGTACCCAATGTCAGCTCAAATCAACAACATCCGCCCGGAATTTGATCGTGAAATCGTTGATAtcgt                                                                     >  2:82279/1‑251 (MQ=255)                                                                                                                                  ccACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAACGTCAGGACAACAAAATTATCCGTTCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGATAAGCGCCAGTAAACCTCTACGAATAACAATAAGGAAACGTACCCAATGTCAGCTCAAATCCACAACATCCGCCCGGAATTTGATCGTGAAATCGTTGATAtcgtc                                                                    >  1:404194/1‑251 (MQ=255)                                                                                                                                                 caccacTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAACGTCAGGACAACAAAATTATCCGTTCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGATAAGCGCCAGTAAACCTCTAc                                                                                                                                                         >  2:123557/1‑151 (MQ=255)                                                                                                                                                 caccacTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAACGTCAGGACAACAAAATTATCCGTTCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTGGTCGCGCTGGATAAGCGCCAGTAAACCTCTAc                                                                                                                                                         <  1:123557/151‑1 (MQ=255)                                                                                                                                                                           aCCGGCTGGGCGGCGCACATTATCGAACAACGTCAGGACAACAAAATTATCCGTTCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGATAAGCGCCAGTAAACCTCTACGAATAACAATAAGGAAACGTACCCAATGTCAGCTCAAATCAACAACATCCGCCCGGAATTTGATCGTGAAATCGTTGATATCGTCGATTACGTCATGAACTACGAAATCAGCTCTAAAGTGGCCTa                           <  2:267561/251‑1 (MQ=255)                                                                                                                                                                                 tGGGCGGCGCACATTATCGAACAACGTCAGGACAACAAAATTATCCGTTCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGATAAGCGCCAGTAAACCTCTACGAATAACAATCAGGAAACGTACCCAATGTCAGCTCAAATCAACAACATCCGCCCGGAATTTGATCGTGAAATCGTTGATATCGTCGATTACGTCATGAACTACGAAATCAGCTCTAAAGTGGCCTACGacac                     >  2:268036/1‑251 (MQ=255)                                                                                                                                                                                           acaTTATCGAACAACGTCAGGACAACAAAATTATCCGTTCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGATAAGCGCCAGTAAACCTCTACGAATAACAATAAGGAAACGTACCCAATGTCAGCTCAAATCAACAACATCCGCCCGGAATTTGATCGTGAAATCGTTGATATCGTCGATTACGTCATGAACTACGAAATCAGCTCTAAAGTGGCCTACGACACCGCACATTAc           >  2:406681/1‑251 (MQ=255)                                                                                                                                                                                                    aacaacGTCAGGACAACAAAATTATCCGTTCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGATAAGCGCCAGTAAACCTCTACGAATAACAATAAGGAAACGTACCCAATGTCAGCTCAAATCAACAACATCCGCCCGGAATTTGATCGTGAAATCGTTGATATCGTCGATTACGTCATGAACTACGAAATCAGCTCTAAAGTGGCCTACGACACCGCACATTACTGCCTGCTc  <  2:77928/251‑1 (MQ=255)                                                                                                                                                                                                                                 |                                                                                                                                                                                                                              CGCAGGAAGGCGGCTCGCTGAAGATGTACAACATCGCCGATCGCCTGGAAACGGTGATGTGGGAGAGCAAAAAGATGTTCCCCAATCTCGACTGGTTCTCTGCTGTTTCCTACAACATGATGGGCGTTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAACGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGATAAGCGCCAGTAAACCTCTACGAATAACAATAAGGAAACGTACCCAATGTCAGCTCAAATCAACAACATCCGCCCGGAATTTGATCGTGAAATCGTTGATATCGTCGATTACGTCATGAACTACGAAATCAGCTCTAAAGTGGCCTACGACACCGCACATTACTGCCTGCTC  >  NC_000913/350895‑351340

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