Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,905,683	G→A	intergenic (‑86/+48)	bglF ← / ← bglG	fused beta‑glucoside‑specific PTS enzymes: IIA component/IIB component/IIC component/transcriptional antiterminator of the bgl operon

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,905,683	0	G	A	100.0%	149.6 / NA	47	intergenic (‑86/+48)	bglF/bglG	fused beta‑glucoside‑specific PTS enzymes: IIA component/IIB component/IIC component/transcriptional antiterminator of the bgl operon
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/28);  total (19/28)

TATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCACGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAA  >  NC_000913/3905549‑3905792                                                                                                                                       |                                                                                                              tATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGAtttt                                                                                                           >  2:456707/1‑139 (MQ=255)    cTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAAACCTTAAAAAAACGGCCTAACCTGAAACAAACAATTTTTTTTAAATTTTgtc                                                                                                        >  2:286788/1‑137 (MQ=255)     tGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGATTTTGccc                                                                                                       >  1:20681/1‑139 (MQ=255)      gcgcCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGATTTTGCCCg                                                                                                      >  2:201488/1‑139 (MQ=255)      gcgcCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGATTTTGCCCg                                                                                                      >  1:674519/1‑139 (MQ=255)        gcCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAAACCTGAAAAAAACGTCCTGCCTTGAATCAAGCAATTTTTTTCAGAtttttatcact                                                                                                    >  2:375220/1‑132 (MQ=255)        gcCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGATTTTGCCCGCt                                                                                                    >  2:799604/1‑139 (MQ=255)        gcCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGATTTTGCCCGCt                                                                                                    >  1:526821/1‑139 (MQ=255)         cccccGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGATTTTGCCCGCtt                                                                                                   <  1:201488/139‑1 (MQ=255)         cccccGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGATTTTGCCCGCtt                                                                                                   <  2:674519/139‑1 (MQ=255)           cccGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGATTTTGCCCGCTTAg                                                                                                 <  1:799604/139‑1 (MQ=255)           cccGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGATTTTGCCCGCTTAg                                                                                                 <  1:456707/139‑1 (MQ=255)             cGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGATTTTGCCCGCTTAGTg                                                                                               <  1:402291/139‑1 (MQ=255)             cGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGATTTTGCCCGCTTAGTg                                                                                               <  1:672774/139‑1 (MQ=255)              gACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGATTTTGCCCGCTTAGTGc                                                                                              <  1:258833/139‑1 (MQ=255)                 tCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGATTTTGCCCGCTTAGTGCGGt                                                                                           <  2:20681/139‑1 (MQ=255)                   cTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGATTTTGCCCGCTTAGTGCGGTaa                                                                                         <  2:200516/139‑1 (MQ=255)                    tGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGATTTTGCCCGCTTAGTGCGGTaa                                                                                         <  2:733476/138‑1 (MQ=255)                    tGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGATTTTGCCCGCTTAGTGCGGTaa                                                                                         >  1:733476/1‑138 (MQ=255)                     gCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGATTTTGCCCGCTTAGTGCGGTAACa                                                                                       <  2:310751/139‑1 (MQ=255)                              tttCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGATTTTGCCCGCTTAGTGCGGTAACAATCCTTTAc                                                                              <  1:514758/139‑1 (MQ=255)                              tttCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGATTTTGCCCGCTTAGTGCGGTAACAATCCTTTAc                                                                              <  2:10779/139‑1 (MQ=255)                                      tAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGATTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGtaat                                                                      >  2:426577/1‑139 (MQ=255)                                      tAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGATTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGtaat                                                                      >  1:230367/1‑139 (MQ=255)                                              tCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGATTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTc                                                              >  1:206548/1‑139 (MQ=255)                                                       gCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGATTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCtt                                                     >  2:356460/1‑139 (MQ=255)                                                             tACCGCTTTGCGGCAAAACTCCAAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGATTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCa                                               <  2:221819/139‑1 (MQ=255)                                                             tACCGCTTTGCGGCAAAACTCCAAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGATTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCa                                               <  1:426577/139‑1 (MQ=255)                                                               ccGCTTTGCGGCAAAACTCCAAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGATTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCGGTGTTCTTTGCGCAcg                                             <  1:434412/139‑1 (MQ=255)                                                               ccGCTTTGCGGCAAAACTCCAAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGATTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCAcg                                             <  1:356460/139‑1 (MQ=255)                                                               ccGCTTTGCGGCAAAACTCCAAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGATTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCAcg                                             <  2:526821/139‑1 (MQ=255)                                                                  cTTTGCGGCAAAACTCCAAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGATTTTGCCCGCTTAGTTCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCAcg                                             >  2:590237/1‑136 (MQ=255)                                                                  cTTTGCGGCAAAACTCCAAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGATTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTc                                                                            >  2:56110/1‑105 (MQ=255)                                                                  cTTTGCGGCAAAACTCCAAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGATTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTc                                                                            <  1:56110/105‑1 (MQ=255)                                                                  cTTTGCGGCAAAACTCCAAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGATTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCAcg                                             <  1:590237/136‑1 (MQ=255)                                                                    ttGCGGCAAAACTCCAAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGATTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCACGCGctc                                        <  2:411132/139‑1 (MQ=255)                                                                       cGGCAAAACTCCAAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGATTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCAcgc                                            >  1:276757/1‑132 (MQ=255)                                                                       cGGCAAAACTCCAAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGATTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCAcgc                                            <  2:276757/132‑1 (MQ=255)                                                                                    aaaaaaaaCCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGATTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCACGCGCTCTATATTTATGGCTaaa                        <  1:18017/139‑1 (MQ=255)                                                                                    aaaaaaaaCCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGATTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCACGCGCTCTATATTTATGGCTaaa                        <  2:206548/139‑1 (MQ=255)                                                                                       aaaaaCCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGATTTTGCCCGCTTAGTGCGGTAACAATCCTTTACt                                                                             >  2:570140/1‑83 (MQ=255)                                                                                       aaaaaCCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGATTTTGCCCGCTTAGTGCGGTAACAATCCTTTACt                                                                             <  1:570140/83‑1 (MQ=255)                                                                                                       cGGCCTGACGTGAATCAAGCAATTTTTTTCAGATTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCAcgcg                                           <  2:425864/101‑1 (MQ=255)                                                                                                       cGGCCTGACGTGAATCAAGCAATTTTTTTCAGATTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCAcgcg                                           >  1:425864/1‑101 (MQ=255)                                                                                                          ccTGACGTGAATCAAGCAATTTTTTTCAGATTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCACGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTg     >  1:207691/1‑136 (MQ=255)                                                                                                          ccTGACGTGAATCAAGCAATTTTTTTCAGATTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCACGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTg     <  2:207691/136‑1 (MQ=255)                                                                                                          ccTGACGTGAATCAAGCAATTTTTTTCAGATTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCACGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGaaa  >  1:109794/1‑139 (MQ=255)                                                                                                                        aGCAATTTTTTTCAGATTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCACGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGt      <  1:763524/121‑1 (MQ=255)                                                                                                                        aGCAATTTTTTTCAGATTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCACGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGt      >  2:763524/1‑121 (MQ=255)                                                                                                                                       |                                                                                                              TATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCACGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAA  >  NC_000913/3905549‑3905792

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