Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,773,602	T→C	V441A (GTT→GCT)	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,773,602	0	T	C	100.0%	108.4 / NA	31	V441A (GTT→GCT)	mtlA	mannitol‑specific PTS enzyme: IIA, IIB and IIC components
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (16/15);  total (16/15)

CAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAG  >  NC_000913/3773522‑3773734                                                                                 |                                                                                                                                     cAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGc                                                                            <  1:139130/139‑1 (MQ=255)  aGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCg                                                                                          >  1:296715/1‑124 (MQ=255)  aGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCg                                                                                          <  2:296715/124‑1 (MQ=255)    cgcgATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCg                                                                                          <  2:989267/122‑1 (MQ=255)    cgcgATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCg                                                                                          >  1:989267/1‑122 (MQ=255)    cgcgATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTg                                                                         >  1:226300/1‑139 (MQ=255)    cgcgATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTg                                                                         >  1:839798/1‑139 (MQ=255)    cgcgATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATACGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTg                                                                         >  2:818115/1‑139 (MQ=255)          caacaaCCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGaa                                                                   <  2:839798/139‑1 (MQ=255)                   gccgccAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGtt                                                          >  1:560052/1‑139 (MQ=255)                   gccgccAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGtt                                                          >  1:477508/1‑139 (MQ=255)                   gccgccAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGtt                                                          >  1:392395/1‑139 (MQ=255)                    ccgccAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTtg                                                         <  2:392395/139‑1 (MQ=255)                                    tCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCacac                                         >  1:933952/1‑139 (MQ=255)                                     cGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGc                                                                                        >  1:523082/1‑91 (MQ=255)                                     cGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGc                                                                                        <  2:523082/91‑1 (MQ=255)                                      gTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCata                                                                                                                     <  1:963670/61‑1 (MQ=255)                                      gTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCata                                                                                                                     >  2:963670/1‑61 (MQ=255)                                      gTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACgg                                       <  2:350742/139‑1 (MQ=255)                                                acCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCg                                                                                          >  1:395521/1‑78 (MQ=255)                                                acCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCg                                                                                          <  2:395521/78‑1 (MQ=255)                                                  cGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAaga                             <  1:793628/137‑1 (MQ=255)                                                  cGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAaga                             >  2:793628/1‑137 (MQ=255)                                                  cGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAgaga                           <  2:651037/139‑1 (MQ=255)                                                  cGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAgaga                           >  1:37790/1‑139 (MQ=255)                                                  cGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAgaga                           <  2:477508/139‑1 (MQ=255)                                                  cGTGACCTGACCGAACGCGCTATACGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAgaga                           <  1:818115/139‑1 (MQ=255)                                                                aCGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGc             >  1:869545/1‑139 (MQ=255)                                                                     cTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACgg                                       <  2:481052/108‑1 (MQ=255)                                                                     cTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACgg                                       >  1:481052/1‑108 (MQ=255)                                                                           gcCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAg  <  2:37790/139‑1 (MQ=255)                                                                                 |                                                                                                                                     CAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAG  >  NC_000913/3773522‑3773734

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