Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,111,968	0	C	T	30.3%	34.0 / 22.6	33	A369V (GCC→GTC)	opgH	osmoregulated periplasmic glucans biosynthesis protein H
Reads supporting (aligned to +/- strand):  ref base C (17/6);  new base T (5/5);  total (22/11)
Fisher's exact test for biased strand distribution p-value = 2.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.40e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GTTATAACCCGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTGCGGGCTGGTGCGCCTGATGGAAGCCAACCCGAACGCCGGGATCATTCAGTCGTCGCCGAAAGCGTCCGGTATGGATACGCTGTATGCGCGCTGTCAGCAGTTCGCGACCCGCGTGTATGGGCCACTGTTTACAGCCGGTTTGCACTTCTGGCAACTTGGCGAGTCGCACTACTGGGGACATAACGCGATTATCCGCGTGAAACCGTTTATCGAGCACTGCGCACTGGCTCCGCTGCCGGGCGAAGGTTCCTTTGCCG  >  NC_000913/1111719‑1112189                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                              gTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTAAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTGCGGGCTGGTGCGCCTGATGGAAGCCAACCCGAACGcc                                                                                                                                                                                                                              >  1:31847/1‑251 (MQ=255)    ataacttggATATCTGCGTCGCAGAGCAACAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTGCGGGCTGGCGCGCCTGATGGAAGCCAACCCGAACGCCggg                                                                                                                                                                                                                           >  2:251402/8‑251 (MQ=255)      aacttggATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTGCGGGCTGGTGCGCCTGATGGAAGCCAACCCGAACGCCGGGAt                                                                                                                                                                                                                         <  1:233158/246‑1 (MQ=255)             atatCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTGCGGGCTGGTGCGCCTGATGGAAGCCAACCCGAACGCCGGGATCATTc                                                                                                                                                                                                                    >  1:220358/1‑249 (MQ=255)                     gTCGCAGAGCAAAAAGCCTGGATGGAGCTTATAGCTGAAGTCGGTGGCGACGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGGGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGCGCGGGCTGGTGCGCCTGATGGAAGCCAACCCGAACGCCGGGATCATTCCGTCGTCGcc                                                                                                                                                                                                          >  1:256398/1‑251 (MQ=255)                           gagCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTGCGGGCTGGTGCGCCTGATGGAAGCCAACCCGAACGTCGGGATCATTCAGTCGTCGCCGAAAGc                                                                                                                                                                                                    >  1:24076/1‑251 (MQ=255)                                    gCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTGCGGGCTGGTGCGCCTGATGGAAGCCAACCCGAACGCCGGGATCATTCAGTCGTCGCCGAAAGCGTCCGGTAt                                                                                                                                                                                           >  1:20867/1‑251 (MQ=255)                                           tggaGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTGCGGGCTGGTGCGCCTGATGGAAGCCAACCCGAACGCCGGGATCATTCAGTCGTCGCCGAAAGCGTCCGGTATGGATACg                                                                                                                                                                                    >  1:266110/1‑251 (MQ=255)                                            ggaGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTGCGGGCTGGTGCGCCTGATGGAAGCCAACCCGAACGCCGGGATCATTCAGTCGTCGCCGAAAGCGTCCGGTATGGATACGc                                                                                                                                                                                   <  2:265268/251‑1 (MQ=255)                                               gCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTGCGGGCTGGTGCGCCTGATGGAAGCCAACCCGAACGTCGGGATCATTCAGTCGTCGCCGAAAGCGTCCGGTATGGATACGCTGt                                                                                                                                                                                <  2:28547/251‑1 (MQ=255)                                                cTTATCGCTGAAGTCGGTGGCGAAGGTCATATTTTCTATCGCCGCCGCCGGCGCCGCGTGAAGCGTAAACGCTGTAATATCGATGACTCATGCCGTCGCTGGGGCAGCCAGTACAGCTACCTGGTGGTGCTGGATACTGACTCGGTAATGACCGGCGATTGTTTGTGAGGGCTGGTGCGCCTGCTGGAAGCCAACCCGAACGTCGGGATCATTCAGTCGGCGCCGCACGCGTCCGCGATGGATACGCTGtt                                                                                                                                                                               >  2:67560/1‑250 (MQ=255)                                                       cTGAAGTCGACGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTACTATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTGCGGGCTGGTGCGCCTGATGGAAGCCAACCCGAACGCCGGGATCATTCAGTCGTCGCCGAAAGCGTCCGGTCTGGATACGCTGTATgcgcgc                                                                                                                                                                        >  2:40258/1‑251 (MQ=255)                                                               ggTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACCTGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTGCGGGCTGGTGCGCCTGATGGAAGCCAACCCGAACGCCGGGATCATTCAGTCGTCGCCGAAAGCGTCCGGTATGGATACGCTGTATGCGCGCTGTcagc                                                                                                                                                                 >  1:267417/1‑250 (MQ=255)                                                               ggTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTGCGGGCTGGTGCGCCTGATGGAAGCCAACCCGAACGCCGGGATCATTCAGTCGTCGCCGAAATCGTCCGGTATGGATACGCTGTATGCGCGCGTTcagca                                                                                                                                                                <  2:40582/251‑1 (MQ=255)                                                                                                    gccgcGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGGTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTGCGGGCTGGTGCGCCTGATGGAAGCCAACCCGAACGTCGGGATCATTCAGTCGTCGCCGAAAGCGTCCGGTATGGATACGCTGTATGCGCGCTGTCAGCAGTTCGCGACCCGCGTGGATGGGCCACTGTTTACAGcc                                                                                                                           <  1:61005/251‑1 (MQ=255)                                                                                                                         gTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTGCGGGCTGGTGCGCCTGATGGAAGCCAACCCGAACGCCGGGATCATTCAGTCGTCGCCGAAAGCGTCCGGTATGGATACGCTGTATGCGCGCTGTCAGCAGTTCGCGACCCGCGTGTATGGGCCACTGTTTACAGCCGGTTTGCACTTCTGGCAACtt                                                                                                      >  1:216265/1‑251 (MQ=255)                                                                                                                              atCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTGCGGGCTGGTGCGCCTGATGGAAGCCAACCCGAACGCCGGGATCATTCAGTCGTCGCCGAAAGCGTCCGGTATGGATACGCTGTATGCGCGCTGTCAGCAGTTCGCGACCCGCGTGTATGGGCCACTGTTTACAGCCGGTTTGCACTTCTGGCAACTTGGCg                                                                                                  >  1:69943/1‑250 (MQ=255)                                                                                                                                  aTGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTGCGGGCTGGTGCGCCTGATGGAAGCCAACCCGAACGTCGGGATCATTCAGTCGTCGCCGAAAGCGTCCGGTATGGATACGCTGTATGCGCGCTGTCAGCAGTTCGCGACCCGCGTGTATGGGCCACTGTTTACAGCCGGTTTGCACTTCTGGCAACTTGGCGAGTCg                                                                                             >  2:83619/1‑251 (MQ=255)                                                                                                                                   tGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTGCGGGCTGGTGCGCCTGATGGAAGCCAACCCGAACGCCGGGATCATTCAGTCGTCGCCGAAAGCGTCCGGTATGGATACGCTGTATGCGCGCTGTCAGCAGTTCGCGACCCGCGTGTATGGGCCACTGTTTACAGCCGGTTTGCACTTCTGGCAACTTGGCGAGTCGc                                                                                            >  1:72463/1‑251 (MQ=255)                                                                                                                                    gACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTGCGGGCTGGTGCGCCTGATGGAAGCCAACCCGAACGTCGGGATCATTCAGTCGTCGCCGAAAGCGTCCGGTATGGATACGCTGTATGCGCGCTGTCAGCAGTTCGCGACCCGCGTGTATGGGCCACTGTTTACAGCCGGTTTGCACTTCTGGCAACTTGGCGAGTCGcc                                                                                           >  2:71230/1‑250 (MQ=255)                                                                                                                                          tGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTGCGGGCTGGTGCGCCTGATGGAAGCCAACCCGAACGCCGGGATCATTCAGTCGTCGCCGAAAGCGTCCGGTATGGATACGCTGTATGCGCGCTGTCAGCAGTTCGCGACCCGCGTGTATGGGCCACTGTTTACAGCCGGTTTGCACTTCTGGCAACTTGGCGAGTCGCACTACTg                                                                                     >  2:81684/1‑251 (MQ=255)                                                                                                                                                   tGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTGCGGGCTGGTGCGCCTGATGGAAGCCAACCCGAACGTCGGGATCATTCAGTCGTCGCCGAAAGCGTCCGGTATGGATACGCTGTATGCGCGCTGTCAGCAGTTCGCGACCCGCGTGTATGGGCCACTGTTTACAGCCGGTTTGCACTTCTGGCAACTTGGCGAGTCGCACTACTGGGGACATa                                                                             >  2:264172/1‑250 (MQ=255)                                                                                                                                                     gggCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTGCGGGCTGGTGCGCCTGATGGAAGCCAACCCGAACGCCGGGATCATTCAGTCGTCGCCGAAAGCGTCCGGTATGGATACGCTGTATGCGCGCTGTCAGCAGTTCGCGACCCGCGTGTATGGGCCACTGTTTACAGCCGGTTTGCACTTCTGGCAACTTGGCGAGTCGCACTACTGGGGACATAAcg                                                                          >  1:119089/1‑251 (MQ=255)                                                                                                                                                            cAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTGCGGGCTGGTGCGCCTGATGGAAGCCAACCCGAACGTCGGGATCATTCAGTCGTCGCCGAAAGCGTCCGGTATGGATACGCTGTATGCGCGCTGTCAGCAGTTCGCGACCCGCGTGTATGGGCCACTGTTTACAGCCGGTTTGCACTTCTGGCAACTTGGCGAGTCGCACTACTGGGGACATAACGCGATTAt                                                                   <  1:32782/251‑1 (MQ=255)                                                                                                                                                                aCAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTGCGGGCTGGTGCGCCTGATGGAAGCCAACCCGAACGCCGGGATCATTCAGTCGTCGCCGAAAGCGTCCGGTATGGATACGCTGTATGCGCGCTGTCAGCAGTTCGCGACCCGCGTGTATGGGCCACTGTTTACAGCCGGTTTGCACTTCTGGCAACTTGGCGAGTCGCACTACTGGGGACATAACGCGATTATCcgc                                                               >  1:94862/1‑251 (MQ=255)                                                                                                                                                                          ggtggtGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTGCGGGCTGGTGCGCCTGATGGAAGCCAACCCGAACGTCGGGATCATTCAGTCGTCGCCGAAAGCGTCCGGTATGGATACGCTGTATGCGCGCTGTCAGCAGTTCGCGACCCGCGTGTATGGGCCACTGTTTACAGCCGGTTTGCACTTCTGGCAACTTGGCGAGTCGCACTACTGGGGACATAACGCGATTATCCGCGTGAAACCGt                                                     <  1:71230/251‑1 (MQ=255)                                                                                                                                                                          ggtggtGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTGCGGGCTGGTGCGCCTGATGGAAGCCAACCCGAACGCCGGGATCATTCAGTCGTCGCCGAAAGCGTCCGGTATGGATACGCTGTATGCGCGCTGTCAGCAGTTCGCGCCCCGCGTGTATTGGCCACTGTTTCCAGCCGGTTTGCACTTCTGGCAACTTGGCGAGTCGCACTACTGGGGACATAACGCGCTTATCCGCGTGACACCGt                                                     >  2:62699/1‑251 (MQ=255)                                                                                                                                                                          ggtggtGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTGCGGGCTGGTGCGCCTGATGGAAGCCAACCCGAACGCCGGGATCATTCAGTCGTCGCCGAAAGCGTCCGGTATGGATACGCTGTATGCGCGCTGTCAGCAGTTCGCGAcc                                                                                                                                                      >  2:10198/1‑154 (MQ=255)                                                                                                                                                                          ggtggtGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTGCGGGCTGGTGCGCCTGATGGAAGCCAACCCGAACGCCGGGATCATTCAGTCGTCGCCGAAAGCGTCCGGTATGGATACGCTGTATGCGCGCTGTCAGCAGTTCGCGAcc                                                                                                                                                      <  1:10198/154‑1 (MQ=255)                                                                                                                                                                                gCTGGATGCTGACTAGGTAATGACCGGTGATTGTTTGTGCGGGCTGGTGCGCCTGATGGAAGCCAACCCGAACGCCGGGATCATTCAGTCGTCGCCGAAAGCGTCCGGTATGGATACGCTGTATGCGCGCTGTCAGCAGTTCGCGACCCGCGTGTATGGGCCACTGTTTACAGCCGGTTTGCACTTCTGGCAACTTGGCGAGTCGCACTACTGGGGACATAACGCGATTATCCGCGTGACACCGTTTATCg                                               >  2:251728/1‑251 (MQ=255)                                                                                                                                                                                     aTGCTGACTCGGTAATGACCGGTGATTGTTTGTGCGGGCTGGTGCGCCTGATGGAAGCCAACCCGAACGTCGGGATCATTCAGTCGTCGCCGAAAGCGTCCGGTATGGATACGCTGTATGCGCGCTGTCAGCAGTTCGCGACCCGCGTGTATGGGCCACTGTTTACAGCCGGTTTGCACTTCTGGCAACTTGGCGAGTCGCACTACTGGGGACATAACGCGATTATCCGCGTGAAACCGTTTATCGAGCAc                                          <  1:68183/251‑1 (MQ=255)                                                                                                                                                                                                       ccGGTTATTGTTTGTGCGGGCTGGTGCGCCTGATGGAAGCCAAACCGAACGCCGGGATCATTCAGTCGTCGCCGAAAGCGTCCGGTAGGGATACGCTGTATGCGCGCTGTCAGCAGTGCGCGACCCGCGTGTATGGGCCACTGTTGACAGCCGGTTTGCACGTCTGGCAACTTGGCGAGTCGCACTACTGTGGACATAACGCGATTATCCGCGTGAAACCGGCTATCGAGCACTGCGCACTGGGTCCGCTg                        <  2:14132/251‑1 (MQ=255)                                                                                                                                                                                                                        gggCTGGTTCGCCTGATGGAAGCCAACCCGAACGCCGGGATCATTCAGTCGTCGCCGAAAGCGTCCGGTTTGGATACGCTGTATGCGCGCTGTCAGCAGTTCTCGCCCCGCGTGTATGGGCCACTGTTTACAGCCGGTTTGCACTTCTGGCAACTTGGCGAGTCGCACTACTGGGGACATAACGCGATTATCCGCGTGAAACCGTTTATCGAGCACTGCGCACTGGCTCCGCTGCCGGGCGAAGGTTCCtt       <  2:233947/251‑1 (MQ=255)                                                                                                                                                                                                                             ggTGCGCCTGATGGAAGCCAACCCGAACGCCGGGATCATTCAGTCGTCGCCGAAAGCGTCCGGTATGGATACGCTGTATGCGCGCTGTCAGCAGTTCGCGACCCGCGTGTATGGGCCACTGTTTACAGCCGGTTTGCACTTCTGGCAACTTGGCGAGTCGCACTACTGGGGACATAACGCGATTATCCGCGTGAAACCGTTTATCGAGCACTGCGCACTGGCTCCGCTGCCGGGCGAAGGTTCCTTTGCCg  >  1:144057/1‑251 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                              GTTATAACCCGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTGCGGGCTGGTGCGCCTGATGGAAGCCAACCCGAACGCCGGGATCATTCAGTCGTCGCCGAAAGCGTCCGGTATGGATACGCTGTATGCGCGCTGTCAGCAGTTCGCGACCCGCGTGTATGGGCCACTGTTTACAGCCGGTTTGCACTTCTGGCAACTTGGCGAGTCGCACTACTGGGGACATAACGCGATTATCCGCGTGAAACCGTTTATCGAGCACTGCGCACTGGCTCCGCTGCCGGGCGAAGGTTCCTTTGCCG  >  NC_000913/1111719‑1112189

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