Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,219,075	A→G	L160L (CTA→CTG)	aceK →	isocitrate dehydrogenase kinase/phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,219,075	0	A	G	100.0%	190.4 / NA	53	L160L (CTA→CTG)	aceK	isocitrate dehydrogenase kinase/phosphatase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (26/27);  total (26/27)

ACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCATTTACAGGTGGCGAACGAACTGTTTTAC  >  NC_000913/4218944‑4219201                                                                                                                                    |                                                                                                                               aCTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTGCCGctgcg                                                                                                                        <  1:291507/140‑1 (MQ=255) aCTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTGCCGctgcg                                                                                                                        <  2:690782/140‑1 (MQ=255)       gAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTGCCGctgcgcctgcg                                                                                                                  >  1:100096/1‑140 (MQ=255)        aGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTGCCGctgcgcctgcgc                                                                                                                 >  2:496993/1‑140 (MQ=255)         gCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTGCCGCTGCGCCTGCGCt                                                                                                                >  1:237884/1‑140 (MQ=255)         gCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTGCCGCTGCGCCTGCGCt                                                                                                                >  2:525133/1‑140 (MQ=255)            gCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTGCCGCTGCGCCTGCGCTGGc                                                                                                             <  1:149197/140‑1 (MQ=255)             cTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACTATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTGCCGCTGCGCCTGCGCTGGCa                                                                                                            >  1:732449/1‑140 (MQ=255)             cTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTGCCGCTGCGCCTGCGCTGGCa                                                                                                            >  1:161799/1‑140 (MQ=255)                        tAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTGCCGCTGCGCCTGCGCTGGCAGAATAAAAGcc                                                                                                 <  2:161799/140‑1 (MQ=255)                         aGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTGCCGCTGCGCCTGCGCTGGCAGAATAAAAGCCg                                                                                                >  1:479844/1‑140 (MQ=255)                          gCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTGCCGCTGCGCCTGCGCTGGCAGAACAAAAGCCGt                                                                                               <  1:496993/140‑1 (MQ=255)                           ctctCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTGCCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTg                                                                                              <  2:237884/140‑1 (MQ=255)                            tctcAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTGCCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGa                                                                                             <  1:777488/140‑1 (MQ=255)                             ctcAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTGCCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGAc                                                                                            <  2:732449/140‑1 (MQ=255)                             ctcAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTGCCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGAc                                                                                            <  2:296756/140‑1 (MQ=255)                                aGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTGCCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATc                                                                                         >  2:500937/1‑140 (MQ=255)                                  ccAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTGCCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCca                                                                                       <  2:747174/140‑1 (MQ=255)                                      agCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTGCCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCcattac                                                                                   >  2:369030/1‑140 (MQ=255)                                          ccgcTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTGCCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGAc                                                                                            >  2:9165/1‑127 (MQ=255)                                          ccgcTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTGCCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGAc                                                                                            <  1:9165/127‑1 (MQ=255)                                                tCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTGCCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGcc                                                                         >  2:246981/1‑140 (MQ=255)                                                 cGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTGCCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCa                                                                        <  2:742362/140‑1 (MQ=255)                                                          ccccGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTGCCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGa                                                               <  1:199151/140‑1 (MQ=255)                                                          ccccGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTGCCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGa                                                               >  2:199151/1‑140 (MQ=255)                                                               cccgcTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTGCCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGc                                                          >  1:787239/1‑140 (MQ=255)                                                               cccgcTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTGCCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGc                                                          <  1:246981/140‑1 (MQ=255)                                                                ccgcTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTGCCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCt                                                         <  1:647932/140‑1 (MQ=255)                                                                ccgcTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTGCCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCt                                                         >  1:104831/1‑140 (MQ=255)                                                                ccgcTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTGCCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCt                                                         <  2:479844/140‑1 (MQ=255)                                                                   cTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTGCCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTggg                                                      <  1:388770/140‑1 (MQ=255)                                                                    tGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTGCCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTgggg                                                     >  1:641386/1‑140 (MQ=255)                                                                    tGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTGCCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTgggg                                                     <  1:94131/140‑1 (MQ=255)                                                                     ggCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTGCCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGga                                                    <  1:369030/140‑1 (MQ=255)                                                                                     cccGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTGCCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGa                                    <  2:787239/140‑1 (MQ=255)                                                                                     cccGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTGCCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGa                                    <  2:641386/140‑1 (MQ=255)                                                                                       cGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTGCCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGaaa                                  <  1:99670/140‑1 (MQ=255)                                                                                       cGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTGCCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGaaa                                  >  1:198071/1‑140 (MQ=255)                                                                                        gATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTGCCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAg                                 <  1:525133/140‑1 (MQ=255)                                                                                                   ggAATCTCTACTGATGCGCGTTATCAGCGACCTGCCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCATTTACAgg                      >  1:581758/1‑140 (MQ=255)                                                                                                          ctACTGATGCGCGTTATCAGCGACCTGCCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCATTTACAGGTGGcgaa               >  1:452254/1‑140 (MQ=255)                                                                                                            aCTGATGCGCGTTATCAGCGACCTGCCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCATTTACAGGTGGcga                <  1:446441/137‑1 (MQ=255)                                                                                                            aCTGATGCGCGTTATCAGCGACCTGCCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCATTTACAGGTGGcga                >  2:446441/1‑137 (MQ=255)                                                                                                             cTGATGCGCGTTATCAGCGACCTGCCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCATTTACAGGTGGcgaacga            <  2:581758/140‑1 (MQ=255)                                                                                                               gATGCGCGTTATCAGCGACCTGCCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCATTTACAGGTGGCgaacgaac          >  1:159549/1‑140 (MQ=255)                                                                                                                  gcgcGTTATCAGCGACCTGCCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCATTTACAGGTGGCGAACGAACTGt       >  1:631153/1‑140 (MQ=255)                                                                                                                  gcgcGTTATCAGCGACCTGCCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCATTTACAGGTGGCGAACGAACTGt       >  2:625521/1‑140 (MQ=255)                                                                                                                   cgcgTTATCAGCGACCTGCCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCATTTACAGGTGGcga                <  2:218985/130‑1 (MQ=255)                                                                                                                   cgcgTTATCAGCGACCTGCCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCATTTACAGGTGGcga                >  1:218985/1‑130 (MQ=255)                                                                                                                   cgcgTTATCAGCGACCTGCCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCATTTACAGGTGGCGAACGAACTGtt      >  2:3211/1‑140 (MQ=255)                                                                                                                   cgcgTTATCAGCGACCTGCCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCATTTACAGGTGGCGAACGAACTGtt      >  2:268583/1‑140 (MQ=255)                                                                                                                     cgTTATCAGCGACCTGCCGCTGCGCCTGCGCTGGCAGAATAAAAGCCATGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCATTTACAGGTGGCGAACGAACTGtttt    >  2:33900/1‑140 (MQ=255)                                                                                                                       ttATCAGCGACCTGCCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCATTTACAGGTGGCGAACGAACTGTTTTAc  <  1:510797/140‑1 (MQ=255)                                                                                                                                    |                                                                                                                               ACTCCCGAGCGGCTTTTTATCTTTAGCTCTCAGCCAGAGCGCCGCTTTCGTACCATTCCCCGCCCGCTGGCGAAAGACTTTCACCCCGATCACGGCTGGGAATCTCTACTGATGCGCGTTATCAGCGACCTACCGCTGCGCCTGCGCTGGCAGAATAAAAGCCGTGACATCCATTACATTATTCGCCATCTGACGGAAACGCTGGGGACAGACAACCTCGCGGAAAGTCATTTACAGGTGGCGAACGAACTGTTTTAC  >  NC_000913/4218944‑4219201

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