Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,218,815	A→G	T74A (ACG→GCG)	aceK →	isocitrate dehydrogenase kinase/phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,218,815	0	A	G	100.0%	237.3 / NA	66	T74A (ACG→GCG)	aceK	isocitrate dehydrogenase kinase/phosphatase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (34/32);  total (34/32)

CCGGTGCGCAGCAGCGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTACGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTA  >  NC_000913/4218681‑4218944                                                                                                                                       |                                                                                                                                  ccGGTGCGCAGCAGCGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGa                                                                                                                               >  1:1048912/1‑139 (MQ=255)                  ttCGAACAGGCCGACGGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTtactac                                                                                                              <  2:1014633/139‑1 (MQ=255)                   ttgAACAGGCCGACGGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGGGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACg                                                                                                             <  2:222279/137‑1 (MQ=255)                      aaCAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACgtgt                                                                                                          >  1:161396/1‑139 (MQ=255)                                             cagcagGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGTGGAcgcgac                                                                                                                         >  2:1255949/1‑99 (MQ=255)                                             cagcagGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGTGGAcgcgac                                                                                                                         <  1:1255949/101‑3 (MQ=255)                                             cagcagGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGt                                                                                   >  1:276537/1‑139 (MQ=255)                                              aggagGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGtt                                                                                  <  2:582802/136‑1 (MQ=255)                                              agcagGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGtt                                                                                  >  2:337909/1‑139 (MQ=255)                                                     gATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGAt                                                                           >  2:234432/1‑139 (MQ=255)                                                           aaaCCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCcg                                                                     <  2:276537/139‑1 (MQ=255)                                                                 tATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGc                                                                                      <  1:1260783/116‑1 (MQ=255)                                                                 tATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGc                                                                                      >  2:1260783/1‑116 (MQ=255)                                                                       tCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGAc                                                                                                                              <  2:281370/70‑1 (MQ=255)                                                                       tCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGAc                                                                                                                              >  1:281370/1‑70 (MQ=255)                                                                        cTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTAcc                                                                       <  2:790715/124‑1 (MQ=255)                                                                        cTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTAcc                                                                       >  1:790715/1‑124 (MQ=255)                                                                         tttACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGAcgcg                                                                                                                           <  1:12459/71‑1 (MQ=255)                                                                         tttACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGAcgcg                                                                                                                           >  2:12459/1‑71 (MQ=255)                                                                           tACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCgg                                                     >  2:849388/1‑139 (MQ=255)                                                                              gATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGgaga                                                  <  2:1048912/139‑1 (MQ=255)                                                                                  atcaCGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGAcgc                                                                                                                            >  2:703296/1‑61 (MQ=255)                                                                                  atcaCGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGAcgc                                                                                                                            <  1:703296/61‑1 (MQ=255)                                                                                   tcaCGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCttt                                             >  2:509433/1‑139 (MQ=255)                                                                                     aCGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCttttt                                           <  2:30224/139‑1 (MQ=255)                                                                                          ggtctggtcGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTtact                                                                                                                <  1:417357/65‑1 (MQ=255)                                                                                          ggtctggtcGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTtact                                                                                                                >  2:417357/1‑65 (MQ=255)                                                                                          ggtctggtcGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACt                                      >  2:251130/1‑139 (MQ=255)                                                                                               ggtcGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCgtg                                 <  1:509433/139‑1 (MQ=255)                                                                                                gtcGTGGAGCAACTGCGTTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACg                                                                                                             >  2:1059678/1‑62 (MQ=255)                                                                                                gtcGTGGAGCAACTGCGTTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACg                                                                                                             <  1:1059678/62‑1 (MQ=255)                                                                                                gtcGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCtt                                                                >  1:1329423/1‑107 (MQ=255)                                                                                                gtcGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCtt                                                                <  2:1329423/107‑1 (MQ=255)                                                                                                gtcGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCgtgt                                >  2:1242358/1‑139 (MQ=255)                                                                                                     ggAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAAc                                       <  2:1199918/127‑1 (MQ=255)                                                                                                     ggAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAAc                                       >  1:1199918/1‑127 (MQ=255)                                                                                                       aGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGGGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCg                         <  1:1129525/139‑1 (MQ=255)                                                                                                         cAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACTCCGCATTTTTACTACGCGTTAGAGAGCGTTACGCCTGGAAGTCGCCGGATTCCCCGCGGTTCGAAATTGCAGAGAGCTTTTTGAACTGCGTGTACTGTCGGt                       <  2:875862/139‑1 (MQ=255)                                                                                                          aaCTGCGCTGCATTACTAACGGCCAAAGTGCGGAcgc                                                                                                                            >  1:704633/1‑37 (MQ=255)                                                                                                          aaCTGCGCTGCATTACTAACGGCCAAAGTGCGGAcgc                                                                                                                            <  2:704633/37‑1 (MQ=255)                                                                                                          aaCTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACt                                      >  2:852640/1‑123 (MQ=255)                                                                                                          aaCTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACt                                      <  1:852640/123‑1 (MQ=255)                                                                                                           aCTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCcg                                                                     >  1:1283271/1‑91 (MQ=255)                                                                                                           aCTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCcg                                                                     <  2:1283271/91‑1 (MQ=255)                                                                                                                      ttACTAACGGCCAAAGTGCGGACGCGGCATTTTta                                                                                                                  <  2:1333621/35‑1 (MQ=255)                                                                                                                      ttACTAACGGCCAAAGTGCGGACGCGGCATTTTta                                                                                                                  >  1:1333621/1‑35 (MQ=255)                                                                                                                      ttACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCg                                                      <  2:1195446/95‑1 (MQ=255)                                                                                                                      ttACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCg                                                      >  1:1195446/1‑95 (MQ=255)                                                                                                                           aaCGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACcgctcgc     <  2:161396/139‑1 (MQ=255)                                                                                                                            aCGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCcg                                                                     <  1:1118929/74‑1 (MQ=255)                                                                                                                            aCGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCcg                                                                     >  2:1118929/1‑74 (MQ=255)                                                                                                                              ggCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGt                                                                                   <  2:1173437/58‑1 (MQ=255)                                                                                                                              ggCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGt                                                                                   >  1:1173437/1‑58 (MQ=255)                                                                                                                              ggCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCcg                                                                     >  1:1035358/1‑72 (MQ=255)                                                                                                                              ggCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCcg                                                                     <  2:1035358/72‑1 (MQ=255)                                                                                                                              ggCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCcg                                                                     <  1:98800/72‑1 (MQ=255)                                                                                                                              ggCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCcg                                                                     >  2:98800/1‑72 (MQ=255)                                                                                                                              ggCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTa  >  2:1062520/1‑139 (MQ=255)                                                                                                                              ggCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTa  >  1:541616/1‑139 (MQ=255)                                                                                                                              ggCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTa  >  1:1250077/1‑139 (MQ=255)                                                                                                                              ggCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGACTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACcgctc       <  1:996546/134‑1 (MQ=255)                                                                                                                              ggCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGACTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACcgctc       >  2:996546/1‑134 (MQ=255)                                                                                                                                 cAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGacca             <  2:451985/125‑1 (MQ=255)                                                                                                                                 cAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGacca             >  1:451985/1‑125 (MQ=255)                                                                                                                                   aaGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGt                       <  1:15245/113‑1 (MQ=255)                                                                                                                                   aaGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACACGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGt                       >  2:15245/1‑113 (MQ=255)                                                                                                                                       |                                                                                                                                  CCGGTGCGCAGCAGCGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTACGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTA  >  NC_000913/4218681‑4218944

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