Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	990,891	0	T	C	52.9%	‑5.3 / 22.3	17	E185G (GAG→GGG)	ssuD	FMNH2‑dependent alkanesulfonate monooxygenase
Reads supporting (aligned to +/- strand):  ref base T (3/5);  new base C (4/5);  total (7/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.20e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CGAATACCGAAACGAATTTTGCGTCCATGCGCGGCAGCTTTCGCCCGCACTTGTTCGATTTTCTCTTTAACCAGTTCCGGCGGTTCGCCCCAGGTGAGGTAGAGATCAACCTGTTCTGCCGCCAGCTCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGTTGAATCGCCGGGAAGAGCAGTTTTGCTCCGCGCACATGAATATGTTTACCGTTGAAATCGACGGTTTC  >  NZ_CP009273/990765‑991009                                                                                                                               |                                                                                                                       cGAATACCGAAACGAATTTTGCGTCCATGCGCGGCAGCTTTCGCCCGCACTTGTTCGATTTTCTCTTTAACCAGTTCCGGCGGTTCGCCCCAGGTGAGGTAGAGATCAACCTGTTCTGCCGCCAGCTCCTGGGCGACAt                                                                                                            <  1:380621/139‑1 (MQ=255)    aTACCGAAACGAATTTTGCGTCCATGCGCGGCAGCTTTCGCCCGCACTTGTTCGATTTTCTCTTTAACCAGTTCCGGCGGTTCGCCCCAGGTGAGGTAGAGATCAACCTGTTCTGCCGCCAGCTCCTGGGCGACATCt                                                                                                          >  2:467805/1‑138 (MQ=255)          aaaCGAATTTTGCGTCCATGCGCGGCAGCTTTCGCCCGCACTTGTTCGATTTTCTCTTTAACCAGTTCCGGCGGTTCGCCCCAGGTGAGGTAGAGATCAACCTGTTCTGCCGCCAGCCCCTGGGCGACATCTGACGATc                                                                                                   >  2:294260/1‑139 (MQ=255)          aaaCGAATTTTGCGTCCATGCGCGGCAGCTTTCGCCCGCACTTGTTCGATTTTCTCTTTAACCAGTTCCGGCGGTTCGCCCCAGGTGAGGTAGAGATCAACCTGTTCTGCCGCCAGCCCCTGGGCGACATCTGACGATc                                                                                                   >  2:226825/1‑139 (MQ=255)           aaCGAATTTTGCGTCCATGCGCGGCAGCTTTCGCCCGCACTTGTTCGATTTTCTCTTTAACCAGTTCCGGCGGTTCGCCCCAGGTGAGGTAGAGATCAACCTGTTCTGCCGCCAGCCCCTGGGCGACATCTGACGATcc                                                                                                  >  1:134000/1‑139 (MQ=255)               aaTTTTGCGTCCATGCGCGGCAGCTTTCGCCCGCACTTGTTCGATTTTCTCTTTAACCAGTTCCGGCGGTTCGCCCCAGGTGAGGTAGAGATCAACCTGTTCTGCCGCCAGCCCCTGGGCGACATCTGACGATCCGCCa                                                                                              <  1:425655/139‑1 (MQ=255)                          cATGCGCGGCAGCTTTCGCCCGCACTTGTTCGATTTTCTCTTTAACCAGTTCCGGCGGTTCGCCCCAGGTGAGGTAGAGATCAACCTGTTCTGCCGCCAGCCCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTg                                                                                    <  1:226825/138‑1 (MQ=255)                                cgGCAGCTTTCGCCCGCACTTGTTCGATTTTCTCTTTAACCAGTTCCGGCGGTTCGCCCCAGGTGAGGTAGAGATCAACCTGTTCTGCCGCCAGCCCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGATa                                                                             <  2:256226/139‑1 (MQ=255)                                              cGCACTTGTTCGATTTTCTCTTTAACCAGTTCCGGCGGTTCGCCCCAGGTGAGGTAGAGATCATCCTGTTCTGCCGCCAGCCCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGTTGAATCg                                                               <  1:294260/139‑1 (MQ=255)                                                                      aCCAGTTCCGGCGGTTCGCCCCAGGTGAGGTAGAGATCAACCTGTTCTGCCGCCAGCTCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGTTGAATCGCCGGGAAGAGCAGTTTTGCTCcgc                                       <  1:609060/139‑1 (MQ=255)                                                                        cAGTTCCGGCGGTTCGCCCCAGGTGAGGTAGAGATCAACCTGTTCTGCCGCCAGCTCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGTTGAATCGCCGGGAAGAGCAGTTTTGCTCCgcgc                                     <  1:562436/139‑1 (MQ=255)                                                                         aGTTCCGGCGGTTCGCCCCAGGTGAGGTAGAGATCAACCTGTTCTGCCGCCAGCTCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGTTGAATCGCCGGGAAGAGCAGTTTTGCTCCGCGca                                    <  2:270534/139‑1 (MQ=255)                                                                            tCCGGCGGTTCGCCCCAGGTGAGGTAGAGATCAACCTGTTCTGCCGCCAGCCCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGTTTAATCGCCGGGAAGAGCAGTTTTGCTCCGCGCACAt                                 >  1:495526/1‑139 (MQ=255)                                                                                gcggTTCGCCCCAGGTGAGGTAGAGATCAACCTGTTCTGCCGCCAGCCCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGTTGAATCGCCGGGAAGAGCAGTTTTGCTCCGCGCACATGAat                             <  2:495526/139‑1 (MQ=255)                                                                                                           cAACCTGTTCTGCCGCCAGCTCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGTTGAATCGCCCGGGAGGGCAGGTTTTCTGCGCGCCCATGAATATGTTTTCCGCTGAAAAAGCCGGtttt  >  2:184193/1‑138 (MQ=255)                                                                                                                   tCTGCCGCCAGCTCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGTTGAATCGCCGGGAAGAGCAGTTTTGCTCCgcgc                                     <  2:435313/96‑1 (MQ=255)                                                                                                                   tCTGCCGCCAGCTCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGTTGAATCGCCGGGAAGAGCAGTTTTGCTCCgcgc                                     >  1:435313/1‑96 (MQ=255)                                                                                                                               |                                                                                                                       CGAATACCGAAACGAATTTTGCGTCCATGCGCGGCAGCTTTCGCCCGCACTTGTTCGATTTTCTCTTTAACCAGTTCCGGCGGTTCGCCCCAGGTGAGGTAGAGATCAACCTGTTCTGCCGCCAGCTCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGTTGAATCGCCGGGAAGAGCAGTTTTGCTCCGCGCACATGAATATGTTTACCGTTGAAATCGACGGTTTC  >  NZ_CP009273/990765‑991009

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