Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	1,303,631	C→T	Q246Y (CAG→TAT)	oppA →	oligopeptide transporter subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	1,303,631	0	C	T	100.0%	200.6 / NA	68	Q246Y (CAG→TAT)	oppA	oligopeptide transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (22/46);  total (22/46)

TAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGT  >  W3110S.gb/1303566‑1303687                                                                  |                                                         tAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTa                                                     <  1:579039/71‑1 (MQ=255) tAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTa                                                     <  1:1714405/71‑1 (MQ=255)   aacgaacgAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAAc                                                   >  1:2949300/1‑71 (MQ=255)   aacgaacgAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAAc                                                   >  1:3313664/1‑71 (MQ=255)   aacgaacgAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAAc                                                   >  1:4220221/1‑71 (MQ=255)    acgaacgaATCGTTCTTGAACGCAGGCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAAcc                                                  <  1:96869/71‑1 (MQ=255)    acgaacgaATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAAcc                                                  <  1:1816859/71‑1 (MQ=255)      gaacgaaTCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCTAAAACCGTTATTAACTATGTAACCTa                                                <  1:218479/71‑1 (MQ=255)      gaacgaaTCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTa                                                <  1:4339830/71‑1 (MQ=255)      gaacgaaTCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTa                                                <  1:313087/71‑1 (MQ=255)      gaacgaaTCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTa                                                <  1:4548665/71‑1 (MQ=255)      gaacgaaTCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTa                                                >  1:3942965/1‑71 (MQ=255)      gaacgaaTCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTa                                                <  1:4594595/71‑1 (MQ=255)      gaacgaaTCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTa                                                <  1:3140804/71‑1 (MQ=255)      gaacgaaTCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTa                                                <  1:89776/71‑1 (MQ=255)      gaacgaaTCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTa                                                <  1:2758332/71‑1 (MQ=255)      gaacgaaTCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTa                                                <  1:2454370/71‑1 (MQ=255)      gaacgaaTCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTa                                                <  1:2362684/71‑1 (MQ=255)      gaacgaaTCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAAATATGTAACCTa                                                <  1:2938012/71‑1 (MQ=255)         ctaaTCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTAt                                               <  1:2783373/67‑1 (MQ=255)           aaTCGTTC‑TCAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGcc                                          <  1:3259904/71‑1 (MQ=255)            aTCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGcc                                          <  1:3507448/71‑1 (MQ=255)            aTCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGcc                                          <  1:4521602/71‑1 (MQ=255)            aTCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGcc                                          <  1:3939997/71‑1 (MQ=255)            aTCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGcc                                          <  1:2353667/71‑1 (MQ=255)            aTCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGcc                                          <  1:293011/71‑1 (MQ=255)             tCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCt                                         <  1:2654723/71‑1 (MQ=255)             tCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCt                                         <  1:1899424/71‑1 (MQ=255)             tCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCt                                         <  1:4211105/71‑1 (MQ=255)             tCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCt                                         <  1:2205424/71‑1 (MQ=255)               gTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTAtt                                              >  1:2605016/1‑64 (MQ=255)               gTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTAt                                       >  1:4361204/1‑71 (MQ=255)               gTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATATAACCTAtt                                              >  1:758159/1‑64 (MQ=255)                  cTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGc                                    >  1:1775636/1‑71 (MQ=255)                     gAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGc                                    <  1:1311591/68‑1 (MQ=255)                     gAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGc                                    <  1:4061581/68‑1 (MQ=255)                     gAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCt                                   <  1:330502/69‑1 (MQ=255)                        cGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGa                              <  1:3915287/71‑1 (MQ=255)                            gCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGt                           >  1:3528490/1‑70 (MQ=255)                              ccGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTAc                        >  1:2357513/1‑71 (MQ=255)                              ccGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTAc                        >  1:4151010/1‑71 (MQ=255)                              ccGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTAc                        >  1:2339589/1‑71 (MQ=255)                               cGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTAcc                       <  1:2141967/71‑1 (MQ=255)                               cGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTAcc                       <  1:4627353/71‑1 (MQ=255)                               cGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTAcc                       <  1:1082670/71‑1 (MQ=255)                               cGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTAcc                       <  1:3754187/71‑1 (MQ=255)                               cGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTAcc                       <  1:3718599/71‑1 (MQ=255)                               cGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTAcc                       <  1:2399733/71‑1 (MQ=255)                               cGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTAcc                       <  1:3140555/71‑1 (MQ=255)                               cGACCTACTGGAACAACGCGAAAACCGTTACTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTAcc                       <  1:1482630/71‑1 (MQ=255)                                 aCCTACTTGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGa                     <  1:3498553/71‑1 (MQ=255)                                      cTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGt                  >  1:2507970/1‑69 (MQ=255)                                      cTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTc                 >  1:2809164/1‑70 (MQ=255)                                      cTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTc                 >  1:1478980/1‑70 (MQ=255)                                      cTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTc                 >  1:300105/1‑70 (MQ=255)                                      cTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTc                 >  1:2678875/1‑70 (MQ=255)                                      cTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTc                 >  1:2529195/1‑70 (MQ=255)                                      cTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTc                 >  1:3814326/1‑70 (MQ=255)                                      cTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCa                <  1:763548/71‑1 (MQ=255)                                        ggAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTTCCTATTGCTTCTGAAGTTACCGATGTCAAc              >  1:1046094/1‑71 (MQ=255)                                           acaacGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGa                     >  1:863010/1‑61 (MQ=255)                                           acaacGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGa                     >  1:259265/1‑61 (MQ=255)                                                 cgAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGt     <  1:2723949/71‑1 (MQ=255)                                                 cgAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGt     <  1:3607408/71‑1 (MQ=255)                                                    aaaCCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAgt  <  1:2486592/71‑1 (MQ=255)                                                     aaCCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAgt  <  1:3154808/70‑1 (MQ=255)                                                       ccGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCg                      <  1:1727236/48‑1 (MQ=39)                                                       ccGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCg                      <  1:92911/48‑1 (MQ=39)                                                                  |                                                         TAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGT  >  W3110S.gb/1303566‑1303687

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