Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	349,825	A→G	intergenic (+253/‑187)	prpB → / → prpC	2‑methylisocitrate lyase/2‑methylcitrate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	349,825	0	A	G	100.0%	112.3 / NA	53	intergenic (+253/‑187)	prpB/prpC	2‑methylisocitrate lyase/2‑methylcitrate synthase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (33/20);  total (33/20)

CGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAA  >  NC_000913/349678‑349971                                                                                                                                                    |                                                                                                                                                   cGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTGCgg                                                                                                                                                 >  1:135984/1‑151 (MQ=14) cGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTGCgg                                                                                                                                                 >  1:607770/1‑151 (MQ=14)    cGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTGCgggg                                                                                                                                               <  1:434473/150‑1 (MQ=14)                 cGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTGCGGGGTGCACTGAACTGTa                                                                                                                                 <  2:53199/151‑1 (MQ=18)                     gCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGTCGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTGCGGGGTGCACTGAACTGTAGGTc                                                                                                                             >  1:289204/1‑151 (MQ=18)                      ccTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTGCGGGGTGCACTGAACTGt                                                                                                                                  >  1:584757/1‑145 (MQ=18)                       cTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTGCGGGGTGCACTGAACTGTAGGTCgg                                                                                                                           >  1:212831/1‑151 (MQ=18)                         gATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTGCGGGGTGCACTGAACTGTAGGTCGGa                                                                                                                          <  2:226048/150‑1 (MQ=18)                              gccgcTGACGTGACTTATCAGGCCTACGGGGAGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTGCGGGGTGCACTGAACTGTAGGTCGGATAAGa                                                                                                                     >  1:623557/1‑150 (MQ=18)                                            ttATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTGCGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCAt                                                                                                      >  2:549338/1‑151 (MQ=14)                                                cAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTGCGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCg                                                                                                   >  1:663827/1‑150 (MQ=14)                                                  ggCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTGCGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGAc                                                                                                 >  2:800105/1‑150 (MQ=14)                                                  ggCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTGCGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGAc                                                                                                 <  2:676424/150‑1 (MQ=14)                                                  ggCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTGCGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGAc                                                                                                 >  1:407905/1‑150 (MQ=14)                                                       aCGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTGCGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCg                                                                                            <  1:102867/150‑1 (MQ=18)                                                         ggggTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTGCGGGGTGCACTGAACTGTAGGTCGGGTAAGACGCATAGCGTCGCATCCGACAACCGATg                                                                                         >  2:322548/1‑151 (MQ=25)                                                               cACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTGCGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGAt                                                                                   <  1:896017/151‑1 (MQ=25)                                                                    aaCTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTGCGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCtgctgccac                                                                              >  2:267248/1‑148 (MQ=255)                                                                    aaCTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTGCGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCAATGCCTGATgccc                                                                               >  2:227207/1‑148 (MQ=255)                                                                     aCTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTGCGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGcgcc                                                                              >  1:853703/1‑150 (MQ=25)                                                                        gTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTGCGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCt                                                                           >  1:700159/1‑150 (MQ=25)                                                                          aGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTGCGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGa                                                                         <  2:846923/150‑1 (MQ=25)                                                                               ggATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTGCGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGAcgcg                                                                     >  1:705749/1‑149 (MQ=25)                                                                               ggATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTGCGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGAc                                                                   <  1:447415/151‑1 (MQ=25)                                                                                    aGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTGCGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATc                                                              >  1:121644/1‑151 (MQ=25)                                                                                      aCGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTGCGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTAt                                                               >  1:58524/1‑148 (MQ=25)                                                                                      aCGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTGCGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATc                                                              >  1:848484/1‑149 (MQ=25)                                                                                       cGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTGCGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAgg                                                           >  1:437486/1‑151 (MQ=25)                                                                                          aTAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTGCGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCt                                                        >  1:603074/1‑151 (MQ=18)                                                                                              cGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTGCGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACACGACTTATCAGGCCTACgg                                                    >  2:67911/1‑151 (MQ=14)                                                                                                 cGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTGCGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGt                                                 <  2:174764/151‑1 (MQ=14)                                                                                                  gCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTGCGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGtg                                                <  1:432079/151‑1 (MQ=14)                                                                                                    aTCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTGCGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGtga                                               >  1:746205/1‑150 (MQ=17)                                                                                                    aTCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTGCGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGtg                                                >  1:811552/1‑149 (MQ=14)                                                                                                       cGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTGCGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGtgaactg                                           <  1:322574/151‑1 (MQ=18)                                                                                                           aaCCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTGCGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTgaactgaact                                       >  1:347885/1‑151 (MQ=18)                                                                                                              cGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTGCGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGt                                     <  1:259019/150‑1 (MQ=18)                                                                                                                 tGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTGCGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAAACGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGtgaac                                             <  1:759842/139‑1 (MQ=21)                                                                                                                 tGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTGCGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAAACGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGtgaac                                             >  2:759792/1‑139 (MQ=21)                                                                                                                    cTGATGCGCCGCTGACGCGACTTATCAGGCCTGCGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCgg                              >  1:218952/1‑151 (MQ=18)                                                                                                                         gcgcCGCTGACGCGACTTATCAGGCCTGCGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATaa                          <  2:465339/150‑1 (MQ=18)                                                                                                                         gcgcCGCTGACGCGACTTATCAGGCCTGCGGGGTGCACTGAACTGTAGGTCGGATAAGACGCAAAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATa                           >  2:230314/1‑149 (MQ=18)                                                                                                                                   cgcgACTTATCAGGCCTGCGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCg               >  2:272984/1‑151 (MQ=18)                                                                                                                                   cgcgACTTATCAGGCCTGCGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCg               <  2:560788/151‑1 (MQ=18)                                                                                                                                   cgcgACTTATCAGGCCTGCGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGCGAATTGAAATGTAGGTCGGATAAGACGCATAGc                >  2:304647/1‑150 (MQ=255)                                                                                                                                     cgACTTATCAGGCCTGCGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTc             <  1:128718/151‑1 (MQ=17)                                                                                                                                      gACTTATCAGGCCTGCGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCg            >  2:4276/1‑151 (MQ=14)                                                                                                                                      gACTTATCAGGCCTGCGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCAGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTc             >  1:674314/1‑150 (MQ=14)                                                                                                                                             cAGGCCTGCGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCg      <  1:509180/150‑1 (MQ=14)                                                                                                                                              gggCCTGCGGGGTGCACTGCAGTGTAGGTCGGATACGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGac    <  2:540506/150‑1 (MQ=9)                                                                                                                                               ggCCTGCGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGaca   <  1:610657/151‑1 (MQ=14)                                                                                                                                               ggCCTGCGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGaca   <  1:623689/151‑1 (MQ=14)                                                                                                                                                 ccTGCGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGacaa  >  2:19062/1‑150 (MQ=14)                                                                                                                                                    |                                                                                                                                                   CGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAA  >  NC_000913/349678‑349971

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