Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	404,111	(G)5→6	coding (409/1116 nt)	yaiC →	diguanylate cyclase, cellulose regualtor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	404,106	1	.	G	100.0%	178.8 / NA	47	M135R (ATG→AGG)	yaiC	diguanylate cyclase, cellulose regualtor
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (29/18);  total (29/18)

CGAGCAGGGCCGTCGATCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGAT‑GGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCA  >  NC_000913/403965‑404249                                                                                                                                               |                                                                                                                                                cGAGCAGGGCCGTCGATCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGAt‑ggtgg                                                                                                                                            <  2:980436/147‑1 (MQ=255)           cGTCGATCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGGCAGGCGGCCCCCg                                                                                                                               <  2:833681/151‑1 (MQ=255)           cGTCGATCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGGCAGGCGGCCCCCg                                                                                                                               <  1:432453/151‑1 (MQ=255)               gATCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGGCAGGCGGCCCCCGTCTg                                                                                                                           >  1:842309/1‑151 (MQ=255)                  ccGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGGCAGGCGGCCCCCGTCTGttt                                                                                                                        >  2:207206/1‑151 (MQ=255)                    gCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGGCAGGCGGCCCCCGTCTGTTTGt                                                                                                                      >  2:791748/1‑151 (MQ=255)                    gCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGGCAGGCGGCCCCCGTCTGTTTGt                                                                                                                      <  1:746231/151‑1 (MQ=255)                        aGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGGCAGGCGGCCCCCGTCTGTTTGTcgc                                                                                                                   <  1:351182/150‑1 (MQ=255)                              gAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTg                                                                                                            <  1:990911/151‑1 (MQ=255)                                    tACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGtt                                                                                                       <  2:765350/150‑1 (MQ=255)                                      cAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGAt                                                                                                    <  2:325118/151‑1 (MQ=255)                                           tAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGAtggtgg                                                                                               >  2:31023/1‑151 (MQ=255)                                                 ccGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCt                                                                                         >  1:198492/1‑151 (MQ=255)                                                 ccGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCt                                                                                         >  2:928132/1‑151 (MQ=255)                                                 ccGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCt                                                                                         >  1:492297/1‑151 (MQ=255)                                                       cAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTg                                                                                   >  1:154249/1‑151 (MQ=255)                                                        aGTATCAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGAGATGTTGATGATTATGTGTCTGAATTTGATGGGGGGCAGGCGGCCCCCGTCTGTTTGTCGAGGGGCTGGTGTTGATGGTGGTTTCCTGCCtt                                                                                    >  2:218540/1‑149 (MQ=255)                                                             tAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCAccc                                                                             >  2:441873/1‑151 (MQ=255)                                                              aGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCt                                                                            >  1:939886/1‑151 (MQ=255)                                                              aGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCt                                                                            <  2:697453/151‑1 (MQ=255)                                                               gCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTc                                                                           >  1:244106/1‑151 (MQ=255)                                                                   gAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGc                                                                       >  1:705696/1‑151 (MQ=255)                                                                      tgtgGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTg                                                                     >  1:809896/1‑150 (MQ=255)                                                                         gggTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACgg                                                                 <  1:722079/151‑1 (MQ=255)                                                                         gggTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACgg                                                                 <  1:341765/151‑1 (MQ=255)                                                                          ggTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACggg                                                                <  1:349585/151‑1 (MQ=255)                                                                                cGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTAc                                                          >  2:236588/1‑151 (MQ=255)                                                                                cGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTAc                                                          >  1:285905/1‑151 (MQ=255)                                                                                  tAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCg                                                        >  2:208657/1‑151 (MQ=255)                                                                                      gggCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCgtgt                                                     <  2:943880/150‑1 (MQ=255)                                                                                      gggCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTAATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTc                                                    >  2:677211/1‑151 (MQ=255)                                                                                          gTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTc                                                >  1:462433/1‑151 (MQ=255)                                                                                          gTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTc                                                <  2:851792/151‑1 (MQ=255)                                                                                             aaCGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAAt                                             >  1:821553/1‑151 (MQ=255)                                                                                                gTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGt                                          >  1:666619/1‑151 (MQ=255)                                                                                                           cACCGCGATGTTGATGA‑TATGTGTCTGAATTTGATGGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGaa                              <  1:749836/151‑1 (MQ=255)                                                                                                             ccGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAAt                             >  2:810528/1‑151 (MQ=255)                                                                                                             ccGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAAt                             <  2:652915/151‑1 (MQ=255)                                                                                                                cgATGTTGATGATTATGTGTCTGAATTTGATGGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATACTGCGCCGCTGGGAtgg                           >  2:190136/1‑150 (MQ=255)                                                                                                                  aTGTTGATGATTATGTGTCTGAATTTGATGGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAAtggtgg                        >  1:650053/1‑151 (MQ=255)                                                                                                                    gTTGATGATTATGTGTCTGAATTTGATGGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGACTGGTGTTGCTGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAAtgg                           >  2:715433/1‑146 (MQ=255)                                                                                                                      tgatgaTTATGTGTCTGAATTTGATGGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGctct                    >  1:725170/1‑151 (MQ=255)                                                                                                                       gatgatTATGTGTCTGAATTTGATGGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGctct                    <  2:678709/150‑1 (MQ=255)                                                                                                                          gatTATGTGTCTGAATTTGATGGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCt                >  2:335850/1‑151 (MQ=255)                                                                                                                                   tCTGAATTTGATGGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGCGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCattat       >  2:730030/1‑151 (MQ=255)                                                                                                                                     tGAATTTGATGGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTg     >  1:952520/1‑151 (MQ=255)                                                                                                                                        aTTTGATGGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCa  <  2:925389/151‑1 (MQ=255)                                                                                                                                        aTTTGATGGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCa  <  2:133185/151‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                CGAGCAGGGCCGTCGATCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGAT‑GGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCA  >  NC_000913/403965‑404249

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