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	96.7%	95.9 / NA	30	N271Y (AAT→TAT)	oppA	oligopeptide ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (15/14);  minor base G (0/1);  total (15/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

TTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACC‑GTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTATGGTTACACTCCACCGT  >  NC_000913/1301751‑1302220                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                     ttCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACC‑GTAGTGGTGAAATCGACATGACGTATAACAGc                                                                                                                                                                                                                              >  1:296206/1‑250 (MQ=255)   cACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGGGCCTATACCTTAAAAGATTGGGTGGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACTTATTTGGCTATTGCTTCTGACGCTACCGCTGTCAACCGC‑ACCAGTATGGGTGCAATCGAATTGACGTATACCAGCATg                                                                                                                                                                                                                           >  2:3806/1‑251 (MQ=255)    aCCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACC‑GTAGTGGTGAAATCGACATGACGTATAACAGCATGc                                                                                                                                                                                                                          <  2:294023/251‑1 (MQ=255)                       aTCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACC‑GTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAg                                                                                                                                                                                                         <  1:186140/249‑1 (MQ=255)                       aTCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACC‑GTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAg                                                                                                                                                                                                         <  1:77388/249‑1 (MQ=255)                       aTCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACC‑GTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGa                                                                                                                                                                                                        >  1:239739/1‑250 (MQ=255)                         cGAGAAATTCGGCGACAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACC‑GTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGc                                                                                                                                                                                                     >  1:302817/1‑251 (MQ=255)                         cGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACC‑GTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAg                                                                                                                                                                                                      <  1:241908/250‑1 (MQ=255)                          gagaAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTGACCGATGTCAACCGCTACC‑GTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCt                                                                                                                                                                                                    <  1:202958/251‑1 (MQ=255)                              aaTTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCATCGGTGCCTATACCTTAAAAGATTGGGGCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTAGTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACC‑GTAGTGGTGAAATCGACAGGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGaaa                                                                                                                                                                                                <  2:285323/251‑1 (MQ=255)                              aaTTCAGCGAAAAATAGCCACAGCCTGGTAATATCGTCACCAACGGTGACTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGAACTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACC‑GTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGaa                                                                                                                                                                                                 >  2:335276/1‑250 (MQ=255)                               aTTCGGCGAAAAATGGACCCAGCCTGGTAATAGCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATGAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACC‑GTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGGTCCAGAATCTGaaaa                                                                                                                                                                                               <  2:54159/251‑1 (MQ=255)                                      gAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACC‑GTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAACAAag                                                                                                                                                                                            >  1:79802/1‑247 (MQ=255)                                             ggACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACC‑GTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGa                                                                                                                                                                                  <  1:399269/250‑1 (MQ=255)                                                 ccAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACC‑GTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAg                                                                                                                                                                             <  2:317498/251‑1 (MQ=255)                                                              atCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGCACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACC‑GTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATcc                                                                                                                                                                >  1:356889/1‑251 (MQ=255)                                                                     ccAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACC‑GTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCTAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAATTTCACGTTTAGCCATACCtg                                                                                                                                                         <  2:348943/251‑1 (MQ=255)                                                                                    ccTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACC‑GTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTAc                                                                                                                                          <  2:472/251‑1 (MQ=255)                                                                                         aaaGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACCGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACC‑GTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAAt                                                                                                                                     >  1:170001/1‑251 (MQ=255)                                                                                                gggTCGTAAAAGTACCACTCGTTCTTGCACGCAGCCCGAACTCCTGGAACAACGCGAAAACCGTTATTCACCAGGTCACGTATTTGCCTATTGCTTCTTAAGTTAACGATGTCAACCGCTACC‑GTAGAGGTGAAACCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGATGACAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTCTTACGAAATTacaaac                                                                                                                              >  2:150093/1‑246 (MQ=255)                                                                                                       aatcgaacgaaTCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCGATTGCTTCTGAAGTTACCGATGTCAACCGCTACC‑GTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTTACGTTGAGCCATACCTGTGGGCTTACTCTGACGAAATTAACATCCAGaaa                                                                                                                        <  1:118412/247‑1 (MQ=255)                                                                                                                           aaCGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCCGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACC‑GTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTACAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGAtgtgcgt                                                                                                   >  1:58843/1‑251 (MQ=255)                                                                                                                            aCGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCGATTTGCCTATTGCTTCTGAAGTTACCGATGACAACCGCTACC‑GTAGGGGTGAAATCGACATGACGGATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGTTCCCGTACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGAtgtgcgtg                                                                                                  <  1:200279/251‑1 (MQ=255)                                                                                                                                                                     aCCAGGTAACATATTTGCCTATTGCTTCTGCAGTTACTGATGTCAACCGCAACC‑GTCGTGGTGAAATCGACATGACGTATAACAGCCTGCCGATCGCATTGTTCCAGAAGCTGAAAAAAGAGATCCAGGACGAAGATCACGTTGATCCATACCTGTGCAATTCATATTACGACATTAAAAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCCCACTCAACATAGGTATCGACCACGATATCAtt                                                         >  2:60674/1‑251 (MQ=255)                                                                                                                                                                      ccAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACC‑GTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGACCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTg                                                        <  2:296206/251‑1 (MQ=255)                                                                                                                                                                                 aTTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACC‑GTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTg                                             >  1:186614/1‑251 (MQ=255)                                                                                                                                                                                      ccTATTGCTTCTGAAGTTACCGATGTCAACCGCTACC‑GTAGTGGTGAAATCGACCTGACGTATAACAGCATGCCGATCGCATTGTTCCAGAAGCTGACAAAAGAGATCCCGGACGAAGTTCCCGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGCAACTAGGTATGGACCGCGATATCATTGTTAATCAAGTTAAAgc                                        >  2:269822/1‑251 (MQ=255)                                                                                                                                                                                           tGCTTCTGAAGTTACCGATGTCAACCGCTACC‑GTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAgg                                   >  1:27741/1‑251 (MQ=255)                                                                                                                                                                                                                        taccgaggtggtgAAACCGACATGACGTAAAACGGCATGCCAATCGAATTGTCCCAAAAGCTAAAAAAAGAGTACCCGGACAAAGTCCAGGTTGATCCATACCTGTGCATTTATTATTAAGAAATAACCAACCAAAAACCGCCATTCACCGATGTTCGTGTGCGTACCGCACTGAAACTAGGTAGGGCCCGGGATATCAGTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTATGGTGACACTc        <  2:186294/243‑1 (MQ=255)                                                                                                                                                                                                                              aGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTATGGTTACACTCCACCGt  >  1:49454/1‑250 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                     TTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACC‑GTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTATGGTTACACTCCACCGT  >  NC_000913/1301751‑1302220

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