Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,301,998	A→G	S273G (AGC→GGC)	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,998	0	A	G	95.5%	71.3 / ‑7.4	22	S273G (AGC→GGC)	oppA	oligopeptide ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (13/8);  total (14/8)
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.96e-01

CGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTATGGTTACACTCCACCGTATACTGATGGCGCAAAATTGAC  >  NC_000913/1301756‑1302242                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     cGGTGCCAAAAGCCGCTATCGAGCAATTAGGCGAAAAAAGAACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGACTCGTTCTTGAACGCAGCCCGACATACTGGAACAACGCGACAACCGTTATTAACCCGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGATACCGTAGTGGTGAAATCGACATGACGAATAACGGCATGCCg                                                                                                                                                                                                                                              >  2:133818/1‑251 (MQ=255)   gTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACGGCATGCCGAt                                                                                                                                                                                                                                            <  2:264784/251‑1 (MQ=255)                                     aaTGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACGGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCg                                                                                                                                                                                                          <  1:133065/251‑1 (MQ=255)                                       tGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACGGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCgg                                                                                                                                                                                                         >  2:73011/1‑250 (MQ=255)                                              aGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACGGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGtt                                                                                                                                                                                                 <  1:126191/251‑1 (MQ=255)                                               gCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACGGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTc                                                                                                                                                                                                >  1:303474/1‑251 (MQ=255)                                                   ggTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGGATCGTTCTTGAACGCAGCCCGACCTACTGGAAGAACGCGAAAACCGTTATTAACCAGGTAACCTATGTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGATACCGTAGGGGTGAAATCGACAGGACGAATAACGGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGTACGATGTTCAc                                                                                                                                                                                              <  2:254892/249‑1 (MQ=255)                                                                     ggTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACGGCATGCCGCTCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCAc                                                                                                                                                                          >  1:183376/1‑251 (MQ=255)                                                                      gTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCTATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATATCGGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACt                                                                                                                                                                         <  1:131035/251‑1 (MQ=255)                                                                                                                            gCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACGGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATgtgcgtgtgcgt                                                                                                                   >  1:40393/1‑251 (MQ=255)                                                                                                                                       tGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACGGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGCAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACtta                                                                                                         >  2:262124/1‑248 (MQ=255)                                                                                                                                                   aaaaCCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACGGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATgg                                                                                             >  2:294539/1‑250 (MQ=255)                                                                                                                                                       acGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACGGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACACGGTATGGACcgc                                                                                        >  2:212017/2‑251 (MQ=255)                                                                                                                                                              taACCAGGTAACCTATTTGCCTATTGCTTCTGCAGTTACCGCTGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACGGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGACATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATc                                                                                  >  2:201120/1‑250 (MQ=255)                                                                                                                                                                                                gTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACGGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGcc                                               <  1:19146/251‑1 (MQ=255)                                                                                                                                                                                                   cccGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACGGCATGCCGATCGAATTGTTCAAGAATCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATGACGAAATTAACAACCATAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCg                                             <  1:354237/249‑1 (MQ=255)                                                                                                                                                                                                           cAACCGCTACCGTAGTAGTGAAATCGACATGACGAGTAACAGCATGCCCATCGAATTGTTACAGAAGCTGAAAAAAGAGCTCACGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAACTTAACAACCAGACACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGAATGGACCGCGATATCATTGTTAATAAAGCGAAAGCGCAGGGCAACATGCCCGCCTATGGtt                                    >  2:39790/1‑251 (MQ=255)                                                                                                                                                                                                                         gtggtgAAATCGACATGACGAATAACGGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTATGGTTACACTCCACCGtat                      >  1:1189/1‑251 (MQ=255)                                                                                                                                                                                                                             tgAAATCGACATGACGAATAACGGCATGCAGATCGAATTGTTCAAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACCTGCCCGCCTATGGTTACACTCCACCGTATACTg                  >  1:260113/1‑251 (MQ=255)                                                                                                                                                                                                                                aaTCGAAATGACGAATAACGGCATGCCGATCGAATGGTTCCAGAAGCTAAAAAAAGAAACCCCGGACGAAGTTACCGTTAATCCAACCCTGTGCACTTACTATACCGAAATTAACAACCAAAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTATGGTTACACTCCACCGTATACTGATg               <  1:228562/251‑1 (MQ=255)                                                                                                                                                                                                                                           cGAATAACGGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTATGGTTACACTCCACAGTATACTGATGGCGCAAAATTg    >  1:143134/1‑251 (MQ=255)                                                                                                                                                                                                                                             aaTAACGGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTATGGTTACACTCCACCGTATACTGATGGCGCAAAATTGAc  >  1:265266/1‑251 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     CGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTATGGTTACACTCCACCGTATACTGATGGCGCAAAATTGAC  >  NC_000913/1301756‑1302242

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