Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,301,919	G→T	Q246Y (CAG→TAT)	oppA →	oligopeptide ABC transporter periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,301,919	0	G	T	100.0%	109.3 / NA	33	Q246Y (CAG→TAT)	oppA	oligopeptide ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (15/18);  total (15/18)

ATTGATGATCACACATTAGAAGTCACCTTAAGTGAACCCGTTCCGTACTTCTATAAATTACTTGTTCACCCATCAACTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATAT  >  NC_000913/1301674‑1302161                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   aTTGATGATCACACATTAGAAGTCACCTTAAGTGAACCCGTTCCGTACTTCTATAAATTACTTGTTCACCCATCAACTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAAc                                                                                                                                                                                                                                               <  2:407600/251‑1 (MQ=255)             ccaTTAGAAGTCACCTTAAGTGAACCCGTTCCGTACTTCTATAAATTACTTGTTCACCCATCAACTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTAt                                                                                                                                                                                                                                   <  1:9418/250‑1 (MQ=255)                     aGTCACCTTAAGTGAACCCGTTCCGTACTTCTATAAATTACTTGTTCACCCATCAACTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTg                                                                                                                                                                                                                           >  1:29827/1‑251 (MQ=255)                            ttAAGTGAACCCGTTCCGTACTTCTATAAATTACTTGTTCACCCATCAACTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTAc                                                                                                                                                                                                                    >  1:143083/1‑251 (MQ=255)                              aaGTGAACCCGTTCCGTACTTCTATAAATTACTTGTTCACCCATCAACTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGACCAACGCGAAAACCGTTATTAACTATGTAACCTATTGGCTTATTgatttgaaggttacgg                                                                                                                                                                                                                  >  1:213745/1‑236 (MQ=255)                                gTGAACCCGTTCCGTACTTCTATACATAACTTGTTCACCCATCAACTTCACCGGTGCCAAAAGCCGCTATCGAGACATTCGGCGAAAAATGTACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGATCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACCCTAAAACCGTTATTAACTATGTAACCTCTTTGCCTATTGCTTCTGAAGTTACCGAt                                                                                                                                                                                                                >  1:208455/1‑251 (MQ=255)                                  gAACCCGTTCCGAACTTCTATAAATTACTTGTTCACACATCAACTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGCACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTACCTATGTAACCTATTTGCCTATTGCTTCTGAAGTTCCCGCTGt                                                                                                                                                                                                              >  2:460854/1‑251 (MQ=255)                                         ttCCGTACTTCTATAAATTACTTGTTCACCCATCAACTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGc                                                                                                                                                                                                       >  2:37730/1‑251 (MQ=255)                                                     ataAATTACTTGTTCACCCATCAACTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAgt                                                                                                                                                                                              <  2:29827/248‑1 (MQ=255)                                                                   ccaccATCAACCTCCCCGGGGCCAAAAACCCCTAACGAGAAATTACGCGAAAAATGGACCCCACCTGGTAATATCGGCACCAACGGTGCCTATTCCCTAAAAGATTGGGGCGGAAACGTACCGATCGTTCCTGAACGCATGCCGACCTACGTGAACAACCCGAAAACCCTTAGTAACTATGGAACCTATTGGCCTAGTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGGGAAATCGACATGAc                                                                                                                                                                             <  2:213745/248‑1 (MQ=255)                                                                          cAACTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGCCCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGAGGAATAAc                                                                                                                                                                      >  2:80149/1‑251 (MQ=255)                                                                                     gTGCCAAAAGCCGCTCTCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGAt                                                                                                                                                           <  2:155288/251‑1 (MQ=255)                                                                                           aaaGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAAtt                                                                                                                                                     >  1:88100/1‑251 (MQ=255)                                                                                             aGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGt                                                                                                                                                   <  2:431066/251‑1 (MQ=255)                                                                                                      cGAGAAATTCGGCGAGAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGc                                                                                                                                          <  1:187566/251‑1 (MQ=255)                                                                                                        agaAATTCGGCGAAAAATGGACCCAGCCTGGTAATAACGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGCCAGGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTg                                                                                                                                        <  2:26875/251‑1 (MQ=255)                                                                                                           aaTTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGaaa                                                                                                                                     <  1:37730/251‑1 (MQ=255)                                                                                                                        aTGGACCCAGCCTGGTAATATCGTCACCATCGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGTCCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTAGTGCGTCTGAAGTTACCGATGTCAACCGCTACCGTAGTTGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGATATCCCgg                                                                                                                        <  1:373094/251‑1 (MQ=255)                                                                                                                            aCCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACATTAAACGATCGGGTCGTAAACCAACGAATCGCTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAACACCGTTATTAACTATGTAACCTCTTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGCCGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACg                                                                                                                     >  1:68598/1‑250 (MQ=255)                                                                                                                                      gTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGtt                                                                                                          <  1:135754/251‑1 (MQ=255)                                                                                                                                      gTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGtt                                                                                                          <  1:145266/251‑1 (MQ=255)                                                                                                                                          tatCGTCCCCAACGGAGCCTATACCTTAAACGTTTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGCAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGCCATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATc                                                                                                      >  2:61747/1‑251 (MQ=255)                                                                                                                                            gcGTCACCAACGGTGCCTGTACCTTAAAAAATGGGGTCGGAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACGGGAACAACGCGAAAACCGTTATTAACTATGGAACCTATTTGCCTATTGCTGCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACAGGACGAATAACAGCATGCCGATCGAAGGGTTCCATGAGCTGAAAAAAGAGATCCCCTACGAATTTCGCGTTGATCCa                                                                                                    <  2:244041/250‑1 (MQ=255)                                                                                                                                                cACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATAc                                                                                                 <  1:63103/250‑1 (MQ=255)                                                                                                                                                   atcaaaaaaatattctctAAAAGATTGGGTCGTAAAGGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCtgt                                                                                             <  2:103825/234‑1 (MQ=255)                                                                                                                                                        gTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACTCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACt                                                                                        <  1:460815/251‑1 (MQ=255)                                                                                                                                                                     aaaaGATTGGGTCGTAAACGAACGAATCGATCTTGACCGCAGCCCGACCTACTGGAACAACGCGAAACCCGTTATTCACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTCCCGAAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGaa                                                                            >  2:285256/1‑250 (MQ=255)                                                                                                                                                                       aaGATGGGGGCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCGTCTGAAGTTACCGATGTCAACCGCTACCGTAGTTGTGAAATCGACATGACGAATAACAGCATGCCGATCGAGTGGTGCCGTGAGCTGAAGAAAGAGGACCCGGACGAGGGTTACGAGGAGCCATACCTGGGCACTTGCTATGACGAAAtt                                                                         <  2:506448/251‑1 (MQ=255)                                                                                                                                                                        aGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGGCGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTa                                                                        <  1:61747/251‑1 (MQ=255)                                                                                                                                                                                    aaaCGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAAc                                                            >  1:189159/1‑251 (MQ=255)                                                                                                                                                                                                                           gAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAAc                     >  1:274391/1‑251 (MQ=255)                                                                                                                                                                                                                            aacaacGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTAACGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGCTCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGGGTGCGTACCGCACTGAAACt                    >  1:341175/1‑251 (MQ=255)                                                                                                                                                                                                                                              attaACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGatat  >  1:186481/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   ATTGATGATCACACATTAGAAGTCACCTTAAGTGAACCCGTTCCGTACTTCTATAAATTACTTGTTCACCCATCAACTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATAT  >  NC_000913/1301674‑1302161

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