Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,143,684	(G)6→7	coding (419/3318 nt)	dgcE →	putative diguanylate cyclase DgcE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,143,678	1	.	G	100.0%	196.9 / NA	53	L138W (TTG→TGG)	dgcE	putative diguanylate cyclase DgcE
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (28/25);  total (28/25)

CCTGGACGACCATCAATATTGTTGAAGCCGTGGTCGGGGCAGTGCTGCTACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTT‑GGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCTGGGATTGTTATTTAAACC  >  NC_000913/2143540‑2143801                                                                                                                                            |                                                                                                                            ccTGGACGACCATCAATATTGTTGAAGCCGTGGTCGGGGCAGTGCTGCTACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTT‑g                                                                                                                            <  2:758636/140‑1 (MQ=255)    ggACGACCATCAATATTGTTGAAGCCGTGGTCGGGGCAGTGCTGCTACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTT‑gggg                                                                                                                         >  2:822753/1‑140 (MQ=255)     gacgacCATCAATATTGTTGAAGCCGTGGTCGGGGCAGTGCTGCTACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTT‑ggggg                                                                                                                        <  2:818858/140‑1 (MQ=255)         acCATCAATATTGTTGAAGCCGTGGTCGGGGCAGTGCTGCTACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTT‑gggggggtg                                                                                                                    <  1:1019896/140‑4 (MQ=255)                  atTGTTGAAGCCGTGGTCGGGGCAGTGCTGCTACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTg                                                                                                            <  1:1098790/140‑1 (MQ=255)                    tgttgAAGCCGTGGTCGGGGCAGTGCTGCTACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTc                                                                                                          >  1:1061927/1‑140 (MQ=255)                    tgttgAAGCCGTGGTCGGGGCAGTGCTGCTACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTc                                                                                                          >  1:1092881/1‑140 (MQ=255)                                         aGTGCTGCTACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGAt                                                                                     <  2:908042/140‑1 (MQ=255)                                         aGTGCTGCTACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGAt                                                                                     >  1:733847/1‑140 (MQ=255)                                           tgctgcTACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATcc                                                                                   >  1:1204195/1‑140 (MQ=255)                                           tgctgcTACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATcc                                                                                   >  1:189037/1‑140 (MQ=255)                                           tgctgcTACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATcc                                                                                   >  2:833434/1‑140 (MQ=255)                                            gctgctACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCct                                                                                  >  1:777403/1‑140 (MQ=255)                                                ctACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCa                                                                              <  2:189037/140‑1 (MQ=255)                                                  aCGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGAt                                                                                     >  2:10334/1‑131 (MQ=255)                                                  aCGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGAt                                                                                     <  1:10334/131‑1 (MQ=255)                                                  aCGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAgg                                                                            <  2:1092881/140‑1 (MQ=255)                                                     tAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCa                                                                         <  2:1061927/140‑1 (MQ=255)                                                       aaTTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCAtt                                                                       <  2:1204032/140‑1 (MQ=255)                                                            ctgcCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGa                                                                  <  1:1271156/140‑1 (MQ=255)                                                             tgcCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGat                                                                 >  1:1183742/1‑140 (MQ=255)                                                             tgcCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGat                                                                 <  2:733847/140‑1 (MQ=255)                                                             tgcCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGat                                                                 >  2:1226450/1‑140 (MQ=255)                                                                  tGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATggg                                                            <  2:1183742/140‑1 (MQ=255)                                                                     tataACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTAc                                                         <  1:1226450/140‑1 (MQ=255)                                                                              ttGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCTCCATTGTTCCGCCTCTGTTGGGGGGGAGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAAt                                                >  2:959736/1‑140 (MQ=255)                                                                              ttGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATGTGGGTACTGTCAGAAt                                                >  2:1175733/1‑140 (MQ=255)                                                                              ttGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAAt                                                >  1:115639/1‑140 (MQ=255)                                                                              ttGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAAt                                                >  1:418072/1‑140 (MQ=255)                                                                              ttGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTAAGGGCATTTTTGATATGGGTACTGTCAGAAt                                                >  2:426282/1‑140 (MQ=255)                                                                               tGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATc                                               >  1:873341/1‑140 (MQ=255)                                                                               tGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATc                                               >  1:1021900/1‑140 (MQ=255)                                                                               tGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATc                                               >  1:1109822/1‑140 (MQ=255)                                                                               tGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATc                                               >  2:413181/1‑140 (MQ=255)                                                                               tGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATc                                               >  1:1281588/1‑140 (MQ=255)                                                                               tGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATc                                               >  2:119526/1‑140 (MQ=255)                                                                                   aaaTCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATc                                           <  2:781982/140‑1 (MQ=255)                                                                                     aTCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCgg                                         <  2:115639/140‑1 (MQ=255)                                                                                     aTCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCgg                                         <  1:1281714/140‑1 (MQ=255)                                                                                     aTCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCGCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGCCGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCgg                                         <  2:1021900/140‑1 (MQ=255)                                                                                          gCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGctc                                    >  1:1097228/1‑140 (MQ=255)                                                                                                        tCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCg                      >  1:893931/1‑140 (MQ=255)                                                                                                        tCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGACCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCa                               >  1:163947/1‑131 (MQ=255)                                                                                                        tCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGACCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCa                               <  2:163947/131‑1 (MQ=255)                                                                                                          tGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGc                        <  1:456671/136‑1 (MQ=255)                                                                                                          tGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGc                        >  2:456671/1‑136 (MQ=255)                                                                                                             cACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCTggg                 <  2:1097228/140‑1 (MQ=255)                                                                                                                tCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCTGGGAtt              <  1:833434/140‑1 (MQ=255)                                                                                                                tCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCTGGGAtt              <  1:119526/140‑1 (MQ=255)                                                                                                                 cGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCTGGGATTg             >  2:631844/1‑140 (MQ=255)                                                                                                                     aGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCTGGGATTGTTAt         <  2:777403/140‑1 (MQ=255)                                                                                                                       cgcCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCTGGGATTGTTAttt       >  2:22155/1‑140 (MQ=255)                                                                                                                       cgcCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCTGGGATTGTTAttt       >  1:1000758/1‑140 (MQ=255)                                                                                                                        gcCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATGTGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCTGGGATTGTTATTTa      <  1:1175733/140‑1 (MQ=255)                                                                                                                          cATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCTGGGATTGTTATTTaaa    >  2:779331/1‑140 (MQ=255)                                                                                                                          cATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCTGGGATTGTTATTTaaa    <  1:224652/140‑1 (MQ=255)                                                                                                                           aTTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCTGGGATTGTTATTTAAAc   <  2:1281588/140‑1 (MQ=255)                                                                                                                            ttGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCTGGGATTGTTATTTAAAcc  <  1:631844/140‑1 (MQ=255)                                                                                                                                            |                                                                                                                            CCTGGACGACCATCAATATTGTTGAAGCCGTGGTCGGGGCAGTGCTGCTACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTT‑GGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCACTGGTGCCGCTGGGATTGTTATTTAAACC  >  NC_000913/2143540‑2143801

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