Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,077,640	+G	intergenic (‑291/+23)	cmtB ← / ← tktA	PTS system mannitol‑specific enzyme IIA component/transketolase 1 isozyme

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,077,640	1	.	G	100.0%	185.1 / NA	72	intergenic (‑291/+23)	cmtB/tktA	PTS system mannitol‑specific enzyme IIA component/transketolase 1 isozyme
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (45/27);  total (45/27)

TGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGC‑ACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGC  >  NC_000913/3077565‑3077709                                                                             |                                                                      tgtgATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGc                                                                        <  1:1789245/76‑1 (MQ=255)    gATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGcgac                                                                     <  2:4760326/76‑4 (MQ=255)       cTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCtt                                                                  >  2:2856211/1‑76 (MQ=255)       cTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCtt                                                                  >  1:2368779/1‑76 (MQ=255)        tACAACACGCCTTATCTATTTCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCttt                                                                 <  2:3894723/76‑1 (MQ=255)        tACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCttt                                                                 >  2:5522775/1‑76 (MQ=255)        tACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCttt                                                                 >  1:3867778/1‑76 (MQ=255)        tACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCttt                                                                 >  2:5628212/1‑76 (MQ=255)        tACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCttt                                                                 <  1:5185783/76‑1 (MQ=255)        gacaacaCGCCTTATCGATTGCTTGGCCGCAAACGGACATATCAAGGGAATAAAAAAGGTCGCCGAAGCGACCttt                                                                 >  1:4387740/2‑76 (MQ=255)         acaacaCGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCtttt                                                                >  2:4195447/1‑76 (MQ=255)          caacaCGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCttttt                                                               >  1:1412072/1‑76 (MQ=255)          caacaCGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCttttt                                                               <  1:953569/76‑1 (MQ=255)             cacGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTAcc                                                            <  2:242632/76‑1 (MQ=255)             cacGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTAcc                                                            >  1:5613353/1‑76 (MQ=255)               cGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGGCGAAGCGGCCGTTTTACCCg                                                          >  1:1413677/1‑76 (MQ=255)                 ccTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGaa                                                        <  1:4624342/76‑1 (MQ=255)                  cTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGaca                                                       >  2:3548700/1‑74 (MQ=255)                  cTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGaaa                                                       <  1:257821/76‑1 (MQ=255)                  cTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGaaa                                                       >  2:3957048/1‑76 (MQ=255)                   ttATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGCAAt                                                      >  2:953646/1‑76 (MQ=255)                   ttATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAAt                                                      >  1:5308020/1‑76 (MQ=255)                    tATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATg                                                     <  1:1473760/76‑1 (MQ=255)                    tATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATg                                                     >  2:1135569/1‑76 (MQ=255)                       cTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTa                                                  <  2:5308019/76‑1 (MQ=255)                        tATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTaa                                                 >  2:611567/1‑76 (MQ=255)                         aTTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAAt                                                >  2:467516/1‑76 (MQ=255)                         aTTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAAt                                                <  1:193247/76‑1 (MQ=255)                          ttGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAAtt                                               <  2:3039663/76‑1 (MQ=255)                          ttGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAAtt                                               <  1:4472822/76‑1 (MQ=255)                            gCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTAc                                             >  2:2859075/1‑76 (MQ=255)                            gCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTAc                                             >  1:4470747/1‑76 (MQ=255)                            gCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTAc                                             <  1:535134/76‑1 (MQ=255)                              ttGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTAcag                                           >  2:343490/1‑76 (MQ=255)                                gTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTAcagca                                         >  2:4724026/1‑76 (MQ=255)                                  ccGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGt                                       <  1:4654870/76‑1 (MQ=255)                                      aaaCGGATATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCtt                                   >  1:3393381/1‑76 (MQ=255)                                      aaaCGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCtt                                   <  1:2746047/76‑1 (MQ=255)                                       aaCGGACATATCAAGGTAATAAAGAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCttt                                  <  2:2726503/76‑1 (MQ=255)                                       aaCGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCttt                                  >  2:3117819/1‑76 (MQ=255)                                        aCGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCtttt                                 <  1:916579/76‑1 (MQ=255)                                          ggACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTAACCGAAATGCTAATTACAGCACTTCTTTTGc                               >  1:3894467/1‑76 (MQ=255)                                          ggACATATCAAGGTAATAAAAAAGGGCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGCTCTTTTGc                               >  2:1623908/1‑76 (MQ=38)                                           gACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCt                              >  2:3653005/1‑76 (MQ=255)                                             cATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCttt                            >  1:4401428/1‑76 (MQ=255)                                              atatCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAAGTACAGCAGTTCTTTTGCtttg                           >  1:2806090/1‑75 (MQ=255)                                               tatCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCg                          >  1:3860155/1‑76 (MQ=255)                                               tatCAAGGTAATAAAAAAGGTCGCCGAAGCGACC‑TTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGc                         <  2:4553115/76‑1 (MQ=255)                                               taacAAGGAAATAAAAAAGGACGCCGAAGCGACCTCTTTACCCGAACTGCTGATTACAGCAGTTCTTTTGCTTTCg                          >  2:3597910/4‑76 (MQ=255)                                                atCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGc                         >  2:4932650/1‑76 (MQ=255)                                                 tCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGca                        <  1:3292307/76‑1 (MQ=255)                                                     ggTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGcaacaa                    <  2:2880707/76‑1 (MQ=255)                                                     ggTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGcaacaa                    >  2:5265557/1‑76 (MQ=255)                                                      gTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCaacaac                   >  2:5526624/1‑76 (MQ=255)                                                      gTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCaacaac                   <  2:3536156/76‑1 (MQ=255)                                                      gTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCaacaac                   <  1:3138293/76‑1 (MQ=255)                                                       taataaAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACg                  <  2:2805535/76‑1 (MQ=255)                                                        aataaAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGt                 >  1:3338955/1‑76 (MQ=255)                                                        aataaAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGt                 >  2:2377116/1‑76 (MQ=255)                                                          taaAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTa               <  2:5198260/76‑1 (MQ=255)                                                            aaaaaGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATc             >  2:40804/1‑76 (MQ=255)                                                            aaaaaGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATc             >  1:4657056/1‑76 (MQ=255)                                                             aaaaGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCa            <  2:3800551/76‑1 (MQ=255)                                                              aaaGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCaacaac                   <  1:1088277/68‑1 (MQ=255)                                                              aaaGGGCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCaacaac                   >  2:1088277/1‑68 (MQ=39)                                                               aaGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTACTTTCGCAACAACGGTATCAAc          >  2:3075595/1‑76 (MQ=255)                                                                 ggTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAg        >  1:1447656/1‑76 (MQ=255)                                                                  gTCGCCGAAGCGACCTTTTTCCCCGACATGCTAATTACAGCAGTTCTTTTGCTTTCTCAACAACGTTATCAACAGt       <  2:452101/76‑1 (MQ=255)                                                                  gTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGt       <  1:1078346/76‑1 (MQ=255)                                                                    cGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGa     >  2:4055670/1‑76 (MQ=255)                                                                     gCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGaa    >  1:1134922/1‑76 (MQ=255)                                                                      ccGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCTGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAg   >  2:2298876/1‑76 (MQ=255)                                                                      ccGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCTGCTCTTTTGCTCTCGCAACCACGTTATCCACAGTGCAg   >  2:2315117/1‑76 (MQ=255)                                                                       cGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGc  >  1:4275722/1‑76 (MQ=255)                                                                             |                                                                      TGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGC‑ACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGC  >  NC_000913/3077565‑3077709

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