Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,093,980 G→A 20.8% D369N (GAC→AAC)  ycdT → diguanylate cyclase, membrane‑anchored

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,093,9800GA20.8% 155.8 / 34.8 77D369N (GAC→AAC) ycdTdiguanylate cyclase, membrane‑anchored
Reads supporting (aligned to +/- strand):  ref base G (24/37);  new base A (8/8);  total (32/45)
Fisher's exact test for biased strand distribution p-value = 5.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01

CAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGATAATCCTGAATACGCTATACCACAAAA  >  NC_000913/1093881‑1094077
                                                                                                   |                                                                                                 
cAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGTGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACgg                                                                                                   <  1:1128481/100‑1 (MQ=255)
cAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACgg                                                                                                   <  1:1006057/100‑1 (MQ=255)
cAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACgg                                                                                                   <  2:573511/100‑1 (MQ=255)
  ggTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGAc                                                                                                 <  1:2335455/100‑1 (MQ=255)
   gTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGAACa                                                                                                <  2:2169839/100‑1 (MQ=255)
     gATAAAAACAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGTGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGAACATc                                                                                              <  2:1784753/100‑1 (MQ=255)
     gATAAAAACAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATc                                                                                              <  1:2023026/100‑1 (MQ=255)
      aTAAAAACAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCg                                                                                             >  2:2503301/1‑100 (MQ=255)
       tAAAAACAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGa                                                                                            >  1:390209/1‑100 (MQ=255)
       tAAAAACAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGa                                                                                            >  2:1188709/1‑100 (MQ=255)
        aaaaaCAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGACGAGTTTGGCGTCTTGCTGACGGACATCGAt                                                                                           >  1:2344264/1‑100 (MQ=255)
         aaaaCAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGCAGAGTTTGGCGTCTTGCTGACGGACATCGATa                                                                                          >  2:858463/1‑100 (MQ=255)
          aaaCAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATAc                                                                                         <  2:1984550/100‑1 (MQ=255)
          aaaCAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATAc                                                                                         <  1:1983011/100‑1 (MQ=255)
          aaaCAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATAc                                                                                         <  1:191754/100‑1 (MQ=255)
           aaCAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACt                                                                                        <  1:1108829/100‑1 (MQ=255)
            aCAGTGGTGAACATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTg                                                                                       <  1:1904340/100‑1 (MQ=255)
              aGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGaa                                                                                     <  2:2207362/100‑1 (MQ=255)
                 ggtgAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGt                                                                                  >  2:1422098/1‑100 (MQ=255)
                    gAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGAACATCGATACTGAACGTGCg                                                                               <  2:1920349/100‑1 (MQ=255)
                    gAATATCATTCGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGTGTGTTGCTGACGGACATCGATACTGAACGTGCg                                                                               <  1:1700952/100‑1 (MQ=255)
                          cATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCt                                                                         <  1:2409030/100‑1 (MQ=255)
                           aTTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTGGGCGGCGAAGAGTATGGCGTCTTGCTGACGGACATCGATACTTACCGTGCGAAAGCtt                                                                        >  1:1262976/1‑100 (MQ=255)
                           aTTGGGAAAAGCAGACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGc                                                                          <  1:2147105/98‑1 (MQ=255)
                            ttGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCttt                                                                       <  2:1610021/100‑1 (MQ=255)
                             tGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTa                                                                      <  2:2165504/100‑1 (MQ=255)
                             tGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTa                                                                      <  1:1163190/100‑1 (MQ=255)
                               ggAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGAGAAAGCTTTAGc                                                                    <  2:1705257/100‑1 (MQ=255)
                                gAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCg                                                                   >  2:1705986/1‑100 (MQ=255)
                                gAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCg                                                                   >  1:1727678/1‑100 (MQ=255)
                                gAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGAACATCGATACTGAACGTGCGAAAGCTTTAGCg                                                                   >  1:640625/1‑100 (MQ=255)
                                gAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGAACATCGATACTGAACGTGCGAAAGCTTTAGCg                                                                   >  1:1018770/1‑100 (MQ=255)
                                 aaaaGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCgg                                                                  <  2:769254/100‑1 (MQ=255)
                                   aaGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGaa                                                                <  1:1018447/100‑1 (MQ=255)
                                      cATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGAACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAgg                                                             >  2:944493/1‑100 (MQ=255)
                                       aTACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGa                                                            >  1:937874/1‑100 (MQ=255)
                                        tACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGAt                                                           >  2:138294/1‑100 (MQ=255)
                                          cGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTc                                                         >  2:2576345/1‑100 (MQ=255)
                                          cGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTc                                                         >  2:978094/1‑100 (MQ=255)
                                           gACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGTAAGCTTTAGCGGAAAGGATTCg                                                        <  1:1505793/100‑1 (MQ=255)
                                           gACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGAACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCg                                                        >  1:751107/1‑100 (MQ=255)
                                            aCCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCgg                                                       >  2:2038610/1‑100 (MQ=255)
                                              cAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGa                                                     <  2:896586/100‑1 (MQ=255)
                                                gATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGAACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGaaa                                                   >  2:1279471/1‑100 (MQ=255)
                                                  tGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGAACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAAt                                                 <  1:1498952/100‑1 (MQ=255)
                                                   gACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATg                                                >  2:768327/1‑100 (MQ=255)
                                                     ccTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGAACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGtt                                              >  2:1039948/1‑100 (MQ=255)
                                                      cTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTg                                             >  1:2575269/1‑100 (MQ=255)
                                                           aGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGt                                        >  1:1713006/1‑100 (MQ=255)
                                                            gcgcgcGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGtt                                       >  2:1697956/1‑100 (MQ=255)
                                                             cgcgcgACGGCGGCGAAGAGTTTGGCGTATTGCTGAAGGACATCGATACTGAACGTGCGAAAGCTTCAGCGGAAAGGATGCGGGAAAAGGTTGAGCGttt                                      >  1:459657/1‑100 (MQ=255)
                                                                    cggcggCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGc                               <  2:1240158/100‑1 (MQ=255)
                                                                    cggcggCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGc                               >  2:184905/1‑100 (MQ=255)
                                                                     ggcggcGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCg                              >  1:2141814/1‑100 (MQ=255)
                                                                       cggcgAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGAt                            <  2:2470385/100‑1 (MQ=255)
                                                                       cggcgAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGAt                            <  1:1188705/100‑1 (MQ=255)
                                                                       cggcgAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGAt                            <  1:2149078/100‑1 (MQ=255)
                                                                       cggcgAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGAt                            <  1:657632/100‑1 (MQ=255)
                                                                       cggcgAAGAGTTTGGCGTCTTGCTGACGAACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGAt                            <  2:443975/100‑1 (MQ=255)
                                                                           gAAGAGTTTGGGGTCTTGCTGACGGACATCGATATTGAACGTGCGAATGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAAGTGGGGATAATc                        <  1:1362255/100‑1 (MQ=255)
                                                                            aaGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGATAATcc                       <  2:1980980/100‑1 (MQ=255)
                                                                            aaGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGATAATcc                       <  1:1352557/100‑1 (MQ=255)
                                                                            aaGAGTTTGGCGTCTTGCTGACGAACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGATAATcc                       <  2:617535/100‑1 (MQ=255)
                                                                             agagTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGATAATCCt                      <  2:1836376/100‑1 (MQ=255)
                                                                             agagTTTGGCGTCTTGCTGACGAACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGATAATCCt                      >  1:2590025/1‑100 (MQ=255)
                                                                              gagTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGATAATCCTg                     <  2:1227334/100‑1 (MQ=255)
                                                                              gagTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGATAATCCTg                     <  1:1521745/100‑1 (MQ=255)
                                                                              gagTTTGGCGTCTTGCTGACGAACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGATAATCCTg                     <  2:2407638/100‑1 (MQ=255)
                                                                               agTTTGGCGTCTTGCTGACGAACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGATAATCCTGa                    <  1:1260534/100‑1 (MQ=255)
                                                                                gTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGATAATCCTGaa                   <  2:1944657/100‑1 (MQ=255)
                                                                                 tttGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGATAATCCTGAAt                  <  2:2557079/100‑1 (MQ=255)
                                                                                     aggtgttGCGGACGGACGCGGGTGCTGGACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGATAATCCTGAATACGc              <  1:2270095/95‑1 (MQ=255)
                                                                                       ggcTTGCTGGCGGACATCGATAGTGAAGGTGCGAATGCTTTAGTTGAAAGGATTCGTCAAAATTTTGAGCGTTTAACTGTCGATAATCCTGAATATGc              <  1:2082859/96‑1 (MQ=255)
                                                                                         cTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGATAATCCTGAATACGCtata          >  1:795462/1‑100 (MQ=255)
                                                                                         cTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGATAATCCTGAATACGCtata          >  2:1204589/1‑100 (MQ=255)
                                                                                             cagACGGACATCGATACTGAACGTGCGAAAGATTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGATAATCCTGAATACGCTATACcac      >  2:1745706/3‑100 (MQ=255)
                                                                                             cTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGATAATCCTGAATACGCTATACcac      <  1:952305/100‑1 (MQ=255)
                                                                                              tGACGAACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGATAATCCTGAATACGCTATACcaca     >  1:1743741/1‑100 (MQ=255)
                                                                                               gACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGATAATCCTGAATACGCTATACCACaa    <  1:2183165/100‑1 (MQ=255)
                                                                                                 cGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGATAATCCTGAATACGCTATACCACaaaa  <  2:1378263/100‑1 (MQ=255)
                                                                                                   |                                                                                                 
CAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGATAATCCTGAATACGCTATACCACAAAA  >  NC_000913/1093881‑1094077

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