Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,627,887 (G)5→6 coding (1052/2244 nt) 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,8821.G97.9% 157.8 / ‑4.0 47coding (1057/2244 nt)yfgFcyclic‑di‑GMP phosphodiesterase, anaerobic
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (27/19);  total (28/19)
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.99e-01

CCGAAAGCTGATAAACATCTTCACCCGGTTCCAGCAAGGGTGACAGCCAGTGAGAAAGTTTTTGCTTGTATTGAATCCGCAGCATGATGCCATAGTTCTTAACCAGCATTTCCATGCCAGGGATGCGTAAATAACAAAGCGCAGACCAGGGGGCATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATAAAAGAGA  >  NC_000913/2627735‑2628025
                                                                                                                                                    |                                                                                                                                                
ccGAAAGCTGATAAACATCTTCACCCGGTTCCAGCAAGGGTGACAGCCAGTGAGAAAGTTTTTGCTTGTATTGAATCCGCAGCATGATGCCATAGTTCTTAACCAGCATTTCCATGCCAGGGATGCGTAAATAACAAAGCGCAGACCAggg                                                                                                                                               <  2:55243/151‑1 (MQ=255)
         gATAAACATCTTCACCCGGTTCCAGCAAGGGTGACAGCCAGTGAGAAAGTTTTTGCTTGTATTGAATCCGCAGCATGATGCCATAGTTCTTAACCAGCATTTCCATGCCAGGGATGCGTAAATAACAAAGCGCAGACCAggggggcatc                                                                                                                                        <  2:63661/149‑6 (MQ=255)
               cATCTTCACCCGGTTCCAGCAAGGGTGACAGCCAGTGAGAAAGTTTTTGCTTGTATTGAATCCGCAGCATGATGCCATAGTTCTTAACCAGCATTTCCATGCCAGGGATGCGTAAATAACAAAGCGCAGACCAGGGGGGCATCACGCAAcg                                                                                                                                >  1:261387/1‑151 (MQ=255)
                aTCTTCACCCGGTTCCAGCAAGGGTGACAGCCAGTGAGAAAGTTTTTGCTTGTATTGAATCCGCAGCATGATGCCATAGTTCTTAACCAGCATTTCCATGCCAGGGATGCGTAAATAACAAAGCGCAGACCAGGGGGGCATCACGCAAcgc                                                                                                                               >  1:231095/1‑151 (MQ=255)
                   ttCACCCGGTTCCAGCAAGGGTGACAGCCAGTGAGAAAGTTTTTGCTTGTATTGAATCCGCAGCATGATGCCATAGTTCTTAACCAGCATTTCCATGCCAGGGATGCGTAAATAACAAAGCGCAGACCAGGGGGGCATCACGCAAcgcgcg                                                                                                                            >  1:293768/1‑151 (MQ=255)
                   ttCACCCGGTTCCAGCAAGGGTGACAGCCAGTGAGAAAGTTTTTGCTTGTATTGAATCCGCAGCATGATGCCATAGTTCTTAACCAGCATTTCCATGCCAGGGATGCGTAAATAACAAAGCGCAGACCAGGGGGGCATCACGCAAcgcgcg                                                                                                                            >  1:53366/1‑151 (MQ=255)
                     cACCCGGTTCCAGCAAGGGTGACAGCCAGTGAGAAAGTTTTTGCTTGTATTGAATCCGCAGCATGATGCCATAGTTCTTAACCAGCATTTCCATGCCAGGGATGCGTAAATAACAAAGCGCAGACCAGGGGGGCATCACGCAACGCGCGAt                                                                                                                          <  2:382576/151‑1 (MQ=255)
                      aCCCGGTTCCAGCAAGGGTGACAGCCAGTGAGAAAGTTTTTGCTTGTATTGAATCCGCAGCATGATGCCATAGTTCTTAACCAGCATTTCCATGCCAGGGATGCGTAAATAACAAAGCGCAGACCAGGGGGGCATCACGCAACGCGCGAtt                                                                                                                         >  1:844074/1‑151 (MQ=255)
                         cGGTTCCAGCAAGGGTGACAGCCAGTGAGAAAGTTTTTGCTTGTATTGAATCCGCAGCATGATGCCATAGTTCTTAACCAGCATTTCCATGCCAGGGATGCGTAAATAACAAAGCGCAGACCAGGGGGGCATCACGCAACGCGCGATTc                                                                                                                        >  2:267055/1‑149 (MQ=255)
                          ggTTCCAGCAAGGGTGACAGCCAGTGAGAAAGTTTTTGCTTGTATTGAATCCGCAGCATGATGCCATAGTTCTTAACCAGCATTTCCATGCCAGGGATGCGTAAATAACAAAGCGCAGACCAGGGGGGCATCACGCAACGCGCGATTCAgg                                                                                                                     >  2:407117/1‑151 (MQ=255)
                                  aaaGGGTGACAGTCAGTGAGAAAGTTTTTGCTTGTATTGAATCCGCAGCATGATGTCATAGTTCTTAACCAGCATTTCCATGCAAGGGATGCGTAAATTACAAAGCGCAGACCAGGGGGGCATCACGCAACGCGCGATTCAGGGCGCGa                                                                                                               >  1:534079/2‑149 (MQ=255)
                                       gTGACAGCCAGTGAGAAAGTTTTTGCTTGTATTGAATCCGCAGCATGATGCCATAGTTCTTAACCAGCATTTCCATGCCAGGGATGCGTAAATAACAAAGCGCAGACCAGGGGGGCATCACGCAACGCGCGATTCAGGGCGCGAACATTTg                                                                                                        >  2:446450/1‑151 (MQ=255)
                                          aCAGCCAGTGAGAAAGTTTTTGCTTGTATTGAATCCGCAGCATGATGCCATAGTTCTTAACCAGCATTTCCATGCCAGGGATGCGTAAATAACAAAGCGCAGACCAGGGGGGCATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCa                                                                                                     <  2:397254/151‑1 (MQ=255)
                                              ccAGTGAGAAAGTTTTTGCTTGTATTGAATCCGCAGCATGATGCCATAGTTCTTAACCAGCATTTCCATGCCAGGGATGCGTAAATAACAAAGCGCAGACCAGGGGGGCATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGAt                                                                                                  >  2:614963/1‑150 (MQ=255)
                                              ccAGTGAGAAAGTTTTTGCTTGTATTGAATCCGCAGCATGATGCCATAGTTCTTAACCAGCATTTCCATGCCAGGGATGCGTAAATAACAAAGCGCAGACCAGGGGGGCATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATg                                                                                                 >  1:820739/1‑151 (MQ=255)
                                                     gaAAGTTTTTGCTTGTATTGAATCCGCAGCATGATGCCATAGTTCTTAACCAGCATTTCCATGCCAGGGATGCGTAAATAACAAAGCGCAGACCAGGGGGGCATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAaccacc                                                                                          <  2:844132/151‑1 (MQ=255)
                                                         gTTTTTGCTTGTATTGAATCCGCAGCATGATGCCATAGTTCTTAACCAGCATTTCCATGCCAGGGATGCGTAAATAACAAAGCGCAGACCAGGGGGGCATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGt                                                                                      <  1:723708/151‑1 (MQ=255)
                                                               gCTTGTATTGAATCCGCAGCATGATGCCATAGTTCTTAACCAGCATTTCCATGCCAGGGATGCGTAAATAACAAAGCGCAGACCAGGGGGGCATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACAt                                                                                >  1:520925/1‑151 (MQ=255)
                                                                           tCCGCAGCATGATGCCATAGTTCTTAACCAGCATTTCCATGCCAGGGATGCGTAAATAACAAAGCGCAGACCAGGGGGGCATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCt                                                                    >  2:117007/1‑151 (MQ=255)
                                                                             cGCAGCATGATGCCATAGTTCTTAACCAGCATTTCCATGCCAGGGATGCGTAAATAACAAAGCGCAGACCAGGGGGGCATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGa                                                                  <  1:722224/151‑1 (MQ=255)
                                                                              gcagcaTGATGCCATAGTTCTTAACCAGCATTTCCATGCCAGGGATGCGTAACTAACAAAGCGCAGACCAGGGGGGCATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGa                                                                  >  2:349581/1‑150 (MQ=255)
                                                                                agcaTGATGCCATAGTTCTTAACCAGCATTTCCATGCCAGGGATGCGTAAATAACAAAGCGCAGACCAGGGGGGCATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATg                                                               >  1:686102/1‑151 (MQ=255)
                                                                                agcaTGATGCCATAGTTCTTAACCAGCATTTCCATGCCAGGGATGCGTAAATAACAAAGCGCAGACCAGGGGGGCATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATg                                                               >  2:138386/1‑151 (MQ=255)
                                                                                   atgatgCCATAGTTCTTAACCAGCATTTCCATGCCAGGGATGCGTAAATAACAAAGCGCAGACCAGGGGGGCATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGt                                                            <  2:730999/151‑1 (MQ=255)
                                                                                     gatgCCATAGTTCTTAACCAGCATTTCCATGCCAGGGATGCGTAAATAACAAAGCGCAGACCAGGGGGGCATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTc                                                           <  2:165188/150‑1 (MQ=255)
                                                                                                tCTTAACCAGCATTTCCATGCCAGGGATGCGTAAATAACAAAGCGCAGACCAGGGGGTCATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGctc                                               >  1:508978/1‑151 (MQ=255)
                                                                                                  ttAACCAGCATTTCCATGCCAGGGATGCGTAAATAACAAAGCGCAGACCAGGGGGGCATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGctctc                                             >  2:487164/1‑151 (MQ=255)
                                                                                                   tAACCAGCATTTCCATGCCAGGGATGCGTAAATAACAAAGCGCAGACCAGGGGGGCATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGctctct                                            >  1:801695/1‑151 (MQ=255)
                                                                                                       cAGCATTTCCAGGCCAGGGATGCGTAAATAACAAAGCGCAGACCAGGGGGGCATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACg                                        >  1:359319/1‑151 (MQ=255)
                                                                                                         gCATTTCCATGCCAGGGATGCGTAAATAACAAAGCGCAGACCAGGGGGGCATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACggg                                      <  1:759045/151‑1 (MQ=255)
                                                                                                                cATGCCAGGGATGCGTAAATAACAAAGCGCAGACCAGGGGGGCATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAg                               <  1:761217/151‑1 (MQ=255)
                                                                                                                 aTGCCAGGGATGCGTAAATAACAAAGCGCAGACCAGGGGGGCATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGt                              >  1:594778/1‑151 (MQ=255)
                                                                                                                          aTGCGTAAATAACAAGGCGCAGACCAGGGGGGCATCACGCACCGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCata                     <  1:36096/151‑1 (MQ=255)
                                                                                                                          aTGCGTAAATAACAAAGCGCAGACCAGGGGGGCATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCata                     <  2:446068/151‑1 (MQ=255)
                                                                                                                          aTGCGTAAATAACAAAGCGCAGACCAGGGGGGCATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCat                      >  2:206959/1‑150 (MQ=255)
                                                                                                                           tGCGTAAATAACAAAGCGCAGACCAGGGGGGCATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCatat                    >  1:33104/1‑151 (MQ=255)
                                                                                                                            gCGTAAATAACAAAGCGCAGACCAGGGGGGCATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCAtata                   <  1:2461/151‑1 (MQ=255)
                                                                                                                              gTAAATAACAAAGCGCAGACCAGGGGGGCATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAAt                 <  2:50453/151‑1 (MQ=255)
                                                                                                                                    aaCAAAGCGCAGACCAGGGGGGCATCACGCAACGCGCGATTCAGGGCGCTAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACa            <  2:743546/150‑1 (MQ=255)
                                                                                                                                    aaCAAAGCGCAGACCAGGGGGGCATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGCACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACa            >  1:32611/1‑150 (MQ=255)
                                                                                                                                    aaCAAAGCGCAGACCAGGGGGGCATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCTGTACAGCCATATAATTGACa            <  2:262204/150‑1 (MQ=255)
                                                                                                                                       aaaGCGCAGACCAGGGGGGCATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATa         <  1:691136/150‑1 (MQ=255)
                                                                                                                                        aaGCGCAGACCAGGGGGGCATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATaa        >  1:53938/1‑150 (MQ=255)
                                                                                                                                          gcgcAGACCAGGGGGGCATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATaaaa      <  2:492909/150‑1 (MQ=255)
                                                                                                                                            gcAGACCAGGGGGGCATCACTCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATAAAagag   <  2:116716/151‑1 (MQ=255)
                                                                                                                                             cggACCAGGGGGGCATCGCGCAGCGCGCAATTCAGGGCGCGAACATTGGCCAGAGGAACAACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGCCTGCCGCACTCCGACGGGTTGCGAGTGCAGCCATATAATTGACATTAAAaga    <  1:832422/147‑1 (MQ=255)
                                                                                                                                             cAGACCAGGGGGGCATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATAAAag     >  2:654464/1‑148 (MQ=255)
                                                                                                                                              aGACCAGGGGGGCATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATAAAAgaga  >  1:709454/1‑150 (MQ=255)
                                                                                                                                              aGACCAGGGGGGCATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATAAAAgaga  >  2:675613/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
CCGAAAGCTGATAAACATCTTCACCCGGTTCCAGCAAGGGTGACAGCCAGTGAGAAAGTTTTTGCTTGTATTGAATCCGCAGCATGATGCCATAGTTCTTAACCAGCATTTCCATGCCAGGGATGCGTAAATAACAAAGCGCAGACCAGGGGGCATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATAAAAGAGA  >  NC_000913/2627735‑2628025

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