Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	1,558,439	A→G	Y518H (TAT→CAT)	dosP ←	oxygen‑sensing cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	1,558,439	0	A	G	100.0%	27.6 / NA	10	Y518H (TAT→CAT)	dosP	oxygen‑sensing cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (1/9);  total (1/9)

TGCGCCTAAAACCAAACGCTCTTTTACCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATACTC  >  NZ_CP009273/1558353‑1558512                                                                                       |                                                                          tGCGCCTAAACCCAAACGCTTTTTTCCCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATc                                                                        <  2:315382/90‑1 (MQ=255) tGCGCCAAAACCCAAACGCTTTTTTCCCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATc                                                                        <  2:350391/90‑1 (MQ=255)                ccGCTTTTTTCCCTTTTCTTTCACCGCCGGGCTGAAGACCTGCCAGCCGTTCCCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGc                                                         <  2:260788/89‑1 (MQ=255)                aCGCTCTTTTACCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGc                                                         <  2:74373/90‑1 (MQ=255)                                     cATCGCGGGGCTAAAGACCGGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATcgcg                                    <  1:139356/90‑1 (MQ=255)                                                 gaagaaCTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGc                                                         <  1:164111/57‑1 (MQ=255)                                                 gaagaaCTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGc                                                         >  2:164111/1‑57 (MQ=255)                                                                  tACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCaa       <  1:413980/90‑1 (MQ=255)                                                                       ccATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATACTc  <  2:478088/90‑1 (MQ=255)                                                                       ccATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGGAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATACTc  <  1:385120/90‑1 (MQ=255)                                                                                       |                                                                          TGCGCCTAAAACCAAACGCTCTTTTACCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATACTC  >  NZ_CP009273/1558353‑1558512

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