Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	1,932,672	A→G	E256G (GAA→GGA)	pyk →	pyruvate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	1,932,672	0	A	G	92.9%	38.7 / ‑3.1	14	E256G (GAA→GGA)	pyk	pyruvate kinase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (6/7);  total (6/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.11e-01

TGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGC  >  NZ_CP009273/1932616‑1932753                                                         |                                                                                  tGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGGAATTGGCGACCCGGAACTGGTCGGCATTCAGaa                                                  <  2:218927/90‑1 (MQ=255) tGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGGAATTGGCGACCCGGAACTGGTCGGCATTCAGaa                                                  <  2:245425/90‑1 (MQ=255) tGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGaa                                                  <  1:312965/90‑1 (MQ=255)             cGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGGAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCcg                                      >  1:207716/1‑90 (MQ=255)             cGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGGAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCcg                                      >  1:523199/1‑90 (MQ=255)             cGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGGAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCcg                                      >  1:55737/1‑90 (MQ=255)                      cGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGGAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGcgcg                             <  2:475470/90‑1 (MQ=255)                         ggTAATGGTTGCACGTGGCGACCTCGGTGTGGGAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCa                          >  1:19295/1‑90 (MQ=255)                         ggTAATGGTTGCACGTGGCGACCTCGGTGTGGGAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCa                          >  2:204220/1‑90 (MQ=255)                                  tGCACGTGGCGACCTCGGTGTGGGAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCg                 <  1:169609/90‑1 (MQ=255)                                    cACGTGGCGACCTCGGTGTGGGAATTGGCGACCCGGAACTGGTCGGCATTCAGaa                                                  <  1:301069/55‑1 (MQ=255)                                    cACGTGGCGACCTCGGTGTGGGAATTGGCGACCCGGAACTGGTCGGCATTCAGaa                                                  >  2:301069/1‑55 (MQ=255)                                                tCGGTGTGGGAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACgg   <  2:186710/90‑1 (MQ=255)                                                 cGGTGTGGGAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGc  <  2:19295/90‑1 (MQ=255)                                                         |                                                                                  TGACATCATCCTCGCCTCTGACGTGGTAATGGTTGCACGTGGCGACCTCGGTGTGGAAATTGGCGACCCGGAACTGGTCGGCATTCAGAAAGCGTTGATCCGTCGTGCGCGTCAGCTAAACCGAGCGGTAATCACGGC  >  NZ_CP009273/1932616‑1932753

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