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

GTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCT  >  NZ_CP009273/1401913‑1402086                                                                                        |                                                                                       gTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTc                                                                                      <  2:394436/90‑1 (MQ=255) gTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTc                                                                                      <  2:2331690/90‑1 (MQ=255)    cGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAg                                                                                   >  1:103746/1‑90 (MQ=255)    cGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAg                                                                                   <  2:256615/90‑1 (MQ=255)     gTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGa                                                                                  >  2:1711956/1‑90 (MQ=255)      tAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGAc                                                                                 <  1:1264093/90‑1 (MQ=255)      tAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGAc                                                                                 <  1:421359/90‑1 (MQ=255)      tAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGAc                                                                                 <  1:1482700/90‑1 (MQ=255)           gACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGAc                                                                            >  1:1803075/1‑90 (MQ=255)           gACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGAc                                                                            >  1:1735634/1‑90 (MQ=255)           gACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGAc                                                                            >  2:1463930/1‑90 (MQ=255)           gACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGAc                                                                            >  1:483123/1‑90 (MQ=255)              tAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTgg                                                                         <  2:1385038/90‑1 (MQ=255)               aaaCGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGt                                                                        <  1:1463930/90‑1 (MQ=255)               aaaCGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGt                                                                        <  1:1431941/90‑1 (MQ=255)               aaaCGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGt                                                                        <  1:121212/90‑1 (MQ=255)                   gTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGt                                                                    <  2:2024387/90‑1 (MQ=255)                    tGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTa                                                                   <  2:2447573/90‑1 (MQ=255)                    tGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTa                                                                   <  1:497124/90‑1 (MQ=255)                    tGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTa                                                                   <  1:426049/90‑1 (MQ=255)                    tGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTa                                                                   <  2:483123/90‑1 (MQ=255)                    tGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTa                                                                   <  1:379157/90‑1 (MQ=255)                    tGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTa                                                                   <  1:872727/90‑1 (MQ=255)                    tGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTa                                                                   <  1:2847012/90‑1 (MQ=255)                    tGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTa                                                                   <  1:1744346/90‑1 (MQ=255)                    tGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTa                                                                   <  1:2407150/90‑1 (MQ=255)                    tGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTa                                                                   <  1:1834293/90‑1 (MQ=255)                     gAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTaa                                                                  >  1:1692466/1‑90 (MQ=255)                     gAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTaa                                                                  >  1:2523604/1‑90 (MQ=255)                     gAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTaa                                                                  >  2:591703/1‑90 (MQ=255)                     gAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTaa                                                                  >  2:1903763/1‑90 (MQ=255)                          cTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCa                                                             <  1:362457/90‑1 (MQ=255)                           ttGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAg                                                                                   <  2:791558/67‑1 (MQ=255)                           ttGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAg                                                                                   >  1:791558/1‑67 (MQ=255)                           ttGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAg                                                            >  1:1820669/1‑90 (MQ=255)                           ttGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAg                                                            <  1:749636/90‑1 (MQ=255)                           ttGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAg                                                            <  2:2560369/90‑1 (MQ=255)                           ttGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAg                                                            <  1:1823083/90‑1 (MQ=255)                           ttGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAg                                                            <  2:2892907/90‑1 (MQ=255)                                  aaCAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAg                                                            <  2:2552457/83‑1 (MQ=255)                                  aaCAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAg                                                            >  1:2552457/1‑83 (MQ=255)                                  aaCAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTg                                                     >  1:1249364/1‑90 (MQ=255)                                  aaCAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTg                                                     >  1:967806/1‑90 (MQ=255)                                  aaCAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTg                                                     >  1:557840/1‑90 (MQ=255)                                      gTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAg                                                            >  2:1367676/1‑79 (MQ=255)                                      gTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAg                                                            <  1:1367676/79‑1 (MQ=255)                                        ttCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTtg                                               >  2:603876/1‑90 (MQ=255)                                           gccgcGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGc                                            <  1:1321317/90‑1 (MQ=255)                                           gccgcGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGc                                            <  2:1735634/90‑1 (MQ=255)                                             cgcgCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGCGt                                          <  1:1711956/90‑1 (MQ=255)                                                        cATAAATTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGCGTAATCATTgaga                               >  1:693791/1‑90 (MQ=255)                                                           aaaTTAACTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATcc                            >  1:2702348/1‑90 (MQ=255)                                                                aaCTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGt                                                 >  2:3019819/1‑64 (MQ=255)                                                                aaCTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGt                                                 <  1:3019819/64‑1 (MQ=255)                                                                aaCTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTa                       >  1:1018691/1‑90 (MQ=255)                                                                aaCTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTa                       >  2:1395303/1‑90 (MQ=255)                                                                aaCTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTa                       >  1:2863525/1‑90 (MQ=255)                                                                aaCTGTTCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAG‑AAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTAt                      >  2:2200862/1‑90 (MQ=255)                                                                      tCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTa                 >  2:471756/1‑90 (MQ=255)                                                                      tCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTa                 >  1:1301224/1‑90 (MQ=255)                                                                      tCCTCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTa                 >  2:731768/1‑90 (MQ=255)                                                                         tCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTa                                                                   >  2:2957678/1‑37 (MQ=255)                                                                         tCAAAGGGCGAAACGGTCCAGACAGGACTGGTGAGTa                                                                   <  1:2957678/37‑1 (MQ=255)                                                                               ggCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAAttttt        <  2:1161638/90‑1 (MQ=255)                                                                                gCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTtgt                                              <  1:1445572/51‑1 (MQ=255)                                                                                gCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTtgt                                              >  2:1445572/1‑51 (MQ=255)                                                                                gCGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTttg       >  1:969452/1‑90 (MQ=255)                                                                                 cGAAACGGTCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTttgt      <  2:2934295/90‑1 (MQ=255)                                                                                     aCGGTCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCt  <  2:681283/90‑1 (MQ=255)                                                                                        |                                                                                       GTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCT  >  NZ_CP009273/1401913‑1402086

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: ‑