Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	1,560,820	T→C	E193G (GAA→GGA)	dosC ←	diguanylate cyclase DosC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	1,560,820	0	T	C	100.0%	28.8 / NA	10	E193G (GAA→GGA)	dosC	diguanylate cyclase DosC
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (7/3);  total (7/3)

AACCACAGGCCAAAATCAGCCTGGCTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTTCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGTTTTCCAGTAACGAGAAGATACGATAGTTTTCATC  >  NZ_CP009273/1560736‑1560902                                                                                     |                                                                                   aaCCACAGGCCAAAATCAGCCTGGCTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTCCCCAt                                                                               >  2:405905/1‑90 (MQ=255)                       ggCTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTCCCCATGAAAGTATTGAGGCTATTTGcc                                                         <  2:258725/90‑1 (MQ=255)                                 cAAACTACTGCCTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTCCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTttc                                               >  2:298223/1‑90 (MQ=255)                                                                        gATAATATCTATTCCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGTTTTCCAGTAACgag                      >  1:326768/1‑76 (MQ=255)                                                                        gATAATATCTATTCCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGTTTTCCAGTAACgag                      <  2:326768/76‑1 (MQ=255)                                                                          taataTCTATTCCCCATGAAGGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGTTTTCCAGTAACGAGAAGATACGATAGtttt      >  2:229784/1‑90 (MQ=255)                                                                          taataTCTATTCCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGTTTTCCAGTAACGAGAAGATACGATAGtttt      >  1:89913/1‑90 (MQ=255)                                                                          taataTCTATTCCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGTTTTCCAGTAACGAGAAGATACGATAGtttt      >  2:398567/1‑90 (MQ=255)                                                                          taataTCTATTCCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTCTCTTCTTCGGCGTTTTCCAGTAACGAGAAGATACGATAGtttt      >  2:49717/1‑90 (MQ=255)                                                                              atCTATTCCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGTTTTCCAGTAACGAGAAGATACGATAGTTTTCATc  <  2:259640/90‑1 (MQ=255)                                                                                     |                                                                                   AACCACAGGCCAAAATCAGCCTGGCTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGATTTTATAGATAATATCTATTTCCCATGAAAGTATTGAGGCTATTTGCCGTTCTTTTTCTTCTTCGGCGTTTTCCAGTAACGAGAAGATACGATAGTTTTCATC  >  NZ_CP009273/1560736‑1560902

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