Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,903,280	G→A	G67S (GGT→AGT)	manY →	mannose‑specific PTS enzyme IIC component

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,903,280	0	G	A	100.0%	145.9 / NA	40	G67S (GGT→AGT)	manY	mannose‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/27);  total (13/27)

ATATGAAAACCGGTATTATTATCGGTGGTACGCTGGAAATGATCGCGCTGGGCTGGATGAACATCGGTGCTGCAGTTGCGCCTGACGCCGCTCTGGCTTCTATCATTTCTACCATTCTGGTTATCGCAGG  >  NC_000913/1903215‑1903344                                                                  |                                                                 atatGAAAACCGGTATTATTATCGGTGGTACGCTGGAAATGATCGCGCTGGGCTGGATGAACATCAGt                                                                <  1:570429/68‑1 (MQ=255) atatGAAAACCGGTATTATTATCGGTGGTACGCTGGAAATGATCGCGCTGGGCTGGATGAACATCAGt                                                                <  1:4183093/68‑1 (MQ=255) atatGAAAACCGGTATTATTATCGGTGGTACGCTGGAAATGATCGCGCTGGGCTGGATGAACATCAGt                                                                <  1:4174851/68‑1 (MQ=255) atatGAAAACCGGTATTATTATCGGTGGTACGCTGGAAATGATCGCGCTGGGCTGGATGAACATCAGt                                                                <  1:2677787/68‑1 (MQ=255) atatGAAAACCGGTATTATTATCGGTGGTACGCTGGAAATGATCGCGCTGGGCTGGATGAACATCAGt                                                                <  1:243268/68‑1 (MQ=255) atatAAAAACCGGTATTAGTATCGGTGGTACGCTGGAAATGAGCGCGCTGGGCTGGATGAACATCAGt                                                                <  1:2186395/68‑1 (MQ=255)       aaaCCGGTATTATTATCGGTGGTACGCTGGAAATGATCGCGCTGGGCTGGATGAACATCAGTGCTGCa                                                          <  1:1069274/68‑1 (MQ=255)                 tattatCGGTGGTACGCTGGAAATGATCGCGCTGGGCTGGATGAACATCAGTGCTGCAGTTGCGCCTg                                                >  1:3459121/1‑68 (MQ=255)                 tattatCGGTGGTACGCTGGAAATGATCGCGCTGGGCTGGATGAACATCAGTGCTGCAGTTGCGCCTg                                                >  1:2773410/1‑68 (MQ=255)                       cGGTGGTACGCTGGAAATGATCGCGCTGGGCTGGATGAACATCAGTGCTGCAGTTGCGCCTGAcgccg                                          >  1:1303306/1‑68 (MQ=255)                        ggtggtACGGTGGAAATGATCGCGCTGGGCTGGATGAACATCAGTGCTGCAGTTGCGCCTGAcgccgc                                         <  1:3107389/68‑1 (MQ=255)                        ggtggtACGCTGGAAATGATCGCGCTGGGCTGGATGAACATCAGTGCTGCAGTTGCGCCTGAcgccgc                                         <  1:1064508/68‑1 (MQ=255)                        ggtggtACGCTGGAAATGATCGCGCTGGGCTGGATGAACATCAGTGCTGCAGTTGCGCCTGAcgccgc                                         <  1:4143121/68‑1 (MQ=255)                        ggtggtACGCTGGAAATGATCGCGCTGGGCTGGATGAACATCAGTGCTGCAGTTGCGCCTGAcgccgc                                         <  1:3385114/68‑1 (MQ=255)                        ggtggtACGCTGGAAATGATCGCGCTGGGCTGGATGAACATCAGTGCTGCAGTTGCGCCTGAcgccgc                                         <  1:3151792/68‑1 (MQ=255)                        ggtggtACGCTGGAAATGATCGCGCTGGGCTGGATGAACATCAGTGCTGCAGTTGCGCCTGAcgccgc                                         <  1:267328/68‑1 (MQ=255)                        ggtggtACGCTGGAAATGATCGCGCTGGGCTGGATGAACATCAGTGCTGCAGTTGCGCCTGAcgccgc                                         <  1:2507547/68‑1 (MQ=255)                        ggtggtACGCTGGAAATGATCGCGCTGGGCTGGATGAACATCAGTGCTGCAGTTGCGCCTGAcgccgc                                         <  1:215917/68‑1 (MQ=255)                        ggtggtACGCTGGAAATGATCGCGCTGGGCTGGATGAACATCAGTGCTGCAGTTGCGCCTGAcgccgc                                         <  1:2587816/68‑1 (MQ=255)                             tACGCTGGAAATGATCGCGCTGGGCTGGATGAACATCAGTGCTGCAGTTGCGCCTGACGCCGCTCTgg                                    <  1:1720022/68‑1 (MQ=255)                              aCGCTGGAAATGATCGCGCTGGGCTGGATGAACATCAGTGCTGCAGTTGCGCCTGACGCCGCTCTGGc                                   <  1:1380075/68‑1 (MQ=255)                               agCTGGAAATGATCGCGCTGGGCTGGATGAACATCAGTGCTGCAGTTGCGCCTGACGCCGCTCTGGCt                                  <  1:2093800/67‑1 (MQ=255)                                gCTGGAAATGATCGCGCTGGGCTGGATGAACATCAGTGCTGCAGTTGCGCCTGACGCCGCTCTGGCtt                                 >  1:1031252/1‑68 (MQ=255)                                gCTGGAAATGATCGCGCTGGGCTGGATGAACATCAGTGCTGCAGTTGCGCCTGACGCCGCTCTGGCtt                                 >  1:2095947/1‑68 (MQ=255)                                gCTGGAAATGATCGCGCTGGGCTGGATGAACATCAGTGCTGCAGTTGCGCCTGACGCCGCTCTGGCtt                                 >  1:3430771/1‑68 (MQ=255)                                  tGGAAATGATCGCGCTGGGCTGGATGAACATCAGTGCTGCAGTTGCGCCTGACGCCGCTCTGGCTTCt                               <  1:62369/68‑1 (MQ=255)                                      aaTGATCGCGCTGGGCTGGATGAACATCAGTGCTGTAGTTGCGCCTGACGCCGCTCTGGCTTCTATCa                           >  1:2754771/1‑68 (MQ=255)                                        tGATCGCGCTGGGCTGGATGAACATCAGTGCTGCAGTTGCGCCTGACGCCGCTCTGGCTTCTTTTAtt                         >  1:3505018/1‑68 (MQ=255)                                        tGATCGCGCTGGGCTGGATGAACATCAGTGCTGCAGTTGCGCCTGACGCCGCTCTGGCTTCTCTCAtt                         <  1:2281067/68‑1 (MQ=255)                                        tGATCGCGCTGGGCTGGATGAACATCAGTGCTGCAGTTGCGCCTGACGCCGCTCTGGCTTCTATCact                         >  1:1949818/1‑66 (MQ=255)                                        tGATCGCGCTGGGCTGGATGAACATCAGTGCTGCAGTTGCGCCTGACGCCGCTCTGGCTTCTATCAtt                         >  1:3701999/1‑68 (MQ=255)                                           tCGCGCTGGGCTGGATGAACATCAGTGCTGCAGTTGCGCCTGACGCCGCTCTGGCTTCTATCATTTCt                      <  1:3500899/68‑1 (MQ=255)                                                ctgggctggATGAACATCAGTGCTGCAGTTGCGCCTGACGCCGCTCTGGCTTCTATCATTTCTACCAt                 <  1:3807504/68‑1 (MQ=255)                                                     cttgATGAACATCAGTGCTGCAGTTGCGCCTGACGCCGCTCTGGCTTCTATCATTTCTACCATTCTgg            <  1:3624434/65‑1 (MQ=255)                                                     ctggATGAACATCAGTGCTGCAGTTGCGCCTGACGCCGCTCTGGCTTCTATCATTTCTACCATTCTgg            <  1:952368/68‑1 (MQ=255)                                                         aTGAACATCAGTGCTGCAGTTGCGCCTGACGCCGCTCTGGCTTCTATCATTTCTACCATTCTGGTTAt        >  1:2714336/1‑68 (MQ=255)                                                          tGAACATCAGTGCTGCAGTTGCGCCTGACGCCGCTCTGGCTTCTATCATTTCTACCATTCTGGTTATc       <  1:4046773/68‑1 (MQ=255)                                                           gAACATCAGTGCTGCAGTTGCGCCTGACGCCGCTCTGGCTTCTATCATTTCTACCATTCTGGTTATCg      >  1:492682/1‑68 (MQ=255)                                                             aCATCAGTGCTGCAGTTGCGCCTGACGCCGCTCTGGCTTCTATCATTTCTACCATTCTGGTTATCGCa    <  1:201308/68‑1 (MQ=255)                                                              cATCAGTGCTGCAGTTGCGCCTGACGCCGCTCTGGCTTCTATCATTTCTACCATTCTGGTTATCGCAg   >  1:2254129/1‑68 (MQ=255)                                                               aTCAGTGCTGCAGTTGCGCCTGACGCCGCTCTGGCTTCTATCATTTCTACCATTCTGGTTATCGCAgg  <  1:3854211/68‑1 (MQ=255)                                                               aTCAGTGCTGCAGTTGCGCCTGACGCCGCTCTGGCTTCTATCATTTCTACCATTCTGGTTATCGCAgg  >  1:326330/1‑68 (MQ=255)                                                                  |                                                                 ATATGAAAACCGGTATTATTATCGGTGGTACGCTGGAAATGATCGCGCTGGGCTGGATGAACATCGGTGCTGCAGTTGCGCCTGACGCCGCTCTGGCTTCTATCATTTCTACCATTCTGGTTATCGCAGG  >  NC_000913/1903215‑1903344

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