Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	1,402,000	A→G	I18T (ATC→ACC)	dgcM ←	diguanylate cyclase DgcM

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	1,402,000	0	A	G	100.0%	38.1 / NA	13	I18T (ATC→ACC)	dgcM	diguanylate cyclase DgcM
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (7/6);  total (7/6)

CGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGT  >  NZ_CP009273/1401916‑1402082                                                                                     |                                                                                   cGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAg                                                                               >  2:333048/1‑90 (MQ=255)        gACTAAGCGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGAc                                                                        >  1:205829/1‑90 (MQ=255)                 tGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTa                                                               <  1:333048/90‑1 (MQ=255)                       cTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCa                                                         >  2:136196/1‑90 (MQ=255)                        ttGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAg                                                        <  1:149732/90‑1 (MQ=255)                        ttGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAg                                                        <  2:342354/90‑1 (MQ=255)                           cATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAGGGt                                                     <  1:156827/90‑1 (MQ=255)                           cATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAGGGt                                                     <  2:205829/90‑1 (MQ=255)                                                     cATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGCGTAATCATTgaga                           >  1:222444/1‑90 (MQ=255)                                                        aaaTTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGtt                                            >  1:88601/1‑70 (MQ=255)                                                        aaaTTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGtt                                            <  2:88601/70‑1 (MQ=255)                                                             aaCTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTa                   >  2:13260/1‑90 (MQ=255)                                                                              cGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTttgt  >  2:29206/1‑90 (MQ=255)                                                                                     |                                                                                   CGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGT  >  NZ_CP009273/1401916‑1402082

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