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%	117.5 / NA	37	*262* (TAG→TAA)	frlD	fructoselysine 6‑kinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (18/19);  total (18/19)

GGAGATTCGTTCATTGCCGGATTCCTTTGCGGCTGGTCTGCGGGGATGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCTTTATTTCAGGGTCGATTATGTCAGCTACGGACCGCTACTCTCATCAACTC  >  NC_000913/3503833‑3504084                                                                                                                        |                                                                                                                                     ggAGATTCGTTCATTGCCGGATTCCTTTGCGGCTGGTCTGCGGGGATGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAAGTATAACGTTGGCGTGAGc                                                                                                                   <  1:141731/139‑1 (MQ=255) ggAGATTCGTTCATTGCCGGATTCCTTTGCGGCTGGTCTGCGGGGATGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAAGTATAACGTTGGCGTGAGc                                                                                                                   >  2:529514/1‑139 (MQ=255)              ttGCCGGATTCCTTTGCGGCTGGTCTGCGGGGATGACATTACCGTAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAAGTATAACGTTGGCGTGAGCATCTTCACGCCaa                                                                                                      <  1:810518/139‑1 (MQ=255)               tGCCGGATTCCTTTGCGGCTGGTCTGCGGGGATGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAAGTATAACGTTGGCGTGAGCATCTTCACGCCAAc                                                                                                     <  1:392368/139‑1 (MQ=255)               tGCCGGATTCCTTTGCGGCTGGTCTGCGGGGATGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAAGTATAACGTTGGCGTGAGCATCTTCACGCCAAc                                                                                                     >  1:512646/1‑139 (MQ=255)               tGCCGGATTCCTTTGCGGCTGGTCTGCGGGGATGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAAGTATAACGTTGGCGTGAGCATCTTCACGCCAAc                                                                                                     >  2:1003013/1‑139 (MQ=255)                         cTTTGCGGCTGGTCTGCGGGGATGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAAGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTAc                                                                                           >  1:188496/1‑139 (MQ=255)                               ggCTGGTCTGCGGGGATGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAAGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCg                                                                                     <  2:512646/139‑1 (MQ=255)                                    gTCTGCGGGGATGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTggtaagta                                                                                                                                   <  2:80079/88‑1 (MQ=255)                                    gTCTGCGGGGATGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTggtaagta                                                                                                                                   >  1:80079/1‑88 (MQ=255)                                    gTCTGCGGGGATGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAAGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGaaaa                                                                                <  2:761275/139‑1 (MQ=255)                                    gTCTGCGGGGATGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAAGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGaaaa                                                                                <  2:1442235/139‑1 (MQ=255)                                      cTGCGGGGATGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAAGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACg                                                                              <  2:188496/139‑1 (MQ=255)                                                gACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAAGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATaac                                                                    <  2:922008/139‑2 (MQ=255)                                                      aCCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAAGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTg                                                                                               <  1:1390570/106‑1 (MQ=255)                                                      aCCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAAGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTg                                                                                               >  2:1390570/1‑106 (MQ=255)                                                                 tAGCGCAGGGAACGGCGTGCGCGGCGAAAACTATTCAGTACCATGGTGCCTGGTAAGTATACCGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCtt                                                   >  2:191568/1‑139 (MQ=255)                                                                  aGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAAGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCttt                                                  >  2:529100/1‑139 (MQ=255)                                                                       aGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAAGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCTTTATTTc                                             >  2:232080/1‑139 (MQ=255)                                                                         ggAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAAGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCt                                                                         <  2:56683/109‑1 (MQ=255)                                                                         ggAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAAGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCt                                                                         >  1:56683/1‑109 (MQ=255)                                                                                   gcgcGGCGAAAACCATTCAGTACCACGGTGCCTGGTAAGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGa                                                                <  1:958077/108‑1 (MQ=255)                                                                                   gcgcGGCGAAAACCATTCAGTACCACGGTGCCTGGTAAGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGa                                                                >  2:958077/1‑108 (MQ=255)                                                                                    cgcgGCGAAAACCATTCAGTACCACGGTGCCTGGTAAGTATAACGTTGGCGTGAGCATCTTCACGCCaa                                                                                                      >  2:380572/1‑69 (MQ=255)                                                                                    cgcgGCGAAAACCATTCAGTACCACGGTGCCTGGTAAGTATAACGTTGGCGTGAGCATCTTCACGCCaa                                                                                                      <  1:380572/69‑1 (MQ=255)                                                                                             aaCCATTCAGTACCACGGTGCCTGGTAAGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCTTTATTTCAGGGTCGATTATGTCAGCTACg                       <  1:529514/139‑1 (MQ=255)                                                                                                   tCAGTACCACGGTGCCTGGTAAGTATAACGTTGGCGTGAGc                                                                                                                   <  1:1315776/41‑1 (MQ=255)                                                                                                   tCAGTACCACGGTGCCTGGTAAGTATAACGTTGGCGTGAGc                                                                                                                   >  2:1315776/1‑41 (MQ=255)                                                                                                     aGTACCACGGTGCCTGGTAAGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAAt                                                                    >  2:142395/1‑86 (MQ=255)                                                                                                     aGTACCACGGTGCCTGGTAAGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAAt                                                                    <  1:142395/86‑1 (MQ=255)                                                                                                          cACGGTGCCTGGTAAGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCTTTATTTCAGGGTCGATTATGTCAGCTACGGACCGCTACtctc          >  2:1464489/1‑139 (MQ=255)                                                                                                            cGGTGCCTGGTAAGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGa                                                                             <  2:1421831/70‑1 (MQ=255)                                                                                                            cGGTGCCTGGTAAGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGa                                                                             >  1:1421831/1‑70 (MQ=255)                                                                                                               tGCCTGGTAAGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCTTTATTTc                                             >  1:620545/1‑99 (MQ=255)                                                                                                               tGCCTGGTAAGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCTTTATTTc                                             <  2:620545/99‑1 (MQ=255)                                                                                                               tGCCTGGTAAGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCTTTATTTCAGGGTCGATTATGTCAGCTACGGACCGCTACTCTCATCaa     <  2:447838/139‑1 (MQ=255)                                                                                                               tGCCTGGTAAGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCTTTATTTCAGGGTCGATTATGTCAGCTACGGACCGCTACTCTCATCaa     >  1:447838/1‑139 (MQ=255)                                                                                                                 ccTGGTAAGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCTTTATTTCAGGGTCGATTATGTCAGCTACGGACCGCTACTCTCATCAAct   >  2:375760/1‑139 (MQ=255)                                                                                                                  cTGGTAAGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCTTTATTTCAGGGTCGATTATGTCAGCTACGGACCGCTACTCTCATCAActc  <  1:375760/139‑1 (MQ=255)                                                                                                                    ggtaagtaTAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCTTTATTTCAGGGTCGATTATGTCAGCTACg                       >  2:742357/1‑116 (MQ=255)                                                                                                                    ggtaagtaTAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCTTTATTTCAGGGTCGATTATGTCAGCTACg                       <  1:742357/116‑1 (MQ=255)                                                                                                                        |                                                                                                                                     GGAGATTCGTTCATTGCCGGATTCCTTTGCGGCTGGTCTGCGGGGATGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCTTTATTTCAGGGTCGATTATGTCAGCTACGGACCGCTACTCTCATCAACTC  >  NC_000913/3503833‑3504084

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