Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,028,160 G→C 16.3% intergenic (‑143/+28) yedQ ← / ← yodC putative membrane‑anchored diguanylate cyclase/uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,028,1600GC16.3% 197.1 / 32.9 74intergenic (‑143/+28)yedQ/yodCputative membrane‑anchored diguanylate cyclase/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (25/37);  new base C (3/9);  total (28/46)
Fisher's exact test for biased strand distribution p-value = 5.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.55e-01

GCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGAC  >  NC_000913/2028063‑2028259
                                                                                                 |                                                                                                   
gCGATAGCGTACCGCAAAACAGCGATAATTTTATAAGGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGAt                                                                                                   <  2:1210626/100‑1 (MQ=255)
 cGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATc                                                                                                  >  1:827032/1‑100 (MQ=255)
 cGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATc                                                                                                  >  2:57624/1‑100 (MQ=255)
 cGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGCGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATc                                                                                                  >  2:1482500/1‑100 (MQ=255)
  gATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATcc                                                                                                 <  2:231042/100‑1 (MQ=255)
    tggcatACCGCAAAACGGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCgg                                                                                               <  2:1977422/95‑1 (MQ=255)
     aGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGc                                                                                              <  1:2017947/100‑1 (MQ=255)
      gCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGcc                                                                                             <  1:943250/100‑1 (MQ=255)
       cGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACATCAGGCCGGATCACATCCGGCCt                                                                                            <  1:2507517/100‑1 (MQ=38)
       cGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCt                                                                                            <  1:649895/100‑1 (MQ=255)
       cGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCt                                                                                            <  1:1308699/100‑1 (MQ=255)
         tACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGa                                                                                          <  2:838847/100‑1 (MQ=255)
         tACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCTgg                                                                                          <  1:1473175/100‑2 (MQ=255)
         tACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCTgg                                                                                          <  1:2103449/100‑2 (MQ=255)
         tACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCTgg                                                                                          <  2:136372/100‑2 (MQ=255)
         tACCGCAAAACAGCGATAATTGTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGa                                                                                          <  2:919857/100‑1 (MQ=255)
          aCCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGAt                                                                                         <  1:2162698/100‑1 (MQ=255)
          aCCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGAt                                                                                         <  2:1212842/100‑1 (MQ=255)
           ccGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATg                                                                                        <  1:2112595/100‑1 (MQ=255)
           ccGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATg                                                                                        >  2:1251475/1‑100 (MQ=255)
           ccGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATg                                                                                        <  2:851382/100‑1 (MQ=255)
           ccGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCTGGTg                                                                                        <  1:1468016/100‑1 (MQ=38)
           ccGCAAAACAGCGATAATGTTATAATGGTGGCACGATAAGACGCGACTATCATCCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCTgg                                                                                          <  2:861522/98‑2 (MQ=255)
            cGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGt                                                                                       <  1:801867/100‑1 (MQ=255)
             gCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCTGgtgtg                                                                                      <  2:635616/100‑2 (MQ=38)
              cAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACGAGGGCGGATGTGATCCGGCCTGATGTTa                                                                                     >  1:2597619/1‑100 (MQ=255)
               aaaaCAACGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTaa                                                                                    <  1:788839/100‑1 (MQ=255)
                aaaCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAAt                                                                                   <  1:427031/100‑1 (MQ=255)
                aaaCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAAt                                                                                   <  1:427006/100‑1 (MQ=255)
                 aaCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATg                                                                                  <  2:609471/100‑1 (MQ=255)
                  aCAGCGATAATTTTATAATGGTGGCATGATAACACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAAt                                                                                   >  2:2119954/1‑98 (MQ=255)
                    aGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGTACACCAGGGCGGATGTGTTCCGGCCTGATGTTAATGAt                                                                                <  2:1000477/99‑1 (MQ=255)
                            ttttATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTca                                                                       >  2:1047799/1‑100 (MQ=255)
                             tttATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTcac                                                                      <  1:2527211/100‑1 (MQ=255)
                             tttATAATGGTGGCATGATAAGACGCGACTATCATACATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTcac                                                                      <  2:158826/100‑1 (MQ=255)
                              ttATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTcaca                                                                     >  2:1678494/1‑100 (MQ=255)
                              ttATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTcaca                                                                     <  1:783755/100‑1 (MQ=255)
                                ataATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACAct                                                                   >  2:2303315/1‑100 (MQ=255)
                                 taATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACActt                                                                  <  2:29017/100‑1 (MQ=255)
                                   aaggtggCATGACAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACCcttc                                                                 <  2:2339982/97‑1 (MQ=255)
                                     ggtggCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACttcttc                                                              <  1:948504/100‑1 (MQ=255)
                                      gtggCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTcc                                                             <  1:513427/100‑1 (MQ=255)
                                        ggCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGc                                                           <  1:1045035/100‑1 (MQ=255)
                                         tcATTATATGATGCGACTATCATTCATTGCAGGCATTAGCACATCAGGCCGTATGTGATCCTTCCTGATCTTAATGTTACCTTTTCCCACTTCTTCCGc                                                           <  2:367772/98‑1 (MQ=255)
                                          cATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAg                                                         <  2:459317/100‑1 (MQ=255)
                                             gATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAAc                                                      <  2:1546983/100‑1 (MQ=255)
                                               tAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGa                                                    >  2:1090334/1‑100 (MQ=255)
                                                 aGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGAct                                                  >  2:1557527/1‑100 (MQ=255)
                                                   aCGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGActcc                                                <  1:1150972/100‑1 (MQ=255)
                                                    cgcgACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGActccc                                               >  2:1129960/1‑100 (MQ=255)
                                                    cgcgACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGActcc                                                <  2:1121516/99‑1 (MQ=255)
                                                               ttcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGaa                                    <  1:1423426/100‑1 (MQ=255)
                                                                tcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAAc                                   <  1:1176937/100‑1 (MQ=255)
                                                                     gCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCt                              >  1:665385/1‑100 (MQ=255)
                                                                     gCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTCCCGCAGAACGACTCCCCTCCCCCGGAACAAGCt                              >  2:298275/1‑100 (MQ=255)
                                                                      cAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTc                             >  2:47325/1‑100 (MQ=255)
                                                                       aGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCg                            <  2:1056806/100‑1 (MQ=255)
                                                                       aGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCg                            >  1:1851740/1‑100 (MQ=255)
                                                                        ggCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCgg                           >  1:852457/1‑100 (MQ=255)
                                                                          cATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGtt                         <  2:999630/100‑1 (MQ=255)
                                                                          cATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGtt                         >  1:1606158/1‑100 (MQ=255)
                                                                            cattaacatcAGGCCGGATCACATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAGCGACTCCCCTCCCCCGGAACAAGCTCGGttt                        >  2:1828752/10‑99 (MQ=25)
                                                                             attaacatcAGGCCGGATCACATCCGGCCTGGTGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCa                      >  2:372688/9‑100 (MQ=25)
                                                                              aGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCAt                     <  1:563293/100‑1 (MQ=255)
                                                                              aGCACACCACGCCGGATGTGATCCGCCCGTATGTTAATGATCCCCTTTCACACTCCTTCCGCAGACCGGCTCCCCTCCCCCGGAACAAGCTCGGTTTCa                      <  2:1715415/99‑1 (MQ=255)
                                                                                cacaCCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGa                   >  1:1453765/1‑100 (MQ=255)
                                                                                     cAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGc              >  1:2225066/1‑100 (MQ=255)
                                                                                      aGGCCGGATCACATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCt             >  2:2143815/1‑100 (MQ=255)
                                                                                       ggCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCtt            >  1:229249/1‑100 (MQ=255)
                                                                                          cGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTccc         >  2:2202056/1‑100 (MQ=255)
                                                                                          cGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTccc         >  1:572372/1‑100 (MQ=255)
                                                                                           ggATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCg        >  1:660386/1‑100 (MQ=255)
                                                                                            gATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCTCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGc       <  1:1839806/100‑1 (MQ=255)
                                                                                            gATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGc       >  2:1591785/1‑100 (MQ=255)
                                                                                              tgtgATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCtt     >  2:17996/1‑100 (MQ=255)
                                                                                               gtgATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTg    <  2:1550780/100‑1 (MQ=255)
                                                                                                 gATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGACCAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGAc  >  1:1155181/1‑100 (MQ=255)
                                                                                                 gATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGAc  >  1:973854/1‑100 (MQ=255)
                                                                                                 gATCCGGCCTGATGGTAATGATCCCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCACGGTTTCATGAAAAGCTTCCAGCTTGAc  >  2:2006057/1‑100 (MQ=255)
                                                                                                 |                                                                                                   
GCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGAC  >  NC_000913/2028063‑2028259

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