Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,627,039 G→A 7.0% Q634* (CAG→TAG)  yfgF ← cyclic‑di‑GMP phosphodiesterase, anaerobic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,627,0390GA 7.0% 178.6 / 3.2 57Q634* (CAG→TAG) yfgFcyclic‑di‑GMP phosphodiesterase, anaerobic
Reads supporting (aligned to +/- strand):  ref base G (23/29);  new base A (2/2);  total (26/31)
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.76e-01

CCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGCAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTGCGCCTGCTTAACATTGGTCAGAGCATTACTTTCGGTGACTTCAAAAATAAGTTGCCACGCTTCAATCTGATATTTAGCCAGCAACTGACTGACTTCAACAGGAAAACGAGCCTGACATACCGAGGTTGG  >  NC_000913/2626895‑2627180
                                                                                                                                                |                                                                                                                                             
ccAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGCAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGAtgctgc                                                                                                                                              >  8:161209/1‑146 (MQ=255)
    aCTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGCAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGAAGCTGCAAGGTg                                                                                                                                        <  3:40193/148‑1 (MQ=255)
     cTATTTGACACAATATTGCTGATAAAACTGCCGTCAATTTTAAGCAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGa                                                                                                                                       <  4:182310/148‑1 (MQ=255)
        tttGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGCAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGAtt                                                                                                                                     <  3:26542/147‑1 (MQ=255)
          tGACACAATATTGCGGATAAAAGTGCCGTCAATTTTAAGCAGATGGGCTTTCACTTTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTg                                                                                                                                   <  1:315863/147‑1 (MQ=255)
          agACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGCAGAGCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTg                                                                                                                                   <  6:138950/146‑1 (MQ=255)
            acacAATATTGCGGATAAAACTGCCGTCAATTTTAAGCAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTgcgc                                                                                                                                <  5:63927/148‑1 (MQ=255)
                 atatTGCGGATAAAACTGCCGTCAATTTTAAGCAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTGCGCCTGCt                                                                                                                           >  6:67131/1‑148 (MQ=255)
                     agaggATAAAACTGCCGTCAATTTTAAGCAGATCAGCATTCACATTTTAAAGACGCGCATAGCATGCGTAGCCCGTGACGAAATCATCAATCGAAATCTCGCAGCATAATTCCTGAAGATGCTTCAAGGTGATTTGCGCCTGCTTaa                                                                                                                        >  2:253515/4‑147 (MQ=255)
                         gATAAAACTGCCGTCAATTTTAAGCAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGTAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTGCGCCTGCTTAACATTg                                                                                                                   >  3:27101/1‑148 (MQ=255)
                           tAAAACTGCCGTCAATTTTAAGCAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTGCGCCTGCTTAACATTGGt                                                                                                                 <  8:22501/148‑1 (MQ=255)
                           tAAAACTGCCGTCAATTTTAAGCAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTGCGCCTGCTTAACATTGGt                                                                                                                 <  2:174037/148‑1 (MQ=255)
                             aaaCTGCCGTCAATTGAAAGCAGCTCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTGCGCCTGCTTAACATTGGTc                                                                                                                <  6:77449/147‑1 (MQ=255)
                               aCTGCCGTCAATTTTAAGCAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTGCGCCTGCTTAACATTGGTCaga                                                                                                             <  5:120998/148‑1 (MQ=255)
                                               aGCAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTGCGCCTGCTTAACATTGGTCAGAGCATTACTTTCGGTGa                                                                                             >  7:103394/1‑148 (MQ=255)
                                                gCAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTGCGCCTGCTTAACATTGGTCAGAGCATTACTTTCGGTGAc                                                                                            >  2:258604/1‑148 (MQ=255)
                                                     tCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTGCGCCTGCTTAACATTGGTCAGAGCATTACTTTCGGTGACTTCa                                                                                        >  8:106496/1‑147 (MQ=255)
                                                          aTTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTGCGCCTGCTTAACATTGGTCAGAGCATTACTTTCGGTGACTTCAAAAATa                                                                                  <  4:105152/148‑1 (MQ=255)
                                                            tCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTGCGCCTGCTTAACATTGGTCAGAGc                                                                                                           <  8:291143/121‑1 (MQ=255)
                                                            tCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTGCGCCTGCTTAACATTGGTCAGAGc                                                                                                           >  7:291143/1‑121 (MQ=255)
                                                            tCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTGCGCCTGCTTAACATTGGTCAGAGCATTACTTTCGGTGACTTCAAAAATAAg                                                                                >  5:306163/1‑148 (MQ=255)
                                                            tCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTGCGCCTGCTTAACATTGGTCAGAGCATTACTTTCGGTGACTTCAAAAATAAg                                                                                >  1:194640/1‑148 (MQ=255)
                                                             cacaTTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTGCGCCTGCTTAACATTGGTCAGAGCATTACTTTCGGTGACTTCAAAAATAAGt                                                                               >  3:162290/1‑148 (MQ=255)
                                                             cacaTTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTACAAGGTGATTTGCGCCTGCTTAACATTGGTCAGAGCATTACTTTCg                                                                                                 >  1:310401/1‑130 (MQ=255)
                                                             cacaTTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTACAAGGTGATTTGCGCCTGCTTAACATTGGTCAGAGCATTACTTTCg                                                                                                 >  1:314886/1‑130 (MQ=255)
                                                             cacaTTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTACAAGGTGATTTGCGCCTGCTTAACATTGGTCAGAGCATTACTTTCg                                                                                                 <  2:314886/130‑1 (MQ=255)
                                                             cacaTTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTACAAGGTGATTTGCGCCTGCTTAACATTGGTCAGAGCATTACTTTCg                                                                                                 <  2:310401/130‑1 (MQ=255)
                                                               caTTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTGCGCCTGCTTAACATTGGTCAGAGCATTACTTTCGGTg                                                                                              <  4:114590/131‑1 (MQ=255)
                                                               caTTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTGCGCCTGCTTAACATTGGTCAGAGCATTACTTTCGGTg                                                                                              >  3:114590/1‑131 (MQ=255)
                                                               caTTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTGCGCCTGCTTAACATTGGTCAGAGCATTACTTTCGGTGACTTCAAAAATAAGTTg                                                                             >  5:290741/1‑148 (MQ=255)
                                                                 tttttAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTTGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTGCGCCTGCTTAACATTGGTCAGAGCATTACTTTCGGTGACTTCAAAAATAAGTTGcc                                                                           <  2:194640/148‑1 (MQ=255)
                                                                      aaGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTGCGCCTGCTTAACATTGGTCAGAGCATTACTTTCGGTGACTTCAAAAATAAGTTGCCACGc                                                                       <  7:106496/147‑1 (MQ=255)
                                                                       aGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGTAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTGCGCCTGCTTAACATTGGTCAGAGCATTACTTTCGGTGACTTCAAAAATAAGTTGCCACGCt                                                                      <  4:27101/147‑1 (MQ=255)
                                                                          cgcgCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTGCGCCTGCTTAACATTGGTCAGAGCATTACTTTCGGTGACTTCAAAAATAAGTTGCCACGCTTCaa                                                                  <  1:258604/148‑1 (MQ=255)
                                                                          cgcgCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTGCGCCTGCTTAACATTGGTCAGAGCATTACTTTCGGTGACTTCAAAAATAAGTTGCCACGCTTCaa                                                                  >  6:69164/1‑148 (MQ=255)
                                                                            cgcATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTGCGCCTGCTTAACATTGGTCAGAGCATTACTTTCGGTGACTTCAAAAATAAGTTGCCACGCTTCAATc                                                                >  5:92571/1‑148 (MQ=255)
                                                                                tAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTGCGCCTGCTTAACATTGGTCAGAGCATTACTTTCGGTGACTTCAAAAATAAGTTGCCACGCTTCAATCTGat                                                            <  1:253515/148‑1 (MQ=255)
                                                                                  gCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTGCGCCTGCTTAACATTGGTCAGAGCATTACTTTCGGTGACTTCAAAAATAAGTTGCCACGCTTCAATCTGatat                                                          >  7:244102/1‑148 (MQ=255)
                                                                                      gCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTGCGCCTGCTTAACATTGGTCAGAGCATTACTTTCGGTGACTTCAAAAATAAGTTGCCACGCTTCAATCTGATAttt                                                        <  5:67131/146‑1 (MQ=255)
                                                                                             cGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTGCGCCTGCTTAACATTGGTCAGAGCATTACTTTCGGTGACTTCAAAAATAAGTTGCCACGCTTCAATCTGATATTTAGCCAGCaa                                               >  5:194160/1‑148 (MQ=255)
                                                                                               gTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTGCGCCTGCTTAACATTGGTCAGAGCATTACTTTCGGTGACTTCAAAAATAAGTTGCCACGCTTCAATCTGATATTTAGCCAGCAact                                             >  8:53211/1‑148 (MQ=255)
                                                                                                tGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTGCGCCTGCTTAACATTGGTCAGAGCATTACTTTCGGTGACTTCAAAAATAAGTTGCCACGCTTCAATCTGATATTTAGCCAGCAactg                                            <  2:254873/148‑1 (MQ=255)
                                                                                                 gCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTGCGCCTGCTTAACATTGGTCAGAGCATTACTTTCGGTg                                                                                              >  1:46057/1‑97 (MQ=255)
                                                                                                 gCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTGCGCCTGCTTAACATTGGTCAGAGCATTACTTTCGGTg                                                                                              <  2:46057/97‑1 (MQ=255)
                                                                                                      aaTCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTGCGCCTGCTTAACATTGGTCAGAGCATTACTTTCGGTGACTTCAAAAATAAGTTGCCACGCTTCAATCTGATATTTAGCCAGCAACtgactgac                                      <  4:162290/148‑1 (MQ=255)
                                                                                                       atcatcAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTGCGCCTGCTTAACATTGGTCAGAGCATTACTTTCGGTGACTTCAAAAATAAGTTGCCACGCTGCAATCTGATATTTAGCCAGCAACTgactgact                                     <  6:290741/148‑1 (MQ=255)
                                                                                                        tcatcaATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTGCGCCTGCTTAACATTGGTCAGAGCATTACTTTCGGTGACTTCAAAAATAAGTTGCCACGCTTCAATCTGATATTTAGCCAGCAACTGACTGACtt                                    >  6:127131/1‑148 (MQ=255)
                                                                                                                     tCTGGCAGCCTAGTTCCTGAAGATGCTGCAAGGTGATTTGCGCCTGCTTAACATTGGTCAGAGCATTACTTTCGGTGACTTCAAAAATa                                                                                  >  6:275012/1‑89 (MQ=255)
                                                                                                                     tCTGGCAGCCTAGTTCCTGAAGATGCTGCAAGGTGATTTGCGCCTGCTTAACATTGGTCAGAGCATTACTTTCGGTGACTTCAAAAATa                                                                                  <  5:275012/89‑1 (MQ=255)
                                                                                                                         gCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTGCGCCTGCTTAACATTGGTCAGAGCATTACTTTCGGTGACTTCAAAAATAAGTTGCCACGCTTCAATCTGATATTTAGCCAGCAACTGACTGACTTCAACAGGAAAACGAg                     <  6:306163/146‑1 (MQ=255)
                                                                                                                            gCCTAATTCCTGAAGATGCTGCAAGGTGATTTGCGCCTGCTTAACATTGGTCAGAGCATTACTTTCGGTGACTTCAAAAATAAGTTGcc                                                                           <  6:3035/89‑1 (MQ=255)
                                                                                                                            gCCTAATTCCTGAAGATGCTGCAAGGTGATTTGCGCCTGCTTAACATTGGTCAGAGCATTACTTTCGGTGACTTCAAAAATAAGTTGcc                                                                           >  5:3035/1‑89 (MQ=255)
                                                                                                                            gCCTAATTCCTGAAGATGCTGCAAGGTGATTTGCGCCTGCTTAACATTGGTCAGAGCATTACTTTCGGTGACTTCAAAAATAAGTTGCCACGCTTCAATCTGATATTTAGCCAGCAACTGACTGACTTCAACAGGAAAACGAGCCTGa                <  5:127131/148‑1 (MQ=255)
                                                                                                                             ccTAATTCCTTAAGATGCTGCAAGGTGATTTGCGCCTGCTTAACATTGGTCAGAGCATTACTTTCGGTGACTTCAAAAATAAGTTGCCACGCTTCAATCTGATATTTAGCCAGCAACTGACTGACTTCAACAGGAAAACGAGCCTGAc               <  7:132132/148‑1 (MQ=255)
                                                                                                                                   tCCTGAAGATGCTGCAAGGTGATTTGCGCCTGCTTAACATTGGTCAGAGCATTACTTTCGGTGACTTCAAAAATAAGTTGCCACGCTTCAATCTGATATTTAGCCAGCAACTGACTGACTTCAACAGGAAAACGAGCCTGACATAcc          >  2:82284/1‑147 (MQ=255)
                                                                                                                                        aaGATGCTGCAAGGTGATTTGCGCCTGCTTAACATTGGTCAGAGCATTACTTTCGGTGACTTCAAAAATAAGTTGCCACGCTTCAATCTGATATTTAGCCAGCAACTGACTGACTTCAACAGGAAAACGAGCCTGACATACCGAGGtt    >  4:77859/1‑148 (MQ=255)
                                                                                                                                          gATGCTGCAAGGTGATTTGCGCCTGCTTAACATTGGTCAGAGCATTACTTTCGGTGACTTCAAAAATAAGTTGCCACGCTTCAATCTGATATTTAGCCAGCAACTGACTGACTTCAACAGGAAAACGAGCCTGACATACCGAGGTTgg  <  6:92571/148‑1 (MQ=255)
                                                                                                                                           aTGCTGCAAGGTGATTTGCGCCTGCTTAACATTGGTCAGAGCATTACTTTCGGTGACTTCAAAAATAAGTTGCCACGCTTCAATCTGATATTTAGCCAGCAACTGACTGACTTCAACAGGAAAACGAGCCTGACATACCGAGGTTg   >  7:73938/1‑146 (MQ=255)
                                                                                                                                                |                                                                                                                                             
CCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGCAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTGCGCCTGCTTAACATTGGTCAGAGCATTACTTTCGGTGACTTCAAAAATAAGTTGCCACGCTTCAATCTGATATTTAGCCAGCAACTGACTGACTTCAACAGGAAAACGAGCCTGACATACCGAGGTTGG  >  NC_000913/2626895‑2627180

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