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

CAGTGCTGCTACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCA  >  NC_000913/2143579‑2143772
                                                                                                    |                                                                                                
cAGTGCTGCTACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGtt                                                                                                   >  1:1363844/1‑100 (MQ=255)
   tgctgcTACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTggg                                                                                               <  1:1081080/100‑1 (MQ=255)
     ctgctACGTAAGTTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTggggg                                                                                            <  1:425426/100‑1 (MQ=255)
       gctACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTggggggg                                                                                          >  1:308406/1‑99 (MQ=255)
            gTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTc                                                                                     <  1:732148/100‑1 (MQ=255)
             tAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCt                                                                                    <  1:894305/100‑1 (MQ=255)
             tAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTgggggggag                                                                                        >  2:1655812/1‑93 (MQ=255)
                 ttGCTGCCTTTGTATTATCTCTTTTAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTgg                                                                                  <  2:1628626/98‑1 (MQ=255)
                 ttGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGtt                                                                                >  2:881755/1‑100 (MQ=255)
                   gctgcCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGt                                                                              <  2:648171/100‑1 (MQ=255)
                    ctgcCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTc                                                                             >  2:269802/1‑100 (MQ=255)
                     tgcCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTcc                                                                            <  1:1353876/100‑1 (MQ=255)
                     tgcCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTcc                                                                            <  2:1358693/100‑1 (MQ=255)
                      gcCGTGGTATAACACAATGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCct                                                                           <  1:1062095/100‑1 (MQ=255)
                        cGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCctgc                                                                         <  2:1363854/100‑1 (MQ=255)
                        cGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCctgc                                                                         >  1:1371699/1‑100 (MQ=255)
                        cGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGTGTTATGGTTGTCctgct                                                                        >  2:1736433/1‑100 (MQ=255)
                         gTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCctgct                                                                        <  1:1611064/100‑1 (MQ=255)
                         gTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCctgct                                                                        >  1:132478/1‑100 (MQ=255)
                         gTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCCCCTCTGTTAGGGAGGTGTTCTGGTTGTCctg                                                                          >  2:757663/1‑98 (MQ=255)
                          tGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCctgctg                                                                       <  1:481387/100‑1 (MQ=255)
                          tGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCctgctg                                                                       >  1:1785668/1‑100 (MQ=255)
                             tataACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACg                                                                    >  2:1153826/1‑100 (MQ=255)
                             tataACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACg                                                                    >  1:1746640/1‑100 (MQ=255)
                                  ccccTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCgg                                                               >  2:136739/1‑100 (MQ=255)
                                    ccTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCGGACGCCCGgag                                                             <  1:1396595/100‑1 (MQ=255)
                                     cTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGgaga                                                            <  2:1200672/100‑1 (MQ=255)
                                      ttGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGAc                                                           <  2:385593/100‑1 (MQ=255)
                                       tGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGAc                                                           <  1:1520740/99‑1 (MQ=255)
                                       tGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACg                                                          >  2:1516455/1‑100 (MQ=255)
                                       tGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACg                                                          <  1:1500824/100‑1 (MQ=255)
                                        gCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGa                                                         >  2:444914/1‑100 (MQ=255)
                                         cAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGAt                                                        >  1:1149329/1‑100 (MQ=255)
                                         cAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGAt                                                        >  1:92048/1‑100 (MQ=255)
                                         cAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGAt                                                        >  2:1112730/1‑100 (MQ=255)
                                           aaaTCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATcc                                                      <  2:1170292/100‑1 (MQ=255)
                                                tGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCa                                                 <  1:444912/100‑1 (MQ=255)
                                                tGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCa                                                 <  2:1214074/100‑1 (MQ=255)
                                                   cTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAggg                                              <  2:92048/100‑1 (MQ=255)
                                                    tGACTGGCTGCGTCTGGCACTCGGCAGCGTCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCTGAGAGGATCGTCTCCggg                                              <  2:177003/99‑1 (MQ=255)
                                                     gACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGATGCCCTGAGAGGATCGTCTCCGGGcg                                            <  1:1563092/100‑2 (MQ=37)
                                                     gACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCa                                            <  2:1559487/100‑1 (MQ=255)
                                                       cTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCAtt                                          >  1:1334373/1‑100 (MQ=255)
                                                          gCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCAttttt                                       >  2:1178409/1‑100 (MQ=255)
                                                          gCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCAttttt                                       <  2:788484/100‑1 (MQ=255)
                                                            tGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGa                                     >  1:1039034/1‑100 (MQ=255)
                                                             gCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGATCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCCCAGGGCATTTTTGat                                    >  1:1173709/1‑100 (MQ=255)
                                                               gTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGatat                                  >  1:6509/1‑100 (MQ=255)
                                                                tCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGatat                                  <  1:1535326/99‑1 (MQ=255)
                                                                tCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATg                                 >  1:407556/1‑100 (MQ=255)
                                                                  tGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATggg                               <  1:837410/100‑1 (MQ=255)
                                                                   ggCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGt                              >  2:1135974/1‑100 (MQ=255)
                                                                   ggCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGt                              <  2:1614688/100‑1 (MQ=255)
                                                                     cACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTAc                            >  2:1169508/1‑100 (MQ=255)
                                                                       cTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACTCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTAc                            <  1:325337/98‑1 (MQ=255)
                                                                         cGGCAGCTCCATTGTTCCGCCTCTGTTGGGGGGTTGTTCTGGTTGTCCTGCTGACGCCAGGAGACGTTCCTCTCAGGGCATTTTTGATATGGGTACTGTc                        <  1:1147025/100‑1 (MQ=255)
                                                                         cGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTc                        <  2:407558/100‑1 (MQ=255)
                                                                          ggCAGCGCCATTGTTCCGCCTCTGTTGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAg                      >  1:657954/1‑100 (MQ=255)
                                                                               cgcCATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTCATATGGGTACTGTCAGAATc                  <  1:1112720/100‑1 (MQ=255)
                                                                                 ccATTGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCa                >  1:782460/1‑100 (MQ=255)
                                                                                     tGTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCgg            <  2:824454/100‑1 (MQ=255)
                                                                                      gTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGTCGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGgc           >  2:1278052/1‑100 (MQ=255)
                                                                                      gTTCCGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGgc           >  2:317327/1‑100 (MQ=255)
                                                                                          cGCCTCTGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGctc       >  2:15786/1‑100 (MQ=255)
                                                                                            ccTCTGTTGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTgg    >  1:23919/1‑100 (MQ=255)
                                                                                             ctctGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTgg    <  2:136625/100‑1 (MQ=255)
                                                                                             ctctGTTGGGGTGGTGTTCTGGTTGTCCTGCTGACGCTCGCAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCGCTgg    <  1:277532/100‑1 (MQ=255)
                                                                                               ctGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCa  >  1:1784477/1‑100 (MQ=255)
                                                                                               ctGTTGGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCa  >  1:1455072/1‑100 (MQ=255)
                                                                                                    |                                                                                                
CAGTGCTGCTACGTAAATTGCTGCCGTGGTATAACCCCTTGCAAAATCTGGCTGACTGGCTGCGTCTGGCACTCGGCAGCGCCATTGTTCCGCCTCTGTTGGGGGGTGTTCTGGTTGTCCTGCTGACGCCCGGAGACGATCCTCTCAGGGCATTTTTGATATGGGTACTGTCAGAATCCATCGGCGCTCTGGCA  >  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: