Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	3,569,316	0	A	C	57.7%	22.2 / 32.8	26	V369G (GTG→GGG)	gntU	gluconate transporter
Reads supporting (aligned to +/- strand):  ref base A (10/1);  new base C (0/15);  total (10/16)
Fisher's exact test for biased strand distribution p-value = 2.07e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.07e-06
Rejected as consensus: Frequency below/above cutoff threshold.

GGCTGACAACAATCGAACCACCAGCGATACAAATCGACAGCGCCGCCATTTGCGCACCGGAGTAGTTCAGTTGTTCAATAACCGGCATCACCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAGCACGAAGC  >  NZ_CP009273/3569227‑3569391                                                                                          |                                                                            ggCTGACAACAATCGAACCACCAGCGATACAAATCGACAGCGCCGCCATTTGCGCACCGGAGTAGTTCAGTTGTTCAATAACCGGCATCa                                                                             >  1:217103/1‑90 (MQ=255)   cTGACAACAATCGAACCACCAGCGATACAAATCGACAGCGCCGCCATTTGCGCACCGGAGTAGTTCAGTTGTTCAATAACCGGCATCAcc                                                                           <  1:796168/90‑1 (MQ=255)            aTCGACCCCCCAGCGATAAAAACCGACAGCGCCGCCATTTGCGCCCCGGAGTAGTTCAGTTGTTCAATAACCGGCATCCCCAGTCCTAcc                                                                  <  2:90405/90‑1 (MQ=255)            aTCAAACCCCCAGCGATACAAATCGACAGCCCCGCCATTTGCGCCCCGGAGTATTTCAGTTGTTCAATAACCGGCATCCCCAGTCCTAcc                                                                  <  2:125232/90‑1 (MQ=255)                 accaccAGCGATACAAATCGACAGCGCCGCCATTTGCGCACCGGAGTAGTTCAGTTGTTCAATAACCGGCATCACCAGTCCTACCGCCGt                                                             >  2:462308/1‑90 (MQ=255)                    accaGCGATACAAATCGACAGCGCCGCCATTTGCGCACCGGAGTAGTTCAGTTGTTCAATAACCGGCATCACCAGTCCTACCGCCGTTaa                                                          >  1:182234/1‑90 (MQ=255)                     caggCGAAAAAAACCGACAGCCCCGCCTTTTGCGCCCCGGAGTAGTTCAGTTTTTCAATACCCGGCATCCCCAGTCCTACCGCCGTTaaa                                                         <  1:263991/87‑1 (MQ=255)                      caGCGATACAAATCAACAGCGCCGCCATTTGGCCCCCGGAGTATTACAGTTCTCAAAAACCCGGCATCCCGAGTCCTACCGCCGTTAAAc                                                        <  2:794059/90‑1 (MQ=255)                          gATACAAATCGACAGCGCCGCCATTTGCGCACCGGAGTAGTTCAGTTGTTCAATAACCGGCATCACCAGTCCTACCGCCGTTAAACAGGc                                                    >  2:103870/1‑90 (MQ=255)                             aCAAATCGACAGCGCCGCCATTTGCGCACCGGAGTAGTTCAGTTGTTCAATAACCGGCATCACCAGTCCTACCGCCGTTAAACAGGCTAc                                                 >  2:22688/1‑90 (MQ=255)                                aaccGACGGCCCCCCCTTTTGCGCCCGGGAGTAGTTCATTTGTTAAATACCCGCCACCCCCAGTCCTACCGCCGTTAAACAGGCTACGGt                                              <  1:74961/87‑1 (MQ=255)                                  tCGACAGCGCCGCCATTTGCGCACCGGAGTAGTTCAGTTGTTCAATAACCGGCATCACCAGTCCTACCGCCGTTAAACAGGCTACggctg                                            >  2:263991/1‑87 (MQ=255)                                     acggcgcCGCCTTTGGCGCCCGGGATTATTTCAGTTGTCAAAAACCCGGCATCCCCATCCCTCCCGCCGTTAAACAGGCTACGGTGGCAg                                         <  1:193975/87‑1 (MQ=255)                                     aCAGCGCCCCCTTTGGCGCCCGGGAGTATTTCGTTTTTTAAAAACCCGCCACCCCCACCCCTACCGCCGTTAAACAGGCTACGGTGGCAg                                         <  1:171762/90‑1 (MQ=255)                                     aCAGCCCCCCCTTTGGCCCCCCGGAGTATTTCATTTTTTAAATACCCGCCACCCCCAGCCCTCCCGCCGTTAAACAGGCTACGGTGGCAg                                         <  1:2792/90‑1 (MQ=255)                                       aGCGCCGCCATTTGCGCCCCGGAGTAGTTCAGTTGTTCAATAACCGGCATCCCCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGaa                                       <  2:123137/90‑1 (MQ=255)                                        gcgcCGCCATTTGCGCACCGGAGTAGTTCAGTTGTTCAATAACCGGCATCACCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAAc                                      >  1:479/1‑90 (MQ=255)                                        gcgcCGCCATTTGCGCACCGGAGTAGTTCAGTTGTTCAATAACCGGCATCACCAGTCCTACCGCCGTCAAACAGGCTACGGTGGCAGAAc                                      >  2:341506/1‑90 (MQ=255)                                           ccgccATTTGCGCACCGGAGTAGTTCAGTTGTTCAATAACCGGCATCACCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCt                                   >  1:459432/1‑90 (MQ=255)                                              ccATTTGCGCACCGGATTAGTTCAGTTGTTAAATAACCGGCATCCCCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGaa                                <  2:30995/90‑1 (MQ=255)                                                      gccccGGAGTATTTCAGTTGTTCAATACCCGGCATCCCCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATgcgc                        <  2:103902/87‑1 (MQ=255)                                                        ccaggAGAAGTTCAGTTGTTCAATACCCGCCACCCCCACCCCCCCCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCAc                      <  1:614184/87‑1 (MQ=255)                                                             agtagtTCAGTTGTTCAATACCCGGCATCCCCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAg                 <  2:100416/90‑1 (MQ=255)                                                                         gttCAATAACCGGCATCACCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAGCACGa     >  1:63719/1‑90 (MQ=255)                                                                         gtcaaaTACCCGGCATCCCCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAGCACGa     <  2:194389/86‑1 (MQ=255)                                                                           tCAATAACCGGCATCACCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAGCACGAAg   >  2:3253/1‑90 (MQ=255)                                                                            aaaTACCCGGCATCCCCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAGCACGAAGc  <  2:182234/89‑1 (MQ=255)                                                                                          |                                                                            GGCTGACAACAATCGAACCACCAGCGATACAAATCGACAGCGCCGCCATTTGCGCACCGGAGTAGTTCAGTTGTTCAATAACCGGCATCACCAGTCCTACCGCCGTTAAACAGGCTACGGTGGCAGAACCCTGAATGATGCGCACTGCAGCTGCCAGCACGAAGC  >  NZ_CP009273/3569227‑3569391

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑