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	93.9%	105.6 / NA	33	Q246Y (CAG→TAT)	oppA	oligopeptide ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base T (17/14);  total (17/16)
Fisher's exact test for biased strand distribution p-value = 2.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

GATGATCACACATTAGAAGTCACCTTAAGTGAACCCGTTCCGTACTTCTATAAATTACTTGTTCACCCATCAACTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGC  >  NC_000913/1301677‑1302156                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                              gatgatCACACATTAGAAGTCACCTTAAGTGAACCCGTTCCGTACTTCTATAAATTACTTGTTCACCCATCAACTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCt                                                                                                                                                                                                                                        >  1:327517/1‑250 (MQ=255)                          tAAGTAAACCCGTTCCGTACTTCTATAAATTACTTGTTCACCCATCAACTTCACCGGTGCCAAAAGCCGCGATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGGAACCTATTGGCCTATTGCTTCTGAAGTTa                                                                                                                                                                                                                <  2:248642/249‑1 (MQ=255)                              tGAACCCGTTCCGTACTTCTATAAATTACTTGTTCACCCATCAACTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATg                                                                                                                                                                                                          >  1:141272/1‑251 (MQ=255)                                  cccGTTCCGTACTTCTATAAATTACTTGTTCACCCATCAAATTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTCATATCGTACCCAACGGTGCCTATACCTTAAAAGCTTGGGTCGTAAACGAACGACTCGTTCTTGAACGCCGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCa                                                                                                                                                                                                       >  1:23126/1‑250 (MQ=255)                                              ttCTATAAATTACTTGTTCACCCATCAACTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGt                                                                                                                                                                                            >  2:121328/1‑249 (MQ=255)                                                    aaaTTACTTGTTCACCCATCACCTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGATCGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATGTGCCTATGGCTTCTGAAGTTACCCATGTCAACCGCTACCGTAGTGGTGa                                                                                                                                                                                    <  1:305867/251‑1 (MQ=255)                                                                        aaCTTCACCGATGCCAAAAGCCGCTATCGAGATATTCCGCGAGAAATGGACCCAGCCTGGTAATCTCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTCACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAAc                                                                                                                                                                 >  2:322847/1‑250 (MQ=255)                                                                            tCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCa                                                                                                                                                             <  2:244351/250‑1 (MQ=255)                                                                                      cAAAAGCCGGTATCGAGAAATGCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATGCCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGGTATTAACTAGGGATCCTATTGGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTGCCTTAGGGGTGAAATCGACATGACAAATAACAGCATGCCGATCGaa                                                                                                                                                  <  1:318307/251‑1 (MQ=255)                                                                                          aGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATTACGAATAACAGCATGCCGATCGAATTGt                                                                                                                                              <  2:327517/251‑1 (MQ=255)                                                                                           gCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGtt                                                                                                                                             <  1:121328/251‑1 (MQ=255)                                                                                              gCTATCGAGAAATTCGGCGAAACATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTcc                                                                                                                                           >  1:70074/1‑250 (MQ=255)                                                                                                     ataaaTTCAGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTAGTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTg                                                                                                                                   <  2:192059/249‑1 (MQ=255)                                                                                                     agaAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTg                                                                                                                                   >  2:83374/1‑251 (MQ=255)                                                                                                               cGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCTACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAgaga                                                                                                                         <  2:232751/251‑1 (MQ=255)                                                                                                               cGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGAGACTATACATTAAAAGATTGGGTCAAAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAAACGTTATTAACTATGTAACCTCTTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAag                                                                                                                            >  1:347768/1‑248 (MQ=255)                                                                                                                   aaaTGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTGACCGATGTCAACCGCGACCGTAGTGGTGAAATCGACATGACGAATAACATCAGGCCGATCGAATTGTCCCAGAAGCTGAAAAAAGAGATccc                                                                                                                     <  2:40143/251‑1 (MQ=255)                                                                                                                       gtaCCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGAc                                                                                                                 <  2:79300/249‑1 (MQ=255)                                                                                                                         aCCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGa                                                                                                               <  1:18503/251‑1 (MQ=255)                                                                                                                                   gTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGtt                                                                                                     <  1:339309/251‑1 (MQ=255)                                                                                                                                    tAATATCGTCACACACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTg                                                                                                    >  1:76420/1‑251 (MQ=255)                                                                                                                                              aCCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCt                                                                                          >  1:270642/1‑251 (MQ=255)                                                                                                                                                  aCGGTGCCTCTACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGc                                                                                      >  1:261703/1‑251 (MQ=255)                                                                                                                                                           ataCCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTa                                                                              <  1:160521/250‑1 (MQ=255)                                                                                                                                                                  aaaaGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGaa                                                                       >  1:256915/1‑250 (MQ=255)                                                                                                                                                                          gggTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAGACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATGGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGTACGAAGGGCACGGTGATCCATACCTGTGCACTTACTATTACGAAATTaacaac                                                              <  1:154966/251‑1 (MQ=255)                                                                                                                                                                                     gaacgaaTCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGcc                                                   <  2:1304/251‑1 (MQ=255)                                                                                                                                                                                         gaaTCGTTCTTGAACGCCGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTCACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGAGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTACCAACCAGAAACCGCCATTc                                               >  1:243416/1‑251 (MQ=255)                                                                                                                                                                                                  ttGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGAtgtg                                      >  2:117139/1‑251 (MQ=255)                                                                                                                                                                                                         cAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGAtgtgcgtgtgc                               >  1:92212/1‑251 (MQ=255)                                                                                                                                                                                                           gCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCCTGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATgtgcgtgtgcgt                             >  1:226625/1‑251 (MQ=255)                                                                                                                                                                                                                          acaacGCAAAAACCGTTATTATCTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCACCCGCTACCGTAGGGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTa              <  2:141272/251‑1 (MQ=255)                                                                                                                                                                                                                                      ccGTTATTAACTATGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGCTCCATACCTGTGCACTTACTATTACGACATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGAACCGCACTGAAACTAGGTATGGACcgc  >  1:168364/1‑251 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                              GATGATCACACATTAGAAGTCACCTTAAGTGAACCCGTTCCGTACTTCTATAAATTACTTGTTCACCCATCAACTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGC  >  NC_000913/1301677‑1302156

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