Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,301,992	A→T	N271Y (AAT→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,992	0	A	T	100.0%	85.6 / NA	28	N271Y (AAT→TAT)	oppA	oligopeptide ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (13/15);  total (13/15)

GCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGC  >  NC_000913/1301760‑1302184                                                                                                                                                                                                                                         |                                                                                                                                                                                                 gCCAAATGCCGCTAGCGAGAAATTCGGCGAAAAATTGACCCAGCCTGGTAATATCGTCACCAACGTTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACTCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTTCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTCCCGTAGTGGTGAAATCTACATGACGTATAACAGCATGCCGATCg                                                                                                                                                                                <  2:211513/251‑1 (MQ=255)      aaGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTg                                                                                                                                                                           <  2:207309/251‑1 (MQ=255)                 gagaAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCt                                                                                                                                                                >  1:5401/1‑251 (MQ=255)                         cGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGaaaaaa                                                                                                                                                         <  2:151246/250‑1 (MQ=255)                                aaaTGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATccc                                                                                                                                                 >  1:87647/1‑251 (MQ=255)                                       cccAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGaa                                                                                                                                          <  2:158265/251‑1 (MQ=255)                                                    tatCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATc                                                                                                                             >  2:329915/1‑251 (MQ=255)                                                       cGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCAt                                                                                                                          >  2:126479/1‑251 (MQ=255)                                                        gTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATc                                                                                                                                                                                 >  2:117205/1‑195 (MQ=255)                                                        gTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATc                                                                                                                                                                                 <  1:117205/195‑1 (MQ=255)                                                           aCCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCt                                                                                                                      <  2:326783/251‑1 (MQ=255)                                                                         taCCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTAtt                                                                                                        >  2:301167/1‑251 (MQ=255)                                                                                                        aTCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATAAACAACCAGAAACCGCCATTCaa                                                                         <  1:126479/251‑1 (MQ=255)                                                                                                                 gAACTCAGCCCGCCCTACTGGAACAACGCGAAAACCGTTATTACCCAGGTAACCTATTTGCCTATTGCTTCGGAAGTTCCCGATTTCAACCGCTACCGTAGTGGTAAAATCGACATGACGTATAACAGCATGTCGATCAAATTGTTCCAGAAGCTGAAAAAACAGATCCCGGACGAAGTTCGCGTTGATCCATACCTGTGCACTTCCTTTTACGAAATTAACAACCATAAACCGCCATTCAACGAtgtgc                                                                 <  2:15323/250‑1 (MQ=255)                                                                                                                          ccGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCAAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACTTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTAc                                                       <  1:364256/251‑1 (MQ=255)                                                                                                                            gACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCg                                                     >  1:20193/1‑251 (MQ=255)                                                                                                                                 aCTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCCCTTACTATTACGAAATTAACACCCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTg                                                <  1:222509/251‑1 (MQ=255)                                                                                                                                    ggAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGaa                                              >  1:362404/1‑250 (MQ=255)                                                                                                                                             cgAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGTTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCtg                                                                                                                     <  2:352347/170‑1 (MQ=255)                                                                                                                                             cgAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGTTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCtg                                                                                                                     >  1:352347/1‑170 (MQ=255)                                                                                                                                             cgAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATg                                    >  1:124654/1‑251 (MQ=255)                                                                                                                                                    cGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACcgcg                             >  2:239112/1‑251 (MQ=255)                                                                                                                                                           aaCCAGGTAACCTATTTGCCTATGGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACGTGTGCACTTACTATTCCGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGGGTGCGTACCGCTCTGAAACTAGGGATGGACCGCGATATCAt                      <  2:171214/251‑1 (MQ=255)                                                                                                                                                                ggTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTta                 <  2:302130/251‑1 (MQ=255)                                                                                                                                                                          ttGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGa        >  2:342580/1‑250 (MQ=255)                                                                                                                                                                            gCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGaaa      <  2:295370/250‑1 (MQ=255)                                                                                                                                                                             ccTATTGCTTCTTAAGTTACGTATGTCCACCGCTACCGTTGTGGAGAAATCGCTATGACGTATAACAGCATGCCGATCCAATTGTTCCAGATCCTGAAAAAAGAGATCCCGGCCGAATTTCACGTTGACCCATACCTGTTCACTTACAATTGCGAAATTAACAACCAGAAACCGCCATTCAACGTGGTGCGTGTTCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTCAAAgc    <  2:110519/251‑1 (MQ=255)                                                                                                                                                                                 tggCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGGGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAGTTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGGGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGAATGGACCGCGATATCATTGTTCATAAAGTGAAAgagc  >  2:195713/3‑246 (MQ=255)                                                                                                                                                                                                                                         |                                                                                                                                                                                                 GCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGC  >  NC_000913/1301760‑1302184

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