Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	352,311	A→G	Y366C (TAC→TGC)	prpD →	2‑methylcitrate dehydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	352,311	0	A	G	100.0%	105.9 / NA	30	Y366C (TAC→TGC)	prpD	2‑methylcitrate dehydratase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (17/13);  total (17/13)

CACGAAGCCTGTATTCGCATCATCGACAAAAAAGGGCCGCTCAATAACCCGGCAGACCGCGATCACTGCATTCAGTACATGGTGGCGATCCCGCTGCTATTCGGGCGCTTAACGGCGGCAGATTACGAGGACAACGTTGCGCAAGATAAACGCATTGACGCCCTGCGCGAGAAGATCAATTGCTTTGAAGATCCGGCATTTACCGCTGACTACCACGACCCGGAAAAACGCGCCATCGCCAATGCCATTA  >  NC_000913/352187‑352436                                                                                                                             |                                                                                                                              cACGAAGCCTGTATTCGCATCATCGACAAAAAAGGGCCGCTCAATAACCCGGCAGACCGCGATCACTGCATTCAGTACATGGTGGCGATCCCGCTGCTATTCGGGCGCTTAACGGCGGCAGATTGCGAGGACAACGTTg                                                                                                                 <  1:124133/139‑1 (MQ=255)    gAAGCCTGTATTCGCATCATCGACAAAAAAGGGCCGCTCAATAACCCGGCAGACCGCGATCACTGCATTCAGTACATGGTGGCGATCCCGCTGCTATTCGGGCGCTTAACGGCGGCAGATTGCGAGGACAACGTTgcgc                                                                                                              >  1:551577/1‑139 (MQ=255)          tGTATTCGCATCATCGACAAAAAAGGGCCGCTCAATAACCCGGCAGACCGCGATCACTGCATTCAGTACATGGTGGCGATCCCGCTGCTATTCGGGCGCTTAACGGCGGCAGATTGCGAGGACAACGTTGCGCAAGATa                                                                                                        <  1:577673/139‑1 (MQ=255)                  catcatCGACAAAAAAGGGCCGCTCAATAACCCGGCAGACCGCGATCACTGCATTCAGTACATGGTGGCGATCCCGCTGCTATTCGGGCGCTTAACGGCGGCAGATTGCGAGGACAACGTTGCGCAAGATAAACGCAtt                                                                                                >  1:525763/1‑139 (MQ=255)                    tcatcGACAAAAAAGGGCCGCTCAATAACCCGGCAGACCGCGATCACTGCATTCAGTACATGGTGGCGATCCCGCTGCTATTCGGGCGCTTAACGGCGGCAGATTGCGAGGACAACGTTGCGCAAGATAAACGCATTGa                                                                                              >  1:744428/1‑139 (MQ=255)                      atcGACAAAAAAGGGCCGCTCAATAACCCGGCAGACCGCGATCACTGCATTCAGTACATGGTGGCGATCCCGCTGCTATTCGGGCGCTTAACGGCGGCAGATTGCGAGGACAACGTTGCGCAAGATAAACGCATTGACg                                                                                            >  2:204766/1‑139 (MQ=255)                        cGACAAAAAAGGGCCGCTCAATAACCCGGCAGACCGCGATCACTGCATTCAGTACATGGTGGCGATCCCGCTGCTATTCGGGCGCTTAACGGCGGCAGATTGCGAGGACAACGTTGCGCAAGATAAACGCATTGACGcc                                                                                          >  1:674021/1‑139 (MQ=255)                         gACAAAAAAGGGCCGCTCAATAACCCGGCAGACCGCGATCACTGCATTCAGTACATGGTGGCGATCCCGCTGCTATTCGGGCGCTTAACGGCGGCAGATTGCGAGGACAACGTTGCGCAAGATAAACGCATTGACGccc                                                                                         <  2:536703/139‑1 (MQ=255)                                 aGGGCCGCTCAATAACCCGGCAGACCGCGATCACTGCATTCAGTACATGGTGGCGATCCCGCTGCTATTCGGGCGCTTAACGGCGGCAGATTGCGAGGACAACGTTGCGCAAGATAAACGCATTGACGCCCTGCGCgag                                                                                 >  1:21214/1‑139 (MQ=255)                                 aGGGCCGCTCAATAACCCGGCAGACCGCGATCACTGCATTCAGTACATGGTGGCGATCCCGCAGCTATTCGGGCGCTTAACGGCGGCAGATTGCGAGGACAACGTTGCGCAAGATAAACGCATTGACGCCCTGCGCgag                                                                                 <  2:551577/139‑1 (MQ=255)                                  gggCCGCTCAATAACCCGGCAGACCGCGATCACTGCATTCAGTACATGGTGGCGATCCCGCTGCTATTCGGGCGCTTAACGGCGGCAGATTGCGAGGACAACGTTGCGCAAGATAAACGCATTGACGCCCTGCGCgaga                                                                                <  2:21214/139‑1 (MQ=255)                                         tCAATAACCCGGCAGACCGCGATCACTGCATTCAGTACATGGTGGCGATCCCGCTGCTATTCGGGCGCTTAACGGCGGCAGATTGCGAGGACAACGTTGCGCAAGATAAACGCATTGACGCCCTGCGCGAGAAGATCaa                                                                         <  2:744428/139‑1 (MQ=255)                                                        aCCGCGATCACTGCATTCAGTACATGGTGGCGATCCCGCTGCTATTCGGGCGCTTAACGGCGGCAGATTGCGAGGACAACGTTGCGCAAGATAAACGCATTGACGCCCTGCGCGAGAAGATCAATTGCTTTGAAGATcc                                                          <  2:674021/139‑1 (MQ=255)                                                          cgcgATCACTGCATTCAGTACATGGTGGCGATCCCGCTGCTATTCGGGCGCTTAACGGCGGCAGATTGCGAGGACAACGTTGCGCAAGATAAACGCATTGACGCCCTGCGCGAGAAGATCAATTGCTTTGAAGATCCgg                                                        >  1:376248/1‑139 (MQ=255)                                                                    gCATTCAGTACATGGTGGCGATCCCGCTGCTATTCGGGCGCTTAACGGCGGCAGATTGCGAGGACAACGTTGCGCAAGATAAACGCATTGACGCCCTGCGCGAGAAGATCAATTGCTTTGAAGATCCGGCATTTACCGc                                              >  1:543628/1‑139 (MQ=255)                                                                           gTACATGGTGGCGATCCCGCTGCTATTCGGGCGCTTAACGGCGGCAGATTGCGAGGACAACGTTGCGCAAGATAAACGCATTGACGCCCTGCGCGAGAAGATCAATTGCTTTGAAGATCCGGCATTTACCGCTGACTAc                                       >  2:300534/1‑139 (MQ=255)                                                                             aCATGGTGGCGATCCCGCTGCTATTCGGGCGCTTAACGGCGGCAGATTGCGAGGACAACGTTGCGCAAGATAAACGCATTGACGCCCTGCGCGAGAAGATCAATTGCTTTGAAGATCCGGCATTTACCGCTGACTACCa                                     <  2:376248/139‑1 (MQ=255)                                                                               aTGGTGGCGATCCCGCTGCTATTCGGGCGCTTAACGGCGGCAGATTGCGAGGACAACGTTGCGCAAGATAAACGCATTGACGCCCTGCGCGAGAAGATCAATTGCTTTGAAGATCCGGCATTTACCGCTGACTACCACg                                   >  2:370873/1‑139 (MQ=255)                                                                                      cGATCCCGCTGCTATTCGGGCGCTTAACGGCGGCAGATTGCGAGGACAACGTTGCGCAAGATAAACGCATTGACGCCCTGCGCGAGAAGATCAATTGCTTTGAAGATCCGGCATTTACCGCTGACTACCACGACCCGGa                            >  1:450750/1‑139 (MQ=255)                                                                                        aTCCCGCTGCTATTCGGGCGCTTAACGGCGGCAGATTGCGAGGACAACGTTGCGCAAGATAAACGCATTGACGCCCTGCGCGAGAAGATCAATTGCTTTGAAGATCCGGCATTTACCGCTGACTACCACGACCCGGaaa                          >  1:376757/1‑139 (MQ=255)                                                                                         tCCCGCTGCTATTCGGGCGCTTAACGGCGGCAGATTGCGAGGACAACGTTGCGCAAGATAAACGCATTGACGCCCTGCGCGAGAAGATCAATTGCTTTGAAGATCCg                                                         <  1:147462/107‑1 (MQ=255)                                                                                         tCCCGCTGCTATTCGGGCGCTTAACGGCGGCAGATTGCGAGGACAACGTTGCGCAAGATAAACGCATTGACGCCCTGCGCGAGAAGATCAATTGCTTTGAAGATCCg                                                         >  2:147462/1‑107 (MQ=255)                                                                                             gctgctATTCGGGCGCTTAACGGCGGCAGATTGCGAGGACAACGTTGCGCAAGATAAACGCATTGACGCCCTGCGCGAGAAGATCAATTGCTTTGAAGATCCGGCATTTACCg                                               >  1:347940/1‑113 (MQ=255)                                                                                             gctgctATTCGGGCGCTTAACGGCGGCAGATTGCGAGGACAACGTTGCGCAAGATAAACGCATTGACGCCCTGCGCGAGAAGATCAATTGCTTTGAAGATCCGGCATTTACCg                                               <  2:347940/113‑1 (MQ=255)                                                                                               tgctATTCGGGCGCTTAACGGCGGCAGATTGCGAGGACAACGTTGCGCAAGATAAACGCATTGACGCCCTGCGCGAGAAGATCAATTGCTTTGAAGATCCGGCATTTACCGCTGACTACCACGACCCGG‑AAAACGCGcc                  <  1:65071/139‑1 (MQ=255)                                                                                                 ctATTCGGGCGCTTAACGGCGGCAGATTGCGAGGACAACGTTGCGCAAGATAAACGCATTGACGCCCTGCGCGAGAAGATCAATTGCTTTGAAGATCCGGCATTTACCGCTGACTACCACGACCCGGAAAAACGCGCCa                 <  1:370873/139‑1 (MQ=255)                                                                                                           gcTTAACGGCGGCAGATTGCGAGGACAACGTTGCGCAAGATAAACGCATTGACGCCCTGCGCGAGAAGATCAATTGCTTTGAAGATCCGGCATTTACCGCTGACTACCACGACCCGGAAAAACGCGCCATCGCCAATGc       >  1:729640/1‑139 (MQ=255)                                                                                                            cTTAACGGCGGCAGATTGCGAGGACAACGTTGCGCAAGATAAACGCATTGACGCCCTGCGCGAGAAGATCAATTGCTTTGAAGATCCGGCATTTACCGCTGACTACCACGACCCGGAAAAACGCGCCATCGCCAATGcc      >  1:119585/1‑139 (MQ=255)                                                                                                              tAACGGCGGCAGATTGCGAGGACAACGTTGCGCAAGATAAACGCATTGACGCCCTGCGCGAGAAGATCAATTGCTTTGAAGATCCGGCATTTACCGCTGACTACCACGACCCGGAAAAACGCGCCATCGCCAATGCCAt    >  1:149958/1‑139 (MQ=255)                                                                                                                aCGGCGGCAGATTGCGAGGACAACGTTGCGCAAGATAAACGCATTGACGCCCTGCGCGAGAAGATCAATTGCTTTGAAGATCCGGCATTTACCGCTGACTACCACGACCCGGAAAAACGCGCCATCGCCAATGCCATTa  <  1:437867/139‑1 (MQ=255)                                                                                                                             |                                                                                                                              CACGAAGCCTGTATTCGCATCATCGACAAAAAAGGGCCGCTCAATAACCCGGCAGACCGCGATCACTGCATTCAGTACATGGTGGCGATCCCGCTGCTATTCGGGCGCTTAACGGCGGCAGATTACGAGGACAACGTTGCGCAAGATAAACGCATTGACGCCCTGCGCGAGAAGATCAATTGCTTTGAAGATCCGGCATTTACCGCTGACTACCACGACCCGGAAAAACGCGCCATCGCCAATGCCATTA  >  NC_000913/352187‑352436

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