Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	108,789	0	T	C	53.8%	5.3 / 12.6	14	S100G (AGC→GGC)	zapD	cell division protein ZapD
Reads supporting (aligned to +/- strand):  ref base T (4/2);  new base C (0/7);  total (4/10)
Fisher's exact test for biased strand distribution p-value = 2.10e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.13e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GATGCTCAGTCGCTGACGCACCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCTTTTAACTGCTGAATTAATGCTTCAATACGGCTCTGGTCCACGCCAGGCACGCCAATCC  >  NZ_CP009273/108706‑108857                                                                                    |                                                                     gATGCTCAGTCGCTGACGCACCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCccgc                                                                <  2:58248/90‑1 (MQ=255)    gCTCAGTCGCTGACGCACCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCcgccgc                                                             >  2:406926/1‑90 (MQ=255)         gTCGCTGACGCACCAGAGCAATAAAACGATCTTCCCGCAAAAATTGCCCGATACGCGGCGCGGAAATTAATCCGCCCCCCGCCGCTTTTa                                                        <  1:361082/90‑1 (MQ=255)                       agagCAATCAAACGATCTTCCCGCAAAATTTGCCCGATACGCGGCGCGAAAATTAATCCGCCCCCCGCCGCTTTTAACTGCTGAATTAAt                                          <  1:427043/90‑1 (MQ=255)                            aaTCAAACGACCTCCACGCAGAAATTGCCCGATTCGCGGCGCGGAAATTAATACGCACCCCGCCGCTTTTAACTGCTGAATTAATGCTTc                                     <  2:182696/90‑1 (MQ=255)                             aTCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCTTTTAACTGCTGAATTAATGCTTCa                                    >  2:290379/1‑90 (MQ=255)                               cAAACAATCTCCACGCAAAAATTCCCCGATACGCGGCGCGGAAATTAATACGCCCCCCGCCGCTTTTAACTGCTGAATTAATGCTTCAAt                                  <  2:400657/90‑1 (MQ=255)                                 aaCGATTTTCACGAAAAAATTCCCCGATACGGGGCGGGAAAATTAATACGCCCCCCGCCGCTTTTAACTGCTGAATTAATGCTTCAATAc                                <  2:300370/90‑1 (MQ=255)                                 aaCGATCTTCACGCAAAAATTGCCCGATCCGCGGCGCGAAAATTAATACGCTCCCCGCCGCTTTTAACTGCTGAATTAATGCTTCAATAc                                <  2:197761/90‑1 (MQ=255)                                 aaCGATCTCCCCGCAAAAATTCCCCGATCCGCGGCGGGAAAATTAATCCGCCCCCCGCCGCTTTTAACTGCTGAATTAATGCTTCAATAc                                <  1:268422/90‑1 (MQ=255)                                           aCGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCTTTTAACTGCTGAATTAATGCTTCAATACGGCTCTGGTc                      >  1:328605/1‑90 (MQ=255)                                               aGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCTTTTAACTGCTGAATTAATGCTTCAATACGGCTCTGGTCCACg                  >  2:171423/1‑90 (MQ=255)                                                               gcggcgcgaAAATTAATCCCCCCCCCGCCGTTTTTAACTGCTGAATTAATGCTTCAATACGGCTCTGGTCCACGCCAGGCACGCCAATcc  <  2:452608/90‑1 (MQ=255)                                                               gcggcgcgaAAATTAATACGCCCCCCGCCGCTTTTAACTGCTGAATTAATGCTTCAATACGGCTCTGGTCCACGCCAGGCACGCCAATcc  <  1:171423/90‑1 (MQ=255)                                                                                    |                                                                     GATGCTCAGTCGCTGACGCACCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCTTTTAACTGCTGAATTAATGCTTCAATACGGCTCTGGTCCACGCCAGGCACGCCAATCC  >  NZ_CP009273/108706‑108857

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