Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,772,302	A→C	K8Q (AAA→CAA)	mtlA →	mannitol‑specific PTS enzyme: IIA, IIB and IIC components

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,772,302	0	A	C	100.0%	202.9 / NA	53	K8Q (AAA→CAA)	mtlA	mannitol‑specific PTS enzyme: IIA, IIB and IIC components
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (22/31);  total (22/31)

ACAGTTAACCGATTCAGTGCCAGATTTCGCAGTATCTACAAGGTCCGGCTACCTCTGCCGCCACATTAACAAAAAACCTCGGGCTTCCAGCCTGCGCGACAGCAAACATAAGAAGGGGTGTTTTTATGTCATCCGATATTAAGATCAAAGTGCAAAGCTTTGGTCGTTTCCTCAGCAACATGGTGATGCCAAATATCGGCGCGTTTATCGCGTGGGGTATCATCACCGCGTTATTTATTCCAACAGGGTGGTTACCGAACGAGACGCTGGCGAAGCTGGTCGGGCCGAT  >  NC_000913/3772156‑3772444                                                                                                                                                   |                                                                                                                                               aCAGTTAACCGATTCAGTGCCAGATTTCGCAGTATCTACAAGGTCCGGCTACCTCTGCCGCCACATTAACAAAAAACCTCGGGCTTCCAGCCTGCGCGACAGCAAACATAAGAAGGGGTGTTTTTATGTCATCCGATATTAAGATCCAAg                                                                                                                                             >  2:924987/1‑150 (MQ=255)     ttAACCGATTCAGTGCCAGATTTCGCAGTATCTACAAGGTCCGGCTACCTCTGCCGCCACATTAACAAAAAACCTCGGGCTTCCAGCCTGCGCGACAGCAAACATAAGAAGGGGTGTTTTTATGTCATCCGATATTAAGATCCAAGTGCa                                                                                                                                         >  1:28294/1‑150 (MQ=255)       aaCCGATTCAGTGCCAGATTTCGCAGTATCTACAAGGTCCGGCTACCTCTGCCGCCACATTAACAAAAAACCTCGGGCTTCCAGCCTGCGCGACAGCAAACATAAGAAGGGGTGTTTTTATGTCATCCGATATTAAGATCCAAGTGCAAAg                                                                                                                                      <  1:17830/151‑1 (MQ=255)               cAGTGCCAGATTTCGCAGTATCTACAAGGTCCGGCTACCTCTGCCGCCACATTAACAAAAAACCTCGGGCTTCCAGCCTGCGCGACAGCAAACATAAGAAGGGGTGTTTTTATGTCATCCGATATTAAGATCCAAGTGCAAAGCTTTGGt                                                                                                                               <  2:882991/150‑1 (MQ=255)                     cAGATTTCGCAGTATCTACAAGGTCCGGCTACCTCTGCCGCCACATTAACAAAAAACCTCGGGCTTCCAGCCTGCGCGACAGCAAACATAAGAAGGGGTGTTTTTATGTCATCCGATATTAAGATCCAAGTGCAAAGCTTTGGTCGTTTcc                                                                                                                        <  2:548981/151‑1 (MQ=255)                       gATTTCGCAGTATCTACAAGGTCCGGCTACCTCTGCCGCCACATTAACAAAAAACCTCGGGCTTCCAGCCTGCGCGACAGCAAACATAAGAAGGGGTGTTTTTATGTCATCCGATATTAAGATCCAAGTGCAAAGCTTTGGTCGTTTCCTc                                                                                                                      <  2:328827/151‑1 (MQ=255)                        aTTTCGCAGTATCTACAAGGTCCGGCTACCTCTGCCGCCACATTAACAAAAAACCTCGGGCTTCCAGCCTGCGCGACAGCAAACATAAGAAGGGGTGTTTTTATGTCATCCGATATTAAGATCCAAGTGCAAAGCTTTGGTCGTTTCCTCa                                                                                                                     <  1:543008/151‑1 (MQ=255)                            cGCAGTATCTACAAGGTCCGGCTACCTCTGCCGCCACATTAACAAAAAACCTCGGGCTTCCAGCCTGCGCGACAGCAAACATAAGAAGGGGTGTTTTTATGTCATCCGATATTAAGATCCAAGTGCAAAGCTTTGGTCGTTTCCTCAGCa                                                                                                                  >  1:61538/1‑150 (MQ=255)                             gCAGTATCTACAAGGTCCGGCTACCTCTGCCGCCACATTAACAAAAAACCTCGGGCTTCCAGCCTGCGCGACAGCAAACATAAGAAGGGGTGTTTTTATGTCATCCGATATTAAGATCCAAGTGCAAAGCTTTGGTCGTTTCCTCAGCAAc                                                                                                                >  1:522197/1‑151 (MQ=255)                                  aTCTACAAGGTCCGGCTACCTCTGCCGCCACATTAACAAAAAACCTCGGGCTTCCAGCCTGCGCGACAGCAAACATAAGAAGGGGTGTTTTTATGTCATCCGATATTAAGATCCAAGTGCAAAGCTTTGGTCGTTTCCTCAGCAACATgg                                                                                                            <  2:799748/150‑1 (MQ=255)                                  aTCTACAAGGTCCGGCTACCTCTGCCGCCACATTAACAAAAAACCTCGGGCTTCCAGCCTGCGCGACAGCAAACATAAGAAGGGGTGTTTTTATGTCATCCGATATTAAGATCCAAGTGCAAAGCTTTGGTCGTTTCCTCAGCAACATgg                                                                                                            <  2:298341/150‑1 (MQ=255)                                  aTCTACAAGGTCCGGCTACCTCTGCCGCCACATTAACAAAAAACCTCGGGCTTCCAGCCTGCGCGACAGCAAACATAAGAAGGGGTGTTTTTATGTCATCCGATATTAAGATCCAAGTGCAAAGCTTTGGTCGTTTCCTCAGCAACATgg                                                                                                            <  2:177438/150‑1 (MQ=255)                                      aCAAGGTCCGGCTACCTCTGCCGCCACATTAACAAAAAACCTCGGGCTTCCAGCCTGCGCGACAGCAAACATAAGAAGGGGTGTTTTTATGTCATCCGATATTAAGATCCAAGTGCAAAGCTTTGGTCGTTTCCTCAGCAACATGGTGAt                                                                                                        <  2:517645/150‑1 (MQ=255)                                              cGGCTACCTCTGCCGCCACATTAACAAAAAACCTCGGGCTTCCAGCCTGCTCGACAGCAAACATAAGAAGGGGTGTTTTTATGTCATCCGATATTAAGATCCAAGTGCAAAGCTTTGGTCGTTTCCTCAGCAACATGGTGATGCCAAatat                                                                                               <  2:677590/151‑1 (MQ=255)                                              cGGCTACCTCTGCCGCCACATTAACAAAAAACCTCGGGCTTCCAGCCTGCGCGACAGCAAACATAAGAAGGGGTGTTTTTATGTCATCCGATATTAAGATCCAAGTGCAAAGCTTTGGTCGTTTCCTCAGCAACATGGTGATGCCAAatat                                                                                               <  1:573995/151‑1 (MQ=255)                                                   aCCTCTGCCGCCACATTAACAAAAAACCTCGGGCTTCCAGCCTGCGCGACAGCAAACATAAGAAGGGGTGTTTTTATGTCATCCGATATTAAGATCCAAGTGCAAAGCTTTGGTCGTTTCCTCAGCAACATGGTGATGCCAAATATCGgc                                                                                           <  1:672380/150‑1 (MQ=255)                                                         gccgccTCATTAACAAAAAACCTCGGGCTTCCAGCCTGCGCGACAGCAAACATAAGAAGGGGTGTTTTTATGTCATCCGATATTAAGATCCAAGTGCAAAGCTTTGGTCGTTTCCTCATCAACATGGTGATGCCAAATATCgg                                                                                            <  2:37388/143‑1 (MQ=255)                                                         gccgccACATTAACAAAAAACCTCGGGCTTCCAGCCTGCGCGACAGCAAACATAAGAAGGGGTGTTTTTATGTCATCCGATATTAAGATCCAAGTGCAAAGCTTTGGTCGTTTCCTCAGCAACATGGTGATGCCAAATATCGGCGCGTTTa                                                                                    <  1:334929/151‑1 (MQ=255)                                                         gccgccACATTAACAAAAAACCTCGGGCTTCCAGCCTGCGCGACAGCAAACATAAGAAGGGGTGTTTTTATGTCATCCGATATTAAGATCCAAGTGCAAAGCTTTGGTCATTTCCTCAGCAACATGGTGATGCCAAATATCGGCGCGttt                                                                                     >  2:129007/1‑150 (MQ=255)                                                                 aTTAACAAAAAACCTCGGGCTTCCAGCCTGCGCGACAGCAAACATAAGAAGGGGTGTTTTTATGTCATCCGATATTAAGATCCAAGTGCAAAGCTTTGGTCGTTTCCTCAGCAACATGGTGATGCCAAATATCGGCGCGTTTATCGCGTgg                                                                            >  1:285624/1‑151 (MQ=255)                                                                    aaCAAAAAACCTCGGGCTTCCAGCCTGCGCGACAGCAAACATAAGAAGGGGTGTTTTTATGTCATCCGATATTAAGATCCAAGTGCAAAGCTTTGGTCGTTTCCTCAGCAACATGGTGATGCCAAATATCGGCGCGTTTATCGCGTGGGGt                                                                         <  2:280093/151‑1 (MQ=255)                                                                      cAAAAAACCTCGGGCTTCCAGCCTGCGCGACAGCAAACATAAGAAGGGGTGTTTTTATGTCATCCGATATTAAGATCCAAGTGCAAAGCTTTGGTCGTTTCCTCAGCAACATGGTGATGCCAAATATCGGCGCGTTTATCGCGTGGGGTat                                                                       >  2:646453/1‑151 (MQ=255)                                                                         aaaaCCTCGGGCTTCCAGCCTGCGCGACAGCAAACATAAGAAGGGGTGTTTTTATGTCATCCGATATTAAGATCCAAGTGCAAAGCTTTGGTCGTTTCCTCAGCAACATGGTGATGCCAAATATCGGCGCGTTTATCGCGTGGGGTatcat                                                                    >  2:214936/1‑151 (MQ=255)                                                                             ccTCGGGCTTCCAGCCTGCGCGACAGCAAACATAAGAAGGGGTGTTTTTATGTCATCCGATATTAAGATCCAAGTGCAAAGCTTTGGTCGTTTCCTCAGCAACATGGTGATGCCAAATATCGGCGCGTTTATCGCGTGGGGTATCATCAcc                                                                <  2:376097/151‑1 (MQ=255)                                                                             ccTCGGGCTTCCAGCCTGCGCGACAGCAAACATAAGAAGGGGTGTTTTTATGTCATCCGATATTAAGATCCAAGTGCAAAGCTTTGGTCGTTTCCTCAGCAACATGGTGATGCCAAATATCGGCGCGTTTATCGCGTGGGGTATCATCAcc                                                                <  2:446169/151‑1 (MQ=255)                                                                                cGGGCTTCCAGCCTGCGCGACAGCAAACATAAGAAGGGGTGTTTTTATGTCATCCGATATTAAGATCCAAGTGCAAAGCTTTGGTCGTTTCCTCAGCAACATGGTGATGCCAAATATCGGCGCGTTTATCGCGTGGGGTATCATCACcgc                                                              <  2:790872/150‑1 (MQ=255)                                                                                  ggCTTCCAGCCTGCGCGACAGCAAACATAAGAAGGGGTGTTTTTATGTCATCCGATATTAAGATCCAAGTGCAAAGCTTTGGTCGTTTCCTCAGCAACATGGTGATGCCAAATATCGGCGCGTTTATCGCGTGGGGTATCATCACCGCGtt                                                           <  2:313966/151‑1 (MQ=255)                                                                                      tCCAGCCTGCGCGACAGCAAACATAAGAAGGGGTGTTTTTATGTCATCCGATATTAAGATCCAAGTGCAAAGCTTTGGTCGTTTCCTCAGCAACATGGTGATGCCAAATATCGGCGCGTTTATCGCGTGGGGTATCATCACCGCGttat                                                         <  2:535938/149‑1 (MQ=255)                                                                                       ccAGCCTGCGCGACAGCAAACATAAGAAGGGGTGTTTTTATGTCATCCGATATTAAGATCCAAGTGCAAAGCTTTGGTCGTTTCCTCAGCAACATGGTGATGCCAAATATCGGCGCGTTTATCGCGTGGGGTATCATCACCGCGttattt                                                       <  1:81556/150‑1 (MQ=255)                                                                                        cAGCCTGCGCGACAGCAAACATAAGAAGGGGTGTTTTTATGTCATCCGATATTAAGATCCAAGTGCAAAGCTTTGGTCGTTTCCTCAGCAACATGGTGATGCCAAATATCGGCGCGTTTATCGCGTGGGGTATCATCACCGCGttatttat                                                     <  1:611273/151‑1 (MQ=255)                                                                                                gccaCAGCAAACATAAGAAGGGGTGTTTTTATGTCATCCGATATTAAGATCCCAGTGCAAAGCTTTGGTCGTTTCCTCAGCAACATGGTGATGCCAAATAACGGGGCGTATATCGAGTGGGGTATAATCACCGCGTCATTTATTTCAAcgg                                             >  2:921221/4‑149 (MQ=255)                                                                                                 cgACAGCAGACATAAGAAGGGGTGTTTTTATGGTATCCTATATTAAGAACCAAGTGGAAAACTTTGGTGTTTTCCTCGACAAATGGGTGGAGCCAAATATGAGCGCTTTTATCGCTTGGGGTGGCAGCACCGGGTATTTTCTTCAAACAg                                             <  2:547804/150‑1 (MQ=255)                                                                                                   aCAGCAAACATAAGAAGGGGTGTTTTTATGTCATCCGATATTAAGATCCAAGTGCAAAGCTTTGGTCGTTTCCTCAGCAACATGGTGATGCCAAATATCGGCGCGTTTATCGCGTGGGGTATCATCACCGCGTTATTTATTCCAACAGggt                                          <  1:150518/151‑1 (MQ=255)                                                                                                       cAAACATAAGAAGGGGTGTTTTTATGTCATCCGATATTAAGATCCAAGTGCAAAGCTTTGGTCGTTTCCTCAGCAACATGGTGATGCCAAATATCGGCGCGTTTATCGCGTGGGGTATCATCACCGCGTTATTTATTCCAACAGGGTGGtt                                      <  2:777682/151‑1 (MQ=255)                                                                                                         aaCATAAGAAGGGGTGTTTTTATGTCATCCGATATTAAGATCCAAGTGCAAAGCTTTGGTCGTTTCCTCAGCAACATGGTGATGCCAAATATCGGCGCGTTTATCGCGTGGGGTATCATCACCGCGTTATTTATTCCAACAGGGTGGTTAc                                    >  1:371398/1‑151 (MQ=255)                                                                                                           cATAAGAAGGGGTGTTTTTATGTCATCCGATATTAAGATCCAAGTGCAAAGCTTTGGTCGTTTCCTCAGCAACATGGTGATGCCAAATATCGGCGCGTTTATCGCGTGGGGTATCATCACCGCGTTATTTATTCCAACAGGGTGGTTACCg                                  <  2:28294/151‑1 (MQ=255)                                                                                                               agaagGGGTGTTTTTATGTCATCCGATATTAAGATCCAAGTGCAAAGCTTTGGTCGTTTCCTCAGCAACATGGTGATGCCAAATATCGGCGCGTTTATCGCGTGGGGTATCATCACCGCGTTATTTATTCCAACAGGGTGGTTACCGAACg                              >  1:799116/1‑151 (MQ=255)                                                                                                                gaagGGGTGTTTTTATGTCATCCGATATTAAGATCCAAGTGCAAAGCTTTGGTCGTTTCCTCAGCAACATGGTGATGCCAAATATCGGCGCGTTTATCGCGTGGGGTATCATCACCGCGTTATTTATTCCAACAGGGTGGTTACCGAACga                             <  2:770360/151‑1 (MQ=255)                                                                                                                gaagGGGTGTTTTTATGTCATCCGATATTAAGATCCAAGTGCAAAGCTTTGGTCGTTTCCTCAGCAACATGGTGATGCCAAATATCGGCGCGTTTATCGCGTGGGGTATCATCACCGCGTTATTTATTCCAACAGGGTGGTTACCGAACg                              >  2:851585/1‑150 (MQ=255)                                                                                                                    gggTGTTTTTATGTCATCCGATATTAAGATCCAAGTGCAAAGCTTTGGTCGTTTCCTCAGCAACATGGTGATGCCAAATATCGGCGCGTTTATCGCGTGGGGTATCATCACCGCGTTATTTATTCCAACAGGGTGGTTACCGAACgaga                           >  2:683759/1‑149 (MQ=255)                                                                                                                    gggTGTTTTTATGTCATCCGATATTAAGATCCAAGTGCAAAGCTTTGGTCGTTTCCTCAGCAACATGGTGATGCCAAATATCGGCGCGTTTATCGCGTGGGGTATCATCACCGCGTTATTTATTCCAACAGGGTGGTTACCGAACGAGACg                         >  2:81349/1‑151 (MQ=255)                                                                                                                      gtgtTTTTATGTCATCCGATATTAAGATCCAAGTGCAAAGCTTTGGTCGTTTCCTCAGCAACATGGTGATGCCAAATATCGGCGCGTTTATCGCGTGGGGTATCATCACCGCGTTATTTATTCCAACAGGGTGGTTACCGAACGAGACGc                        >  2:427280/1‑150 (MQ=255)                                                                                                                      gtgtTTTTATGTCATCCGATATTAAGATCCAAGTGCAAAGCTTTGGTCGTTTCCTCAGCAACATGGTGATGCCAAATATCGGCGCGTTTATCGCGTGGGGTATCATCACCGCGTTATTTATTCCAACAGGGTGGTTACCGAACGAGACGCt                       >  2:311557/1‑151 (MQ=255)                                                                                                                      gtgtTTTTATGTCATCCGATATTAAGATCCAAGTGCAAAGCTTTGGTCGTTTCCTCAGCAACATGGTGATGCCAAATATCGGCGCGTTTATCGCGTGGGGTATCATCACCGCGTTATTTATTCCAACAGGGTGGTTACCGAACGAGACGCt                       <  2:87463/151‑1 (MQ=255)                                                                                                                                gTCATCCGATATTAAGATCCAAGTGCAAAGCTTTGGTCGTTTCCTCAGCAACATGGTGATGCCAAATATCGGCGCGTTTATCGCGTGGGGTATCATCACCGCGTTATTTATTCCAACAGGGTGGTTACCGAACGAGACGCTGGCGAAGCt              <  2:367099/150‑1 (MQ=255)                                                                                                                                   aTCCGATATTAAGATCCAAGTGCAAAGCTTTGGTCGTTTCCTCAGCAACATGGTGATGCCAAATATCGGCGCGTTTATCGCGTGGGGTATCATCACCGCGTTATTTATTCCAACAGGGTGGTTACCGAACGAGACGCTGGCGAAGCTGGTc          >  1:725994/1‑151 (MQ=255)                                                                                                                                     ccGATATTAAGATCCAAGTGCAAAGCTTTGGTCGTTTCCTCAGCAACATGGTGATGCCAAATATCGGCGCGTTTATCGCGTGGGGTATCATCACCGCGTTATTTATTCCAACAGGGTGGTTACCGAACGAGACGCTGGCGAAGCTGGTCgg        <  2:752758/151‑1 (MQ=255)                                                                                                                                       gATATTAAGATCCAAGTGCAAAGCTTTGGTCGTTTCCTCAGCAACATGGTGATGCCAAATATCGGCGCGTTTATCGCGTGGGGTATCATCACCGCGTTATTTATTCCAACAGGGTGGTTACCGAACGAGACGCTGGCGAAGCTGGTCGGGc      >  2:842935/1‑151 (MQ=255)                                                                                                                                         tatTAAGATCCAAGTGCAAAGCTTTGGTCGTTTCCTCAGCAACATGGTGATGCCAAATATCGGCGCGTTTATCGCGTGGGGTATCATCACCGCGTTATTTATTCCAACAGGGTGGTTACCGAACGAGACGCTGGCGAAGCTGGTCGGGCCg    >  2:5608/1‑151 (MQ=255)                                                                                                                                         tatTAAGATCCAAGTGCAAAGCTTTGGTCGTTTCCTCAGCAACATGGTGATGCCAAATATCGGCGCGTTTATCGCGTGGGGCATCATCACCGCGTTATTTATTCCAACAGGGTGGTTACCGAACGAGACGCTGGCGAAGCTGGTCGGGCCg    >  2:360545/1‑151 (MQ=255)                                                                                                                                          atTAAGATCCAAGTGCAAAGCTTTGGTCGTTTCCTCAGCAACATGGTGATGCCAAATATCGGCGCGTTTATCGCGTGGGGTATCATCACCGCGTTATTTATTCCAACAGGGTGGTTACCGAACGAGACGCTGGCGAAGCTGGTCGGGCCg    <  1:85387/150‑1 (MQ=255)                                                                                                                                           ttAAGATCCAAGTGCAAAGCTTTGGTCGTTTCCTCAGCAACATGGTGATGCCAAATATCGGCGCGTTTATCGCGTGGGGTATCATCACCGCGTTATTTATTCCAACAGGGTGGTTACCGAACGAGACGCTGGCGAAGCTGGTCGGGCCgat  >  1:82087/1‑151 (MQ=255)                                                                                                                                           ttAAGATCCAAGTGCAAAGCTTTGGTCGTTTCCTCAGCAACATGGTGATGCCAAATATCGGCGCGTTTATCGCGTGGGGTATCATCACCGCGTTATTTATTCCAACAGGGTGGTTACCGAACGAGACGCTGGCGAAGCTGGTCGGGCCgat  >  2:731950/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                               ACAGTTAACCGATTCAGTGCCAGATTTCGCAGTATCTACAAGGTCCGGCTACCTCTGCCGCCACATTAACAAAAAACCTCGGGCTTCCAGCCTGCGCGACAGCAAACATAAGAAGGGGTGTTTTTATGTCATCCGATATTAAGATCAAAGTGCAAAGCTTTGGTCGTTTCCTCAGCAACATGGTGATGCCAAATATCGGCGCGTTTATCGCGTGGGGTATCATCACCGCGTTATTTATTCCAACAGGGTGGTTACCGAACGAGACGCTGGCGAAGCTGGTCGGGCCGAT  >  NC_000913/3772156‑3772444

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