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

CGCAGCAGCGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTACGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTA  >  NC_000913/4218687‑4218944                                                                                                                                 |                                                                                                                                  cgcAGCAGCGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATAATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGc                                                                                                                         >  1:1198514/1‑139 (MQ=255)     gcagcGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCAttt                                                                                                                     >  2:467015/1‑139 (MQ=255)         cGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAGAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATT‑‑‑‑ACGGCCAAAGTGCGGACGCGGCATTTTTACTACg                                                                                                             <  1:629925/139‑1 (MQ=255)           tttCGCACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAATAACCGTATCCGGCTTTACGATCATCCCGTTGGTCTGGTCGTATAGCAACTGCGCTGCATTACTAACGGTCAAAGTTCGGACCGGGCATTTTtaata                                                                                                               >  1:1133519/1‑139 (MQ=255)                  cAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTa                                                                                                        <  2:77896/139‑1 (MQ=255)                     gCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAag                                                                                                     <  2:1121489/139‑1 (MQ=255)                     gCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAag                                                                                                     >  1:1121489/1‑139 (MQ=255)                               aTGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTacac                                                                                           >  1:250322/1‑139 (MQ=255)                                tatatgagacagggCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGAGTTAAAGAGCATTACAcc                                                                                          <  2:803310/127‑1 (MQ=255)                                       cagcagGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGt                                                                                   >  2:854344/1‑139 (MQ=255)                                                       aCCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGAGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGAATTTTTACTACGTGTTAAAGAGCATTACACGCGGCTGTTGccctgtctcttata                                                                   >  1:803310/1‑127 (MQ=255)                                                               cATCTTTACGATAATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCgaga                                                           <  2:1198514/139‑1 (MQ=255)                                                                          tcatcaCGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGGGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGc                                                <  1:303389/139‑1 (MQ=255)                                                                           catcaCGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCt                                               >  2:660146/1‑139 (MQ=255)                                                                            atcaCGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACg                                                                                                             <  1:831549/76‑1 (MQ=255)                                                                            atcaCGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACg                                                                                                             >  2:831549/1‑76 (MQ=255)                                                                            atcaCGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACAcc                                                                                          >  2:267234/1‑95 (MQ=255)                                                                            atcaCGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACAcc                                                                                          <  1:267234/95‑1 (MQ=255)                                                                            atcaCGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCtt                                              >  1:671771/1‑139 (MQ=255)                                                                            atcaCGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCtt                                              >  2:1028274/1‑139 (MQ=255)                                                                            atcaCGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTAC‑‑GTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCtttt                                            >  1:484539/1‑139 (MQ=255)                                                                            atcaCGTTGGTCTGGTCGAGGATCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGTCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCtt                                              >  1:678215/1‑139 (MQ=255)                                                                             tcaCGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAag                                                                                                     >  1:606904/1‑83 (MQ=255)                                                                             tcaCGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAag                                                                                                     <  2:606904/83‑1 (MQ=255)                                                                             tcaCGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCttt                                             >  1:1206094/1‑139 (MQ=255)                                                                                 gTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTa                                         <  1:660146/139‑1 (MQ=255)                                                                                 gTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTa                                         <  1:854344/139‑1 (MQ=255)                                                                                    ggtctggtcGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCgg                                                                                                                          >  2:433696/1‑55 (MQ=255)                                                                                    ggtctggtcGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCgg                                                                                                                          <  1:433696/55‑1 (MQ=255)                                                                                      tctggtcGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACgtg                                                                                                           <  2:96136/68‑1 (MQ=255)                                                                                      tctggtcGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACgtg                                                                                                           >  1:96136/1‑68 (MQ=255)                                                                                         ggtcGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACTCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCgtg                                 <  2:1133519/139‑1 (MQ=255)                                                                                               ggAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACgtgt                                                                                                          >  1:1065224/1‑60 (MQ=255)                                                                                               ggAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACgtgt                                                                                                          <  2:1065224/60‑1 (MQ=255)                                                                                               ggAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACgtg                                                                                                           >  1:504241/1‑59 (MQ=255)                                                                                               ggAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACgtg                                                                                                           <  2:504241/59‑1 (MQ=255)                                                                                                gAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCa                                                                                                                        >  2:1089777/1‑45 (MQ=255)                                                                                                gAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCa                                                                                                                        <  1:1089777/45‑1 (MQ=255)                                                                                                gAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTc                          <  2:346996/139‑1 (MQ=255)                                                                                                gAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTc                          <  2:250322/139‑1 (MQ=255)                                                                                                    aaCTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCAGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGtt                      >  1:126349/1‑139 (MQ=255)                                                                                                     aCTGCGCTGCATTACTAACGGCCAAAGTGCGGAcgc                                                                                                                            >  2:716743/1‑36 (MQ=255)                                                                                                           cTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGac               <  1:1028274/139‑1 (MQ=255)                                                                                                              cATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCcg                                                                     >  2:428505/1‑82 (MQ=255)                                                                                                              cATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCcg                                                                     <  1:428505/82‑1 (MQ=255)                                                                                                                      aCGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACcgctcgct    <  1:1093051/139‑1 (MQ=255)                                                                                                                        ggCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTa  >  1:528070/1‑139 (MQ=255)                                                                                                                        ggCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTa  >  2:459365/1‑139 (MQ=255)                                                                                                                                 |                                                                                                                                  CGCAGCAGCGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTACGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTA  >  NC_000913/4218687‑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: ‑