Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	2,176,396	C→T	M6I (ATG→ATA)	gatC ←	galactitol‑specific enzyme IIC component of PTS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,176,396	0	C	T	100.0%	333.2 / NA	98	M6I (ATG→ATA)	gatC	galactitol‑specific enzyme IIC component of PTS
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (66/32);  total (66/32)

AAATAATGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATAACGCATGACTTCTGAAAACATAGGTCACCCCTGTAAGATAGTCAGAATTTTATTTTGTA  >  W3110S.gb/2176341‑2176451                                                        |                                                        aaaTAATGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATAACGTATGAc                                                    >  1:650295/1‑61 (MQ=255) aaaTAATGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATAACGTATGACt                                                   >  1:4583967/1‑62 (MQ=255) aaaTAATGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATAACGTATGACt                                                   >  1:3254017/1‑62 (MQ=255) aaaTAATGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATAACGTATGACt                                                   >  1:4624583/1‑62 (MQ=255) aaaTAATGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATAACGTATGACt                                                   >  1:485568/1‑62 (MQ=255) aaaTAATGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATAACGTATGACt                                                   >  1:1813176/1‑62 (MQ=255)          tgaCAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATAACGTATGACTTCTGaaaa                                           >  1:3569121/1‑61 (MQ=255)          tgaCAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATAACGTATGACTTCTGaaaa                                           >  1:1795326/1‑61 (MQ=255)          tgaCAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATAACGTATGACTTCTGaaaa                                           >  1:509050/1‑61 (MQ=255)          tgaCAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATAACGTATGACTTCTGaaaa                                           >  1:1548493/1‑61 (MQ=255)          tgaCAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATAACGTATGACTTCTGaaaa                                           >  1:1405535/1‑61 (MQ=255)          tgaCAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATAACGTATGACTTCTGaaaa                                           >  1:409139/1‑61 (MQ=255)          tgaCAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATAACGTATGACTTCTGaaaa                                           >  1:3183132/1‑61 (MQ=255)            aCAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATAACGTATGACTTCTGAAAACAt                                        <  1:4006185/62‑1 (MQ=255)            aCAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATAACGTATGACTTCTGAAAACAt                                        <  1:92357/62‑1 (MQ=255)                         tCACCGTAGGGCCGAGGTCGAGAATATAACGTATGACTTCTGaa                                             <  1:2591679/44‑1 (MQ=255)                         tCACCGTAGGGCCGAGGTCGAGAATATAACGTATGACTTCTGaa                                             <  1:2789341/44‑1 (MQ=255)                         tCACCGTAGGGCCGAGGTCGAGAATATAACGTATGACTTCTGaa                                             <  1:2923247/44‑1 (MQ=255)                         tCACCGTAGGGCCGAGGTCGAGAATATAACGTATGACTTCTGaa                                             <  1:3850004/44‑1 (MQ=255)                         tCACCGTAGGGCCGAGGTCGAGAATATAACGTATGACTTCTGaa                                             <  1:1965445/44‑1 (MQ=255)                         tCACCGTAGGGCCGAGGTCGAGAATATAACGTATGACTTCTGaa                                             <  1:1211153/44‑1 (MQ=255)                         tCACCGTAGGGCCGAGGTCGAGAATATAACGTATGACTTCTGaa                                             <  1:1387747/44‑1 (MQ=255)                         tCACCGTAGGGCCGAGGTCGAGAATATAACGTATGACTTCTGaa                                             <  1:1204762/44‑1 (MQ=255)                         tCACCGTAGGGCCGAGGTCGAGAATATAACGTATGACTTCTGAAAACATAGGTCACCCCt                             >  1:3873649/1‑60 (MQ=255)                         tCACCGTAGGGCCGAGGTCGAGAATATAACGTATGACTTCTGAAAACATAGGTCACCCCt                             >  1:1169475/1‑60 (MQ=255)                         tCACCGTAGGGCCGAGGTCGAGAATATAACGTATGACTTCTGAAAACATAGGTCACCCCt                             >  1:2244978/1‑60 (MQ=255)                         tCACAGTAGGGCCGAGGTCGAGAATATAACGTATGACTTCTGAAAACATAGGTCACCCCt                             >  1:3447278/1‑60 (MQ=255)                           aCCGTAGGGCCGAGGTCGAGAATATAACGTATGACTTCTGAAAACATAGGTCAcc                                <  1:3471744/55‑1 (MQ=255)                           aCCGTAGGGCCGAGGTCGAGAATATAACGTATGACTTCTGAAAACATAGGTCAcc                                <  1:1713501/55‑1 (MQ=255)                           aCCGTAGGGCCGAGGTCGAGAATATAACGTATGACTTCTGAAAACATAGGTCAcc                                <  1:1449048/55‑1 (MQ=255)                           aCCGTAGGGCCGAGGTCGAGAATATAACGTATGACTTCTGAAAACATAGGTCAcc                                <  1:975800/55‑1 (MQ=255)                           aCCGTAGGGCCGAGGTCGAGAATATAACGTATGACTTCTGAAAACATAGGTCAcc                                <  1:2266914/55‑1 (MQ=255)                                aGGGCCGAGGTCGAGAATATAACGTATGACTTCTGAAAACATAGGTCACCCCTGTAAGATa                     >  1:817258/1‑61 (MQ=255)                                aGGGCCGAGGTCGAGAATATAACGTATGACTTCTGAAAACATAGGTCACCCCTGTAAGATa                     >  1:2274146/1‑61 (MQ=255)                                aGGGCCGAGGTCGAGAATATAACGTATGACTTCTGAAAACATAGGTCACCCCTGTAAGATa                     >  1:3190543/1‑61 (MQ=255)                                aGGGCCGAGGTCGAGAATATAACGTATGACTTCTGAAAACATAGGTCACCCCTGTAAGATa                     >  1:1405440/1‑61 (MQ=255)                                aGGGCCGAGGTCGAGAATATAACGTATGACTTCTGAAAACATAGGTCACCCCTGTAAGATa                     >  1:43288/1‑61 (MQ=255)                                  ggCCGAGGTCGAGAATATAACGTATGACTTCTGaaaa                                           >  1:3843520/1‑37 (MQ=255)                                  ggCCGAGGTCGAGAATATAACGTATGACTTCTGaaaa                                           >  1:328543/1‑37 (MQ=255)                                  ggCCGAGGTCGAGAATATAACGTATGACTTCTGaaaa                                           >  1:3364439/1‑37 (MQ=255)                                  ggCCGAGGTCGAGAATATAACGTATGACTTCTGaaaa                                           >  1:3858870/1‑37 (MQ=255)                                  ggCCGAGGTCGAGAATATAACGTATGACTTCTGaaaa                                           >  1:3696286/1‑37 (MQ=255)                                  ggCCGAGGTCGAGAATATAACGTATGACTTCTGaaaa                                           >  1:1612608/1‑37 (MQ=255)                                  ggCCGAGGTCGAGAATATAACGTATGACTTCTGaaaa                                           >  1:394314/1‑37 (MQ=255)                                  ggCCGAGGTCGAGAATATAACGTATGACTTCTGaaaa                                           >  1:2685881/1‑37 (MQ=255)                                  ggCCGAGGTCGAGAATATAACGTATGACTTCTGaaaa                                           >  1:638460/1‑37 (MQ=255)                                  ggCCGAGGTCGAGAATATAACGTATGACTTCTGaaaa                                           >  1:2453682/1‑37 (MQ=255)                                  ggCCGAGGTCGAGAATATAACGTATGACTTCTGaaaa                                           >  1:2253310/1‑37 (MQ=255)                                  ggCCGAGGTCGAGAATATAACGTATGACTTCTGaaaa                                           >  1:488654/1‑37 (MQ=255)                                  ggCCGAGGTCGAGAATATAACGTATGACTTCTGaaaa                                           >  1:1687241/1‑37 (MQ=255)                                  ggCCGAGGTCGAGAATATAACGTATGACTTCTGaaaa                                           >  1:4585775/1‑37 (MQ=255)                                  ggCCGAGGTCGAGAATATAACGTATGACTTCTGAAAACa                                         >  1:3855057/1‑39 (MQ=255)                                  ggCCGAGGTCGAGAATATAACGTATGACTTCTGAAAACa                                         >  1:4064805/1‑39 (MQ=255)                                  ggCCGAGGTCGAGAATATAACGTATGACTTCTGAAAACa                                         >  1:4144110/1‑39 (MQ=255)                                  ggCCGAGGTCGAGAATATAACGTATGACTTCTGAAAACa                                         >  1:1019235/1‑39 (MQ=255)                                  ggCCGAGGTCGAGAATATAACGTATGACTTCTGAAAACa                                         >  1:4388141/1‑39 (MQ=255)                                  ggCCGAGGTCGAGAATATAACGTATGACTTCTGAAAACa                                         >  1:983236/1‑39 (MQ=255)                                  ggCCGAGGTCGAGAATATAACGTATGACTTCTGAAAACa                                         >  1:1009255/1‑39 (MQ=255)                                  ggCCGAGGTCGAGAATATAACGTATGACTTCTGAAAACa                                         >  1:1349179/1‑39 (MQ=255)                                  ggCCGAGGTCGAGAATATAACGTATGACTTCTGAAAACa                                         >  1:1661018/1‑39 (MQ=255)                                  ggCCGAGGTCGAGAATATAACGTATGACTTCTGAAAACa                                         >  1:1768118/1‑39 (MQ=255)                                  ggCCGAGGTCGAGAATATAACGTATGACTTCTGAAAACa                                         >  1:2210405/1‑39 (MQ=255)                                  ggCCGAGGTCGAGAATATAACGTATGACTTCTGAAAACATAGGTCACCCCTGTaa                         >  1:3503941/1‑55 (MQ=255)                                  ggCCGAGGTCGAGAATATAACGTATGACTTCTGAAAACATAGGTCACCCCTGTaa                         >  1:46139/1‑55 (MQ=255)                                  ggCCGAGGTCGAGAATATAACGTATGACTTCTGAAAACATAGGTCACCCCTGTaa                         >  1:4023292/1‑55 (MQ=255)                                  ggCCGAGGTCGAGAATATAACGTATGACTTCTGAAAACATAGGTCACCCCTGTaa                         >  1:3903329/1‑55 (MQ=255)                                  ggCCGAGGTCGAGAATATAACGTATGACTTCTGAAAACATAGGTCACCCCTGTaa                         >  1:2715499/1‑55 (MQ=255)                                  ggCCGAGGTCGAGAATATAACGTATGACTTCTGAAAACATAGGTCACCCCTGTaa                         >  1:3417158/1‑55 (MQ=255)                                  ggCCGAGGTCGAGAATATAACGTATGACTTCTGAAAACATAGGTCACCCCTGTaa                         >  1:3803388/1‑55 (MQ=255)                                  ggCCGAGGTCGAGAATATAACGTATGACTTCTGAAAACATAGGTCACCCCTGTAAGATa                     >  1:349604/1‑59 (MQ=255)                                  ggCCGAGGTCGAGAATATAACGTATGACTTCTGAAAACATAGGTCACCCCTGTAAGATa                     >  1:3489263/1‑59 (MQ=255)                                  ggCCGAGGTCGAGAATATAACGTATGACTTCTGAAAACATAGGTCACCCCTGTAAGATa                     >  1:832883/1‑59 (MQ=255)                                  ggCCGAGGTCGAGAATATAACGTATGACTTCTGAAAACATAGGTCACCCCTGTAAGATa                     >  1:1684785/1‑59 (MQ=255)                                  ggCCGAGGTCGAGAATATAACGTATGACTTCTGAAAACATAGGTCACCCCTGTAAGATa                     >  1:4517641/1‑59 (MQ=255)                                  ggCCGAGGTCGAGAATATAACGTATGACTTCTGAAAACATAGGTCACCCCTGTAAGATa                     >  1:1895808/1‑59 (MQ=255)                                  ggCCGAGGTCGAGAATATAACGTATGACTTCTGAAAACATAGGTCACCCCTGTAAGATa                     >  1:2185796/1‑59 (MQ=255)                                  ggCCGAGGTCGAGAATATAACGTATGACTTCTGAAAACATAGGTCACCCCTGTAAGATa                     >  1:1024528/1‑59 (MQ=255)                                  ggCCGAGGTCGAGAATATAACGTATGACTTCTGAAAACATAGGTCACCCCTGTAAGATa                     >  1:2835075/1‑59 (MQ=255)                                  ggCCGAGGTCGAGAATATAACGTATGACTTCTGAAAACATAGGTCACACATGTaa                         >  1:1583342/1‑55 (MQ=255)                                         gTCGAGAATATAACGTATGACTTCTGAAAACATATGTAACCCCTGTAAGATAGTCAGAattt           <  1:1290304/62‑1 (MQ=255)                                         gTCGAGAATATAACGTATGACTTCTGAAAACATAGGTCACCCCTGTAAGATAGTCAGAattt           <  1:2116868/62‑1 (MQ=255)                                          tCGAGAATATAACGTATGACTTCTGAAAACATAGGTCACCCCTGTAAGATAGTCAGAatttt          <  1:4215160/62‑1 (MQ=255)                                          tCGAGAATATAACGTATGACTTCTGAAAACATAGGTCACCCCTGTAAGATAGTCAGAatttt          <  1:4616434/62‑1 (MQ=255)                                          tCGAGAATATAACGTATGACTTCTGAAAACATAGGTCACCCCTGTAAGATAGTCAGAatttt          <  1:2543072/62‑1 (MQ=255)                                          tCGAGAATATAACGTATGACTTCTGAAAACATAGGTCACCCCTGTAAGATAGTCAGAatttt          <  1:2693111/62‑1 (MQ=255)                                          tCGAGAATATAACGTATGACTTCTGAAAACATAGGTCACCCCTGTAAGATAGTCAGAatttt          <  1:1397483/62‑1 (MQ=255)                                          tCGAGAATATAACGTATGACTTCTGAAAACATAGGTCACCCCTGTAAGATAGTCAGAattt           <  1:926967/61‑1 (MQ=255)                                            gagaATATAACGTATGACTTCTGAAAACATAGGTCACCCCTGTAAGATAGTCAGAatttt          <  1:1163977/60‑1 (MQ=255)                                            gagaATATAACGTATGACTTCTGAAAACATAGGTCACCCCTGTAAGATAGTCAGAatttt          <  1:1576347/60‑1 (MQ=255)                                             agaATATAACGTATGACTTCTGAAAACATAGGTCACCCCTGTAAGATAGTCAGAattttatt       >  1:1612512/1‑62 (MQ=255)                                             agaATATAACGTATGACTTCTGAAAACATAGGTCACCCCTGTAAGATAGTCAGAattttat        >  1:853665/1‑61 (MQ=255)                                                  ataACGTATGACTTCTGAAAACATAGGTCACCCCTGTAAGATAGTCAGAATTTTATTTTGTa  <  1:1606776/62‑1 (MQ=255)                                                  ataACGTATGACTTCTGAAAACATAGGTCACCCCTGTAAGATAGTCAGAATTTTATTTTGTa  <  1:4684203/62‑1 (MQ=255)                                                  ataACGTATGACTTCTGAAAACATAGGTCACCCCTGTAAGATAGTCAGAATTTTATTTTGTa  <  1:4681803/62‑1 (MQ=255)                                                  ataACGTATGACTTCTGAAAACATAGGTCACCCCTGTAAGATAGTCAGAATTTTATTTTGTa  <  1:4681357/62‑1 (MQ=255)                                                  ataACGTATGACTTCTGAAAACATAGGTCACCCCTGTAAGATAGTCAGAATTTTATTTTGTa  <  1:3191071/62‑1 (MQ=255)                                                  ataACGTATGACTTCTGAAAACATAGGTCACCCCTGTAAGATAGTCAGAATTTTATTTTGTa  <  1:2975173/62‑1 (MQ=255)                                                  ataACGTATGACTTCTGAAAACATAGGTCACCCCTGTAAGATAGTCAGAATTTTATTTTGTa  <  1:2829873/62‑1 (MQ=255)                                                        |                                                        AAATAATGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATAACGCATGACTTCTGAAAACATAGGTCACCCCTGTAAGATAGTCAGAATTTTATTTTGTA  >  W3110S.gb/2176341‑2176451

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