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.6 / ‑3.3	30	N271Y (AAT→TAT)	oppA	oligopeptide ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base T (15/14);  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 = 2.11e-01

CACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTATGGTTACACTCCACCGTATACTGATGGCGCAAAAT  >  NC_000913/1301753‑1302238                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                       cACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCAt                                                                                                                                                                                                                                              <  1:426896/250‑1 (MQ=255)    cGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGGCACCAACGGGGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTGCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATGGCGTCTGAAGTTGCCGATGTCAACCGCGACCTTAGTGGTGAAATCGACATGGCGTATAACAGCATGCCg                                                                                                                                                                                                                                          <  2:450879/251‑1 (MQ=255)            aaaGCCGCTATCGAGAAATTCGGCGAAAAATTGACCCAGCCTGGTTATATCGTCACCAACGGTGCCGATACCTTAAAATATTGGGTCGTAAACGAACGCATCGTTCTTGAACGCAGCCCGTCCTACTGGACCAACGCGAAAACCTTTATTAACCAGGTAAACTATTTGTCTATTGCTTCTGATGTTACCGATGTCAACCGCTACCGTAGTGGTTAACTCGACATGACGTATAACAGCATGTCGATCGAAtt                                                                                                                                                                                                                                  <  2:350413/251‑1 (MQ=255)                          gaAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTg                                                                                                                                                                                                                     >  1:336036/1‑250 (MQ=255)                                cGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCATAAGCTGAAAAAag                                                                                                                                                                                                              <  2:393701/251‑1 (MQ=255)                                         aTGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCgg                                                                                                                                                                                                     >  1:243496/1‑251 (MQ=255)                                                     tGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACg                                                                                                                                                                                         >  2:327209/1‑251 (MQ=255)                                                      ggTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTAGTAACCAGGTAACCTATTTGCCTATTGCGTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACg                                                                                                                                                                                         <  1:457864/250‑1 (MQ=255)                                                       gTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGtt                                                                                                                                                                                       >  2:128116/1‑251 (MQ=255)                                                          atatCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGAt                                                                                                                                                                                    <  2:336036/251‑1 (MQ=255)                                                            atCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATc                                                                                                                                                                                   <  1:447240/250‑1 (MQ=255)                                                                       cGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGATAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGAACGAAGTTCACGTTGATCCATACCTGTGCa                                                                                                                                                                       <  2:352150/251‑1 (MQ=255)                                                                                          gATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTa                                                                                                                                                     >  2:148948/1‑250 (MQ=255)                                                                                           aTTGGGTCGTAAACGAACGCATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCAGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTaac                                                                                                                                                   >  1:177263/1‑251 (MQ=255)                                                                                             tGGGTCGTAAACGAACGAAGCGTTCTTGAACGCAACCCGACCTACTGGAACAACGCGAAAACCGTTATTATCCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCATAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTaaca                                                                                                                                                  <  1:61097/250‑1 (MQ=255)                                                                                                       acgaacgaATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTaaca                                                                                                                                                  <  1:24778/240‑1 (MQ=255)                                                                                                       acgaacgaATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTaaca                                                                                                                                                  >  2:24778/1‑240 (MQ=255)                                                                                                        cgaacgaaTCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGc                                                                                                                                      >  1:65394/1‑251 (MQ=255)                                                                                                         gaacgaaTCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGcc                                                                                                                                     <  1:70294/251‑1 (MQ=255)                                                                                                          aacgaaTCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGATACCGATGTCAACCGCTACCGTCGTGGTGAAATCGACATGACGTCTAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGACATTAACAACCAGAAACCGcc                                                                                                                                     >  2:462034/1‑250 (MQ=255)                                                                                                                  gTTCTTGAACGCAGCCCTACCTACTTGAACAACGCGAAAACCGTTATTAACCAGGTAACTGATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACAGGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAATAGATCCCGGACGAAGTTCACGTTGATCCATACCTTTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGa                                                                                                                            <  2:408575/251‑1 (MQ=255)                                                                                                                                   gACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGACGTTAACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACTAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCg                                                                                                           >  1:303117/1‑251 (MQ=255)                                                                                                                                                          ccGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGCTGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATAGAATTGTTCCAGAAGCTGACAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGGGCACTTACTATTACGAAATTAACAACCCGAAACCGCCATTCAACGATGTGCGTGTGCGTCCCGCACTGAAACTAGGTATGGACcgc                                                                                    >  1:244779/1‑251 (MQ=255)                                                                                                                                                          ccGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACcgc                                                                                    <  1:327209/251‑1 (MQ=255)                                                                                                                                                               attaACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGatat                                                                               <  2:197061/251‑1 (MQ=255)                                                                                                                                                                                          gCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAggg                                                    >  2:98318/1‑251 (MQ=255)                                                                                                                                                                                          gCTTCGGAAGTTACCGATGGCAACCGCTACCGTATGGGGGAAATCTATACGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGTTCCCGGACGAAGTGCGCGTTGACCCATGCCTGTGCACTTACGATTACAAAATTAACAACCAGAAACCGCCATTCAACGATGTGGGTGTGCGTACCGCACTGAAACCAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAggg                                                    <  1:273417/251‑1 (MQ=255)                                                                                                                                                                                                                         gTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTATGGTTACACTCCACCg                     >  1:293482/1‑251 (MQ=255)                                                                                                                                                                                                                              ggtgAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTATGGTTACACTCCACCGTATAc                >  1:125970/1‑251 (MQ=255)                                                                                                                                                                                                                                            gACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAGACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTATGGTTACACTCCACCGTATACTGATGGCGCAAAAt  >  1:276112/1‑251 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                       CACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTATGGTTACACTCCACCGTATACTGATGGCGCAAAAT  >  NC_000913/1301753‑1302238

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: ‑