Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,503,952	G→A	*262* (TAG→TAA)	frlD →	fructoselysine 6‑kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,503,952	0	G	A	100.0%	104.7 / NA	33	*262* (TAG→TAA)	frlD	fructoselysine 6‑kinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/20);  total (13/20)

CGACACCATGGGTGCCGGAGATTCGTTCATTGCCGGATTCCTTTGCGGCTGGTCTGCGGGGATGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCTTTATTTCAGGGTCGATTATGTCAGCTACGGACCGCTACTCTCATCAACT  >  NC_000913/3503817‑3504083                                                                                                                                        |                                                                                                                                    cGACACCATGGGTGCCGGAGATTCGTTCATTGCCGGATTCCTTTGCGGCTGGTCTGCGGGGATGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTggtaagta                                                                                                                                  <  2:137223/139‑1 (MQ=255)         tGGGTGCCGGAGATTCGTTCATTGCCGGATTCCTTTGCGGCTGGTCTGCGGGGATGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAAGTATAACGTTg                                                                                                                          <  1:110439/139‑1 (MQ=255)            gTGCCGGAGATTCGTTCATTGCCGGATTCCTTTGCGGCTGGTCTGCGGGGATGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAAGTATAACGTTGGCg                                                                                                                       <  2:626481/139‑1 (MQ=255)               ccGGAGATTCGTTCATTGCCGGATTCCTTTGCGGCTGGTCTGCGGGGATGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAAGTATAACGTTGGCGTGa                                                                                                                    <  1:454601/139‑1 (MQ=255)                 ggAGATTCGTTCATTGCCGGATTCCTTTGCGGCTGGTCTGCGGGGATGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAAGTATAACGTTGGCGTGAGc                                                                                                                  >  2:473490/1‑139 (MQ=255)                          ttCATTGCCGGATTCCTTTGCGGCTGGTCTGCGGGGATGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAAGTATAACGTTGGCGTGAGCATCTTCACg                                                                                                         <  1:473490/139‑1 (MQ=255)                            cATTGCCGGATTCCTTTGCGGCTGGTCTGCGGGGATGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAAGTATAACGTTGGCGTGAGCATCTTCACGcc                                                                                                       <  1:707213/139‑1 (MQ=255)                               tGCCGGATTCCTTTGCGGCTGGTCTGCGGGGATGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAAGTATAACGTTGGCGTGAGCATCTTCACGCCAAc                                                                                                    <  1:325637/139‑1 (MQ=255)                                   ggATTCCTTTGCGGCTGGTCTGCGGGGATGACATTACCGCAGGCGATAGCGCTGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAAGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGc                                                                                                <  1:584545/139‑1 (MQ=255)                                     aTTCCTTTGCGGCTGGTCTGCGGGGATGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAAGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTg                                                                                              >  1:419223/1‑139 (MQ=255)                                           ttGCGGCTGGTCTGCGGGGATGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAAGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACtt                                                                                        >  1:754832/1‑139 (MQ=255)                                                gCTGGTCTGCGGGGATGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAAGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCgg                                                                                   >  2:43955/1‑139 (MQ=255)                                                      cTGCGGGGATGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAAGTATAATGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGaaaa                                                                               <  1:298507/137‑1 (MQ=255)                                                      cTGCGGGGATGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAAGTATAATGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGaaaa                                                                               >  2:298507/1‑137 (MQ=255)                                                      cTGCGGGGATGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAAGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACg                                                                             <  2:419223/139‑1 (MQ=255)                                                                gACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGCACCACGGTGCCTGGTAAGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAAt                                                                   <  1:409428/139‑1 (MQ=255)                                                                                 tAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAAGTATAACGTTGGCGTGAGCATCtt                                                                                                             <  1:48975/80‑1 (MQ=255)                                                                                 tAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAAGTATAACGTTGGCGTGAGCATCtt                                                                                                             >  2:48975/1‑80 (MQ=255)                                                                                            aCGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAAGTATAACGTTGGCGTGAGCATCTTCACGCCAACg                                                                                                   <  2:42526/79‑1 (MQ=255)                                                                                            aCGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAAGTATAACGTTGGCGTGAGCATCTTCACGCCAACg                                                                                                   >  1:42526/1‑79 (MQ=255)                                                                                            aCGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAAGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCTTTATTTCAGGGt                                       <  1:43955/139‑1 (MQ=255)                                                                                                        gcgAAAACCATTCAGTACCACGGTGCCTGGTAAGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACt                                                                                         >  1:351951/1‑77 (MQ=255)                                                                                                        gcgAAAACCATTCAGTACCACGGTGCCTGGTAAGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACt                                                                                         <  2:351951/77‑1 (MQ=255)                                                                                                        gcgAAAACCATTCAGTACCACGGTGCCTGGTAAGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCTTTATTTCAGGGTCGATTATGTCAg                           >  1:58559/1‑139 (MQ=255)                                                                                                                  ttCAGTACCACGGTGCCTGGTAAGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCTTTATTTCAGGGTCGATTATGTCAGCTACGGACCg                 <  2:58559/139‑1 (MQ=255)                                                                                                                     aGTACCACGGTGCCTGGTAAGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAattca                                                       <  1:502620/98‑1 (MQ=255)                                                                                                                     aGTACCACGGTGCCTGGTAAGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAattca                                                       >  2:502620/1‑98 (MQ=255)                                                                                                                     aGTACCACGGTGCCTGGTAAGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCTTTATTTCAGGGTCGATTATGTCAGCTACGGACCGCTa              <  1:203932/139‑1 (MQ=255)                                                                                                                       tACCACGGTGCCTGGTAAGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCTTTATTTCAGGGTCGATTATGTCAGCTACGGACCGCTAct            <  2:662948/139‑1 (MQ=255)                                                                                                                               tGCCTGGTAAGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCTTTATTTCAGGGTCGATTATGTCAGCTACGGa                    >  2:86125/1‑123 (MQ=255)                                                                                                                               tGCCTGGTAAGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCTTTATTTCAGGGTCGATTATGTCAGCTACGGa                    <  1:86125/123‑1 (MQ=255)                                                                                                                                 ccTGGTAAGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCTTTATTTCAGGGTCGATTATGTCAGCTa                        >  2:133661/1‑117 (MQ=255)                                                                                                                                 ccTGGTAAGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCTTTATTTCAGGGTCGATTATGTCAGCTa                        <  1:133661/117‑1 (MQ=255)                                                                                                                                 ccTGGTAAGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCTTTATTTCAGGGTCGATTATGTCAGCTACGGACCGCTACTCTCATCAAct  >  2:112925/1‑139 (MQ=255)                                                                                                                                        |                                                                                                                                    CGACACCATGGGTGCCGGAGATTCGTTCATTGCCGGATTCCTTTGCGGCTGGTCTGCGGGGATGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCTTTATTTCAGGGTCGATTATGTCAGCTACGGACCGCTACTCTCATCAACT  >  NC_000913/3503817‑3504083

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