Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,905,639	(A)9→8	intergenic (‑42/+92)	bglF ← / ← bglG	beta‑glucoside specific PTS enzyme II/BglG kinase/BglG phosphatase/transcriptional antiterminator BglG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,905,631	0	A	.	100.0%	198.9 / NA	47	intergenic (‑34/+100)	bglF/bglG	beta‑glucoside specific PTS enzyme II/BglG kinase/BglG phosphatase/transcriptional antiterminator BglG
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (17/30);  total (17/30)

TTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCACGCGC  >  NC_000913/3905501‑3905752                                                                                                                                   |                                                                                                                          ttaatttaaAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCAGTCATAACGTGCCCTCGACCGCTTTTCGGCGAAACTCCaaa                                                                                                                         >  2:132678/1‑133 (MQ=255)         aaaCGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCaaaaaaaacct                                                                                                                 >  1:50471/1‑130 (MQ=255)                aTCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGaaaaaa                                                                                                         <  2:260920/133‑1 (MQ=255)                     gTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGc                                                                                                    <  1:512686/133‑1 (MQ=255)                      ttGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGcc                                                                                                   >  2:370753/1‑133 (MQ=255)                      ttGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGc                                                                                                    <  2:84608/132‑1 (MQ=255)                      ttGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGc                                                                                                    >  1:84608/1‑132 (MQ=255)                       tGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCt                                                                                                  <  2:285317/133‑1 (MQ=255)                          gcAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGAc                                                                                               <  1:376157/133‑1 (MQ=255)                                   cAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCaa                                                                                      >  2:29577/1‑133 (MQ=255)                                        tCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑‑AAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAAtt                                                                                >  2:471746/1‑133 (MQ=255)                                          acaATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAAttt                                                                               >  1:500225/1‑133 (MQ=255)                                            aaTGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAAttttt                                                                             >  1:361973/1‑133 (MQ=255)                                                      gcgcCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGtttt                                                                   >  2:42644/1‑133 (MQ=255)                                                        gcCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGc                                                                 >  2:313734/1‑133 (MQ=255)                                                        gcCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGc                                                                 >  1:440556/1‑133 (MQ=255)                                                          ccccGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGccc                                                               >  2:361274/1‑133 (MQ=255)                                                             cGACTCCTGCTACTATTTTTCTGGCTAACTCCTTCATAACTCGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCt                                                            <  2:207549/133‑1 (MQ=255)                                                             cGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCt                                                            <  1:311912/133‑1 (MQ=255)                                                             cGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCt                                                            <  1:370753/133‑1 (MQ=255)                                                                      cGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAAATTTTTTCAGGTTTTGCCCGCTTAGTGCGGt                                                   >  1:207549/1‑133 (MQ=255)                                                                        aCTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTaa                                                 <  2:308928/133‑1 (MQ=255)                                                                          tATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACa                                               <  1:361274/133‑1 (MQ=255)                                                                                  tGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTa                                       <  2:50471/133‑1 (MQ=255)                                                                                  tGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTa                                       <  2:359175/133‑1 (MQ=255)                                                                                    gCTAACTCCTGCATATCTTGCCCTATCCCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACt                                     <  2:237139/133‑1 (MQ=255)                                                                                    gCTAACTCCGTCATATCTTGCCCTCTACCGATTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACt                                     <  1:132678/133‑1 (MQ=255)                                                                                    gCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACt                                     <  1:29577/133‑1 (MQ=255)                                                                                     cTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTc                                    <  2:254386/133‑1 (MQ=255)                                                                                       aaCTCCGTCATAACTTGCCCTCTACAGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAg                                  >  1:60783/1‑133 (MQ=255)                                                                                        aCTCCGTCATTACTTTCCCTCTACGGCTGTGCTGCAAAACTTC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGt                                 <  2:175632/133‑1 (MQ=255)                                                                                          tCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGtaa                               <  2:25396/133‑1 (MQ=255)                                                                                              tCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCTATTACTCTGTaataat                           >  1:175632/1‑133 (MQ=255)                                                                                                   aCTTGCCCTCTCCCGGTTTGCGGCGAATTTCC‑TAAAAAAACATGAAAAAACCGGTCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTc                      <  2:175266/133‑1 (MQ=255)                                                                                                    cTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGtaa                               >  2:89546/1‑123 (MQ=255)                                                                                                    cTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGtaa                               <  1:89546/123‑1 (MQ=255)                                                                                                       gCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTc                      >  2:31360/1‑129 (MQ=255)                                                                                                       gCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTc                      <  1:31360/129‑1 (MQ=255)                                                                                                                  cTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCa       <  2:60783/133‑1 (MQ=255)                                                                                                                  cTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCa       >  1:107557/1‑133 (MQ=255)                                                                                                                  cTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCa       >  2:278826/1‑133 (MQ=255)                                                                                                                   tttGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCAc      <  1:278826/133‑1 (MQ=255)                                                                                                                   tttGCGGCAAAACTCC‑‑AAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCAcg     <  1:471746/133‑1 (MQ=255)                                                                                                                    ttGCGGCAAAACTCC‑CAACAAATCCTGAAACAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCAcg     <  1:248168/133‑1 (MQ=255)                                                                                                                    ttGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCAcg     <  2:107557/133‑1 (MQ=255)                                                                                                                    ttGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCAcg     <  2:500225/133‑1 (MQ=255)                                                                                                                      gCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCAcgcg   <  1:412671/133‑1 (MQ=255)                                                                                                                      gCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCAcgcg   <  2:306926/133‑1 (MQ=255)                                                                                                                      gCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCAcgcg   <  1:269504/133‑1 (MQ=255)                                                                                                                      gCGGCAAAACACC‑AACAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCAcgcg   <  2:218543/133‑1 (MQ=255)                                                                                                                       cGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCACgcgc  <  2:302410/133‑1 (MQ=255)                                                                                                                                   |                                                                                                                          TTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCACGCGC  >  NC_000913/3905501‑3905752

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