Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,567,472 G→C 8.6% intergenic (‑332/+32) dosC ← / ← yddW diguanylate cyclase, cold‑ and stationary phase‑induced oxygen‑dependent biofilm regulator/liprotein, glycosyl hydrolase homolog

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,567,4720GC 8.6% 191.2 / 8.0 58intergenic (‑332/+32)dosC/yddWdiguanylate cyclase, cold‑ and stationary phase‑induced oxygen‑dependent biofilm regulator/liprotein, glycosyl hydrolase homolog
Reads supporting (aligned to +/- strand):  ref base G (28/25);  new base C (3/2);  total (31/27)
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.49e-01

GTGACAGCGAATATCGTATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCGGGCACAGCATCGTTAAGA  >  NC_000913/1567328‑1567618
                                                                                                                                                |                                                                                                                                                  
gTGACAGCGAATATCGTATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCAcc                                                                                                                                                >  1:725463/1‑149 (MQ=255)
   aCAGCGAATATCGTATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGc                                                                                                                                             <  3:464952/149‑1 (MQ=255)
   aCAGCGAATATCGTATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGc                                                                                                                                             >  1:444440/1‑149 (MQ=255)
     aGCGAATATCGTATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCtt                                                                                                                                           <  1:443529/149‑1 (MQ=255)
     aGCGAATATCGTATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGCTTTACCCACCGGCtt                                                                                                                                           <  4:476748/149‑1 (MQ=255)
       cGAATATCGTATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGt                                                                                                                                         <  2:157140/149‑1 (MQ=255)
             tCGTATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACaaa                                                                                                                                   <  1:274332/149‑1 (MQ=255)
               gTATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGt                                                                                                                                 >  3:527341/1‑149 (MQ=255)
                tATGGTTTTTCATATTCATACATTTTTATGAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGACGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTa                                                                                                                                >  2:598016/1‑149 (MQ=255)
                  tGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAg                                                                                                                              >  2:481326/1‑149 (MQ=255)
                   ggTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGa                                                                                                                             <  1:752848/149‑1 (MQ=255)
                     tttttCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAAt                                                                                                                           >  4:433558/1‑149 (MQ=255)
                                 ataCATTTGTATTAGGGGTTTATGGCGGATTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACt                                                                                                               >  4:114719/1‑149 (MQ=255)
                                  taCATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTg                                                                                                              >  3:631305/1‑149 (MQ=255)
                                    cATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGt                                                                                                                                         >  2:387783/1‑120 (MQ=255)
                                    cATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGt                                                                                                                                         <  1:387783/120‑1 (MQ=255)
                                     aTTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGccc                                                                                                           <  1:654277/149‑1 (MQ=255)
                                     aTTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGccc                                                                                                           >  3:695086/1‑149 (MQ=255)
                                           attaGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGa                                                                                                     <  1:481326/149‑1 (MQ=255)
                                               gggATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCAGCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACttt                                                                                                 <  4:723103/149‑1 (MQ=255)
                                                  aTTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCa                                                                                              <  3:499445/149‑1 (MQ=255)
                                                           tGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGa                                                                                     >  2:142816/1‑149 (MQ=255)
                                                           tGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGCTTTACCCACCGGCTTGTTAc                                                                                                                                      >  4:498070/1‑100 (MQ=37)
                                                           tGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGCTTTACCCACCGGCTTGTTAc                                                                                                                                      <  3:498070/100‑1 (MQ=37)
                                                             tttAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGAcc                                                                                   >  4:358499/1‑149 (MQ=255)
                                                                      gtgtGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGa                                                                          >  1:258557/1‑149 (MQ=255)
                                                                      gtgtGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGa                                                                          <  3:114719/149‑1 (MQ=255)
                                                                      gtgtGGTTAATTTGACCTCAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGa                                                                          <  1:598016/149‑1 (MQ=255)
                                                                      atctgGTTAATTTGGATTTAGTAAACATGATTATTAGGTGGATAGTTTAAGGGGGTTACAAACCGGTTGGGAAAGCCCCGGCTTGTTAAAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGa                                                                          <  1:28574/146‑1 (MQ=255)
                                                                           gTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCtt                                                                                                                                           <  4:463172/79‑1 (MQ=255)
                                                                           gTTAATTTGACTTAAGTAAGCATGATTAGTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCtt                                                                                                                                           >  3:463172/1‑79 (MQ=255)
                                                                            ttAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGc                                                                    >  3:690046/1‑149 (MQ=255)
                                                                             tAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCg                                                                   <  1:534810/149‑1 (MQ=255)
                                                                             tAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCg                                                                   <  3:120732/149‑1 (MQ=255)
                                                                                tttGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGtt                                                                >  3:460216/1‑149 (MQ=255)
                                                                                         agtaagCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGAt                                                       <  4:105194/149‑1 (MQ=255)
                                                                                              gCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTc                                                           >  3:438983/1‑140 (MQ=255)
                                                                                              gCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTcc                                                  >  3:14001/1‑149 (MQ=255)
                                                                                              gCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACGTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTc                                                           <  4:438983/140‑1 (MQ=255)
                                                                                                aTGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTc                                                >  3:211179/1‑149 (MQ=255)
                                                                                                      attaGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCTAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGaa                                          >  3:716656/1‑149 (MQ=255)
                                                                                                           tGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGa                                     >  3:499143/1‑149 (MQ=255)
                                                                                                               aTAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGt                                 >  4:701659/1‑149 (MQ=255)
                                                                                                                   tttAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCg                             >  1:387535/1‑149 (MQ=255)
                                                                                                                       agagGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTaa                         >  3:520978/1‑149 (MQ=255)
                                                                                                                        gagGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAAt                        <  3:165606/149‑1 (MQ=255)
                                                                                                                             tAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCgg                   <  1:142816/149‑1 (MQ=255)
                                                                                                                             tAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCgg                   <  2:258557/149‑1 (MQ=255)
                                                                                                                             tAACAAGCCGGTGCTTTACCCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGAATTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCgg                   >  4:410556/1‑149 (MQ=255)
                                                                                                                              aaCAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCggg                  <  4:520978/149‑1 (MQ=255)
                                                                                                                              aaCAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTAACGGAACAAGATGGTGCCGCTAATTTCggg                  >  1:426706/1‑149 (MQ=255)
                                                                                                                                 aaGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCGGGcac               <  4:460216/149‑1 (MQ=255)
                                                                                                                                 aaGCCGGTGCTTTACCCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCGGGcac               >  3:418472/1‑149 (MQ=255)
                                                                                                                                  aGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCGGGcaca              <  4:14001/149‑1 (MQ=255)
                                                                                                                                         gggTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAg                                                                         >  3:557678/1‑83 (MQ=255)
                                                                                                                                         gggTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAg                                                                         <  4:557678/83‑1 (MQ=255)
                                                                                                                                            tAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCGGGCACAGCATCGTTaa    >  3:144332/1‑149 (MQ=255)
                                                                                                                                              aaGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCGGGCACAGCATCGTTAAGa  >  4:417258/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
GTGACAGCGAATATCGTATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCGGGCACAGCATCGTTAAGA  >  NC_000913/1567328‑1567618

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑