Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,146,293	C→T	R1010C (CGC→TGC)	yegE →	putative diguanylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,146,293	0	C	T	100.0%	204.1 / NA	54	R1010C (CGC→TGC)	yegE	putative diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (27/27);  total (27/27)

AATGATCTGCTTGAGCAGCTGGAAAATAGCCCTCTACCACCACGGTTATTACATCTGATTATTCCGGCTGAAGCGATTTTAGATCACGCAGAAAGCGTGCAAAAACTGCGGCTGGCGGGATGTCGGATAGTGCTCAGCCAGGTGGGCCGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCGACTCG  >  NC_000913/2146146‑2146434                                                                                                                                                    |                                                                                                                                              aaTGATCTGCTTGAGCAGCTGGAAAATAGCCCTCTACCACCACGGTTATTACATCTGATTATTCCGGCTGAAGCGATTTTAGATCACGCAGAAAGCGTGCAAAAACTGCGGCTGGCGGGATGTCGGATAGTGCTCAGCCAGGTGGGCtgcg                                                                                                                                            <  2:109669/151‑1 (MQ=255)  aTGATCTGCTTGAGCAGCTGGAAAATAGCCCTCTACCACCACGGTTATTACATCTGATTATTCCGGCTGAAGCGATTTTAGATCACGCAGAAAGCGTGCAAAAACTGCGGCTGGCGGGATGTCGGATAGTGCTCAGCCAGGTGGGCTGCGa                                                                                                                                           <  1:790042/151‑1 (MQ=255)        tGCTTGAGCAGCTGGAAAATAGCCCTCTACCACCACGGTTATTACATCTGATTATTCCGGCTGAAGCGATTTTAGATCACGCAGAAAGCGTGCAAAAACTGCGGCTGGCGGGATGTCGGATAGTGCTCAGCCAGGTGGGCTGCGATCTGc                                                                                                                                      >  2:902525/1‑150 (MQ=255)         gCTTGAGCAGCTGGAAAATAGCCCTCTACCACCACGGTTATTACATCTGATTATTCCGGCTGAAGCGATTTTAGATCACGCAGAAAGCGTGCAAAAACTGCGGCTGGCGGGATGTCGGATAGTGCTCAGCCAGGTGGGCTGCGATCt                                                                                                                                        <  1:831234/147‑1 (MQ=255)         gCTTGAGCAGCTGGAAAATAGCCCTCTACCACCACGGTTATTACATCTGATTATTCCGGCTGAAGCGATTTTAGATCACGCAGAAAGCGTGCAAAAACTGCGGCTGGCGGGATGTCGGATAGTGCTCAGCCAGGTGGGCTGCGATCTGCaa                                                                                                                                    <  2:598120/151‑1 (MQ=255)             gAGCAGCTGGAAAATAGCCCTCTACCACCACGGTTATTACATCTGATTATTCCGGCTGAAGCGATTTTAGATCACGCAGAAAGCGTGCAAAAACTGCGGCTGGCGGGATGTCGGATAGTGCTCAGCCAGGTGGGCTGCGATCTGCAAATCt                                                                                                                                >  2:578469/1‑151 (MQ=255)                cagcTGGAAAATAGCCCTCTACCACCACGGTTATTACATCTGATTATTCCGGCTGAAGCGATTTTAGATCACGCAGAAAGCGTGCAAAAACTGCGGCTGGCGGGATGTTGGATAGTGCTCAGCCAGGTGGGCTGCGATCTGCAAATCtt                                                                                                                               <  2:666299/149‑1 (MQ=255)                   cTGGAAAATAGCCCTCTACCACCACGGTTATTACATCTGATTATTCCGGCTGAAGCGATTTTAGATCACGCAGAAAGCGTGCAAAAACTGCGGCTGGCGGGATGTCGGATAGTGCTCAGCCAGGTGGGCTGCGATCTGCAAATCTTCAACt                                                                                                                          >  1:13450/1‑151 (MQ=255)                       aaaaTAGCCCTCTACCACCACGGTTATTACATCTGATTATTCCGGCTGAAGCGATTTTAGATCACGCAGAAAGCGTGCAAAAACTGCGGCTGGCGGGATGTCGGATAGTGCTCAGCCAGGTGGGCTGCGATCTGCAAATCTTCAACTCGCt                                                                                                                      <  1:729065/151‑1 (MQ=255)                       aaaaTAGCCCTCTACCACCACGGTTATTACATCTGATTATTCCGGCTGAAGCGATTTTAGATCACGCAGAAAGCGTGCAAAAACTGCGGCTGGCGGGATGTCGGATAGTGCTCAGCCAGGTGGGCTGCGATCTGCAAATCTTCAACTCGCt                                                                                                                      <  1:353564/151‑1 (MQ=255)                        aaaTAGCCCTCTACCACCACGGTTATTACATCTGATTATTCCGGCTGAAGCGATTTTAGATCACGCAGAAAGCGTGCAAAAACTGCGGCTGGCGGGATGTCGGATAGTGCTCAGCCACGTGGGCTGCGATCTGCAAATCTTCAACTCGCTg                                                                                                                     >  1:875644/1‑151 (MQ=255)                               ccTCTACCACCACGGTTATTACATCTGATTATTCCGGCTGAAGCGATTTTAGATCACGCAGAAAGCGTGCAAAAACTGCGGCTGGCGGGATGTCGGATAGTGCTCAGCCAGGTGGGCTGCGATCTGCAAATCTTCAACTCGCTGAAAGCg                                                                                                               >  1:120844/1‑150 (MQ=255)                                ctctACCACCACGGTTATTACATCTGATTATTCCGGCTGAAGCGATTTTAGATCACGCAGAAAGCGTGCAAAAACTGCGGCTGGCGGGATGTCGGATAGTGCTCAGCCAGGTGGGCTGCGATCTGCAAATCTTCAACTCGCTGAAAGCGa                                                                                                              >  2:487434/1‑150 (MQ=255)                                  ctACCACCACGGTTATTACATCTGATTATTCCGGCTGAAGCGATTTTAGATCACGCAGAAAGCGTGCAAAAACTGCGGCTGGCGGGATGTCGGATAGTGCTCAGCCAGGTGGGCTGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAat                                                                                                            >  2:786131/1‑150 (MQ=255)                                    accaccACGGTTATTACATCTGATTATTCCGGCTGAAGCGATTTTAGATCACGCAGAAAGCGTGCAAAAACTGCGGCTGGCGGGATGTCGGATAGTGCTCAGCCAGGTGGGCTGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATg                                                                                                         >  2:167621/1‑151 (MQ=255)                                         cacGGTTATTACATCTGATTATTCCGGCTGAAGCGATTTTAGATCACGCAGAAAGCGTGCAAAAACTGCGGCTGGCGGGATGTCGGATAGTGCTCAGCCAGGTGGGCTGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGc                                                                                                       >  1:55457/1‑148 (MQ=255)                                         cacGGTTATTACATCTGATTATTCCGGCTGAAGCGATTTTAGATCACGCAGAAAGCGTGCAAAAACTGCGGCTGGCGGGATGTCGGATAGTGCTCAGCCAGGTGGGCTGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGa                                                                                                    <  1:349695/151‑1 (MQ=255)                                              ttattaCATCTGATTATTCCGGCTGAAGCGATTTTAGATCACGCAGAAAGCGTGCAAAAACTGCGGCTGGCGGGATGTCGGATAGTGCTCAGCCAGGTGGGCTGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTAcc                                                                                               >  2:433624/1‑151 (MQ=255)                                                attaCATCTGATTATTCCGGCTGAAGCGATTTTAGATCACGCAGAAAGCGTGCAAAAACTGCGGCTGGCGGGATGTCGGATAGTGCTCAGCCAGGTGGGCTGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTg                                                                                             >  2:708441/1‑151 (MQ=255)                                                 ttaCATCTGATTATTCCGGCTGAAGCGATTTTAGATCACGCAGAAAGCGTGCAAAAACTGCGGCTGGCGGGATGTCGGATAGTGCTCAGCCAGGTGGGCTGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGc                                                                                            >  1:215427/1‑151 (MQ=255)                                                   aCATCTGATTATTCCGGCTGAAGCGATTTTAGATCACGCAGAAAGCGTGCAAAAACTGCGGCTGGCGGGATGTCGGATAGTGCTCAGCCAGGTGGGCTGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTa                                                                                          <  2:334570/151‑1 (MQ=255)                                                             attCCGGCTGAAGCGATTTTAGATCACGCAGAAAGCGTGCAAAAACTGCGGCTGGCGGGATGTCGGATAGTGCTCAGCCAGGTGGGCTGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTg                                                                                >  1:808383/1‑151 (MQ=255)                                                               tCCGGCTGAAGCGATTTTAGATCACGCAGAAAGCGTGCAAAAACTGCGGCTGGCGGGATGTCGGATAGTGCTCAGCCAGGTGGGCTGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGa                                                                               <  1:282248/150‑1 (MQ=255)                                                                      gAAGCGATTTTAGATCACGCAGAAAGCGTGCAAAAACTGCGGCTGGCGGGATGTCGGATAGTGCTCAGCCAGGTGGGCTGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATgc                                                                        <  1:641243/150‑1 (MQ=255)                                                                       aaGCGATTTTAGATCACGCAGAAAGCGTGCAAAAACTGCGGCTGGCGGGATGTCGGATAGTGCTCAGCCAGGTGGGCTGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATgcgc                                                                      >  1:685434/1‑151 (MQ=255)                                                                           gATTTTAGATCACGCAGAAAGCGTGCAAAAACTGCGGCTGGCGGGATGTCGGATAGTGCTCAGCCAGGTGGGCTGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATg                                                                         <  1:520471/144‑1 (MQ=255)                                                                           gATTTTAGATCACGCAGAAAGCGTGCAAAAACTGCGGCTGGCGGGATGTCGGATAGTGCTCAGCCAGGTGGGCTGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAAc                                                                  <  1:320345/151‑1 (MQ=255)                                                                            aTTTTAGATCACGCAGAAAGCGTGCAAAAACTGCGGCTGGCGGGATGTCGGATAGTGCTCAGCCAGGTGGGCTGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACg                                                                 >  2:883315/1‑151 (MQ=255)                                                                                   aTCACGCAGAAAGCGTGCAAAAACTGCGGCTGGCGGGATGTCGGATAGTGCTCAGCCAGGTGGGCTGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAggg                                                          <  2:754725/151‑1 (MQ=255)                                                                                     cACGCAGAAAGCGTGCAAAAACTGCGGCTGGCGGGATGTCGGATAGTGCTCAGCCAGGTGGGCTGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTa                                                        <  1:182113/151‑1 (MQ=255)                                                                                         cAGAAAGCGTGCAAAAACTGCGGCTGGCGGGATGTCGGATAGTGCTCAGCCAGGTGGGCTGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAAttt                                                    <  2:546798/151‑1 (MQ=255)                                                                                         cAGAAAGCGTGCAAAAACTGCGGCTGGCGGGATGTCGGATAGTGCTCAGCCAGGTGGGCTGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAAttt                                                    <  2:8549/151‑1 (MQ=255)                                                                                             aaGCGTGCAAAAACTGCGGCTGGCGGGATGTCGGATAGTGCTCAGCCAGGTGGGCTGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTgatg                                                >  1:755431/1‑151 (MQ=255)                                                                                             aaGCGTGCAAAAACTGCGGCTGGCGGGATGTCGGATAGTGCTCAGCCAGGTGGGCTGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTgatg                                                >  2:726829/1‑151 (MQ=255)                                                                                             aaGCGTGCAAAAACTGCGGCTGGCGGGATGTCGGATAGTGCTCAGCCAGGTGGGCTGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTgatg                                                <  1:305354/151‑1 (MQ=255)                                                                                                cGTGCAAAAACTGCGGCTGGCGGGATGTCGGATAGTGCTCAGCCAGGTGGGCTGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTgatgga                                              >  2:10979/1‑150 (MQ=255)                                                                                                   gCAAAAACTGCGGCTGGCGGGATGTCGGATAGTGCTCAGCCAGGTGGGCTGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTgatggatg                                            >  2:706660/1‑149 (MQ=255)                                                                                                      aaaaCTGCGGCTGGCGGGATGTCGGATAGTGCTCAGCCAGGTGGGCTGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATg                                       >  1:600096/1‑151 (MQ=255)                                                                                                        aaCTGCGGCTGGCGGGATGTCGGATAGTGCTCAGCCAGGTGGGCTGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCt                                     <  1:871568/151‑1 (MQ=255)                                                                                                         aCTGCGGCTGGCGGGATGTCGGATAGTGCTCAGCCAGGTGGGCTGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTg                                    >  1:752087/1‑151 (MQ=255)                                                                                                            gCGGCTGGCGGGATGTCGGATAGTGCTCAGCCAGGTGGGCTGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGa                                   <  2:808318/149‑1 (MQ=255)                                                                                                             cGGCTGGCGGGATGTCGGATAGTGCTCAGCCAGGTGGGCTGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTa                                <  1:88110/151‑1 (MQ=255)                                                                                                             cGGCTGGCGGGATGTCGGATAGTGCTCAGCCAGGTGGGCTGCGATCTGCAAATCTTCAACTCGCTGAAAGCAAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTa                                <  1:685662/151‑1 (MQ=255)                                                                                                              ggCTGGCGGGATGTCGGATAGTGCTCAGCCAGGTGGGCTGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTAc                               <  1:352788/151‑1 (MQ=255)                                                                                                              ggCTGGCGGGATGTCGGATAGTGCTCAGCCAGGTGGGCTGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTAc                               >  1:349339/1‑151 (MQ=255)                                                                                                               gCTGGCGGGATGTCGGATAGTGCTCAGCCAGGTGGGCTGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACg                              <  2:178757/151‑1 (MQ=255)                                                                                                                cTGGCGGGATGTCGGATAGTGCTCAGCCAGGTGGGCTGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACg                              >  1:133654/1‑150 (MQ=255)                                                                                                                 tGGCGGGATGTCGGATAGTGCTCAGCCAGGTGGGCTGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGat                            >  1:792799/1‑151 (MQ=255)                                                                                                                  ggCGGGATGTCGGATAGTGCTCAGCCAGGTGGGCTGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGatt                           <  2:681921/151‑1 (MQ=255)                                                                                                                                 aGTGCTCAGCCAGGTGGGCTGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCt            <  2:349313/151‑1 (MQ=255)                                                                                                                                   tGCTCAGCCAGGTGGGCTGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTc           >  1:125279/1‑150 (MQ=255)                                                                                                                                       cAGCCAGGTGGGCTGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCg       <  2:774101/150‑1 (MQ=255)                                                                                                                                        aGCCAGGTGGGCTGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCGAc     <  2:738340/151‑1 (MQ=255)                                                                                                                                          ccAGGTGGGCTGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCGACTc   >  2:581538/1‑151 (MQ=255)                                                                                                                                           cAGGTGGGCTGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCGACTCg  >  1:422823/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                              AATGATCTGCTTGAGCAGCTGGAAAATAGCCCTCTACCACCACGGTTATTACATCTGATTATTCCGGCTGAAGCGATTTTAGATCACGCAGAAAGCGTGCAAAAACTGCGGCTGGCGGGATGTCGGATAGTGCTCAGCCAGGTGGGCCGCGATCTGCAAATCTTCAACTCGCTGAAAGCGAATATGGCAGATTACCTGCTACTTGATGGTGAGTTATGCGCCAACGTGCAGGGTAATTTGATGGATGAGATGCTGATTACGATTATTCAGGGGCACGCTCAGCGACTCG  >  NC_000913/2146146‑2146434

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