Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,872,853	G→C	G271G (GGG→GGC)	dgcJ →	putative diguanylate cyclase DgcJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,872,853	0	G	C	100.0%	164.5 / NA	42	G271G (GGG→GGC)	dgcJ	putative diguanylate cyclase DgcJ
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (22/20);  total (22/20)

GTCGCAGGGCAGCTAAAAGGGATTGTGCTGCTGGATATTAACAAAAACAATCTGCGGAATATCTTTTATACTCATGACCGCCCTCTCCTCTGGCGTTTTCTCAATGTCACGCTAACCGATACCGATTCGGGGCGCGACATTATCATCAACCAGAGCGAAGATAATCTGTTCCAGTATGTCAGTTACGTCCATGACTTACCGGGCGGCATTCGTGTCTCGTTATCCATTGATATTCTTTACTTTATCACGTCTTCGTGGAAAAGCGT  >  NC_000913/1872722‑1872987                                                                                                                                    |                                                                                                                                       gTCGCAGGGCAGCTAAAAGGGATTGTGCTGCTGGATATTAACAAAAACAATCTGCGGAATATCTTTTATACTCATGACCGCCCTCTCCTCTGGCGTTTTCTCAATGTCACGCTAACCGATACCGATTCGGGCCGCGACAt                                                                                                                                >  1:15724/1‑140 (MQ=255) gTCGCAGGGCAGCTAAAAGGGATTGTGCTGCTGGATATTAACAAAAACAATCTGCGGAATATCTTTTATACTCATGACCGCCCTCTCCTCTGGCGTTTTCTCAATGTCACGCTAACCGATACCGATTCGGGCCGCGACAt                                                                                                                                <  2:706830/140‑1 (MQ=255) gTCGCAGGGCAGCTAAAAGGGATTGTGCTGCTGGATATTAACAAAAACAATCTGCGGAATATCTTTTATACTCATGACCGCCCTCTCCTCTGGCGTTTTCTCAATGTCACGCTAACCGATACCGATTCGGGCCGCGACAt                                                                                                                                >  1:241656/1‑140 (MQ=255) gTCGCAGGGCAGCTAAAAGGGATTGTGCTGCTGGATATTAACAAAAACAATCTGCGGAATATCTTTTATACTCATGACCGCCCTCTCCTCTGGCGTTTTCTCAATGTCACGCTAACCGATACCGATTCGGGCCGCGACAt                                                                                                                                <  1:345566/140‑1 (MQ=255) gTCGCAGGGCAGCTAAAAGGGATTGTGCTGCTGGATATTAACAAAAACAATCTGCGGAATATCTTTTATACTCATGACCGCCCTCTCCTCTGGCGTTTTCTCAATGTCACGCTAACCGATACCGATTCGGGCCGCGACAt                                                                                                                                <  2:7934/140‑1 (MQ=255) gTCGCAGGGCAGCTAAAAGGGATTGTGCTGCTGGATATTAACAAAAACAATCTGCGGAATATCTTTTATACTCATGACCGCCCTCTCCTCTGGCGTTTTCTCAATGTCACGCTAACCGATACCGATTCGGGCCGCGACAt                                                                                                                                >  1:59404/1‑140 (MQ=255)   cGCAGGGCAGCTAAAAGGGATTGTGCTGCTGGATATTAACAAAAACAATCTGCGGAATATCTTTTATACTCATGACCGCCCTCTCCTCTGGCGTTTTCTCAATGTCACGCTAACCGATACCGATTCGGGCCGCGACAt                                                                                                                                >  2:936026/1‑138 (MQ=255)   cGCAGGGCAGCTAAAAGGGATTGTGCTGCTGGATATTAACAAAAACAATCTGCGGAATATCTTTTATACTCATGACCGCCCTCTCCTCTGGCGTTTTCTCAATGTCACGCTAACCGATACCGATTCGGGCCGCGACAt                                                                                                                                <  1:936026/138‑1 (MQ=255)           aGCTAAAAGGGATTGTGCTGCTGGATATTAACAAAAACAATCTGCGGAATATCTTTTATACTCATGACCGCCCTCTCCTCTGGCGTTTTCTCAATGTCACGCTAACCGATACCGATTCGGGCCGCGACATTATCATCAAc                                                                                                                      <  1:451713/140‑1 (MQ=255)            gCTAAAAGGGATTGTGCTGCTGGATATTAACAAAAACAATCTGCGGAATATCTTTTATACTCATGACCGCCCTCTCCTCTGGCGTTTTCTCAATGTCACGCTAACCGATACCGATTCGGGCCGCGACATTATCATCAAcc                                                                                                                     >  2:893098/1‑140 (MQ=255)                       ttGTGCTGCTGGATATTAACAAAAACAATCTGCGGAATATCTTTTATACTCATGACCGCCCTCTCCTCTGGCGTTTTCTCAATGTCACGCTAACCGATACCGATTCGGGCCGCGACATTATCATCAACCAGAGCGAAGa                                                                                                           >  1:1005205/1‑139 (MQ=255)                       ttGTGCTGCTGGATATTAACAAAAACAATCTGCGGAATATCTTTTATACTCATGACCGCCCTCTCCTCTGGCGTTTTCTCAATGTCACGCTAACCGATACCGATTCGGGCCGCGACATTATCATCAACCAGAGCGAAGa                                                                                                           <  2:1005205/139‑1 (MQ=255)                        tgtgCTGCTGGATATTAACAAAAACAATCTGCGGAATATCTTTTATACTCATGACCGCCCTCTCCTCTGGCGTTTTCTCAATGTCACGCTAACCGATACCGATTCGGGCCGCGACATTATCATCAACCAGAGCGAAGATa                                                                                                         <  2:451537/140‑1 (MQ=255)                        tgtgCTGCTGGATATTAACAAAAACAATCTGCGGAATATCTTTTATACTCATGACCGCCCTCTCCTCTGGCGTTTTCTCAATGTCACGCTAACCGATACCGATTCGGGCCGCGACATTATCATCAACCAGAGCGAAGATa                                                                                                         <  2:241656/140‑1 (MQ=255)                          tgctgcTGGATATTAACAAAAACAATCTGCGGAATATCTTTTATACTCATGACCGCCCTCTCCTCTGGCGTTTTCTCAATGTCACGCTAACCGATACCGATTCGGGCCGCGACATTATCATCAACTAGAGCGAAGATAAt                                                                                                       >  1:112149/1‑140 (MQ=255)                          tgctgcTGGATATTAACAAAAACAATCTGCGGAATATCTTTTATACTCATGACCGCCCTCTCCTCTGGCGTTTTCTCAATGTCACGCTAACCGATACCGATTCGGGCCGCGACATTATCATCAACCAGAGCGAAGATAAt                                                                                                       <  1:631648/140‑1 (MQ=255)                          tgctgcTGGATATTAACAAAAACAATCTGCGGAATATCTTTTATACTCATGACCGCCCTCTCCTCTGGCGTTTTCTCAATGTCACGCTAACCGATACCGATTCGGGCCGCGACATTATCATCAACCAGAGCGAAGATAAt                                                                                                       <  1:893098/140‑1 (MQ=255)                              gctgGATATTAACAAAAACAATCTGCGGAATATCTTTTATACTCATGACCGCCCTCTCCTCTGGCGTTTTCTCAATGTCACGCTAACCGATACCGATTCGGGCCGCGACATTATCATCAACCAGAGCGAAGATAATCTGt                                                                                                   <  2:986025/140‑1 (MQ=255)                                tgGATATTAACAAAAACAATCTGCGGAATATCTTTTATACTCATGACCGCCCTCTCCTCTGGCGTTTTCTCAATGTCACGCTAACCGATACCGATTCGGGCCGCGACATTATCATCAACCAGAGCGAAGATAATCTGTTc                                                                                                 <  1:164512/140‑1 (MQ=255)                                    tatTAACAAAAACAATCTGCGGAATATCTTTTATACTCATGACCGCCCTCTCCTCTGGCGTTTTCTCAATGTCACGCTAACCGATACCGATTCGGGCCGCGACATTATCATCAACCAGAGCGAAGATAATCTGTTCCAGt                                                                                             >  2:816053/1‑140 (MQ=255)                                       tAACAAAAACAATCTGCGGAATATCTTTTATACTCATGACCGCCCTCTCCTCTGGCGTTTTCTCAATGTCACGCTAACCGATACCGATTCGGGCCGCGACATTATCATCAACTAGAGCGAAGATAATCTGTTCCAGTATg                                                                                          <  2:112149/140‑1 (MQ=255)                                            aaaaCAATCTGCGGAATATCTTTTATACTCATGACCGCCCTCTCCTCTGGCGTTTTCTCAATGTCACGCTAACCGATACCGATTCGGGCCGCGACATTATCATCAACCAGAGCGAAGATAATCTGt                                                                                                   >  2:1014868/1‑126 (MQ=255)                                            aaaaCAATCTGCGGAATATCTTTTATACTCATGACCGCCCTCTCCTCTGGCGTTTTCTCAATGTCACGCTAACCGATACCGATTCGGGCCGCGACATTATCATCAACCAGAGCGAAGATAATCTGt                                                                                                   <  1:1014868/126‑1 (MQ=255)                                                               cTTTTATACTCATGACCGCCCTCTCCTCTGGCGTTTTCTCAATGTCACGCTAACCGATACCGATTCGGGCCGCGACATTATCATCAACCAGAGCGAAGATAATCTGTTCCAGTATGTCAGTTACGTCCATGACTTACCgg                                                                  >  1:76684/1‑140 (MQ=255)                                                                      aCTCATGACCGCCCTCTCCTCTGGCGTTTTCTCAATGTCACGCTAACCGATACCGATTCGGGCCGCGACATTATCATCAACCAGAGCGAAGATAATCTGTTCCAGTATGTCAGTTACGTCCATGACTTACCGGGCGGCAt                                                           >  1:798947/1‑140 (MQ=255)                                                                            gACCGCCCTCTCCTCTGGCGTTTTCTCAATGTCACGCTAACCGATACCGATTCGGGCCGCGACATTATCATCAACCAGAGCGAAGATAATCTGTTCCAGTATGTCAGTTACGTCCATGACTTACCGGGCGGCATTCgtgt                                                     >  2:872271/1‑140 (MQ=255)                                                                                 cccTCTCCTCTGGCGTTTTCTCAATGTCACGCTAACCGATACCGATTCGGGCCGCGACATTATCATCAACCAGAGCGAAGATAATCTGTTCCAGTATGTCAGTTACGTCCATGACTTACCGGGCGGCATTCgtg                                                      >  2:405164/1‑134 (MQ=255)                                                                                 cccTCTCCTCTGGCGTTTTCTCAATGTCACGCTAACCGATACCGATTCGGGCCGCGACATTATCATCAACCAGAGCGAAGATAATCTGTTCCAGTATGTCAGTTACGTCCATGACTTACCGGGCGGCATTCgtg                                                      <  1:405164/134‑1 (MQ=255)                                                                                 cccTCTCCTCTGGCGTTTTCTCAATGTCACGCTAACCGATACCGATTCGGGCCGCGACATTATCATCAACCAGAGCGAAGATAATCTGTTCCAGTATGTCAGTTACGTCCATGACTTACCGGGCGGCATTCGTGTCTCGt                                                >  2:518073/1‑140 (MQ=255)                                                                                    tctcctctGGCGTTTTCTCAATGTCACGCTAACCGATACCGATTCGGGCCGCGACATTATCATCAACCAGAGCGAAGATAATCTGTTCCAGTATGTCAGTTACGTCCATGACTTACCGGGCGGCATTCGTGTCTCGTTAt                                             >  2:498334/1‑140 (MQ=255)                                                                                       cctctGGCGTTTTCTCAATGTCACGCTAACCGATACCGATTCGGGCCGCGACATTATCATCAACCAGAGCGAAGATAATCTGTTCCAGTATGTCAGTTACGTCCATGACTTACCGGGCGGCATTCGTGTCTCGTTATCCa                                          <  1:518073/140‑1 (MQ=255)                                                                                        ctctGGCGTTTTCTCAATGTCACGCTAACCGATACCGATTCGGGCCGCGACATTATCATCAACCAGAGCGAAGATAATCTGTTCCAGTATGTCAGTTACGTCCATGACTTACCGGGCGGCATTCGTGTCTCGTTATCCAt                                         <  2:335640/140‑1 (MQ=255)                                                                                           tGGCGTTTTCTCAATGTCACGCTAACCGATACCGATTCGGGCCGCGACATTATCATCAACCAGAGCGAAGATAATCTGTTCCAGTATGTCAGTTACGTCCATGACTTACCGGGCGGCATTCGTGTCTCGTTATCCATTg                                       <  1:777549/139‑1 (MQ=255)                                                                                           tGGCGTTTTCTCAATGTCACGCTAACCGATACCGATTCGGGCCGCGACATTATCATCAACCAGAGCGAAGATAATCTGTTCCAGTATGTCAGTTACGTCCATGACTTACCGGGCGGCATTCGTGTCTCGTTATCCATTg                                       >  2:777549/1‑139 (MQ=255)                                                                                           tGGCGTTTTCTCAATGTCACGCTAACCGATACCGATTCGGGCCGCGACATTATCATCAACCAGAGCGAAGATAATCTGTTCCAGTATGTCAGTTACGTCCATGACTTACCGGGCGGCATTCGTGTCTCGTTATCCATTGa                                      >  1:108679/1‑140 (MQ=255)                                                                                            ggCGTTTTCTCAATGTCACGCTAACCGATACCGATTCGGGCCGCGACATTATCATCAACCAGAGCGAAGATAATCTGTTCCAGTATGTCAGTTACGTCCATGACTTACCGGGCGGCATTCGTGTCTCGTTATCCATTg                                       >  2:666066/1‑138 (MQ=255)                                                                                            ggCGTTTTCTCAATGTCACGCTAACCGATACCGATTCGGGCCGCGACATTATCATCAACCAGAGCGAAGATAATCTGTTCCAGTATGTCAGTTACGTCCATGACTTACCGGGCGGCATTCGTGTCTCGTTATCCATTg                                       <  1:666066/138‑1 (MQ=255)                                                                                                  ttCTCAATGTCACGCTAACCGATACCGATTCGGGCCGCGACATTATCATCAACCAGAGCGAAGATAATCTGTTCCAGTATGTCAGTTACGTCCATGACTTACCGGGCGGCATTCGTGTCTCGTTATCCATTGATATTctt                               <  1:872271/140‑1 (MQ=255)                                                                                                            cACGCTAACCGATACCGATTCGGGCCGCGACATTATCATCAACCAGAGCGAAGATAATCTGTTCCAGTATGTCAGTTACGTCCATGACTTACCGGGCGGCATTCGTGTCTCGTTATCCATTGATATTCTTTACTTTATCa                     >  2:742852/1‑140 (MQ=255)                                                                                                             aCGCTAACCGATACCGATTCGGGCCGCGACATTATCATCAACCAGAGCGAAGATAATCTGTTCCAGTATGTCAGTTACGTCCATGACTTACCGGGCGGCATTCGTGTCTCGTTATCCATTGATATTCTTTACTTTATCAc                    >  2:10250/1‑140 (MQ=255)                                                                                                                              ttCGGGCCGCGACATTATCATCAACCAGAGCGAAGATAATCTGTTCCAGTATGTCAGTTACGTCCATGACTTACCGGGCGGCATTCGTGTCTCGTTATCCATTGATATTCTTTACTTTATCACGTCTTCGTGGAAAAGCg   >  2:212113/1‑140 (MQ=255)                                                                                                                               tCGGGCCGCGACATTATCATCAACCAGAGCGAAGATAATCTGTTCCAGTATGTCAGTTACGTCCATGACTTACCGGGCGGCATTCGTGTCTCGTTATCCATTGATATTCTTTACTTTATCACGTCTTCGTGGAAAAGCgt  >  2:647782/1‑140 (MQ=255)                                                                                                                                    |                                                                                                                                       GTCGCAGGGCAGCTAAAAGGGATTGTGCTGCTGGATATTAACAAAAACAATCTGCGGAATATCTTTTATACTCATGACCGCCCTCTCCTCTGGCGTTTTCTCAATGTCACGCTAACCGATACCGATTCGGGGCGCGACATTATCATCAACCAGAGCGAAGATAATCTGTTCCAGTATGTCAGTTACGTCCATGACTTACCGGGCGGCATTCGTGTCTCGTTATCCATTGATATTCTTTACTTTATCACGTCTTCGTGGAAAAGCGT  >  NC_000913/1872722‑1872987

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