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%	146.3 / NA	42	I18T (ATC→ACC)	dgcM	diguanylate cyclase DgcM
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (14/28);  total (14/28)

GTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTA  >  NZ_CP009273/1401913‑1402065                                                                                        |                                                                  gTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTc                                                                 <  1:1067795/90‑1 (MQ=255) gTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTc                                                                 <  2:342977/90‑1 (MQ=255) gTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTc                                                                 <  1:791519/90‑1 (MQ=255) gTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTc                                                                 <  1:967029/90‑1 (MQ=255)      tAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGAc                                                            >  2:639404/1‑90 (MQ=255)      tAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGAc                                                            <  2:276920/90‑1 (MQ=255)      tAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGAc                                                            >  2:836303/1‑90 (MQ=255)         cTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAgg                                                         <  1:816580/90‑1 (MQ=255)           gACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGAc                                                       >  2:760597/1‑90 (MQ=255)           gACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGAc                                                       >  2:325082/1‑90 (MQ=255)              tAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTgg                                                    <  1:191579/90‑1 (MQ=255)                   gTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGt                                               <  1:670174/90‑1 (MQ=255)                    tGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTa                                              <  2:216971/90‑1 (MQ=255)                    tGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTa                                              <  1:760597/90‑1 (MQ=255)                    tGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTa                                              <  1:301604/90‑1 (MQ=255)                    tGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTa                                              <  2:738237/90‑1 (MQ=255)                     gAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTaa                                             >  1:71783/1‑90 (MQ=255)                     gAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTaa                                             >  1:704187/1‑90 (MQ=255)                           ttGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAg                                       <  2:803920/90‑1 (MQ=255)                           ttGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAg                                       >  2:879857/1‑90 (MQ=255)                           ttGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAg                                       <  1:836303/90‑1 (MQ=255)                           ttGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAg                                       <  1:78554/90‑1 (MQ=255)                                   aCAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGa                               <  1:1175431/90‑1 (MQ=255)                                     aGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAg                             <  1:894480/90‑1 (MQ=255)                                      gTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGt                            >  1:664814/1‑90 (MQ=255)                                           gccgcGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGc                       <  1:1155842/90‑1 (MQ=255)                                             cgcgCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGCGt                     <  2:764043/90‑1 (MQ=255)                                             cgcgCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGCGt                     <  2:71783/90‑1 (MQ=255)                                             cgcgCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGCGt                     <  1:335027/90‑1 (MQ=255)                                              gcgcTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGCGTa                    <  1:325082/90‑1 (MQ=255)                                                         aTAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGAt         <  2:798913/90‑1 (MQ=255)                                                                aaCTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTa                                              >  2:381007/1‑46 (MQ=255)                                                                aaCTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTa                                              <  1:381007/46‑1 (MQ=255)                                                                aaCTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAg                                       <  1:1170917/53‑1 (MQ=255)                                                                aaCTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAg                                       >  2:1170917/1‑53 (MQ=255)                                                                aaCTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTa  >  2:236976/1‑90 (MQ=255)                                                                     ttCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGtt                           <  1:197825/60‑1 (MQ=255)                                                                     ttCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGtt                           >  2:197825/1‑60 (MQ=255)                                                                      tCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTcc                                         >  2:912777/1‑45 (MQ=255)                                                                      tCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTcc                                         <  1:912777/45‑1 (MQ=255)                                                                              gggCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGc                       >  2:1022766/1‑55 (MQ=255)                                                                              gggCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGc                       <  1:1022766/55‑1 (MQ=255)                                                                                        |                                                                  GTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTA  >  NZ_CP009273/1401913‑1402065

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