Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,871,162	A→G	T260A (ACA→GCA)	cdgI →	putative c‑di‑GMP binding protein CdgI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,871,162	0	A	G	100.0%	198.0 / NA	54	T260A (ACA→GCA)	cdgI	putative c‑di‑GMP binding protein CdgI
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (26/28);  total (26/28)

AGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCACACGATTAGCCAGCGAACTATGGCCGTTAATAGCATTATTATGTCTGGCATCGCTATGCTGCAACTTACTTCTACTGACTCTGGATGAGTATAATTATACTATCTGGTATATCAGTCGCGGGATTGAAGTTTCC  >  NC_000913/1871027‑1871296                                                                                                                                        |                                                                                                                                       aGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGcgca                                                                                                                                      >  1:1025205/1‑138 (MQ=255) aGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGcgca                                                                                                                                      <  2:1025205/138‑1 (MQ=255)    aGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCGCACGAtt                                                                                                                                 <  1:472872/140‑1 (MQ=255)      ttATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCGCACGATTAg                                                                                                                               >  2:360368/1‑140 (MQ=255)       tataACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCGCACGATTAGc                                                                                                                              >  2:791755/1‑140 (MQ=255)              cTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCGCACGATTAGCCAGCGaa                                                                                                                       >  1:90951/1‑140 (MQ=255)                           gTATGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCGCACGATTAGCCAGCGAACTATGGCCGTTaa                                                                                                          >  2:186195/1‑140 (MQ=255)                                tcgtcgATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCGCACGATTAGCCAGCGAACTATGGCCGTTAATAGCa                                                                                                     >  2:188891/1‑140 (MQ=255)                                 cgtcgATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCGCACGATTAGCCAGCGAACTATGGCCGTTAATAGCat                                                                                                    <  1:186195/140‑1 (MQ=255)                                  gtcgATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCGCACGATTAGCCAGCGAACTATGGCCGTTAATAGCatt                                                                                                   >  1:923789/1‑140 (MQ=255)                                   tcgATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCGCACGATTAGCCAGCGAACTa                                                                                                                    >  2:602253/1‑122 (MQ=255)                                   tcgATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCGCACGATTAGCCAGCGAACTa                                                                                                                    <  1:602253/122‑1 (MQ=255)                                    cgATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCGCACGATTAGCCAGCGAACTATGGCCGTTAATAGCattat                                                                                                 >  1:158593/1‑140 (MQ=255)                                    cgATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCGCACGATTAGCCAGCGAACTATGGCCGTTAATAGCa                                                                                                     >  1:77278/1‑136 (MQ=255)                                    cgATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCGCACGATTAGCCAGCGAACTATGGCCGTTAATAGCa                                                                                                     <  2:77278/136‑1 (MQ=255)                                      aTTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCGCACGATTAGCCAGCGAACTATGGCCGTTAATAg                                                                                                       <  1:940019/132‑1 (MQ=255)                                      aTTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCGCACGATTAGCCAGCGAACTATGGCCGTTAATAg                                                                                                       >  2:940019/1‑132 (MQ=255)                                          cTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCGCACGATTAGCCAGCGAACTATGGCCGTTAATAGCATTATTATGTc                                                                                           <  1:134075/140‑1 (MQ=255)                                             tCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCGCACGATTAGCCAGCGAACTATGGCCGTTAATAGCATTATTATGTCTgg                                                                                        >  1:390919/1‑140 (MQ=255)                                             tCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCGCACGATTAGCCAGCGAACTATGGCCGTTAATAGCATTATTATGTCTgg                                                                                        >  2:154701/1‑140 (MQ=255)                                             tCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCGCACGATTAGCCAGCGAACTATGGCCGTTAATAGCATTATTATGTCTgg                                                                                        >  1:950280/1‑140 (MQ=255)                                               cGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCGCACGATTAGCCAGCGAACTATGGCCGTTAATAGCATTATTATGTCTGGCa                                                                                      <  1:767760/140‑1 (MQ=255)                                                ggACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCGCACGATTAGCCAGCGAACTATGGCCGTTAATAGCATTATTATGTCTGGCAt                                                                                     <  1:360368/140‑1 (MQ=255)                                                  aCAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCGCACGATTAGCCAGCGAACTATGGCCGTTAATAGCATTATTATGTCTGGCATCg                                                                                   <  2:950280/140‑1 (MQ=255)                                                  aCAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCGCACGATTAGCCAGCGAACTATGGCCGTTAATAGCATTATTATGTCTGGCATCg                                                                                   <  1:188891/140‑1 (MQ=255)                                                     aCCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCGCACGATTAGCCAGCGAACTATCGCCGTTAATAGCATTATTATGTCTGGCATCGCTa                                                                                <  2:923789/140‑1 (MQ=255)                                                          aCTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCGCACGATTAGCCAGCGAACTATGGCCGTTAATAGCATTATTATGTCTGGCATCGCTAtgctg                                                                           <  2:166555/140‑1 (MQ=255)                                                          aCTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCGCACGATTAGCCAGCGAACTATGGCCGTTAATAGCATTATTATGTCTGGCATCGCTAtgctg                                                                           <  2:390919/140‑1 (MQ=255)                                                              cGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCGCACGATTAGCCAGCGAACTATGGCCGTTAATAGCATTATTATGTCTGGc                                                                                       <  1:503300/124‑1 (MQ=255)                                                              cGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCGCACGATTAGCCAGCGAACTATGGCCGTTAATAGCATTATTATGTCTGGc                                                                                       >  2:503300/1‑124 (MQ=255)                                                              cGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCGCACGATTAGCCAGCGAACTATGGCCGTTAATAGCATTATTATGTCTGGCATCGCT‑TGCTGCAact                                                                      >  2:1027094/1‑140 (MQ=255)                                                               gACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCGCACGATTAGCCAGCGAACTATGGCCGTTAATAGCATTATTATGTCTGGCATTGCTATGCTGCAact                                                                      <  2:90951/140‑1 (MQ=255)                                                                aCCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCGCACGATTAGCCAGCGAACTATGGCCGTTAATAGCATTATTATGTCTGGCATCGCTATGCTGCAactt                                                                     >  1:833987/1‑140 (MQ=255)                                                                 ccTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCGCACGATTAGCCAGCGAACTATGGCCGTTAATAGCATTATTATGTCTGGCATCGCTATGCTGCAactta                                                                    >  1:921040/1‑140 (MQ=255)                                                                        aaTCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCGCACGATTAGCCAGCGAACTATGGCCGTTAATAGCATTATTATGTCTGGCATCGCTATGCTGCAACTTACTTCTAc                                                             <  1:781343/140‑1 (MQ=255)                                                                          tCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCGCACGATTAGCCAGCGAACTATGGCCGTTAATAGCATTATTATGTCTGGCATCGCTATGCTGCAACTTACTTCTACt                                                            <  2:1024385/139‑1 (MQ=255)                                                                          tCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCGCACGATTAGCCAGCGAACTATGGCCGTTAATAGCATTATTATGTCTGGCATCGCTATGCTGCAACTTACTTCTACt                                                            >  1:1024385/1‑139 (MQ=255)                                                                                 taCAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCGCACGATTAGCCAGCGAACTATGGCCGTTAATAGCATTATTATGTCTGGCATCGCTATGCTGCAACTTACTTCTACTGACTCTgg                                                    >  1:363131/1‑140 (MQ=255)                                                                                  aCAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCGCACGATTAGCCAGCGAACTATGGCCGTTAATAGCATTATTATGTCTGGCATCGCT‑TGCTGCAACTTACTTCTACTGACTCTGGAt                                                  <  1:1027094/140‑1 (MQ=255)                                                                                   cAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCGCACGATTAGCCAGCGAACTATGGCCGTTAATAGCATTATTATGTCTGGCATCGCTATGCTGCAACTTACTTCTACTGACTCTGGAt                                                  <  1:221836/140‑1 (MQ=255)                                                                                      aaaTATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCGCACGATTAGCCAGCGAACTATGGCCGTTAATAGCATTATTATGTCTGGCATCGCTATGCTGCAACTTACTTCTACTGACTCTGGATGAg                                               <  2:282318/140‑1 (MQ=255)                                                                                        atatTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCGCACGATTAGCCAGCGAACTATGGCCGTTAATAGCATTATTATGTCTGGCATCGCTATGCTGCAACTTACTTCTACTGACTCTGGATGAg                                               <  2:782413/138‑1 (MQ=255)                                                                                        atatTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCGCACGATTAGCCAGCGAACTATGGCCGTTAATAGCATTATTATGTCTGGCATCGCTATGCTGCAACTTACTTCTACTGACTCTGGATGAg                                               >  1:782413/1‑138 (MQ=255)                                                                                        atatTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCGCACGATTAGCCAGCGAACTATGGCCGTTAATAGCATTATTATGTCTGGCATCGCTATGCTGCAACTTACTTCTACTGACTCTGGATGAGta                                             >  1:170569/1‑140 (MQ=255)                                                                                           ttGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCGCACGATTAGCCAGCGAACTATGGCCGTTAATAGCATTATTATGTCTGGCATCGCTATGCTGCAACTTACTTCTACTGACTCTGGATGAGTATaa                                          >  2:91779/1‑140 (MQ=255)                                                                                              gTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCGCACGATTAGCCAGCGAACTATGGCCGTTAATAGCATTATTATGTCTGGCATCGCTATGCTGCAACTTACTTCTACTGACTCTGGATGAGTATAATta                                       <  2:363131/140‑1 (MQ=255)                                                                                                ttGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCGCACGATTAGCCAGCGAACTATGGCCGTTAATAGCATTATTATGTCTGGCATCGCTATGCTGCAACTTACTTCTACTGACTCTGGATGAGTATAATtata                                     >  2:881586/1‑140 (MQ=255)                                                                                                  gTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCGCACGATTAGCCAGCGAACTATGGCCGTTAATAGCATTATTATGTCTGGCATCGCTATGCTGCAACTTACTTCTACTGACTCTGGATGAGTATAATTATACt                                   <  1:91779/140‑1 (MQ=255)                                                                                                  gTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCGCACGATTAGCCAGCGAACTATGGCCGTTAATAGCATTATTATGTCTGGCATCGCTATGCTGCAACTTACTTCTACTGACTCTGGATGAGTATAATTATACt                                   <  2:26190/140‑1 (MQ=255)                                                                                                   tCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCGCACGATTAGCCAGCGAACTATGGCCGTTAATAGCATTATTATGTCTGGCATCGCTATGCTGCAACTTACTTCTACTGACTCTGGATGAGTATAATTATACTa                                  <  2:833987/140‑1 (MQ=255)                                                                                                   tCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCGCACGATTAGCCAGCGAACTATGGCCGTTAATAGCATTATTATGTCTGGCATCGCTATGCTGCAACTTACTTCTACTGACTCTGGATGAGTATAATTATACTa                                  <  2:607723/140‑1 (MQ=255)                                                                                                   tCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCGCACGATTAGCCAGCGAACTATGGCCGTTAATAGCATTATTATGTCTGGCATCGCTATGCTGCAACTTACTTCTACTGACTCTGGATGAGTATAATTATACTa                                  <  2:921040/140‑1 (MQ=255)                                                                                                     tgtgGGCATTTTTACTGTTCTTTATTATCATGCGCGCACGATTAGCCAGCGAACTATGGCCGTTAATAGCATTATTATGTCTGGCATCGCTATGCTGCAACTTACTTCTACTGACTCTGGATGAGTATAATTATACTATc                                <  1:154701/140‑1 (MQ=255)                                                                                                                           tattatCATGCGCGCACGATTAGCCAGCGAACTATGGCCGTTAATAGCATTATTATGTCTGGCATCGCTATGCTGCAACTTACTTCTACTGACTCTGGATGAGTATAATTATACTATCTGGTATATCAGTCGCGGGATTg          >  2:83312/1‑140 (MQ=255)                                                                                                                                   tGCGCGCACGATTAGCCAGCGAACTATGGCCGTTAATAGCATTATTATGTCTGGCATCGCTATGCTGCAACTTACTTCTACTGACTCTGGATGAGTATAATTATACTATCTGGTATATCAGTCGCGGGATTGAAGTTTcc  >  2:397141/1‑140 (MQ=255)                                                                                                                                   tGCGCGCACGATTAGCCAGCGAACTATGGCCGTTAATAGCATTATTATGTCTGGCATCGCTATGCTGCAACTTACTTCTACTGACTCTGGATGAGTATAATTATACTATCTGGTATATCAGTCGCGGGATTGAAGTTTcc  >  1:335665/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                       AGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCACACGATTAGCCAGCGAACTATGGCCGTTAATAGCATTATTATGTCTGGCATCGCTATGCTGCAACTTACTTCTACTGACTCTGGATGAGTATAATTATACTATCTGGTATATCAGTCGCGGGATTGAAGTTTCC  >  NC_000913/1871027‑1871296

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