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

CCGAAAGCTGATAAACATCTTCACCCGGTTCCAGCAAGGGTGACAGCCAGTGAGAAAGTTTTTGCTTGTATTGAATCCGCAGCATGATGCCATAGTTCTTAACCAGCATTTCCATGCCAGGGATGCGTAAATAACAAAGCGCAGACCA‑GGGGG‑CATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATAAAAGAGA  >  NC_000913/2627735‑2628025
                                                                                                                                                    |                                                                                                                                                
ccGAAAGCTGATAAACATCTTCACCCGGTTCCAGCAAGGGTGACAGCCAGTGAGAAAGTTTTTGCTTGTATTGAATCCGCAGCATGATGCCATAGTTCTTAACCAGCATTTCCATGCCAGGGATGCGTAAATAACAAAGCGCAGACCA‑ggg                                                                                                                                               <  2:55243/151‑1 (MQ=255)
         gATAAACATCTTCACCCGGTTCCAGCAAGGGTGACAGCCAGTGAGAAAGTTTTTGCTTGTATTGAATCCGCAGCATGATGCCATAGTTCTTAACCAGCATTTCCATGCCAGGGATGCGTAAATAACAAAGCGCAGACCA‑ggggggcatc                                                                                                                                        <  2:63661/149‑6 (MQ=255)
               cATCTTCACCCGGTTCCAGCAAGGGTGACAGCCAGTGAGAAAGTTTTTGCTTGTATTGAATCCGCAGCATGATGCCATAGTTCTTAACCAGCATTTCCATGCCAGGGATGCGTAAATAACAAAGCGCAGACCAGGGGGG‑CATCACGCAAcg                                                                                                                                >  1:261387/1‑151 (MQ=255)
                aTCTTCACCCGGTTCCAGCAAGGGTGACAGCCAGTGAGAAAGTTTTTGCTTGTATTGAATCCGCAGCATGATGCCATAGTTCTTAACCAGCATTTCCATGCCAGGGATGCGTAAATAACAAAGCGCAGACCAGGGGGG‑CATCACGCAAcgc                                                                                                                               >  1:231095/1‑151 (MQ=255)
                   ttCACCCGGTTCCAGCAAGGGTGACAGCCAGTGAGAAAGTTTTTGCTTGTATTGAATCCGCAGCATGATGCCATAGTTCTTAACCAGCATTTCCATGCCAGGGATGCGTAAATAACAAAGCGCAGACCAGGGGGG‑CATCACGCAAcgcgcg                                                                                                                            >  1:293768/1‑151 (MQ=255)
                   ttCACCCGGTTCCAGCAAGGGTGACAGCCAGTGAGAAAGTTTTTGCTTGTATTGAATCCGCAGCATGATGCCATAGTTCTTAACCAGCATTTCCATGCCAGGGATGCGTAAATAACAAAGCGCAGACCAGGGGGG‑CATCACGCAAcgcgcg                                                                                                                            >  1:53366/1‑151 (MQ=255)
                     cACCCGGTTCCAGCAAGGGTGACAGCCAGTGAGAAAGTTTTTGCTTGTATTGAATCCGCAGCATGATGCCATAGTTCTTAACCAGCATTTCCATGCCAGGGATGCGTAAATAACAAAGCGCAGACCAGGGGGG‑CATCACGCAACGCGCGAt                                                                                                                          <  2:382576/151‑1 (MQ=255)
                      aCCCGGTTCCAGCAAGGGTGACAGCCAGTGAGAAAGTTTTTGCTTGTATTGAATCCGCAGCATGATGCCATAGTTCTTAACCAGCATTTCCATGCCAGGGATGCGTAAATAACAAAGCGCAGACCAGGGGGG‑CATCACGCAACGCGCGAtt                                                                                                                         >  1:844074/1‑151 (MQ=255)
                         cGGTTCCAGCAAGGGTGACAGCCAGTGAGAAAGTTTTTGCTTGTATTGAATCCGCAGCATGATGCCATAGTTCTTAACCAGCATTTCCATGCCAGGGATGCGTAAATAACAAAGCGCAGACCAGGGGGG‑CATCACGCAACGCGCGATTc                                                                                                                        >  2:267055/1‑149 (MQ=255)
                          ggTTCCAGCAAGGGTGACAGCCAGTGAGAAAGTTTTTGCTTGTATTGAATCCGCAGCATGATGCCATAGTTCTTAACCAGCATTTCCATGCCAGGGATGCGTAAATAACAAAGCGCAGACCAGGGGGG‑CATCACGCAACGCGCGATTCAgg                                                                                                                     >  2:407117/1‑151 (MQ=255)
                                  aaaGGGTGACAGTCAGTGAGAAAGTTTTTGCTTGTATTGAATCCGCAGCATGATGTCATAGTTCTTAACCAGCATTTCCATGCAAGGGATGCGTAAATTACAAAGCGCAGACCAGGGGGG‑CATCACGCAACGCGCGATTCAGGGCGCGa                                                                                                               >  1:534079/2‑149 (MQ=255)
                                       gTGACAGCCAGTGAGAAAGTTTTTGCTTGTATTGAATCCGCAGCATGATGCCATAGTTCTTAACCAGCATTTCCATGCCAGGGATGCGTAAATAACAAAGCGCAGACCAGGGGGG‑CATCACGCAACGCGCGATTCAGGGCGCGAACATTTg                                                                                                        >  2:446450/1‑151 (MQ=255)
                                          aCAGCCAGTGAGAAAGTTTTTGCTTGTATTGAATCCGCAGCATGATGCCATAGTTCTTAACCAGCATTTCCATGCCAGGGATGCGTAAATAACAAAGCGCAGACCAGGGGGG‑CATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCa                                                                                                     <  2:397254/151‑1 (MQ=255)
                                              ccAGTGAGAAAGTTTTTGCTTGTATTGAATCCGCAGCATGATGCCATAGTTCTTAACCAGCATTTCCATGCCAGGGATGCGTAAATAACAAAGCGCAGACCAGGGGGG‑CATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGAt                                                                                                  >  2:614963/1‑150 (MQ=255)
                                              ccAGTGAGAAAGTTTTTGCTTGTATTGAATCCGCAGCATGATGCCATAGTTCTTAACCAGCATTTCCATGCCAGGGATGCGTAAATAACAAAGCGCAGACCAGGGGGG‑CATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATg                                                                                                 >  1:820739/1‑151 (MQ=255)
                                                     gaAAGTTTTTGCTTGTATTGAATCCGCAGCATGATGCCATAGTTCTTAACCAGCATTTCCATGCCAGGGATGCGTAAATAACAAAGCGCAGACCAGGGGGG‑CATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAaccacc                                                                                          <  2:844132/151‑1 (MQ=255)
                                                         gTTTTTGCTTGTATTGAATCCGCAGCATGATGCCATAGTTCTTAACCAGCATTTCCATGCCAGGGATGCGTAAATAACAAAGCGCAGACCAGGGGGG‑CATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGt                                                                                      <  1:723708/151‑1 (MQ=255)
                                                               gCTTGTATTGAATCCGCAGCATGATGCCATAGTTCTTAACCAGCATTTCCATGCCAGGGATGCGTAAATAACAAAGCGCAGACCAGGGGGG‑CATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACAt                                                                                >  1:520925/1‑151 (MQ=255)
                                                                           tCCGCAGCATGATGCCATAGTTCTTAACCAGCATTTCCATGCCAGGGATGCGTAAATAACAAAGCGCAGACCAGGGGGG‑CATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCt                                                                    >  2:117007/1‑151 (MQ=255)
                                                                             cGCAGCATGATGCCATAGTTCTTAACCAGCATTTCCATGCCAGGGATGCGTAAATAACAAAGCGCAGACCAGGGGGG‑CATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGa                                                                  <  1:722224/151‑1 (MQ=255)
                                                                              gcagcaTGATGCCATAGTTCTTAACCAGCATTTCCATGCCAGGGATGCGTAACTAACAAAGCGCAGACCAGGGGGG‑CATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGa                                                                  >  2:349581/1‑150 (MQ=255)
                                                                                agcaTGATGCCATAGTTCTTAACCAGCATTTCCATGCCAGGGATGCGTAAATAACAAAGCGCAGACCAGGGGGG‑CATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATg                                                               >  1:686102/1‑151 (MQ=255)
                                                                                agcaTGATGCCATAGTTCTTAACCAGCATTTCCATGCCAGGGATGCGTAAATAACAAAGCGCAGACCAGGGGGG‑CATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATg                                                               >  2:138386/1‑151 (MQ=255)
                                                                                   atgatgCCATAGTTCTTAACCAGCATTTCCATGCCAGGGATGCGTAAATAACAAAGCGCAGACCAGGGGGG‑CATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGt                                                            <  2:730999/151‑1 (MQ=255)
                                                                                     gatgCCATAGTTCTTAACCAGCATTTCCATGCCAGGGATGCGTAAATAACAAAGCGCAGACCAGGGGGG‑CATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTc                                                           <  2:165188/150‑1 (MQ=255)
                                                                                                tCTTAACCAGCATTTCCATGCCAGGGATGCGTAAATAACAAAGCGCAGACCA‑GGGGGTCATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGctc                                               >  1:508978/1‑151 (MQ=255)
                                                                                                  ttAACCAGCATTTCCATGCCAGGGATGCGTAAATAACAAAGCGCAGACCAGGGGGG‑CATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGctctc                                             >  2:487164/1‑151 (MQ=255)
                                                                                                   tAACCAGCATTTCCATGCCAGGGATGCGTAAATAACAAAGCGCAGACCAGGGGGG‑CATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGctctct                                            >  1:801695/1‑151 (MQ=255)
                                                                                                       cAGCATTTCCAGGCCAGGGATGCGTAAATAACAAAGCGCAGACCAGGGGGG‑CATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACg                                        >  1:359319/1‑151 (MQ=255)
                                                                                                         gCATTTCCATGCCAGGGATGCGTAAATAACAAAGCGCAGACCAGGGGGG‑CATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACggg                                      <  1:759045/151‑1 (MQ=255)
                                                                                                                cATGCCAGGGATGCGTAAATAACAAAGCGCAGACCAGGGGGG‑CATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAg                               <  1:761217/151‑1 (MQ=255)
                                                                                                                 aTGCCAGGGATGCGTAAATAACAAAGCGCAGACCAGGGGGG‑CATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGt                              >  1:594778/1‑151 (MQ=255)
                                                                                                                          aTGCGTAAATAACAAGGCGCAGACCAGGGGGG‑CATCACGCACCGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCata                     <  1:36096/151‑1 (MQ=255)
                                                                                                                          aTGCGTAAATAACAAAGCGCAGACCAGGGGGG‑CATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCata                     <  2:446068/151‑1 (MQ=255)
                                                                                                                          aTGCGTAAATAACAAAGCGCAGACCAGGGGGG‑CATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCat                      >  2:206959/1‑150 (MQ=255)
                                                                                                                           tGCGTAAATAACAAAGCGCAGACCAGGGGGG‑CATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCatat                    >  1:33104/1‑151 (MQ=255)
                                                                                                                            gCGTAAATAACAAAGCGCAGACCAGGGGGG‑CATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCAtata                   <  1:2461/151‑1 (MQ=255)
                                                                                                                              gTAAATAACAAAGCGCAGACCAGGGGGG‑CATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAAt                 <  2:50453/151‑1 (MQ=255)
                                                                                                                                    aaCAAAGCGCAGACCAGGGGGG‑CATCACGCAACGCGCGATTCAGGGCGCTAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACa            <  2:743546/150‑1 (MQ=255)
                                                                                                                                    aaCAAAGCGCAGACCAGGGGGG‑CATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGCACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACa            >  1:32611/1‑150 (MQ=255)
                                                                                                                                    aaCAAAGCGCAGACCAGGGGGG‑CATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCTGTACAGCCATATAATTGACa            <  2:262204/150‑1 (MQ=255)
                                                                                                                                       aaaGCGCAGACCAGGGGGG‑CATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATa         <  1:691136/150‑1 (MQ=255)
                                                                                                                                        aaGCGCAGACCAGGGGGG‑CATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATaa        >  1:53938/1‑150 (MQ=255)
                                                                                                                                          gcgcAGACCAGGGGGG‑CATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATaaaa      <  2:492909/150‑1 (MQ=255)
                                                                                                                                            gcAGACCAGGGGGG‑CATCACTCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATAAAagag   <  2:116716/151‑1 (MQ=255)
                                                                                                                                             cggACCAGGGGGG‑CATCGCGCAGCGCGCAATTCAGGGCGCGAACATTGGCCAGAGGAACAACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGCCTGCCGCACTCCGACGGGTTGCGAGTGCAGCCATATAATTGACATTAAAaga    <  1:832422/147‑1 (MQ=255)
                                                                                                                                             cAGACCAGGGGGG‑CATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATAAAag     >  2:654464/1‑148 (MQ=255)
                                                                                                                                              aGACCAGGGGGG‑CATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATAAAAgaga  >  1:709454/1‑150 (MQ=255)
                                                                                                                                              aGACCAGGGGGG‑CATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATAAAAgaga  >  2:675613/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
CCGAAAGCTGATAAACATCTTCACCCGGTTCCAGCAAGGGTGACAGCCAGTGAGAAAGTTTTTGCTTGTATTGAATCCGCAGCATGATGCCATAGTTCTTAACCAGCATTTCCATGCCAGGGATGCGTAAATAACAAAGCGCAGACCA‑GGGGG‑CATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTCCACATACGCAAGCCGCTGAATGCGTCTGACTACCGCTCTCTGACGGGTTGCCAGTACAGCCATATAATTGACAATAAAAGAGA  >  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: ‑