Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	351,040	A→G	T343T (ACA→ACG)	prpC →	2‑methylcitrate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	351,040	0	A	G	97.5%	145.3 / NA	40	T343T (ACA→ACG)	prpC	2‑methylcitrate synthase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (17/22);  total (17/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GGCGGCTCGCTGAAGATGTACAACATCGCCGATCGCCTGGAAACGGTGATGTGGGAGAGCAAAAAGATGTTCCCCAATCTCGACTGGTTCTCTGCTGTTTCCTACAACATGATGGGCGTTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAACGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGATAAGCGCCAGTAA  >  NC_000913/350903‑351181                                                                                                                                          |                                                                                                                                              ggcggcTCGCTGAAGATGTACAACATCGCCGATCGCCTGGAAACGGTGATGTGGGAGAGCAAAAAGATGTTCCCCAATCTCGACTGGTTCTCTGCTGTTTCCTACAACATGATGGGCGTTCCCACCGAGATGTTCACGCCACTGTTTGTTa                                                                                                                                  <  1:917486/151‑1 (MQ=255)    ggcTCGCTGAAGATGTACAACATCGCCGATCGCCTGGAAACGGTGATGTGGGAGAGCAAAAAGATGTTCCCCAATCTCGACTGGTTCTCTGCTGTTTCCTACAACATGATGGGCGTTCCCACCGAGATGTTCACGCCACTGTTTGTTATCg                                                                                                                               <  1:206458/151‑1 (MQ=255)       tCGCTGAAGATGTACAACATCGCCGATCGCCTGGAAACGGTGATGTGGGAGAGCAAAAAGATGTTCCCCAATCTCGACTGGTTCTCTGCTGTTTCCTACAACATGATGGGCGTTCCCACCGAGATGTTCACGCCACTGTTTGTTATCGccc                                                                                                                            >  2:922213/1‑151 (MQ=255)           tGAAGATGTACAACATCGCCGATCGCCTGGAAACGGTGATGTGGGAGAGCAAAAAGATGTTCCCCAATCTCGACTGGTTCTCTGCTGTTTCCTACAACATGATGGGCGTTCCCACCGAGATGTTCACGCCACTGTTTGTTATCGCCCGCGt                                                                                                                        >  1:867798/1‑151 (MQ=255)                aTGTACAACATCGCCGATCGCCTGGAAACGGTGATGTGGGAGAGCAAAAAGATGTTCCCCAATCTCGACTGGTTCTCTGCTGTTTCCTACAACATGATGGGCGTTCCCACCGAGATGTTCACGCCACTGTTTGTTATCGCCCGCGTCACCg                                                                                                                   >  1:283613/1‑151 (MQ=255)                      aacaTCGCCGATCGCCTGGAAACGGTGATGTGGGAGAGCAAAAAGATGTTCCCCAATCTCGACTGGTTCTCTGCTGTTTCCTACAACATGATGGGCGTTCCCACCGAGATGTTCACGCCACTGTTTGTTATCGCCCGCGTCACCGGCTggg                                                                                                             >  1:438841/1‑151 (MQ=255)                         aTCGCCGATCGCCTGGAAACGGTGATGTGGGAGAGCAAAAAGATGTTCCCCAATCTCGACTGGTTCTCTGCTGTTTCCTACAACATGATGGGCGTCCCCACCGAGATGTTCACGCCACTGTTTGTTATCGCCCGCGTCACCGGCTGggc                                                                                                            <  2:980728/149‑1 (MQ=255)                           cGCCGATCGCCTGGAAACGGTGATGTGGGAGAGCAAAAAGATGTTCCCCAATCTCGACTGGTTCTCTGCTGTTTCCTACAACATGATGGGCGTTCCCACCGGGATGTTCACGCCACTGTTTGTTATCGCCCGCGTCACCGGCTGGgcggcg                                                                                                        <  2:932932/151‑1 (MQ=255)                           cGCCGATCGCCTGGAAACGGTGATGTGGGAGAGCAAAAAGATGTTCCCCAATCTCGACTGGTTCTCTGCTGTTTCCTACAACATGATGGGCGTTCCCACCGAGATGTTCACGCCACTGTTTGTTATCGCCCGCGTCACCGGCTGGgcggcg                                                                                                        <  2:831254/151‑1 (MQ=255)                                     cTGGAAACGGTGATGTGGGAGAGCAAAAAGATGTTCCCCAATCTCGACTGGTTCTCTGCTGTTTCCTACAACATGATGGGCGTTCCCACCGAGATGTTCACGCCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCg                                                                                              >  1:337141/1‑151 (MQ=255)                                     cTGGAAACGGTGATGTGGGAGAGCAAAAAGATGTTCCCCAATCTCGACTGGTTCTCTGCTGTTTCCTACAACATGATGGGCGTTCCCACCGAGATGTTCACGCCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCg                                                                                              >  1:913707/1‑151 (MQ=255)                                        gAAACGGTGATGTGGGAGAGCAAAAAGATGTTCCCCAATCTCGACTGGTTCTCTGCTGTTTCCTACAACATGATGGGCGTTCCCACCGAGATGTTCACGCCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGaac                                                                                           <  2:877088/151‑1 (MQ=255)                                         aaaGGGTGATGTGGGAGGGCAAAAAGATGTTCCCCAATCTCGACTGGTTCTCTGCTGTTTCCTACAACATGATGGGCGTTCCCCCCGAGATGTTCACGCCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGaac                                                                                           <  1:540146/150‑1 (MQ=255)                                            cGGTGATGTGGGAGAGCAAAAAGATGTTCCCCAATCTCGACTGGTTCTCTGCTGTTTCCTACAACATGATGGGCGTTCCCACCGAGATGTTCACGCCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAACg                                                                                       >  1:745185/1‑151 (MQ=255)                                                 aTGTGGGAGAGCAAAAAGATGTTCCCCAATCTCGACTGGTTCTCTGCTGTTTCCTACAACATGATGGGCGTTCCCACCGAGATGTTCACGCCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAACGTCAgg                                                                                  >  2:1024253/1‑151 (MQ=255)                                                       gagaGCAAAAAGATGTTCCCCAATCTCGACTGGTTCTCTGCTGTTTCCTACAACATGATGGGCGTTCCCACCGAGATGTTCACGCCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAACGTCAGGacaac                                                                             >  1:316673/1‑150 (MQ=255)                                                        agagCAAAAAGATGTTCCCCAATCTCGACTGGTTCTCTGCTGTTTCCTACAACATGATGGGCGTTCCCACCGAGATGTTCACGCCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAACGTCAGGAcaacaa                                                                           <  1:983377/151‑1 (MQ=255)                                                         gagCAAAAAGATGTTCCCCAATCTCGACTGGTTCTCTGCTGTTTCCTACAACATGATGGGCGTTCCCACCGAGATGTTCACGCCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAACGTCAGGACAACaaa                                                                          <  2:337116/151‑1 (MQ=255)                                                              aaaaGATGTTCCCCAATCTCGACTGGTTCTCTGCTGTTTCCTACAACATGATGGGCGTTCCCACCGAGATGTTCACGCCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAACGTCAGGACAACAAAATTAt                                                                     >  2:452338/1‑151 (MQ=255)                                                              aaaaGATGTTCCCCAATCTCGACTGGTTCTCTGCTGTTTCCTACAACATGATGGGCGTTCCCACCGAGATGTTCACGCCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAACGTCAGGACAACAAAATTAt                                                                     >  2:939300/1‑151 (MQ=255)                                                                   aTGTTCCCCAATCTCGACTGGTTCTCTGCTGTTTCCTACAACATGATGGGCGTTCCCACCGAGATGTTCACGCCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAACGTCAGGACAACAAAATTATCCGTc                                                                <  1:265454/151‑1 (MQ=255)                                                                      ttCCCCAATCTCGACTGGTTCTCTGCTGTTTCCTACAACATGATGGGCGTTCCCACCGAGATGTTCACGCCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAACGTCAGGACAACAAAATTATCCGTCCt                                                              <  2:371510/150‑1 (MQ=255)                                                                         cccAATCTCGACTGGTTCTCTGCTGTTTCCTACAACATGATGGGCGTTCCCACCGAGATGTTCACGCCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAACGTCAGGACAACAAAATTATCCGTCCTTCCg                                                          >  1:998344/1‑151 (MQ=255)                                                                         cccAATCTCGACTGGTTCTCTGCTGTTTCCTACAACATGATGGGCGTTCCCACCGAGATGTTCACGCCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAACGTCAGGACAACAAAATTATCCGTCCTTCCg                                                          <  2:487470/151‑1 (MQ=255)                                                                          ccAATCTCGACTGGTTCTCTGCTGTTTCCTACAACATGATGGGCGTTCCCACCGAGATGTTCACGCCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAACGTCAGGACAACAAAATTATCCGTCCTTCCGc                                                         <  1:780847/151‑1 (MQ=255)                                                                                  gtagggTTCGCTGCTGTTTCCTACAACATGATGGGCGGTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAACGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATg                                                 <  2:207962/147‑1 (MQ=255)                                                                                   aCTGGTTCTCTGCTGTTTCCTACAACATGATGGGCGTTCCCACCGAGATGTTCACGCCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAACGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGt                                                <  1:86417/151‑1 (MQ=255)                                                                                      ggTTCTCTGCTGTTTCCTACAACATGATGGGCGTTCCCACCGAGATGTTCACGCCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAACGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTgg                                             <  1:286865/151‑1 (MQ=255)                                                                                                   ttCCTACAACATGATGGGCGTTCCCACCGAGATGTTCACGCCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAACGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGc                                >  2:703844/1‑151 (MQ=255)                                                                                                        acaacaTGATGGGCGTTCCCCCCGAGATGTTCACGCCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAACGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCg                             <  2:973390/149‑1 (MQ=255)                                                                                                                  gggCGTTCCCACCGAGATGTTCACGCCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAACGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCgcgc                   >  2:672514/1‑149 (MQ=255)                                                                                                                      gTTCCCACCGAGATGTTCACGCCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAACGTCAGGACAACAAAATTATCCGTCCtt                                                             >  1:90172/1‑103 (MQ=255)                                                                                                                       ttCCCACCGAGATGTTCACGCCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAACGTCAGGACAACAAAATTATCCGTCCtt                                                             <  2:90170/102‑1 (MQ=255)                                                                                                                         cccACCGAGATGTTCACGCCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAACGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGATaa            <  1:375875/149‑1 (MQ=255)                                                                                                                         cccACCGAGATGTTCACGCCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAACGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGATAAg           <  2:276653/150‑1 (MQ=255)                                                                                                                             ccGAGATGTTCACGCCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAACGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGATAAGCGcc       <  2:782994/150‑1 (MQ=255)                                                                                                                              cGAGATGTTCACGCCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAACGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGATAAGCGCCAg     <  2:1019250/151‑1 (MQ=255)                                                                                                                                agaTGTTCACGCCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAACGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGATAAGCGCCAGt    >  1:66704/1‑150 (MQ=255)                                                                                                                                 gaTGTTCACGCCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAACGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGATAAGCGCCAGTaa  <  2:402394/151‑1 (MQ=255)                                                                                                                                  aTGTTCACGCCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAACGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGATAAGCGCCAGTaa  >  1:139688/1‑150 (MQ=255)                                                                                                                                          |                                                                                                                                              GGCGGCTCGCTGAAGATGTACAACATCGCCGATCGCCTGGAAACGGTGATGTGGGAGAGCAAAAAGATGTTCCCCAATCTCGACTGGTTCTCTGCTGTTTCCTACAACATGATGGGCGTTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAACGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGATAAGCGCCAGTAA  >  NC_000913/350903‑351181

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