Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,904,500	A→G	G366G (GGT→GGC)	bglF ←	beta‑glucoside specific PTS enzyme II/BglG kinase/BglG phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,904,500	0	A	G	100.0%	158.1 / NA	45	G366G (GGT→GGC)	bglF	beta‑glucoside specific PTS enzyme II/BglG kinase/BglG phosphatase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (27/18);  total (27/18)

TCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGATACCAAACAGACTCGTCAACGCCGCTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTCGCAGAGGAAGACGCCGAGCGCCGCCCCGACCTGCGCCATAATGGCGGGCATTAACAGCGGGATCATG  >  NC_000913/3904391‑3904626                                                                                                              |                                                                                                                               tCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGATGCCAAACAGACTCGTCAAcgccgc                                                                                                         >  1:51831/1‑133 (MQ=255) tCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGATGCCAAACAGACTCGTCAAcgccgc                                                                                                         >  2:32231/1‑133 (MQ=255)                         tggCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGATGCCAAACAGACTCGTCAACGCCGCTGATCCCGCCACCACTTTTTTCTg                                                                                 >  2:380943/1‑133 (MQ=255)                             ccccAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGATGCCAAACAGACTCGTCAACGCCGCTGATCCCGCCACCACTTTTTTCTGCGCa                                                                             >  2:334547/1‑133 (MQ=255)                              cccAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGATGCCAAACAGACTCGTCAACGCCGCTGATCCCGCCACCACTTTTTTCTGCGCAt                                                                            <  1:380943/133‑1 (MQ=255)                                          cTGATACAGGCGATAACAAAGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGATGCCAAACAGACTCGTCAACGCCGCTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTCGCaga                                                                <  1:120316/133‑1 (MQ=255)                                           tGATACAGGCGATAACAAAGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGATGCCAAACAGACTCGTCAACGCCGCTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTCGCagag                                                               >  1:491672/1‑133 (MQ=255)                                             aTACAGGCGATAACAAAGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGATGCCAAACAGACTCGTCAACGCCGCTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTCGCAGAGGa                                                             >  2:53266/1‑133 (MQ=255)                                             aTACAGGCGATAACAAAGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGATGCCAAACAGACTCGTCAACGCCACTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTCGCAGAGGa                                                             >  2:569099/1‑133 (MQ=255)                                                   gCGATAACAAAGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGATGCCAAACAGACTCGTCAACGCCGCTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTCGCAGAGGAAGACGc                                                       >  2:374857/1‑133 (MQ=255)                                                       tAACAAAGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGATGCCAAACAGACTCGTCAACGCCGCTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTCGCAGAGGAAGACGCCGAg                                                   >  2:606816/1‑133 (MQ=255)                                                       tAACAAAGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGATGCCAAACAGACTCGTCAACGCCGCTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTCGCAGAGGAAGACGCCGAg                                                   >  2:77427/1‑133 (MQ=255)                                                       tAACAAAGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGATGCCAAACAGACTCGTCAACGCCGCTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTCGCAGAGGAAGACGCCGAg                                                   >  1:137760/1‑133 (MQ=255)                                                       tAACAAAGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGCGATGCCAAACAGACTCGTCAACTCCGCTGATCCCGCCACCACTTTTTTCTGCGCATCGAGTTTGCTGAGGAAGACGCCGAg                                                   >  1:202561/1‑133 (MQ=255)                                                         aCAAAGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGATGCCAAACAGACTCGTCAACGCCGCTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTCGCAGAGGAAGACGCCGAgcg                                                 <  1:334547/133‑1 (MQ=255)                                                         aCAAAGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGATGCCAAACAGACTCGTCAACGCCGCTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTCGCAGAGGAAGACGCCGAgcg                                                 <  2:343013/133‑1 (MQ=255)                                                         aCAAAGGGCTACTTATGCAGCATGTTGACGCCATATACCGCTGGTTCGGTGATGCCAAACAGACTCGTCAACGCCGCTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTCGCAGAGGAAGACGCCGAgcg                                                 <  1:207463/133‑1 (MQ=255)                                                                 gTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGATGCCAAACAGACTCGTCAACGCCGCTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTCGCAGAGGAAGACGCCGAGCGCCGCCCCg                                         >  2:141815/1‑133 (MQ=255)                                                                 gTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGATGCCAAACAGACTCGTCAACGCCGCTGATCCCGCCACCACTTTTTTCTGCGCATCGAGTTCGCAGAGGAAGACCCCGATCGCCGCCCCg                                         >  2:198484/1‑133 (MQ=255)                                                                  tACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGATGCCAAACAGACTCGTCAACGCCGCTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTCGCAGAGGAAGACGCCGAGCGCCGCCCCGa                                        >  1:354112/1‑133 (MQ=255)                                                                    cTTACGCGTCAGGTTAACCCCAAATACCGCAGGTTCGGTGATGCCAAACAGACTCGTCAACGCCGCTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTCGCAGAGGAAGACGCCGAGCGCCGCCCCGAcc                                      <  1:211777/133‑1 (MQ=255)                                                                    cTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGATGCCAAACAGACTCGTCAACGCCGCTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTCGCAGAGGAAGACGCCGAgcg                                                 >  1:512106/1‑122 (MQ=255)                                                                    cTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGATGCCAAACAGACTCGTCAACGCCGCTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTCGCAGAGGAAGACGCCGAgcg                                                 <  2:512106/122‑1 (MQ=255)                                                                      tACGCGGCCGGTTGATGCAATACACCGCTGTTTCGGTGATGCCAAACAGCCTCGTCAACGCCGCTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTCGCAGAGGAAGACGCCGAGCGCCGCCCCGACCTg                                    <  1:198484/133‑1 (MQ=255)                                                                      tACGCGGCAGGTTGAGGCCATATAGCGCTGGTTCGGTGATGTCAAACAGACTCGTCAACGCCGCTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTCGCAGAGGAAGACGCCGAGCGCCGCCCCGACCTg                                    <  2:264003/133‑1 (MQ=255)                                                                      tACGCGGCAGGTTGAAGCTATATACCGCTGGTTCGGTGATGCCAAACAGACTCGTCAACGCCGCTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTCGCAGAGGAAGACGCCGAGCGCCGCCCCGACCTg                                    <  1:257696/133‑1 (MQ=255)                                                                         gcgGCAGGTTGACGCCATATACCGCTGGTTCGGTGATGCCAAACAGACTCGTCAACGCCGTTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTCGCAGAGGAAGACGCCGAGCGCCGCCCCGACCTgcgc                                 >  2:475476/1‑133 (MQ=255)                                                                         gcgGCAGGTTGACGCCATATACCGCTGGTTCGGTGATGCCAAACAGACTCGTCAACGCCGCTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTCGCAGAGGAAGACGCCGAGCGCCGCCCCGACCTgcgc                                 >  2:85586/1‑133 (MQ=255)                                                                          cgGCAGGTTGACGCCATATACCGCTGGTTCGGTGATGCCAAACAGACTCGTCAACGCCGCTGATCCCGCCACCACTTTTTTCTGCGCATCGAGTTCGCAGGGGAAGACGCCGAGCGCGGCCCCGACCTGCGcc                                >  2:215486/1‑133 (MQ=255)                                                                           ggCAGGTTGACGCCATATACCGCTGGTTCGGTGATGCCAAACAGACTCGTCAACGCCGCTGATCCCGCCACCACTTTTTTCTGCGCATCGCTTTCGCAGAGGAAAACGGCGAGCGCCGCCCCGACCTGCGCCa                               >  2:211777/1‑133 (MQ=255)                                                                           ggCAGGTTGACGCCATATACCGCTGGTTCGGTGATGCCAAACAGACTCGTCAACGCCGCTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTCGCAGAGGAAGACGCCGAGCGCCGCCCCGACCTGCGCCa                               >  1:481065/1‑133 (MQ=255)                                                                                gTTGACGCCATATACCGCTGGTTCGGTGATGCCAAACAGACTCGTCAACGCCGCTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTCGCAGAGGAAGACGCCGAGCGCCGCTCCGACCTGCACCATAATg                          >  1:264003/1‑133 (MQ=255)                                                                                 ttGACGCCATATACCGCTGGTTCGGTGATGCCAAACAGACTCGTCAACGCCACTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTCGCAGAGGAAGACGCCGAGCGCCGCCCCGACCTGCGCCATAATgg                         <  1:569099/133‑1 (MQ=255)                                                                                   gACGCCATATACCGCTGGTTCGGTGATGCCAAACAGACTCGTCAACGCCGCTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTCGCAGAGGAAGACGCCGAGCGCCGCCCCGACCTGCGCCATAATGGCg                       <  1:85586/133‑1 (MQ=255)                                                                                              ccGCTGGTCCGGTGATGCCAAACAGACTCGCCAACGCCGCTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTCGCAGAGGAAGACGCCGAGCGCCGCCCCGACCTGCGCCATAATGGCGGGCATTAACAg            <  1:163825/133‑1 (MQ=255)                                                                                               cGCTGGTTCGGTGATGCCAAACAGACTCGTCAACGCCGCTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTCGCAGAGGAAGACGCCGCGCGCCGCCACGACCTGAGCCATAATAGCGGGCATTAACAGc           >  2:163825/1‑133 (MQ=255)                                                                                                gCTGGTTCGGTGATGCCAAACAGACTCGTCAACGCCGCTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTCGCAGAGGAAGACGCCGAGCGCCGCCCCGACCTGCGCCATAATGGCGGGCATTAACAGCg          <  1:32231/133‑1 (MQ=255)                                                                                                    gTTCGGTGATGCCAAACAGACTCGTCAACGCCGCTGATCCCGCCACCCCTTTTTTCTGCGCATCGCGTTCGCAGAGGAAGACGCCGAGCGCCGCCCCGACCTGCGCCATAATGGCGGGCATTAACAGCg          <  1:39867/129‑1 (MQ=255)                                                                                                    gTTCGGTGATGCCAAACAGACTCGTCAACGCCGCTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTCGCAGAGGAAGACGCCGAGCGCCGCCCCGACCTGCGCCATAATGGCGGGCATTAACAGCg          >  2:39867/1‑129 (MQ=255)                                                                                                    gTTCGGTGATGCCAAACAGACTCGTCAACGCCGCTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTCGCAGAGGAAGACGCCGAGCGCCGCCCCGACCTGCGCCATAATGGCGGGCATTAACAGCGGGAt      >  2:497200/1‑133 (MQ=255)                                                                                                       cGGTGATGCCAAACAGACTCGTCAACGCCGCTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTCGCAGAGGAAGACGCCGAGCGCCGCCCCGACCTGCGCCATAATGGCGGGCATTAACAGCGGGATCAt   >  2:72418/1‑133 (MQ=255)                                                                                                       cGGTGATGCCAAACAGACTCGTCAACGCCGCTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTCGCAGAGGAAGACGCCGAGCGCCGCCCCGACCTGCGCCATAATGGCGGGCATTAACAGCGGGATCAt   >  1:54002/1‑133 (MQ=255)                                                                                                        ggTGATGCCAAACAGACTCGTCAACGCCGTTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTCGCAGAGGAAGACGCCGAGCGCCGCCCCGACCTGCGCCATAATGGCGGGCATTAACAGCGGGATCATg  <  1:475476/133‑1 (MQ=255)                                                                                                        ggTGATGCCAAACAGACTCGTCAACGCCGCTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTCGCAGAGGAAGACGCCGAGCGCCGCCCCGACCTGCGCCATAATGGCGGGCATTAACAGCGGGATCATg  <  2:54002/133‑1 (MQ=255)                                                                                                        ggTGATGCCAAACAGACTCGTCAACGCCGCTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTCGCAGAGGAAGACGCCGAGCGCCGCCCCGACCTGCGCCATAATGGCGGGCATTAACAGCGGGATCATg  <  2:354112/133‑1 (MQ=255)                                                                                                              |                                                                                                                               TCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGATACCAAACAGACTCGTCAACGCCGCTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTCGCAGAGGAAGACGCCGAGCGCCGCCCCGACCTGCGCCATAATGGCGGGCATTAACAGCGGGATCATG  >  NC_000913/3904391‑3904626

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: ‑