Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,143,684 (G)6→7 coding (419/3318 nt) yegE → putative diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,143,6781.G92.1% 206.4 / 13.0 64coding (413/3318 nt)yegEputative diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base . (3/2);  new base G (27/31);  total (31/33)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

CAGTGCTGCTACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTT‑GGGG‑GG‑TGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCA  >  NC_000913/2143579‑2143772
                                                                                                    |                                                                                                
cAGTGCTGCTACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGtt                                                                                                   >  1:1363844/1‑100 (MQ=255)
   tgctgcTACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTT‑ggg                                                                                               <  1:1081080/100‑1 (MQ=255)
     ctgctACGTAAGTTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTT‑gggg‑g                                                                                            <  1:425426/100‑1 (MQ=255)
       gctACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTT‑gggg‑ggg                                                                                          >  1:308406/1‑99 (MQ=255)
            gTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGG‑GG‑TGTTc                                                                                     <  1:732148/100‑1 (MQ=255)
             tAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGG‑GG‑TGTTCt                                                                                    <  1:894305/100‑1 (MQ=255)
             tAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTT‑gggg‑gggag                                                                                        >  2:1655812/1‑93 (MQ=255)
                 ttGCTGCCTTTGTATTATCTCTTTTAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGG‑GG‑TGTTCTgg                                                                                  <  2:1628626/98‑1 (MQ=255)
                 ttGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGG‑GG‑TGTTCTGGtt                                                                                >  2:881755/1‑100 (MQ=255)
                   gctgcCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGG‑GG‑TGTTCTGGTTGt                                                                              <  2:648171/100‑1 (MQ=255)
                    ctgcCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGG‑GG‑TGTTCTGGTTGTc                                                                             >  2:269802/1‑100 (MQ=255)
                     tgcCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGG‑GG‑TGTTCTGGTTGTcc                                                                            <  1:1353876/100‑1 (MQ=255)
                     tgcCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGG‑GG‑TGTTCTGGTTGTcc                                                                            <  2:1358693/100‑1 (MQ=255)
                      gcCGTGGTATAACACAATGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGG‑GG‑TGTTCTGGTTGTCct                                                                           <  1:1062095/100‑1 (MQ=255)
                        cGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGG‑GG‑TGTTCTGGTTGTCctgc                                                                         <  2:1363854/100‑1 (MQ=255)
                        cGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGG‑GG‑TGTTCTGGTTGTCctgc                                                                         >  1:1371699/1‑100 (MQ=255)
                        cGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTT‑GGGG‑GG‑TGTTATGGTTGTCctgct                                                                        >  2:1736433/1‑100 (MQ=255)
                         gTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGG‑GG‑TGTTCTGGTTGTCctgct                                                                        <  1:1611064/100‑1 (MQ=255)
                         gTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGG‑GG‑TGTTCTGGTTGTCctgct                                                                        >  1:132478/1‑100 (MQ=255)
                         gTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCCCCTCTGTTAGGGA‑GG‑TGTTCTGGTTGTCctg                                                                          >  2:757663/1‑98 (MQ=255)
                          tGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGG‑GG‑TGTTCTGGTTGTCctgctg                                                                       <  1:481387/100‑1 (MQ=255)
                          tGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGG‑GG‑TGTTCTGGTTGTCctgctg                                                                       >  1:1785668/1‑100 (MQ=255)
                             tataACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGG‑GG‑TGTTCTGGTTGTCCTGCTGACg                                                                    >  2:1153826/1‑100 (MQ=255)
                             tataACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGG‑GG‑TGTTCTGGTTGTCCTGCTGACg                                                                    >  1:1746640/1‑100 (MQ=255)
                                  ccccTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGG‑GG‑TGTTCTGGTTGTCCTGCTGACGCCCgg                                                               >  2:136739/1‑100 (MQ=255)
                                    ccTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGG‑GG‑TGTTCTGGTTGTCCTGCGGACGCCCGgag                                                             <  1:1396595/100‑1 (MQ=255)
                                     cTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGG‑GG‑TGTTCTGGTTGTCCTGCTGACGCCCGgaga                                                            <  2:1200672/100‑1 (MQ=255)
                                      ttGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGG‑GG‑TGTTCTGGTTGTCCTGCTGACGCCCGGAGAc                                                           <  2:385593/100‑1 (MQ=255)
                                       tGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGG‑GG‑TGTTCTGGTTGTCCTGCTGACGCCCGGAGAc                                                           <  1:1520740/99‑1 (MQ=255)
                                       tGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGG‑GG‑TGTTCTGGTTGTCCTGCTGACGCCCGGAGACg                                                          >  2:1516455/1‑100 (MQ=255)
                                       tGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGG‑GG‑TGTTCTGGTTGTCCTGCTGACGCCCGGAGACg                                                          <  1:1500824/100‑1 (MQ=255)
                                        gCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGG‑GG‑TGTTCTGGTTGTCCTGCTGACGCCCGGAGACGa                                                         >  2:444914/1‑100 (MQ=255)
                                         cAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGG‑GG‑TGTTCTGGTTGTCCTGCTGACGCCCGGAGACGAt                                                        >  1:1149329/1‑100 (MQ=255)
                                         cAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGG‑GG‑TGTTCTGGTTGTCCTGCTGACGCCCGGAGACGAt                                                        >  1:92048/1‑100 (MQ=255)
                                         cAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGG‑GG‑TGTTCTGGTTGTCCTGCTGACGCCCGGAGACGAt                                                        >  2:1112730/1‑100 (MQ=255)
                                           aaaTCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGG‑GG‑TGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATcc                                                      <  2:1170292/100‑1 (MQ=255)
                                                tGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGG‑GG‑TGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCa                                                 <  1:444912/100‑1 (MQ=255)
                                                tGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGG‑GG‑TGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCa                                                 <  2:1214074/100‑1 (MQ=255)
                                                   cTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGG‑GG‑TGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAggg                                              <  2:92048/100‑1 (MQ=255)
                                                    tGACTGGCTGCGTCTGGCACTCGGCAGCGTCATTGTTCCGCCTCTGTTGGGGG‑GG‑TGTTCTGGTTGTCCTGCTGACGCCCTGAGAGGATCGTCTCCggg                                              <  2:177003/99‑1 (MQ=255)
                                                     gACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGG‑GG‑TGTTCTGGTTGTCCTGCTGATGCCCTGAGAGGATCGTCTCCGGGcg                                            <  1:1563092/100‑2 (MQ=37)
                                                     gACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGG‑GG‑TGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCa                                            <  2:1559487/100‑1 (MQ=255)
                                                       cTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGG‑GG‑TGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCAtt                                          >  1:1334373/1‑100 (MQ=255)
                                                          gCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGG‑GG‑TGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCAttttt                                       >  2:1178409/1‑100 (MQ=255)
                                                          gCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGG‑GG‑TGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCAttttt                                       <  2:788484/100‑1 (MQ=255)
                                                            tGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGG‑GG‑TGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGa                                     >  1:1039034/1‑100 (MQ=255)
                                                             gCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGG‑GG‑TGATCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCCCAGGGCATTTTTGat                                    >  1:1173709/1‑100 (MQ=255)
                                                               gTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGG‑GG‑TGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGatat                                  >  1:6509/1‑100 (MQ=255)
                                                                tCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGG‑GG‑TGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGatat                                  <  1:1535326/99‑1 (MQ=255)
                                                                tCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGG‑GG‑TGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATg                                 >  1:407556/1‑100 (MQ=255)
                                                                  tGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGG‑GG‑TGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATggg                               <  1:837410/100‑1 (MQ=255)
                                                                   ggCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGG‑GG‑TGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGt                              >  2:1135974/1‑100 (MQ=255)
                                                                   ggCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGG‑GG‑TGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGt                              <  2:1614688/100‑1 (MQ=255)
                                                                     cACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGG‑GG‑TGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTAc                            >  2:1169508/1‑100 (MQ=255)
                                                                       cTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGG‑GG‑TGTTCTGGTTGTCCTGCTGACTCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTAc                            <  1:325337/98‑1 (MQ=255)
                                                                         cGGCAGCTCCATTGTTCCGCCTCTGTT‑GGGG‑GGTTGTTCTGGTTGTCCTGCTGACGCCAGGAGACGTTCCTCTCAGGGCATTTTTGATATGGGTACTGTc                        <  1:1147025/100‑1 (MQ=255)
                                                                         cGGCAGCGCCATTGTTCCGCCTCTGTTGGGGG‑GG‑TGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTc                        <  2:407558/100‑1 (MQ=255)
                                                                          ggCAGCGCCATTGTTCCGCCTCTGTT‑GGGG‑GG‑TGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAg                      >  1:657954/1‑100 (MQ=255)
                                                                               cgcCATTGTTCCGCCTCTGTTGGGGG‑GG‑TGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTCATATGGGTACTGTCAGAATc                  <  1:1112720/100‑1 (MQ=255)
                                                                                 ccATTGTTCCGCCTCTGTTGGGGG‑GG‑TGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCa                >  1:782460/1‑100 (MQ=255)
                                                                                     tGTTCCGCCTCTGTTGGGGG‑GG‑TGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCgg            <  2:824454/100‑1 (MQ=255)
                                                                                      gTTCCGCCTCTGTTGGGGG‑GG‑TGTTCTGGTTGTCCTGCTGACGCCCGGAGTCGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGgc           >  2:1278052/1‑100 (MQ=255)
                                                                                      gTTCCGCCTCTGTTGGGGG‑GG‑TGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGgc           >  2:317327/1‑100 (MQ=255)
                                                                                          cGCCTCTGTTGGGGG‑GG‑TGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGctc       >  2:15786/1‑100 (MQ=255)
                                                                                            ccTCTGTT‑GGGG‑GG‑TGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTgg    >  1:23919/1‑100 (MQ=255)
                                                                                             ctctGTTGGGGG‑GG‑TGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTgg    <  2:136625/100‑1 (MQ=255)
                                                                                             ctctGTT‑GGGGTGG‑TGTTCTGGTTGTCCTGCTGACGCTCGCAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCGCTgg    <  1:277532/100‑1 (MQ=255)
                                                                                               ctGTTGGGGG‑GG‑TGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCa  >  1:1784477/1‑100 (MQ=255)
                                                                                               ctGTTGGGGG‑GG‑TGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCa  >  1:1455072/1‑100 (MQ=255)
                                                                                                    |                                                                                                
CAGTGCTGCTACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTT‑GGGG‑GG‑TGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCA  >  NC_000913/2143579‑2143772

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