Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,873,384	A→G	K448K (AAA→AAG)	dgcJ →	putative diguanylate cyclase DgcJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,873,384	0	A	G	100.0%	188.9 / NA	53	K448K (AAA→AAG)	dgcJ	putative diguanylate cyclase DgcJ
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (25/28);  total (25/28)

CGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACAAA  >  NC_000913/1873248‑1873516                                                                                                                                         |                                                                                                                                     cGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGt                                                                                                                                   <  1:412836/140‑1 (MQ=255) agTTATTAGCTCAGGCGCTTAAGCAGTCGATACGTAAAACTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATAATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGt                                                                                                                                   <  2:333986/139‑1 (MQ=255)     attaGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGc                                                                                                                               >  2:723874/1‑140 (MQ=255)             aGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATc                                                                                                                       >  2:477101/1‑140 (MQ=255)              ggCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATcg                                                                                                                      >  1:483597/1‑140 (MQ=255)              ggCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATAATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGAGTCTGCAACatatgt                                                                                                                      >  1:333986/1‑138 (MQ=255)               gCGATTAAACAGTCGATTCGTATAACTGATTATGCCAACCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATcgc                                                                                                                     <  1:340187/140‑1 (MQ=255)               gCGATTAAACAGTCGATCCGTACAAGTGATTTTGCTATCAGACTCGGTCGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATcgc                                                                                                                     <  2:194155/140‑1 (MQ=255)                 gATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCgcgc                                                                                                                   <  1:689336/140‑1 (MQ=255)                   ttAAACAGTCGATTCGTGAAAGTGATTATGCAATCCGACCCGGTGGCGTTGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGcgccg                                                                                                                 <  2:188778/140‑1 (MQ=255)                   ttAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGcgccg                                                                                                                 <  2:615855/140‑1 (MQ=255)                    tAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGcgccgc                                                                                                                <  1:81258/140‑1 (MQ=255)                           tCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAaga                                                                                                         >  1:419302/1‑140 (MQ=255)                           tCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTAGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGa                                                                                                             >  2:70045/1‑136 (MQ=255)                           tCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTAGATTCGACGCCGCAAATTGCAGCACAACTGCCTAAACGTATCGAAAAGCGTCCGCAACATATCGCGCCGCAGa                                                                                                             <  1:70045/136‑1 (MQ=255)                             gATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCATCTCCAACATATCCCGCCGCAGAAAgaga                                                                                                       >  1:313966/1‑140 (MQ=255)                                           aTTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTcc                                                                                         >  1:580030/1‑140 (MQ=255)                                            ttATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCg                                                                                        >  2:457376/1‑140 (MQ=255)                                                 ccATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGACTGCAACATATCACGCGGCAGAAAGGGATCTGCTTCTCTTCCGGTAtt                                                                                   >  2:216607/1‑140 (MQ=255)                                                   aTCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCGGCAACATAGCGCTCCGAAGAAAGAGGTCGGTTTCTCTTCCGGTATTTa                                                                                 >  2:258874/1‑140 (MQ=255)                                                    tCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTtc                                                                                               >  2:165678/1‑125 (MQ=255)                                                    tCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTtc                                                                                               <  1:165678/125‑1 (MQ=255)                                                    tCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTAc                                                                                >  1:628379/1‑140 (MQ=255)                                                    tCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTAc                                                                                >  2:706929/1‑140 (MQ=255)                                                    tCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTCCAACATATGGCGGCCCAGAAAGAGAGCGGCGTCTCTTCCCGTATTTAc                                                                                >  2:358644/1‑140 (MQ=255)                                                       gACTCGGTGGCGATGGATTCTGTATGATTGTTGTCGACTCGACGCCGCAAATTGCAGCACAACCGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTAcgcg                                                                             <  1:243693/140‑1 (MQ=255)                                                        aCTCGGTGGCGATGAGTTCTGCATCATTCTTGTCGATACGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGa                                                                            <  1:338369/140‑1 (MQ=255)                                                         cTCGGTGGGGATGAATGCTGCACCATTCTTCTCCATTCGGCGCCGCAATTTGCAGCACAACTGCCTGAACGTACCGAAAAGCGTCTGCAACACATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGAt                                                                           <  2:357990/140‑1 (MQ=255)                                                         cTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGAt                                                                           >  2:557116/1‑140 (MQ=255)                                                         cTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGAt                                                                           >  2:696311/1‑140 (MQ=255)                                                           cGGTGGCGATGAATTTTGCATCCTTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATga                                                                         <  1:317095/140‑1 (MQ=255)                                                           cGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATga                                                                         <  2:628379/140‑1 (MQ=255)                                                            ggTGGCGAGGAATTCTCTATCACTCTTGTCCTTTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATgaa                                                                        <  2:355818/140‑1 (MQ=255)                                                                     gAATTCTGCATCATTCTTGTCGCTTCGACACCGCAAACTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGa                                                               <  1:195500/140‑1 (MQ=255)                                                                      aaTTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGAt                                                              <  1:464813/140‑1 (MQ=255)                                                                        ttCTGCATCATTCTTCTCGATTCGAGGCCGCAAATTGGAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCACAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATAc                                                            <  1:358644/140‑1 (MQ=255)                                                                              atcatTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACa                                                      >  1:448677/1‑140 (MQ=255)                                                                               tcatTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACat                                                     <  2:526704/140‑1 (MQ=255)                                                                                catTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACatg                                                    >  2:178367/1‑140 (MQ=255)                                                                                 atTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACatga                                                   <  1:723874/140‑1 (MQ=255)                                                                                       gTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGta                                             >  1:558462/1‑140 (MQ=255)                                                                                          gATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATaa                                          >  1:632945/1‑140 (MQ=255)                                                                                             tCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGc                                       >  2:573499/1‑140 (MQ=255)                                                                                                    cgcAAATTGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGACACGTTACATGATGCGCATAAAGCTTCCGAt                                <  2:448677/140‑1 (MQ=255)                                                                                                          ttGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGt                          <  2:525427/140‑1 (MQ=255)                                                                                                          ttGCAGCACAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGt                          <  1:557116/140‑1 (MQ=255)                                                                                                              agcaCAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTtat                      <  1:477101/140‑1 (MQ=255)                                                                                                              agcaCAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTtat                      <  2:419302/140‑1 (MQ=255)                                                                                                              agcaCAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTtat                      >  2:76219/1‑140 (MQ=255)                                                                                                               gcaCAACTGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTtata                     <  1:550281/140‑1 (MQ=255)                                                                                                                      tGCCTGAACGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAAt              <  1:76219/140‑1 (MQ=255)                                                                                                                            aaCGTATCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACAGGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAg        >  2:486313/1‑140 (MQ=255)                                                                                                                                  tCGAAAAGCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGAGGCGTATAACGCTGCCGGTGAGCGTTTATATGTCAATAAGCAGAACaaa  >  2:249137/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                     CGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACAAA  >  NC_000913/1873248‑1873516

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