Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	349,847	T→C	intergenic (+275/‑165)	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,847	0	T	C	100.0%	73.9 / NA	29	intergenic (+275/‑165)	prpB/prpC	2‑methylisocitrate lyase/2‑methylcitrate synthase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (15/14);  total (15/14)

TGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCT  >  NC_000913/349738‑349977                                                                                                              |                                                                                                                                   tGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCGCATCCGAc                                                                                                       >  2:799701/1‑139 (MQ=17) tGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCGCATCCGAc                                                                                                       >  2:599394/1‑139 (MQ=17) tGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCGCATCCGAc                                                                                                       <  1:901277/139‑1 (MQ=18)   cACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCGGCTGACGCGACTTATCAGGCCAACGGGGGGCACTGAACGGTAGGCCGGATAGGAAGCATAGCGTCGCATCCGACaa                                                                                                     >  2:147495/1‑139 (MQ=11)             tAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTg                                                                                           >  1:1037822/1‑139 (MQ=255)                  cGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATgcg                                                                                      <  2:1010386/139‑1 (MQ=255)                     aTAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGcgccg                                                                                   >  2:650990/1‑139 (MQ=38)                     aTAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGcgccg                                                                                   >  2:306275/1‑139 (MQ=38)                         gACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGa                                                                               <  1:599394/139‑1 (MQ=255)                               tAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCg                                                                                                               <  1:110909/101‑1 (MQ=18)                               tAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCg                                                                                                               >  2:110909/1‑101 (MQ=17)                                   gTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTAt                                                                     <  2:1086336/139‑1 (MQ=17)                                             aCAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACg                                                           <  1:306275/139‑1 (MQ=17)                                               aaCCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACggg                                                         <  2:298910/139‑1 (MQ=18)                                                    aTGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGtgaa                                                    >  1:353496/1‑139 (MQ=21)                                                    aTGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCACTGACGCGACTTATCAGGCCTACGGGGtgaa                                                    <  2:1037822/139‑1 (MQ=16)                                                           aTGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCg                                                                                                               >  2:59740/1‑73 (MQ=255)                                                           aTGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCg                                                                                                               <  1:59740/73‑1 (MQ=21)                                                                ccgcTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAgg                                        >  1:293825/1‑139 (MQ=21)                                                                  gcTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTc                                      <  2:280021/139‑1 (MQ=21)                                                                    tGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCgg                                    <  2:353496/139‑1 (MQ=25)                                                                      aCGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGAt                                  <  2:203962/139‑1 (MQ=25)                                                                                      cTACGGGGTGCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCgg                                    >  2:209643/1‑121 (MQ=17)                                                                                      cTACGGGGTGCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCgg                                    <  1:209643/121‑1 (MQ=17)                                                                                      cTACGGGGTGCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCg                  <  1:25265/139‑1 (MQ=18)                                                                                           gggTGCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATcc             >  2:233503/1‑139 (MQ=18)                                                                                              tGCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGac          >  1:324974/1‑139 (MQ=18)                                                                                                    gAACTGTAGGCCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCTTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAAt    >  2:843257/1‑139 (MQ=21)                                                                                                      aCTGTAGGCCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAAtct  >  1:492376/1‑139 (MQ=37)                                                                                                              |                                                                                                                                   TGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCT  >  NC_000913/349738‑349977

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