Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,877,900	G→A	M62I (ATG→ATA)	dgcP →	diguanylate cyclase DgcP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,877,900	0	G	A	100.0%	161.3 / NA	46	M62I (ATG→ATA)	dgcP	diguanylate cyclase DgcP
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (30/16);  total (30/16)

GTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGTCTCGTGGGATTACTCGTTATGCAAACGCGCCATTGATGAAAACTGCTTTTTCAGCGATGAAGTCCCCGACCGTTGGGGTGACTGTATTGCGGCACGCAATCTTGGCATCACC  >  NC_000913/1877765‑1878035                                                                                                                                        |                                                                                                                                        gTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATATACa                                                                                                                                     <  2:40242/140‑1 (MQ=255)      tGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATATACATCCCg                                                                                                                                >  2:700523/1‑140 (MQ=255)      tGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATATACATCCCg                                                                                                                                <  2:38460/140‑1 (MQ=255)         aaaGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGCGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATATACATCCCagag                                                                                                                             <  2:947616/140‑1 (MQ=255)            gTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATATACATCCCGGAGAAt                                                                                                                          >  2:534200/1‑140 (MQ=255)             tCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATATACATCCCGGAGAAtt                                                                                                                         >  1:299130/1‑140 (MQ=255)                ggTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATATACATCCCGGAGAATTTTa                                                                                                                      >  2:559885/1‑140 (MQ=255)                      aCAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATATACATCCCGGAGAATTTTACCGtct                                                                                                                >  1:278373/1‑140 (MQ=255)                              tGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATATACATCCCGGAGAATTTTACCGTCTCGTGGGAt                                                                                                        <  1:534200/140‑1 (MQ=255)                                gagaTGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATATACATCCCGGAGAATTTTACCGTCTCGTGGGATTa                                                                                                      >  1:250039/1‑140 (MQ=255)                                     gCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATATACATCCCGGAGAATTTTACCGTCTCGTGGGATTACTCGt                                                                                                 <  1:214598/140‑1 (MQ=255)                                     gCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATATACATCCCGGAGAATTTTACCGTCTCGTGGGATTACTCGt                                                                                                 >  2:546555/1‑140 (MQ=255)                                     gCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATATACATCCCGGAGAATTTTACCGTCTCGTGGGATTACTCGt                                                                                                 >  2:214598/1‑140 (MQ=255)                                     gCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATATACATCCCGGAGAATTTTACCGTCTCGTGGGATTACTCGt                                                                                                 >  1:291572/1‑140 (MQ=255)                                     gCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATATACATCCCGGAGAATTTTACCGTCTCGTGGGATTACTCGt                                                                                                 <  1:546555/140‑1 (MQ=255)                                        ggAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATATACATCCCGGAGAATTTTACCGTCTCGTGGGATTACTCGTTAt                                                                                              >  2:775416/1‑140 (MQ=255)                                               gTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATATACATCCCGGAGAATTTTACCGTCTCGTGGGATTACTCGTTATGCAAAcg                                                                                       <  2:278373/140‑1 (MQ=255)                                                tCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATATACATCCCGGAGAATTTTACCGTCTCGTGGGATTACTCGTTATGCAAAcgc                                                                                      >  2:785454/1‑140 (MQ=255)                                                  aCTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATATACATCCCGGAGAATTTTACCGTCTCGTGGGATTACTCGTTATGCAAACgcgc                                                                                    >  2:696802/1‑140 (MQ=255)                                                         tGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATATACATCCCGGAGAATTTTACCGTCTCGTGGGATTACTCGTTATGCAAACGCGCCATtgat                                                                             >  2:738686/1‑140 (MQ=255)                                                              tCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATATACATCCCGGAGAATTTTACCGTCTCGTGGGATTACTCGTTATGCAAACGCGCCATTGATGaaaa                                                                        >  2:887427/1‑140 (MQ=255)                                                              tCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATATACATCCCGGAGAATTTTACCGTCTCGTGGGATTACTCGTTATGCAAACGCGCCATTGATGaaaa                                                                        >  1:225834/1‑140 (MQ=255)                                                              tCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATATACATCCCGGAGAATTTTACCGTCTCGTGGGATTACTCGTTATGCAAACGCGCCATTGATGaaaa                                                                        >  1:351145/1‑140 (MQ=255)                                                                aaCCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATATACATCCCGGAGAATTTTACCGTCTCGTGGGATTACTCGTTATGCAAACGCGCCATTGATGAAAACt                                                                      >  2:614841/1‑140 (MQ=255)                                                                       ctGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATATACATCCCGGAGAATTTTACCGTCTCGTGGGATTACTCGTTATGCAAACGCGCCATTGATGAAAACTGCttttt                                                               <  1:696802/140‑1 (MQ=255)                                                                            cAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATATACATCCCGGAGAATTTTACCGTCTCGTGGGATTACTCGTTATGCAAACGCGCCATTGATGAAAACTGCTTTTTCAGc                                                           <  1:160196/139‑1 (MQ=255)                                                                            cAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATATACATCCCGGAGAATTTTACCGTCTCGTGGGATTACTCGTTATGCAAACGCGCCATTGATGAAAACTGCTTTTTCAGc                                                           >  2:160196/1‑139 (MQ=255)                                                                               aGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATATACATCCCGGAGAATTTTACCGTCTCGTGGGATTACTCGTTATGCAAACGCGCCATTGATGAAAACTGCTTTTTCAGCGATg                                                       <  1:795653/140‑1 (MQ=255)                                                                                     tGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATATACATCCCGGAGAATTTTACCGTCTCGTGGGATTACTCGTTATGCAAACGCGCCATTGATGAAAACTGCTTTTTCAGCGATGAAGTcc                                                 >  2:573818/1‑140 (MQ=255)                                                                                       tCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATATACATCCCGGAGAATTTTACCGTCTCGTGGGATTACTCGTTATGCAAACGCGCCATTGATGAAAACTGCTTTTTCAGCGATGAAGTcccc                                               <  1:785454/140‑1 (MQ=255)                                                                                       tCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATATACATCCCGGAGAATTTTACCGTCTCGTGGGATTACTCGTTATGCAAACGCGCCATTGATGAAAACTGCTTTTTCAGCGATGAAGTcccc                                               <  2:1002715/140‑1 (MQ=255)                                                                                       tCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATATACATCCCGGAGAATTTTACCGTCTCGTGGGATTACTCGTTATGCAAACGCGCCATTGATGAAAACTGCTTTTTCAGCGATGAAGTcccc                                               <  2:622074/140‑1 (MQ=255)                                                                                         gAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATATACATCCCGGAGAATTTTACCGTCTCGTGGGATTACTCGTTATGCAAACGCGCCATTGATGAAAACTGCTTTTTCAGCGATGAAGTCCCCGa                                             >  2:887356/1‑140 (MQ=255)                                                                                             cgcgcCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATATACATCCCGGAGAATTTTACCGTCTCGTGGGATTACTCGTTATGCAAACGCGCCATTGATGAAAACTGCTTTTTCAGCGATGAAGTCCCCGACCGt                                         >  1:625680/1‑140 (MQ=255)                                                                                                     cagcaTATTATGTTTGCCCGTAACAGCCAGAAAATATACATCCCGGAGAATTTTACCGTCTCGTGGGATTACTCGTTATGCAAACGCGCCATTGATGAAAACTGCTTTTTCAGCGATGAAGTCCCCGACCGTTGGGGTGa                                 <  1:738686/140‑1 (MQ=255)                                                                                                                    gCCCGTAACAGCCAGAAAATATACATCCCGGAGAATTTTACCGTCTCGTGGGATTACTCGTTATGCAAACGCGCCATTGATGAAAACTGCTTTTTCAGCGATGAAGTCCGCGACCGTTGGGGTGACTGTATTGCGGCACg                  >  2:640202/1‑140 (MQ=255)                                                                                                                    gCCCGTAACAGCCAGAAAATATACATCCCGGAGAATTTTACCGTCTCGTGGGATTACTCGTTATGCAAACGCGCCATTGATGAAAACTGCTTTTTCAGCGATGAAGTCCCCGACCGTTGGGGTGACTGTATTGCGGCACg                  >  2:268913/1‑140 (MQ=255)                                                                                                                    gCCCGTAACAGCCAGAAAATATACATCCCGGAGAATTTTACCGTCTCGTGGGATTACTCGTTATGCAAACGCGCCATTGATGAAAACTGCTTTTTCAGCGATGAAGTCCCCGACCGTTGGGGTGACTGTATTGCGGCACg                  >  2:197642/1‑140 (MQ=255)                                                                                                                    gCCCGTAACAGCCAGAAAATATACATCCCGGAGAATTTTACCGTCTCGTGGGATTACTCGTTATGCAAACGCGCCATTGATGAAAACTGCTTTTTCAGCGATGAAGTCCCCGACCGTTGGGGTGACTGTATTGCGGCACg                  >  1:384336/1‑140 (MQ=255)                                                                                                                    gCCCGTAACAGCCAGAAAATATACATCCCGGAGAATTTTACCGTCTCGTGGGATTACTCGTTATGCAAACGCGCCATTGATGAAAACTGCTTTTTCAGCGATGAAGTCCCCGACCGTTGGGGTGACTGTATTGCGGCACg                  >  2:449055/1‑140 (MQ=255)                                                                                                                       cGTAACAGCCAGAAAATATACATCCCGGAGAATTTTACCGTCTCGTGGGATTACTCGTTATGCAAACGCGCCATTGATGAAAACTGCTTTTTCAGCGATGAAGTCCCCGACCGTTGGGGTGACTGTATTGCGGCACGCaa               >  2:732187/1‑140 (MQ=255)                                                                                                                        gTAACAGCCAGAAAATATACATCCCGGAGAATTTTACCGTCTCGTGGGATTACTCGTTATGCAAACGCGCCATTGATGAAAACTGCTTTTTCAGCGATGAAGTCCCCGACCGTTGGGGTGACTGTATTGCGGCACGCAAt              >  2:179113/1‑140 (MQ=255)                                                                                                                           aCAGCCAGAAAATATACATCCCGGAGAATTTTACCGTCTCGTGGGATTACTCGTTATGCAAACGCGCCATTGATGAAAACTGCTTTTTCAGCGATGAAGTCCCCGACCGTTGGGGTGACTGTATTGCGGCACGCAATCtt           >  1:450299/1‑140 (MQ=255)                                                                                                                               ccAGAAAATATACATCCCGGAGAATTTTACCGTCTCGTGGGATTACTCGTTATGCAAACGCGCCATTGATGAAAACTGCTTTTTCAGCGATGAAGTCCCCGACCGTTGGGGTGACTGTATTGCGGCACGCAATCTTGGCa       >  1:468014/1‑140 (MQ=255)                                                                                                                                  gAAAATATACATCCCGGAGAATTTTACCGTCTCGTGGGATTACTCGTTATGCAAACGCGCCATTGATGAAAACTGCTTTTTCAGCGATGAAGTCCCCGACCGTTGGGGTGACTGTATTGCGGCACGCAATCTTGGCATca    <  1:775416/140‑1 (MQ=255)                                                                                                                                    aaaTATACATCCCGGAGAATTTTACCGTCTCGTGGGATTACTCGTTATGCAAACGCGCCATTGATGAAAACTGCTTTTTCAGCGATGAAGTCCCCGACCGTTGGGGTGACTGTATTGCGGCACGCAATCTTGGCATcacc  <  1:732187/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                        GTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGTCTCGTGGGATTACTCGTTATGCAAACGCGCCATTGATGAAAACTGCTTTTTCAGCGATGAAGTCCCCGACCGTTGGGGTGACTGTATTGCGGCACGCAATCTTGGCATCACC  >  NC_000913/1877765‑1878035

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