Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	1,673,542	A→G	K287E (AAG→GAG)	ydgH →	DUF1471 family protein YdgH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	1,673,542	0	A	G	88.9%	11.2 / ‑3.6	9	K287E (AAG→GAG)	ydgH	DUF1471 family protein YdgH
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (1/7);  total (2/7)
Fisher's exact test for biased strand distribution p-value = 2.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.99e-01

GGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATC  >  NZ_CP009273/1673460‑1673617                                                                                   |                                                                            ggTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAgg                                                                      >  1:22836/1‑90 (MQ=255)     ccGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCt                                                                  <  2:37771/90‑1 (MQ=255)      cGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTa                                                                 <  2:219998/90‑1 (MQ=255)       gTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTaa                                                                <  1:175094/90‑1 (MQ=255)                                       gCAACATGCCCGAAGTCTCCTATCAGGTTGCGAAAGGGGCCGGGGAGAAGGGTGCTAAGTATTACCACATCCCCCGCCAGTGGCAGGAAc                                <  2:91363/90‑1 (MQ=255)                                         aaCATGACCGAAGTCCCCTATCGGGTTGGAAAGCGTGCCGCGGAGAAAGGTGCTAAGTTTTCCCACACCCCCCGCCAGTGGCAGGAACGt                              <  2:170550/90‑1 (MQ=255)                                         aaCAGGCCCGAAGTCTCCTATCGGGTTGCGAAACGGGCCGGGGAGAAGGGTGCTAAGTATTCCCACACCCCCCGCCAGTGGCAGGAACGt                              <  2:153234/90‑1 (MQ=255)                                                                   tGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGa    >  1:284539/1‑90 (MQ=255)                                                                     cGAAACGTGCCGCGGAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATc  <  1:143645/90‑1 (MQ=255)                                                                                   |                                                                            GGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATC  >  NZ_CP009273/1673460‑1673617

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