Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,904,521	C→T	A359A (GCG→GCA)	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,521	0	C	T	100.0%	136.1 / NA	38	A359A (GCG→GCA)	bglF	beta‑glucoside specific PTS enzyme II/BglG kinase/BglG phosphatase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (20/18);  total (20/18)

ATAACAAAGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGATACCAAACAGACTCGTCAACGCCGCTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTCGCAGAGGAAGACGCCGAGCGCCGCCCCGACCTGCGCCAT  >  NC_000913/3904444‑3904598                                                                              |                                                                              aTAACAAAGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGATACCAAACAGACTCGTCAACgctgct                                                                            <  1:762938/81‑1 (MQ=255)  tAACAAAGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGATACCAAACAGACTCGTCAACgctgct                                                                            >  2:762983/1‑80 (MQ=255)       aaGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGATACCAAACAGACTCGTCAACGCTGCTGATccc                                                                      <  1:762939/81‑1 (MQ=255)          gggTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGATACCAAACAGACTCGTCAACGCTGCTGATCCCGcc                                                                   <  1:762940/81‑1 (MQ=255)          gggTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGATACCAAACAGACTCGTCAACGCTGCTGATCCCGcc                                                                   <  1:762941/81‑1 (MQ=255)            gTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGATACCAAACAGACTCGTCAACGCTGCTGATCCCGccac                                                                 >  1:763011/1‑81 (MQ=255)            gTACTTACGCGGCAGGTTGACGCCATATACAGCTGGTTCGGTGATACCAAACAGACTCGTCAACGCTGCTGATCCCGccac                                                                 <  1:762942/81‑1 (MQ=255)                  aCGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGATACCAAACAGACTCGTCAACGCTGCTGATCCCGCCACCACttt                                                           >  2:763003/1‑81 (MQ=255)                    gcgGCAGGTTGACGCTATATACCGCTGGTTCGGTGATACCAAACAGACTCGTCAACGCTGCTGATCCCGCCACCACttttt                                                         >  2:762986/1‑81 (MQ=255)                      ggCAGGTTGACGCCATATACCGCTGGTTCGGTGATACCAAACAGACTCGTCAACGCTGCTGATCCCGCCACCACTTTTTTc                                                       >  2:763026/1‑81 (MQ=255)                       gCAGGTTGACGCCATATACCGCTGGTTCGGTGATACCAAACAGACTCGTCAACGCTGCTGATCCCGCCACCACTTTTTTCt                                                      >  2:762990/1‑81 (MQ=255)                           gTTGACGCCATATACCGCTGGTTCGGTGATACCAAACAGACTCGTCAACGCTGCTGATCCCGCCACCACTTTTTTCTgcgc                                                  <  2:762943/81‑1 (MQ=255)                           gTTGACGCCATATACCGCTGGTTCGGTGATACCAAACAGACTCGTCAACGCTGCTGATCCCGCCACCACTTTTTTCTgcg                                                   >  2:762998/1‑80 (MQ=255)                               aCGCCATATACCGCTGGTTCGGTGATACCAAACAGACTCGTCAACGCTGCTGATCCCGCCACCACTTTTTTCTGCGCATcg                                              >  1:762984/1‑81 (MQ=255)                                  ccATATACCGCTGGTTCGGTGATACCAAACAGACTCGTCAACGCTGCTGATCCCGCCACCACTTTTTTCTGCGCATCGCGt                                           <  1:762944/81‑1 (MQ=255)                                   cATATACCGCTGGTTCGGTGATACCAAACAGACTCGTCAACGCTGCTGATCCCGCCACCACTTTTTTCTGCGCATCGCGtt                                          >  2:762996/1‑81 (MQ=255)                                    atatACCGCTGGTTCGGTGATACCAAACAGACTCGTCAACGCTGCTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTc                                         >  2:762994/1‑81 (MQ=255)                                        aCCGCTGGTTCGGTGATACCAAACAGACTCGTCAACGCTGCTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTCg                                        >  1:762995/1‑78 (MQ=255)                                        aCCGCTGGTTCGGTGATACCAAACAGACTCGTCAACGCTGCTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTCg                                        >  2:763001/1‑78 (MQ=255)                                         ccGCTGGTTCGGTGATACCAAACAGACTCGTCAACGCTGCTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTCGCag                                     <  2:762945/80‑1 (MQ=255)                                          cGCTGGTTCGGTGATACCAAACAGACTCGTCAACGCTGCTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTCGCagag                                   >  2:762988/1‑81 (MQ=255)                                               gTTCGGTGATACCAAACAGACTCGTCAACGCTGCTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTCGCAGAGGaa                                <  2:762947/79‑1 (MQ=255)                                               gTTCGGTGATACCAAACAGACTCGTCAACGCTGCTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTCGCAGAGGAAg                               <  1:762946/80‑1 (MQ=255)                                               gTTCGGTGATACCAAACAGACTCGTCAACGCTGCTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTCGCAGAGGAAGa                              >  1:762997/1‑81 (MQ=255)                                                ttCGGTGATACCAAACAGACTCGTCAACGCTGCTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTCGCAGAGGAAGAc                             <  2:762948/81‑1 (MQ=255)                                                 tCGGTGATACCAAACAGACTCGTCAACGCTGCTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTCGCAGAGGAAGACg                            >  1:763013/1‑81 (MQ=255)                                                  cGGTGATACCAAACAGACTCGTCAACGCTGCTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTCGCAGAGGAAGACGc                           >  2:763016/1‑81 (MQ=255)                                                   ggTGATACCAAACAGACTCGTCAACGCTGCTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTCGCAGAGGAAGACGcc                          <  1:762949/81‑1 (MQ=255)                                                    gTGATACCAAACAGACTCGTCAACGCTGCTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTCGCAGAGGAAGACGCCg                         <  1:762950/81‑1 (MQ=255)                                                    gTGATACCAAACAGACTCGTCAACGCTGCTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTCGCAGAGGAAGACGCCg                         <  2:762951/81‑1 (MQ=255)                                                       aTACCAAACAGACTCGTCAACGCTGCTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTCGCAGAGGAAGACGCCGAgc                      >  1:763002/1‑81 (MQ=255)                                                          ccAAACAGACTCGTCAACGCTGCTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTCGCAGAGGAAGACGCCGAGcgcc                   <  2:762952/81‑1 (MQ=255)                                                               cAGACTCGTCAACGCTGCTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTCGCAGAGGAAGACGCCGAGCGCCGcccc              <  2:762953/81‑1 (MQ=255)                                                               cAGACTCGTCAACGCTGCTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTCGCAGAGGAAGACGCCGAGCGCCGcccc              <  2:762954/81‑1 (MQ=255)                                                                aGACTCGTCAACGCTGCTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTCGCAGAGGAAGACGCCGAGCGCCGcccc              >  1:763000/1‑80 (MQ=255)                                                                aGACTCGTCAACGCTGCTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTCGCAGAGGAAGACGCCGAGCGCCGCCCCg             <  1:762955/81‑1 (MQ=255)                                                                      gTCAACGCTGCTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTCGCAGAGGAAGACGCCGAGCGCCGCCCCGACCTgc       >  1:763004/1‑81 (MQ=255)                                                                      gTCAACGCTGCTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTCGCAGAGGAAGACGCCGAGCGCCGCCCCGACCTgc       >  2:763012/1‑81 (MQ=255)                                                                      gTCAACGCTGCTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTCGCAGAGGAAGACGCCGAGCGCCGCCCCGACCTgc       <  2:762956/81‑1 (MQ=255)                                                                         aaCGCTGCTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTCGCAGAGGAAGACGCCGAGCGCCGCCCCGACCTgcgc     >  1:763009/1‑80 (MQ=255)                                                                           cgctgcTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTCGCAGAGGAAGACGCCGAGCGCCGCCCCGACCTGCGCCAt  <  2:762957/81‑1 (MQ=255)                                                                              |                                                                              ATAACAAAGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGATACCAAACAGACTCGTCAACGCCGCTGATCCCGCCACCACTTTTTTCTGCGCATCGCGTTCGCAGAGGAAGACGCCGAGCGCCGCCCCGACCTGCGCCAT  >  NC_000913/3904444‑3904598

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