Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,028,159 T→A 20.6% intergenic (‑142/+29) 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,1590TA20.6% 155.7 / 37.8 73intergenic (‑142/+29)yedQ/yodCputative membrane‑anchored diguanylate cyclase/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (31/26);  new base A (7/8);  total (38/35)
Fisher's exact test for biased strand distribution p-value = 7.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.24e-01

CTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGA  >  NC_000913/2028060‑2028258
                                                                                                   |                                                                                                   
cTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGAtgt                                                                                                     >  2:2292883/1‑100 (MQ=255)
cTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTACAGGCATTAGCACACCAGGCCGGAtgt                                                                                                     >  1:1339523/1‑100 (MQ=255)
 tGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGAtgtg                                                                                                    >  1:1880721/1‑100 (MQ=255)
    cGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATc                                                                                                 <  2:845536/100‑1 (MQ=255)
     gATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATcc                                                                                                <  2:525799/100‑1 (MQ=255)
       tAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCgg                                                                                              >  2:680456/1‑100 (MQ=255)
         gCGTACTGTAAAACAGCGATAATTTTATTGTGGTGGCATGATAAGACGCGACTATCATTCATTGCGGCCATTAGCACACCAGGCCGGATGTGATCCGGcc                                                                                            <  1:1107945/100‑1 (MQ=255)
         gCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGcc                                                                                            <  2:1622627/100‑1 (MQ=255)
         gCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGcc                                                                                            <  2:120123/100‑1 (MQ=255)
         gCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGcc                                                                                            <  2:383478/100‑1 (MQ=255)
          cGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCt                                                                                           <  2:1057515/100‑1 (MQ=255)
          cGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCt                                                                                           <  2:2033427/100‑1 (MQ=255)
          cGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCt                                                                                           <  2:2572210/100‑1 (MQ=255)
          cGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCt                                                                                           <  2:2071738/100‑1 (MQ=255)
           gTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACATCAGGCCGGATCACATCCGGCCt                                                                                           <  1:2188068/99‑1 (MQ=38)
           gTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCTg                                                                                          <  2:445724/100‑1 (MQ=255)
              ccGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATg                                                                                       >  1:995473/1‑100 (MQ=255)
              ccGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATg                                                                                       >  2:840620/1‑100 (MQ=255)
              ccGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATg                                                                                       >  2:49214/1‑100 (MQ=255)
              ccGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATg                                                                                       >  2:127384/1‑100 (MQ=255)
              ccGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATg                                                                                       <  2:871160/100‑1 (MQ=255)
              ccGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCTGGTg                                                                                       <  2:64645/100‑1 (MQ=38)
               cGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCTGGTg                                                                                       <  1:1476666/99‑1 (MQ=38)
                gCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCTGgtgtg                                                                                     <  1:275416/100‑2 (MQ=38)
                   aaaCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAAt                                                                                  <  1:2263620/100‑1 (MQ=255)
                   aaaCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAAt                                                                                  >  2:1141016/1‑100 (MQ=255)
                   aaaCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAAt                                                                                  <  1:76155/100‑1 (MQ=255)
                    aaCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATg                                                                                 <  1:423915/100‑1 (MQ=255)
                        gCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATAc                                                                             <  1:289768/100‑1 (MQ=255)
                          gATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATAccc                                                                           >  2:222117/1‑100 (MQ=255)
                               ttttATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTca                                                                      >  1:520708/1‑100 (MQ=255)
                                 ttATAATTGTGGCATGATAAGACGCGACTATCATTCATGGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTcaca                                                                    <  1:2143251/100‑1 (MQ=255)
                                  tataATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCacac                                                                   >  1:2356724/1‑100 (MQ=255)
                                   ataATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACAct                                                                  <  1:2468584/100‑1 (MQ=255)
                                    taATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACActt                                                                 >  2:2532790/1‑100 (MQ=255)
                                    taATGGTAGCATAATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTTATCCGGCCTGATGTTAATGATACCCTTTCACActt                                                                 <  2:1200944/100‑1 (MQ=255)
                                           ggCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGc                                                          >  1:1613609/1‑100 (MQ=255)
                                            gCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCa                                                         <  1:1211549/100‑1 (MQ=255)
                                             cATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAg                                                        >  2:2543629/1‑100 (MQ=255)
                                              aTGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGa                                                       >  1:1222072/1‑100 (MQ=255)
                                                gATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAAc                                                     <  2:1222076/100‑1 (MQ=255)
                                                    aaacacaaCTATCATTAATTGCAGGCATTAGCACACAAGGCCAGGTGTGATCCGCCCTGATGTTAATCATACCCTTTCACACGTCTTCCGCAGAACTAc                                                  >  1:2425299/8‑99 (MQ=255)
                                                    aGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGAct                                                 >  1:541307/1‑100 (MQ=255)
                                                       cgcgACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGActccc                                              >  1:69062/1‑100 (MQ=255)
                                                         cgACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGActcccct                                            <  2:2463895/100‑1 (MQ=255)
                                                             tATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACtcccctcccc                                        <  2:944460/100‑1 (MQ=255)
                                                              aTCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTccccc                                       <  1:1175714/100‑1 (MQ=255)
                                                              aTCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTccccc                                       >  2:1953077/1‑100 (MQ=255)
                                                              aTCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTccccc                                       >  2:2509310/1‑100 (MQ=255)
                                                               tcattcatTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCg                                      <  1:898232/100‑1 (MQ=255)
                                                                 attcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGGCCCACTCCCCTccacc                                       >  1:2394367/1‑97 (MQ=255)
                                                                 attcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGa                                    >  1:1183476/1‑100 (MQ=255)
                                                                     attGCATGCAATAGTACACCAGGCCTGATGTGAGCCGGCCTGATGTTAATTATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAAc                                  <  1:1835583/98‑1 (MQ=255)
                                                                      ttGCAGGCATTAGCACACCAGGCCGGATAGGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAg                               <  2:1624850/100‑1 (MQ=255)
                                                                       tGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGc                              >  2:1916771/1‑100 (MQ=255)
                                                                       tGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGc                              >  2:2027754/1‑100 (MQ=255)
                                                                        gCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCt                             <  2:1906923/100‑1 (MQ=255)
                                                                         cAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTc                            >  2:2515005/1‑100 (MQ=255)
                                                                           ggCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCgg                          >  2:1590346/1‑100 (MQ=255)
                                                                           ggCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCgg                          >  2:1264766/1‑100 (MQ=255)
                                                                           ggCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCgg                          <  2:829638/100‑1 (MQ=255)
                                                                           ggCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCgg                          <  2:827168/100‑1 (MQ=255)
                                                                               cattaacatcAGGCCGGATCACATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCACCAGGAAAAAGCTCGGtt                        >  2:1120519/10‑98 (MQ=21)
                                                                                 ttaacatcAGGCCGGATCACATCCGGCCTGGTGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCAt                    >  2:1646667/8‑100 (MQ=25)
                                                                                   cacaCCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGa                  >  2:290769/1‑100 (MQ=255)
                                                                                   aacatCAGGCCGGATCACATCCGGCCTGGTGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGa                  >  2:309271/2‑100 (MQ=37)
                                                                                   aacatCAGGCCGGATCACATCCGGCCTGGTGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGa                  >  1:1190115/2‑100 (MQ=37)
                                                                                    acatCAGGCCGGATCACATCCGGCCTGGTGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGaa                 >  2:1551993/1‑100 (MQ=37)
                                                                                    acatCAGGCCGGATCACATCCGGCCTGGTGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGaa                 >  2:736515/1‑100 (MQ=37)
                                                                                     cacCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGaaa                >  2:1217046/1‑100 (MQ=255)
                                                                                       tcAGGCCGGATCACATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGaaaa               >  1:2065691/2‑99 (MQ=255)
                                                                                       ccAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGaaaa               >  1:1427148/1‑99 (MQ=255)
                                                                                         aGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCt            >  1:441954/1‑100 (MQ=255)
                                                                                          ggCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCTTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCtt           >  2:2260192/1‑100 (MQ=255)
                                                                                           gCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTc          >  1:1781353/1‑100 (MQ=255)
                                                                                                aTGTGAGCCAGCCTGATGTTAATGATACCCTTGCTCTCTTCCTCCGCAGAACGACACCCCGCCCCCGGTACAAGCTCGGTTTCAGGAAAAGCTTCCCGc      <  1:91976/99‑1 (MQ=255)
                                                                                                 tgtgATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCtt    <  1:1916761/100‑1 (MQ=255)
                                                                                                  gtgCTCCGGCCTGATGTTAATGCTACCCTTTCACACTTCTTCCGCAGACCGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTg   <  1:2524396/100‑1 (MQ=255)
                                                                                                   tgATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACCCCCCTCCCCCGCAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGa  >  1:1467394/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
CTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGA  >  NC_000913/2028060‑2028258

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: