Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,508,659	G→A	G256G (GGC→GGT)	glk ←	glucokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,508,659	0	G	A	100.0%	104.2 / NA	33	G256G (GGC→GGT)	glk	glucokinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (18/15);  total (18/15)

ACGATGAGATACACCGGAATATCATGGACATATTCTTTAAAGCGCCCTTTATCTTCAAATGCGGCACGGAAACCGGAGGCTTTGAAGAACTCAAGGAAGCGCGGCACGATACCGCCCGCAATAAACACGCCGCCAAATGTCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGGCCCATAATGACGCAAAACAGCGACAATGCGCGGCGGCAATCGGTGCAGCTGTCAGCCAGCGCGCGTTCGGTAATATCTTTTGGCTTGAG  >  NC_000913/2508528‑2508790                                                                                                                                    |                                                                                                                                    aCGATGAGATACACCGGAATATCATGGACATATTCTTTAAAGCGCCCTTTATCTTCAAATGCGGCACGGAAACCGGAGGCTTTGAAGAACTCAAGGAAGCGCGGCACGATACCGCCCGCAATAAACACGCCACCAAATg                                                                                                                              <  2:471991/139‑1 (MQ=255)     tGAGATACACCGGAATATCATGGACATATTCTTTAAAGCGCCCTTTATCTTCAAATGCGGCACGGAAACCGGAGGCTTTGAAGAACTCAAGGAAGCGCGGCACGATACCTCCCGCAATAAACACGCCACCAAATGTccc                                                                                                                          >  1:768473/1‑139 (MQ=255)     tGAGATACACCGGAATATCATGGACATATTCTTTAAAGCGCCCTTTATCTTCAAATGCGGCACGGAAACCGGAGGCTTTGAAGAACTCAAGGAAGCGCGGCACGATACCGCCCGCAATAAACACGCCACCAAATGTccc                                                                                                                          >  1:892300/1‑139 (MQ=255)               cGGAATATCATGGACATATTCTTTAAAGCGCCCTTTATCTTCAAATGCGGCACGGAAACCGGAGGCTTTGAAGAACTCAAGGAAGCGCGGCACGATACCGCCCGCAATAAACACGCCACCAAATGTCCCGAGATTGAgc                                                                                                                >  2:980881/1‑139 (MQ=255)                 gAATATCATGGACATATTCTTTAAAGCGCCCTTTATCTTCAAATGCGGCACGGAAACCGGAGGCTTTGAAGAACTCAAGGAAGCGCGGCACGATACCGCCCGCAATAAACACGCCACCAAATGTCCCGAGATTGAgcgc                                                                                                              >  2:708580/1‑139 (MQ=255)                       cATGGACATATTCTTTAAAGCGCCCTTTATCTTCAAATGCGGCACGGAAACCGGAGGCTTTGAAGAACTCAAGGAAGCGCGGCACGATACCGCCCGCAATAAACACGCCACCAAATGTCCCGAGATTGAGCGCCAGAtt                                                                                                        >  2:5526/1‑139 (MQ=255)                         tGGACATATTCTTTAAAGCGCCCTTTATCTTCAAATGCGGCACGGAAACCGGAGGCTTTGAAGAACTCAAGGAAGCGCGGCACGATACCGCCCGCAATAAACACGCCACCAAATGTCCCGAGATTGAGCGCCAGATTgc                                                                                                      >  2:636790/1‑139 (MQ=255)                                 ttCTTTAAAGCGCCCTTTATCTTCAAATGCGGCACGGAAACCGGAGGCTTTGAAGAACTCAAGGAAGCGCGGCACGATACCGCCCGCAATAAACACGCCACCAAATGTCCCGAGATTGAGCGCCAGATTGCCGCCaaaa                                                                                              <  1:143478/139‑1 (MQ=255)                                       aaaGCGCCCTTTATCTTCAAATGCGGCACGGAAACCGGAGGCTTTGAAGAACTCAAGGAAGCGCGGCACGATACCGCCCGCAATAAACACGCCACCaaa                                                                                                                                <  1:147993/99‑1 (MQ=255)                                       aaaGCGCCCTTTATCTTCAAATGCGGCACGGAAACCGGAGGCTTTGAAGAACTCAAGGAAGCGCGGCACGATACCGCCCGCAATAAACACGCCACCaaa                                                                                                                                >  2:147993/1‑99 (MQ=255)                                           cgcCCTTTATCTTCAAATGCGGCACGGAAACCGGAGGCTTTGAAGAACTCAAGGAAGCGCGGCACGATACCGCCCGCAATAAACACGCCACCAAATGTCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGGCCCATat                                                                                    >  2:397845/1‑138 (MQ=255)                                              ccTTTATCTTCAAATGCGGCACGGAAACCGGAGGCTTTGAAGAACTCAAGGAAGCGCGGCACGATACCGCCCGCAATAAACACGCCACCAAATGTCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGGCCCATAATGa                                                                                 >  1:60660/1‑139 (MQ=255)                                              ccTTTATCTTCAAATGCGGCACGGAAACCGGAGGCTTTGAAGAACTCAAGGAAGCGCGGCACGATACCGCCCGCAATAAACACGCCACCAAATGTCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGGCCCATAATGa                                                                                 >  1:621614/1‑139 (MQ=255)                                                        cAAATGCGGCACGGAAACCGGAGGCTTTGAAGAACTCAAGGAAGCGCGGCACGATACCGCCCGCAATAAACACGCCACCAAATGTCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGGCCCATATTGACGCAAAACAg                                                                       <  1:397845/139‑1 (MQ=255)                                                        cAAATGCGGCACGGAAACCGGAGGCTTTGAAGAACTCAAGGAAGCGCGGCACGATACCGCCCGCAATAAACACGCCACCAAATGTCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGGCCCATAATGACGCAAAACAg                                                                       >  2:874035/1‑139 (MQ=255)                                                        cAAATGCGGCACGGAAACCGGAGGCTTTGAAGAACTCAAGGAAGCGCGGCACGATACCGCCCGCAATAAACACGCCACCAAATGTCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGGCCCATAATGACGCAAAACAg                                                                       >  1:451383/1‑139 (MQ=255)                                                          aaTGCGGCACGGAAACCGGAGGCTTTGAAGAACTCAAGGAAGCGCGGCACGATACCGCCCGCAATAAACACGCCACCAAATGTCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGGCCCATAATGACGCAAAACAGCg                                                                     >  2:653930/1‑139 (MQ=255)                                                                           ggAGGCTTTGAAGAACTCAAGGAAGCGCGGCACGATACCGCCCGCAATAAACACGCCACCAAATGTCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGGCCCATAATGACGCAAAACAGCGACAATGCGCGGCGGCaa                                                    <  2:621614/139‑1 (MQ=255)                                                                                 tttGAAGAACTCAAGGAAGCGCGGCACGATACCGCCCGCAATAAACACGCCACCAAATGTCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGGCCCATAATGACGCAAAACAGCGACAATGCGCGGCGGCAATCGGTg                                              >  2:727880/1‑139 (MQ=255)                                                                                  ttGAAGAACTCAAGGAAGCGCGGCACGATACCGCCCGCAATAAACACGCCACCAAATGTCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGGCCCATAATGACGCAAAACAGCGACAATGCGCGGCGGCAATCGGTGc                                             <  1:5526/139‑1 (MQ=255)                                                                                        aaCTCAAGGAAGCGCGGCACGATACCGCCCGCAATAAACACGCCACCAAATGTCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGGCCCATAATGACGCAAAACAGCGACAATGCGCGGCGGCAATCGGTGCAGCTGt                                       <  1:653930/139‑1 (MQ=255)                                                                                                  aGCGCGGCACGATACCGCCCGCAATAAACACGCCACCAAATGTCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGGCCCATAATGACGCAAAACAGCGACAATGCGCGGCGGCAATCGGTGCAGCTGTCAGCCAgcgc                             >  2:1007395/1‑139 (MQ=255)                                                                                                     gcgGCACGATGCCTCCCGCAATAGACACGCTACCACATGGCCCGAAATTGAGCGCCAGATAGCCGCCAAAACGGCCCATAATGCCGCAAAACAGCGACAATGCGCGACGGCAATGGGTTCAGCTGTCAGCCAGcgcgcg                          <  2:768473/139‑1 (MQ=255)                                                                                                     gcgGCACGATACCGCCCGCAATAAACACGCCACCAAATGTCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGGCCCATAATGACGCAAAACAGCGACAATGCGCGGCGGCAATCGGTGCAGCTGTCAGCCAGcgcgcg                          <  1:1007395/139‑1 (MQ=255)                                                                                                      cgGCACGATACCGCCCGCAATAAACACGCCACCAAATGTCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGGCCCATAATGACGCAAAACAGCGACAATGCGCGGCGGCAATCGGTGCAGCTGTCAGCCAGCGCGCGt                         <  2:60660/139‑1 (MQ=255)                                                                                                       ggCACGATACCGCCCGCAATAAACACGCCACCAAATGTCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGGCCCATAATGACGCAAAACAGCGACAATGCGCGGCGGCAATCGGTGCAGCTGTCAGCCAGCGCGCGtt                        <  1:636790/139‑1 (MQ=255)                                                                                                       ggCACGATACCGCCCGCAATAAACACGCCACCAAATGTCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGGCCCATAATGACGCAAAACAGCGACAATGCGCGGCGGCAATCGGTGCAGCTGTCAGCCAGCGCGCGtt                        <  1:274662/139‑1 (MQ=255)                                                                                                        gCACGATACCGCCCGCAATAAACACGCCACCAAATGTCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGGCCCATAATGACGCAAAACAGCGACAATGCGCGGCGGCAAt                                                   >  2:1075832/1‑111 (MQ=255)                                                                                                        gCACGATACCGCCCGCAATAAACACGCCACCAAATGTCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGGCCCATAATGACGCAAAACAGCGACAATGCGCGGCGGCAAt                                                   <  1:1075832/111‑1 (MQ=255)                                                                                                                   cccgcAATAAACACGCCACCAAATGTCCCGAGATTGAGCGCCAGAtt                                                                                                        <  2:149712/47‑1 (MQ=255)                                                                                                                   cccgcAATAAACACGCCACCAAATGTCCCGAGATTGAGCGCCAGAtt                                                                                                        >  1:149712/1‑47 (MQ=255)                                                                                                                     cgcAATAAACACGCCACCAAATGTCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGGCCCATAATGACGCAAAACAGCGACAATGAGCGGCGGCAAACGGTGCAGCTGTCAGTCAGCGCGCGTTCGGTAATATCtttt          >  2:82134/1‑139 (MQ=255)                                                                                                                             acacGCCACCAAATGTCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGGCCCATAATGACGCAAAACAGCGACAATGCGCGGCGGCAATCGGTGCAGCTGTCAGCCAGCGCGCGTTCGGTAATATCTTTTGGCTTgag  <  1:727880/139‑1 (MQ=255)                                                                                                                                    |                                                                                                                                    ACGATGAGATACACCGGAATATCATGGACATATTCTTTAAAGCGCCCTTTATCTTCAAATGCGGCACGGAAACCGGAGGCTTTGAAGAACTCAAGGAAGCGCGGCACGATACCGCCCGCAATAAACACGCCGCCAAATGTCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGGCCCATAATGACGCAAAACAGCGACAATGCGCGGCGGCAATCGGTGCAGCTGTCAGCCAGCGCGCGTTCGGTAATATCTTTTGGCTTGAG  >  NC_000913/2508528‑2508790

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