Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,685,370	G→A	S83S (AGC→AGT)	pdeK ←	putative c‑di‑GMP phosphodiesterase PdeK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,685,370	0	G	A	100.0%	191.1 / NA	54	S83S (AGC→AGT)	pdeK	putative c‑di‑GMP phosphodiesterase PdeK
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (26/28);  total (26/28)

CGTTCGAGCGAGTAAACGCCCAGCGAGATTTGAACCGGTAGCTCAAACAGGCGAGTAACCATTACCGGCACCGGGCGCTCTGGAATAAAACTTTTACGCAGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGGCTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGACAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGAT  >  NC_000913/3685233‑3685506                                                                                                                                          |                                                                                                                                         cGTTCGAGCGAGTAAACGCCCAGCGAGATTTGAACCGGTAGCTCAAACAGGCGAGTAACCATTACCGGCACCGGGCGCTCTGGAATAAAACTTTTACGCAGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGACt                                                                                                                                        >  1:81366/1‑140 (MQ=255)             tAAACGCCCAGCGAGATTTGAACCGGTAGCTCAAACAGGCGAGTAACCATTACCGGCACCGGGCGCTCTGGAATAAAACTTTTACGCAGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGACTGACCACGCCTGc                                                                                                                            >  2:1219692/1‑140 (MQ=255)             tAAACGCCCAGCGAGATTTGAACCGGTAGCTCAAACAGGCGAGTAACCATTACCGGCACCGGGCGCTCTGGAATAAAACTTTTACGCAGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGACTGACCACGCCTGc                                                                                                                            >  1:849984/1‑140 (MQ=255)             tAAACGCCCAGCGAGATTTGAACCGGTAGCTCAAACAGGCGAGTAACCATTACCGGCACCGGGCGCTCTGGAATAAAACTTTTACGCAGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGACTGACCACGCCTGc                                                                                                                            >  1:679774/1‑140 (MQ=255)                  gCCCAGCGAGATTTGAACCGGTAGCTCAAACAGGCGAGTAACCATTACCGGCACCGGGCGCTCTGGAATAAAACTTTTACGCAGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGACTGACCACGCCTGCCGGtt                                                                                                                       <  1:74400/140‑1 (MQ=255)                  gCCCAGCGAGATTTGAACCGGTAGCTCAAACAGGCGAGTAACCATTACCGACACCGGGCGCTCTGGAATAAAACTTTTACGCAGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGACTGACCACGCCTGCCGGtt                                                                                                                       <  1:1171180/140‑1 (MQ=255)                     cAGCGAGATTTGAACCGGTAGCTCAAACAGGCGAGTAACCATTACCGGCACCGGGCGCTCTGGAATAAAACTTTTACGCAGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGACTGACCACGCCTGCCGGTTTaa                                                                                                                    >  1:1086164/1‑140 (MQ=255)                      aGCGAGATTTGAACCGGTAGCTCAAACAGGCGAGTAACCATTACCGGCACCGGGCGCTCTGGAATAAAACTTTTACGCAGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGACTGACCACGCCTGCCGGTTTAAt                                                                                                                   >  2:551940/1‑140 (MQ=255)                        cGAGATTTGAACCGGTAGCTC‑AACAGGCGAGTAACCATTACCGGCACCGGGCGCTCTGGAATAAAACTTTTACGCAGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGACTGACCACGCCTGCCGGTTTAATGCt                                                                                                                <  1:845531/140‑1 (MQ=255)                         gagaTTTGAACCGGTAGCTCAAACAGGCGAGTAACCATTACCGGCACCGGGCGCTCTGGAATAAAACTTTTACGCAGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGACTGACCACGCCTGCCGGTTTAATGCt                                                                                                                <  1:328493/140‑1 (MQ=255)                         gagaTTTGAACCGGTAGCTCAAACAGGCGAGTAACCATTACCGGCACCGGGCGCTCTGGAATAAAACTTTTACGCAGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGACTGACCACGCCTGCCGGTTTAATGCt                                                                                                                <  2:612149/140‑1 (MQ=255)                            aTTTGAACCGGTAGCTCAAACAGGCGAGTAACCATTACCGGCACCGGGCGCTCTGGAATAAAACTTTTACGCAGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGACTGACCACGCCTGCCGGTTTAATGCTGGc                                                                                                             <  1:743958/140‑1 (MQ=255)                                 aaCCGGTAGCTCAAACAGGCGAGTAACCATTACCGGCACCGGGCGCTCTGGAATAAAACTTTTACGCAGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGACTGACCACGCCTGCCGGTTTAATGCTGGCAAGaa                                                                                                        >  1:253615/1‑140 (MQ=255)                                      gTAGCTCAAACAGGCGAGTAACCATTACCGGCACCGGGCGCTCTGGAATAAAACTTTTACGCAGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGACTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCt                                                                                                   >  1:944345/1‑140 (MQ=255)                                      gTAGCTCAAACAGGCGAGTAACCATTACCGGCACCGGGCGCTCTGGAATAAAACTTTTACGCAGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGACTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCt                                                                                                   <  1:660792/140‑1 (MQ=255)                                      gTAGCTCAAACAGGCGAGTAACCATTACCGGCACCGGGCGCTCTGGAATAAAACTTTTACGCAGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGACTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCt                                                                                                   <  1:618747/140‑1 (MQ=255)                                      gTAGCTCAAACAGGCGAGTAACCATTACCGGCACCGGGCGCTCTGGAATAAAACTTTTACGCAGCGCCCGGAACTGGTTAGGCAGCACTACATCGGCACGACTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCt                                                                                                   >  1:267995/1‑140 (MQ=255)                                        aGCTCAAACAGGCGAGTAACCATTACCGGCACCGGGCGCTCTGGAATAAAACTTTTACGCAGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGACTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTc                                                                                                 <  1:1219692/140‑1 (MQ=255)                                           tCAAACAGGCGAGTAACCATTACCGGCACCGGGCGCTCTGGAATAAAACTTTTACGCAGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGACTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGc                                                                                              >  2:831840/1‑140 (MQ=255)                                               aCAGGCGAGTAACCATTACCGGCACCGGGCGCTCTGGAATAAAACTTTTACGCAGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGACTGACCACGCCTGCCGGTTTAATGC‑GGCAAGAATAGCTTCCGCTTCgg                                                                                         <  1:592236/140‑1 (MQ=255)                                                cAGGCGAGTAACCATTACCGGCACCGGGCGCTCTGGAATAAAACTTTTACGCAGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGACTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCgg                                                                                         <  1:763461/140‑1 (MQ=255)                                                        tAACCATTACCGGCACCGGGCGCTCTGGAATAAAACTTTTACGCAGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGACTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCg                                                                                 <  2:679774/140‑1 (MQ=255)                                                        tAACCATTACCGGCACCGGGCGCTCTGGAATAAAACTTTTACGCAGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGACTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCg                                                                                 <  1:1010308/140‑1 (MQ=255)                                                        tAACCATTACCGGCACCGGGCGCTCTGGAATAAAACTTTTACGCAGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGACTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCg                                                                                 <  2:81366/140‑1 (MQ=255)                                                        tAACCATTACCGGCACCGGGCGCTCTGGAATAAAACTTTTACGCAGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGACTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCg                                                                                 <  2:849984/140‑1 (MQ=255)                                                                    gCACCGGGCGCTCTGGAATAAAACTTTTACGCAGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGACTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATa                                                                     >  1:300296/1‑140 (MQ=255)                                                                     cACCGGGCGCTCTGGAATAAAACTTTTACGCAGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGACTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAg                                                                    <  1:728284/140‑1 (MQ=255)                                                                     cACCGGGCGCTCTGGAATAAAACTTTTACGCAGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGACTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAg                                                                    >  2:462388/1‑140 (MQ=255)                                                                      aCCGGGCGCTCTGGAATAAAACTTTTACGCAGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGACTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGc                                                                   >  2:631475/1‑140 (MQ=255)                                                                         gggCGCTCTGGAATAAAACTTTTACGCAGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGACTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGa                                                                >  1:1209447/1‑140 (MQ=255)                                                                          ggCGCTCTGGAATAAAACTTTTACGCAGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGACTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCt                                                                  >  2:945284/1‑137 (MQ=255)                                                                          ggCGCTCTGGAATAAAACTTTTACGCAGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGACTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCt                                                                  <  1:945284/137‑1 (MQ=255)                                                                                  ggAATAAAACTTTTACGCAGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGACTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAggg                                                       <  1:1101308/140‑1 (MQ=255)                                                                                             tttACGCAGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGACTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGa                                               <  1:1204824/137‑1 (MQ=255)                                                                                             tttACGCAGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGACTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGa                                               >  2:1204824/1‑137 (MQ=255)                                                                                             tttACGCAGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGACTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGACag                                            >  1:906939/1‑140 (MQ=255)                                                                                              ttACGCAGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGACTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGACaga                                           >  1:629937/1‑140 (MQ=255)                                                                                                      cgcCTGGAACTGGTTAGGCAGCACTACATCGGCACGACTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGACAGAGCGAGCaa                                   >  2:288390/1‑140 (MQ=255)                                                                                                      cgcCTGGAACTGGTTAGGCAGCACTACATCGGCACGACTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGACAGAGCGAGCaa                                   <  1:462388/140‑1 (MQ=255)                                                                                                              aCTGGTTAGGCAGCACTACATCGGCACGACTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGACAGAGCGAGCAATGCTTTcc                           <  2:944345/140‑1 (MQ=255)                                                                                                                tGGTTAGGCAGCACTACATCGGCACGACTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGACAGAGCGAGCAATGCTTTCCAg                         >  1:936289/1‑140 (MQ=255)                                                                                                                tGGTTAGGCAGCACTACATCGGCACGACTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGACAGAGCGAGCAATGCTTTCCAg                         <  1:831840/140‑1 (MQ=255)                                                                                                                tGGTTAGGCAGCACTACATCGGCACGACTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGACAGAGCGAGCAATGCTTTCCAg                         <  1:288390/140‑1 (MQ=255)                                                                                                                  gTTAGGCAGCACTACATCGGCACGACTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGACAGAGCGAGCAATGCTTTCCAGtt                       >  2:1075934/1‑140 (MQ=255)                                                                                                                     aGGCAGCACTACATCGGCACGACTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGACAGAGCGAGCAATGCTTTCCAGTTGCg                    <  2:603062/140‑1 (MQ=255)                                                                                                                       gcagcaCTACATCGGCACGACTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGACAGAGCGAGCAATGCTTTCCAGTTGCGTa                  >  1:1115222/1‑140 (MQ=255)                                                                                                                        cagcaCTACATCGGCACGACTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGACAGAGCGAGCAATGCTTTCCAGTTGCGTAg                 <  2:86213/140‑1 (MQ=255)                                                                                                                          gcaCTACATCGGCACGACTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGACAGAGCGAGCAATGCTTTCCAGTTGCGTAGcc               >  2:244378/1‑140 (MQ=255)                                                                                                                            aCTACATCGGCACGACTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGACAGAGCAAGCAATGCTTTCCAGTTGCGTAGCCgt             >  1:613851/1‑140 (MQ=255)                                                                                                                              tACATCGGCACGTCTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGACAGGGCGAGCAATGCTTTCCAGTTGCGTAGCCgtgt           <  2:267995/140‑1 (MQ=255)                                                                                                                                   cGGCACGACTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGACAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAg      >  1:645847/1‑140 (MQ=255)                                                                                                                                     gCACGACTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGACAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCg    >  2:967396/1‑140 (MQ=255)                                                                                                                                      cACGACTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGACAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGa   <  2:253615/140‑1 (MQ=255)                                                                                                                                       aCGACTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGACAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGAt  <  1:631475/140‑1 (MQ=255)                                                                                                                                       aCGACTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGACAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGAt  <  1:1122738/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                         CGTTCGAGCGAGTAAACGCCCAGCGAGATTTGAACCGGTAGCTCAAACAGGCGAGTAACCATTACCGGCACCGGGCGCTCTGGAATAAAACTTTTACGCAGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGGCTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGACAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGAT  >  NC_000913/3685233‑3685506

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