Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	3,679,399	0	G	T	57.1%	4.9 / 14.8	14	N1045K (AAC→AAA)	bcsC	cellulose synthase complex outer membrane protein BcsC
Reads supporting (aligned to +/- strand):  ref base G (4/2);  new base T (0/8);  total (4/10)
Fisher's exact test for biased strand distribution p-value = 1.50e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.55e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

AGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGAC  >  NZ_CP009273/3679319‑3679457                                                                                 |                                                           aGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTg                                                   <  1:292823/90‑1 (MQ=255)      catcaGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTc                                              >  2:395308/1‑90 (MQ=255)                       gCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTa                             >  2:526634/1‑90 (MQ=255)                             tttgGGTGGGTGAATGCGCCCACGAGCCAGACCCCCCCAGCCCCCACGACCATTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCaa                       <  1:319522/90‑1 (MQ=255)                                 ggTGCGTGAATGCGCCCACGAGCCAGCCGCCCCCAGCTCCCCCGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGa                   <  1:175024/90‑1 (MQ=255)                                   tGCGTGAATGCGCCCACGAGCCAGACGCCCCCAGCCCCCACGACCATTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACa                 <  1:178475/90‑1 (MQ=255)                                      gTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGt              >  1:258074/1‑90 (MQ=255)                                      gTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGt              >  1:75063/1‑90 (MQ=255)                                             cgcccACGAGCCACCCGCCCCCACCCCCCACGCCCTTTTTTCCGTCCGCCCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTg       <  2:258074/87‑1 (MQ=255)                                              ccccACGAGCCAGCCGCCCCCCCCCCCCACGACCATTTTTCCGTGCCCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGc      <  1:504060/88‑1 (MQ=255)                                                  aCGCGCCAGCCCCCCCCACCCCCCACGCCCATTTTTCCGGGCCCCCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGAc  <  1:370362/90‑1 (MQ=255)                                                  aCGCGCCACACCCCCCCAGCCCCCCCGCCCTTTTTTCCGTGCCCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGAc  <  2:15869/90‑1 (MQ=255)                                                  aCGAGCCAGCCGCCCCCAGCTCCCCCGACCATTTTCCCGGGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGAc  <  1:155982/90‑1 (MQ=255)                                                  aCGAGCCAGACCCCCCCGCCCCCCCCGACCTTTTTTCCGGGCCCCCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCCGCGGAc  <  2:6537/90‑1 (MQ=255)                                                                                 |                                                           AGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGAC  >  NZ_CP009273/3679319‑3679457

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