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%	134.4 / NA	39	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 (20/19);  total (20/19)

TCGCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCC  >  NC_000913/3773502‑3773724                                                                                                     |                                                                                                                           tCGCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCa                                                                                               <  1:1250308/130‑1 (MQ=255) tCGCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCa                                                                                               >  2:1250308/1‑130 (MQ=255)      gATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGc                                                                                 >  2:506989/1‑139 (MQ=255)       aTTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCg                                                                                <  2:149454/139‑1 (MQ=255)       aTTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCg                                                                                <  1:823119/139‑1 (MQ=255)           cTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCt                                                                            >  2:795432/1‑139 (MQ=255)            tGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTg                                                                           <  2:555283/139‑1 (MQ=255)                 cTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCg                                                                                <  2:982025/129‑1 (MQ=255)                 cTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTaca                                                                       >  2:410562/1‑138 (MQ=255)                 cTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTaca                                                                       <  1:410562/138‑1 (MQ=255)                 cTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATACGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCg                                                                                >  1:982025/1‑129 (MQ=255)                   aaCAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACTGAACGCGCTATGCGCCAGGCTCCGCAGGc                                                                                                                   <  2:889589/92‑1 (MQ=255)                   aaCAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACTGAACGCGCTATGCGCCAGGCTCCGCAGGc                                                                                                                   >  1:889589/1‑92 (MQ=255)                          cgATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCa                                                                                               <  1:1369611/105‑1 (MQ=255)                          cgATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCa                                                                                               >  2:1369611/1‑105 (MQ=255)                          cgATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCg                                                                                >  2:162803/1‑120 (MQ=255)                          cgATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCg                                                                                <  1:162803/120‑1 (MQ=255)                                       gccgccAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGc                                                                                                                       >  1:459400/1‑68 (MQ=255)                                       gccgccAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGc                                                                                                                       <  2:459400/68‑1 (MQ=255)                                       gccgccAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGcac                                                                                                                 <  2:1018353/74‑1 (MQ=255)                                       gccgccAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGcac                                                                                                                 >  1:1018353/1‑74 (MQ=255)                                       gccgccAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCg                                                                                <  2:601607/107‑1 (MQ=255)                                       gccgccAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCg                                                                                >  1:601607/1‑107 (MQ=255)                                               aTGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCa                                        >  1:332574/1‑139 (MQ=255)                                                         cGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCt                                                                                        <  2:1119047/81‑1 (MQ=255)                                                         cGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCt                                                                                        >  1:1119047/1‑81 (MQ=255)                                                         cGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGc                                                                                 >  2:651693/1‑88 (MQ=255)                                                         cGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGc                                                                                 <  1:651693/88‑1 (MQ=255)                                                         cGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCt                                                                            <  2:1097551/93‑1 (MQ=255)                                                         cGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCt                                                                            >  1:1097551/1‑93 (MQ=255)                                                          gTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCAGAGCATATTTCGCTGACCAACTTCCTCGACAGc                                                                                 >  1:900083/1‑87 (MQ=255)                                                          gTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGc                                                                                 <  2:900083/87‑1 (MQ=255)                                                            catcaCTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCa                           >  2:1066866/1‑139 (MQ=255)                                                                    acCGTGACCTGACTGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAaga                   >  2:535775/1‑139 (MQ=255)                                                                    acCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAaga                   <  1:403048/139‑1 (MQ=255)                                                                      cGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAgaga                 <  1:795432/139‑1 (MQ=255)                                                                          acctgaccGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACTGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGt             >  2:814281/1‑139 (MQ=255)                                                                                ccGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGa       <  2:87166/139‑1 (MQ=255)                                                                                     cgcgCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGcc  >  2:827442/1‑139 (MQ=255)                                                                                                     |                                                                                                                           TCGCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCC  >  NC_000913/3773502‑3773724

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