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%	212.7 / NA	59	T74A (ACG→GCG)	aceK	isocitrate dehydrogenase kinase/phosphatase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (30/29);  total (30/29)

GCGCAGCAGCGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTACGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCC‑GGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTC  >  NC_000913/4218686‑4218947                                                                                                                                  |                                                                                                                                      gcgcAGCAGCGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCgg                                                                                                                              <  2:1377587/139‑1 (MQ=255)       cagcGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCAtttt                                                                                                                        >  2:1541288/1‑139 (MQ=255)                     ggCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTa                                                                                                            <  1:664754/137‑1 (MQ=255)                     ggCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTa                                                                                                            >  2:664754/1‑137 (MQ=255)                             ggCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTa                                                                                                  <  1:215395/139‑1 (MQ=255)                                                gATGAAAAACCGTATCCATCTTTACGTGCATAACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGATGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCC‑GGAt                                                                              <  1:1036850/139‑1 (MQ=255)                                                              tCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCC‑GGATTACCCGCGCTTCga                                                                >  2:966916/1‑139 (MQ=255)                                                                   cTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTg                                                                                     <  2:718080/114‑1 (MQ=255)                                                                   cTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTg                                                                                     >  1:718080/1‑114 (MQ=255)                                                                           tcatcaCGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTtact                                                                                                                    <  1:66143/75‑1 (MQ=255)                                                                           tcatcaCGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTtact                                                                                                                    >  2:66143/1‑75 (MQ=255)                                                                            catcaCGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTaaa                                                                                                          >  2:1575956/1‑84 (MQ=255)                                                                            catcaCGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTaaa                                                                                                          <  1:1575956/84‑1 (MQ=255)                                                                             atcaCGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTtact                                                                                                                    <  2:1646423/73‑1 (MQ=255)                                                                             atcaCGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTtact                                                                                                                    >  1:1646423/1‑73 (MQ=255)                                                                             atcaCGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACgtg                                                                                                               <  1:1372863/78‑1 (MQ=255)                                                                             atcaCGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACgtg                                                                                                               >  2:1372863/1‑78 (MQ=255)                                                                                  gTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCC‑GGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTa                                            <  1:966916/139‑1 (MQ=255)                                                                                     ggtctggtcGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCC‑GGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACt                                         <  2:199759/139‑1 (MQ=255)                                                                                      gtctggtcGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCC‑GGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACt                                         <  1:1127118/138‑1 (MQ=255)                                                                                      gtctggtcGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCC‑GGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACt                                         >  2:1127118/1‑138 (MQ=255)                                                                                       tctggtcGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCC‑gg                                                                                <  2:65402/98‑1 (MQ=255)                                                                                       tctggtcGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCC‑gg                                                                                >  1:65402/1‑98 (MQ=255)                                                                                       tctggtcGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCC‑GGATTAcc                                                                          <  1:878125/104‑1 (MQ=255)                                                                                       tctggtcGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCC‑GGATTAcc                                                                          >  2:878125/1‑104 (MQ=255)                                                                                         tggtcGTGGAGCAACTGCGCTGCATTACTAACGGACAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAagag                                                                                                       >  1:850561/1‑74 (MQ=255)                                                                                         tggtcGTGGAGCAACTGCGCTGCATTACTAACGGACAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAagag                                                                                                       <  2:850561/74‑1 (MQ=255)                                                                                           gtcGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCC‑GGATTACCCGCGCTTCGAGATTGCGgaga                                                     >  1:756455/1‑121 (MQ=255)                                                                                           gtcGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCC‑GGATTACCCGCGCTTCGAGATTGCGgaga                                                     <  2:756455/121‑1 (MQ=255)                                                                                             cGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCC‑GGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTAc                                 <  1:1541288/139‑1 (MQ=255)                                                                                                ggAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCgg                                                                                                                              >  1:1189201/1‑44 (MQ=255)                                                                                                ggAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCgg                                                                                                                              <  2:1189201/44‑1 (MQ=255)                                                                                                ggAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCC‑GGATTACCcg                                                                        >  1:1077060/1‑97 (MQ=255)                                                                                                ggAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCC‑GGATTACCcg                                                                        <  2:1077060/97‑1 (MQ=255)                                                                                                ggAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCC‑GGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGt                              >  2:1401593/1‑139 (MQ=255)                                                                                                      aCTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCgg                                                                                           >  2:1476992/1‑73 (MQ=255)                                                                                                      aCTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCgg                                                                                           <  1:1476992/73‑1 (MQ=255)                                                                                                         gcgcTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCC‑GGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTAttt                     >  1:537779/1‑139 (MQ=255)                                                                                                          cgcTGCATTACTAACGACCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCC‑GGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTg                    >  2:1361511/1‑139 (MQ=255)                                                                                                            cTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCC‑GGATTACCCGCGCTTCgaga                                                              >  2:145260/1‑95 (MQ=255)                                                                                                            cTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCC‑GGATTACCCGCGCTTCgaga                                                              <  1:145260/95‑1 (MQ=255)                                                                                                            cTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCC‑GGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGac                  >  1:480708/1‑139 (MQ=255)                                                                                                               cATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTTCTACGTGTTAAAGAGCATTACACCCGGCTGTTGCC‑GGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGaccac               >  1:1654984/1‑139 (MQ=255)                                                                                                                  tACTAACGACCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCC‑GGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACcgc            <  1:1361511/139‑1 (MQ=255)                                                                                                                    cTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCC‑GGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACcgctc          <  2:130982/139‑1 (MQ=255)                                                                                                                    cTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCC‑GGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACcgctc          >  1:130982/1‑139 (MQ=255)                                                                                                                       aCGGCCAAAGTGCGGACGCGGCATTTTTTCTACGTGTTAAAGAGCATTACACCCGGCTGTTGCC‑GGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACcgctcgct       <  2:1654984/139‑1 (MQ=255)                                                                                                                         ggCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACcg             >  1:629567/1‑132 (MQ=255)                                                                                                                         ggCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACcg             <  2:629567/132‑1 (MQ=255)                                                                                                                         ggCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCC‑GGATTACCCGCGCTTCg                                                                 >  1:1339618/1‑79 (MQ=255)                                                                                                                         ggCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCC‑GGATTACCCGCGCTTCg                                                                 >  2:1518182/1‑79 (MQ=255)                                                                                                                         ggCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCC‑GGATTACCCGCGCTTCg                                                                 <  1:1518182/79‑1 (MQ=255)                                                                                                                         ggCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCC‑GGATTACCCGCGCTTCg                                                                 <  2:1339618/79‑1 (MQ=255)                                                                                                                         ggCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCC‑GGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGaccac               >  2:1276808/1‑129 (MQ=255)                                                                                                                         ggCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCC‑GGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGaccac               <  1:1276808/129‑1 (MQ=255)                                                                                                                         ggCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCC‑GGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTa     >  1:1128449/1‑139 (MQ=255)                                                                                                                            cAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCC‑GGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTc  >  1:1380609/1‑139 (MQ=255)                                                                                                                              aaGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCC‑GGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTc                             >  2:698482/1‑110 (MQ=255)                                                                                                                              aaGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCC‑GGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTc                             <  1:698482/110‑1 (MQ=255)                                                                                                                                  |                                                                                                                                      GCGCAGCAGCGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTACGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCC‑GGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTC  >  NC_000913/4218686‑4218947

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