Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,965,105 C→T 7.8% F439F (TTC→TTT)  yegE → phosphodiesterase domain‑containing diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,965,1050CT 7.8% 147.1 / 2.7 53F439F (TTC→TTT) yegEphosphodiesterase domain‑containing diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base C (25/22);  new base T (2/2);  total (28/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

GGCAAGCAACCCACGATGCGCTGACCGGGCTGGCCAACCGCCGCGAGTTCGAATACCGC‑C‑TGGAGCAGGCGCTCAACGACCTGGCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCGGCACATCTGCGCGGTGCTGCAATCCGGCCTGCGTGAGGGCGATACGCTGGCGCGGCTGGGCGGCGATGAATTTGGCGTGCTGCTGGAAAACTGCCC  >  NC_002947/1964972‑1965251
                                                                                                                                       |                                                                                                                                                  
ggCAAGCAACCCACGATGCGCTGACCGGGCTGGCCAACCGCCGCGAGTTCGAATACCGC‑C‑TGGAGCAGGCGCTCAACGACCTGGCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGac                                                                                                                                     >  8:58325/1‑149 (MQ=255)
 gcaagcaaCCCACGATGCGCTGACCGGGCTGGCCAACCGCCGCGAGTTCGAATACCGC‑C‑TGGAGCAGGCGCTCAACGACGTGGCGCGCCTGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGaca                                                                                                                                    <  4:287694/149‑1 (MQ=255)
  caagcaaCCCACGATGCGCTGACCGGGCTGGCCAACCGCCGCGAGTTCGAATACCGC‑C‑TGGAGCAGGCGCTCAACGACCTGTCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCTACGacac                                                                                                                                   <  5:246926/149‑1 (MQ=255)
   aagcaaCCCACGATGCGCTGACCGGGCTGGCCAACCGCCGCGAGTTCGAATACCGC‑C‑TGGAGCAGGCGCTCAACGACCTGGCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGCTTGAACTGGTcga                                                                                                                                         <  4:325492/142‑3 (MQ=255)
   aagcaaCCCACGATGCGCTGACCGGGCTGGCCAACCGCCGCGAGTTCGAATACCGC‑C‑TGGAGCAGGCGCTCAACGACCTGGCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGCTTGAACTGGTcga                                                                                                                                         >  3:325492/1‑140 (MQ=255)
       aaCCCACGATGCGCTGACCGGGCTGGCCAACCGCCGCGAGTTCGAATACCGC‑C‑TGGAGCAGGCGCTCAACGACCTGGCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCtgtg                                                                                                                              <  5:11982/149‑1 (MQ=255)
       aaCCCACGATGCGCTGACCGGGCTGGCCAACCGCCGCGAGTTCGAATACCGC‑C‑TGGAGCAGGCGCTCAACGACCTGGCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGGCAACGACACCtgtg                                                                                                                              >  6:267278/1‑149 (MQ=255)
           caggttGCGCTTACCTGGTGGGCCAACCGCCGGTATTTGTAATACCGCTCAAGTAGGAGGCGCTCAACGACCTGGCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCTGTGGc                                                                                                                            <  1:111800/144‑1 (MQ=255)
              gATGCGCTGACCGGGCTGGCCAACCGCCGCGAGTTCGAATACCGC‑C‑TGGAGCAGGCGCTCAACGACCTGGCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGCTTGAACTGGTcgag                                                                                                                                        >  2:238984/1‑129 (MQ=255)
              gATGCGCTGACCGGGCTGGCCAACCGCCGCGAGTTCGAATACCGC‑C‑TGGAGCAGGCGCTCAACGACCTGGCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGCTTGAACTGGTcgag                                                                                                                                        <  1:238984/132‑4 (MQ=255)
                         cGGGCTGGCCAACCGCCGCGAGTTCGAATACCGC‑C‑TGGAGCAGGCGCTCAACGACCTGGCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGAt                                                                                                             >  4:69498/1‑148 (MQ=255)
                          gggCTGGCCAACCGCCGCGAGTTCGAATACGGC‑C‑TGGAGCAGGCGCTCAACGACCTGGCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGa                                                                                                           <  2:304064/149‑1 (MQ=255)
                          gggCTGGCCAACCGCCGCGAGTTCGAATACCGC‑C‑TGGAGCAGGCGCTCAACGACCTGGCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGa                                                                                                           <  3:69498/149‑1 (MQ=255)
                          gggCTGGCCAACCGCCGCGAGTTCGAATACCGC‑C‑TGGAGCAGGCGCTCAACGACCTGGCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGa                                                                                                           <  7:206543/149‑1 (MQ=255)
                            gCTGGCCAACCGCCGCGAGTTCGAATACCGC‑C‑TGGAGCAGGCGCTCAACGACCTGGCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGt                                                                                                         <  7:58325/149‑1 (MQ=255)
                                gCCAACCGCCGCGAGTTCGAATACCGC‑C‑TGGAGCAGGCGCTCAACGACCTGGCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTtgc                                                                                                      <  1:273943/148‑1 (MQ=255)
                                 ccAACCGGCGCGAGTTCGTTTACCGC‑C‑TGGAGCAGGCGCTCAACGACCTGGCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTGGACCTGTTGAAGCTGGTCAACGACACCTGTTGCCACGCTACCGGCGATGAGTtgctg                                                                                                    <  4:99563/149‑1 (MQ=255)
                                     ccgccgCGAGTTCGAATACCGC‑C‑TGGAGCAGGCGCTCAACGACCTGGCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCgg                                                                                                 >  4:178409/1‑148 (MQ=255)
                                           cgAGTTCGAATACCGC‑C‑TGGAGCAGGCGCTCAACGACCTGGCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCtgtg                                                                                                                              >  2:273850/1‑113 (MQ=255)
                                           cgAGTTCGAATACCGC‑C‑TGGAGCAGGCGCTCAACGACCTGGCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCtgtg                                                                                                                              <  1:273850/113‑1 (MQ=255)
                                             aGTTCGAATACCGC‑C‑TGGAGCAGGCGCTCAACGACCTGGCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCGGCACATCTgc                                                                                        <  3:178409/149‑1 (MQ=255)
                                                    aTACCGC‑C‑TGGAGCAGGCGCTCAACGACCTGGCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCGGCACATCTGCGCGGtgc                                                                                 <  3:149270/149‑1 (MQ=255)
                                                       ccGC‑C‑TGGAGCAGGCGCTCAACGACCTGGCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCGGCACATCTgcgc                                                                                      >  7:323622/1‑141 (MQ=255)
                                                       ccGC‑C‑TGGAGCAGGCGCTCAACGACCTGGCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCGGCACATCTgcgc                                                                                      <  8:323622/141‑1 (MQ=255)
                                                       ccGC‑C‑TGGAGCAGGCGCTCAACGACCTGGCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCGGCACATCTGCGCGGtgctg                                                                               >  6:304460/1‑148 (MQ=255)
                                                                gAGCAGGCGCTCAACGACCTGGCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCGGCACATCTGCGCGGTGCTGCAATCCgg                                                                       >  7:181261/1‑149 (MQ=255)
                                                                     ggCGCTCAACGACCTGGCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCGGCACATCTGCGCGGTGCTGCAATCCGGCCTGc                                                                  <  8:181261/149‑1 (MQ=255)
                                                                       cgcTCACCGACCTGGCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTg                                                                                                                                              <  2:284299/71‑1 (MQ=255)
                                                                                  cTGGCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTTCTGCGGCACATCTGCGCGGTGCTGCAATCCGGCCTGCGTGAGGGCGATAc                                                     >  6:227844/1‑149 (MQ=255)
                                                                                  cTGGCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCGGCACAt                                                                                            <  6:58768/110‑1 (MQ=255)
                                                                                  cTGGCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCGGCACAt                                                                                            >  5:58768/1‑110 (MQ=255)
                                                                                            aGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCGGCACATCTGCGCGGTGCTGCAATCCGGCCTGCGTGAGGGCGATACGCTGGCGCgg                                           >  3:214901/1‑149 (MQ=255)
                                                                                                cGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCGGCACATCTGCGCGGTGCTGCAATCCGGCCTGCGTGAGGGCGATACGCTGGCGCGGCTgg                                       >  2:177998/1‑149 (MQ=255)
                                                                                                  ggcggcATTCGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCGGCACATCTGCGCGGTGCTGCAATCCGGCCTGCGTGAGGGCGATACGCTGGCGCGGCTGggc                                     >  1:120461/1‑149 (MQ=255)
                                                                                                        aTTCGCTGATGTTCCTCGACCTCGACCAGTTAAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCGGCACATATGCGCGGTGCTGCAATCCGGCCTGCGTGAGGGCGATACGCTGGCGCGGCTGGGCGGCGAt                               >  3:133779/1‑149 (MQ=255)
                                                                                                           cGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGAt                                                                                                             <  5:255793/68‑1 (MQ=255)
                                                                                                           cGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGAt                                                                                                             >  6:255793/1‑68 (MQ=255)
                                                                                                            gCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCGGCACATCTGCGCGGTGCTGCAATCCGGCCTGCGTGAGGGCGATACGCTGGCGCGGCTGGGCGGCGATg                              >  4:213605/1‑146 (MQ=255)
                                                                                                            gCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCGGCACATCTGCGCGGTGCTGCAATCCGGCCTGCGTGAGGGCGATACGCTGGCGCGGCTGGGCGGCGATg                              <  3:213605/146‑1 (MQ=255)
                                                                                                                aTGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCGGCACATCTGCGCGGTGCTGCAATCCGGCCTGCGTGAGGGCGATACGCTGGCGCGGCTGGGCGGCGATGAATTTg                        <  2:120461/148‑1 (MQ=255)
                                                                                                                  gTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCGGCACATCTGCGCGGTGCTGCAATCCGGCCTGCGTGAGGGCGATACGCTGGCGCGGCTGGGCGGCGATGAATTTGGCg                     >  1:129848/1‑149 (MQ=255)
                                                                                                                   ttCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCGGCACATCTGCGCGGTGCTGCAATCCGGCCTGCGTGAGGGCGATACGCTGGCGCGGCTGGGCGGCGATGAATTTGGCg                     <  5:227844/148‑1 (MQ=255)
                                                                                                                     cctcgacctcgaCCAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCGGCACTTCTGCGCGGTGCTGCAATCCGGCCTGCGTGAGGGCGATACGCTGGCGCGGCTGGGCGGCGATGAATTTGGCGtgc                  >  4:9102/1‑149 (MQ=255)
                                                                                                                      ctcgacctcgacCAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCGGCACATCTGCGCGGTGCTGCAATCCGGCCTGCGTGAGGGCGATACGCTGGCGCGGCTGGGCGGCGATGAATTTGGCGtgct                 >  3:335599/1‑149 (MQ=255)
                                                                                                                       tcgacctcgaccAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCGGCACATCTGCGCGGTGCTGCAATCCGGCCTGCGTGAGGGCGATACGCTGGCGCGGCTGGGCGGCGATGAATTTGGCGtgctg                >  3:79073/1‑149 (MQ=255)
                                                                                                                         gacctcgaccAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCGGCACATCTGCGCGGTGCTGCAATCCGGCCTGCGTGAGGGCGATACGCTGGCGCGGCTGGGCGGCGATGAATTTGGCGTgctgct              >  6:83523/1‑149 (MQ=255)
                                                                                                                            ctcgcccAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCGGCACATCTGCGCGGTGCTGCAATCCGGCCTGCGTGAGGGCGATACGCTGGCGCGGCTGGGCGGCGATGAATTTGGCGTGCTGCTgg            >  3:305924/1‑148 (MQ=255)
                                                                                                                                  cAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCGGCACATCTGCGCGGTGCTGCAATCCGGCCTGCGTGAGGGCGATACGCTGGCGCGGCTGGGCGGCGATGAATTTGGCGTGCTGCTGGAAAACTg     >  6:41281/1‑149 (MQ=255)
                                                                                                                                  cAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCGGCACATCTGCGCGGTGCTGCAATCCGGCCTGCGTGAGGGCGATACGCTGGCGCGGCTGGGCGGCGATGAATTTGGCGTGCTGCTGGAAAACTg     >  6:36115/1‑149 (MQ=255)
                                                                                                                                   aGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCGGCACATCTGCGCGGTGCTGCAATCCGGCCTGCGTGAGGGCGATACGCTGGCGCGGCTGGGCGGCGATGAATTTGGCGTGCTGCTGGAAAACTgc    <  4:79073/149‑1 (MQ=255)
                                                                                                                                     ttCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCGGCACATCTGCGCGGTGCTGCAATCCGGCCTGCGTGAGGGCGATACGCTGGCGCGGCTGGGCGGCGATGAATTTGGc                      >  2:218275/1‑129 (MQ=255)
                                                                                                                                     ttCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCGGCACATCTGCGCGGTGCTGCAATCCGGCCTGCGTGAGGGCGATACGCTGGCGCGGCTGGGCGGCGATGAATTTGGc                      <  1:218275/129‑1 (MQ=255)
                                                                                                                                     ttCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCGGCACATCTGCGCGGTGCTGCAATCCGGCCTGCGTGAGGGCGATACGCTGGCGCGGCTGGGCGGCGATGAATTTGGCGTGCTGCTGGAAAACTgccc  >  3:215123/1‑149 (MQ=255)
                                                                                                                                       |                                                                                                                                                  
GGCAAGCAACCCACGATGCGCTGACCGGGCTGGCCAACCGCCGCGAGTTCGAATACCGC‑C‑TGGAGCAGGCGCTCAACGACCTGGCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCGGCACATCTGCGCGGTGCTGCAATCCGGCCTGCGTGAGGGCGATACGCTGGCGCGGCTGGGCGGCGATGAATTTGGCGTGCTGCTGGAAAACTGCCC  >  NC_002947/1964972‑1965251

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