Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,611,964 T→C 19.7% intergenic (‑110/+2) yneF ← / ← yneG putative membrane‑bound diguanylate cyclase/DUF4186 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,611,9640TC19.7% 107.2 / 24.0 56intergenic (‑110/+2)yneF/yneGputative membrane‑bound diguanylate cyclase/DUF4186 family protein
Reads supporting (aligned to +/- strand):  ref base T (23/22);  new base C (6/5);  total (29/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 = 7.51e-01

CCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTAT  >  NC_000913/1611866‑1612063
                                                                                                  |                                                                                                   
ccAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCACATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATtc                                                                                                    >  1:1402537/1‑100 (MQ=255)
ccAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAAtt                                                                                                     >  2:1137533/1‑99 (MQ=255)
 cAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATtct                                                                                                   >  2:2297680/1‑100 (MQ=255)
        tCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCt                                                                                            >  2:1226539/1‑100 (MQ=255)
         cAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTg                                                                                           >  2:220878/1‑100 (MQ=255)
         cAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATCCTCATGGCTg                                                                                           >  2:137338/1‑100 (MQ=255)
         cAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATCCTCATGGCTg                                                                                           >  2:2244307/1‑100 (MQ=255)
          aaTACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGa                                                                                          >  1:2124454/1‑100 (MQ=255)
            tACTGACCGAACCAATATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGAtt                                                                                        <  1:2364468/100‑1 (MQ=255)
             aCTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATCCTCATGGCTGATTc                                                                                       >  2:1104488/1‑100 (MQ=255)
              cTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATCCTCATGGCTGATTCa                                                                                      >  1:840745/1‑100 (MQ=255)
                gACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCCTGGCTGATTCATc                                                                                    >  1:2114362/1‑100 (MQ=255)
                gACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTGGTCGGTTTCCGCAATTCTCATGGCTGCTTCATc                                                                                    >  2:594641/1‑100 (MQ=255)
                  ccGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTg                                                                                  <  2:2124430/100‑1 (MQ=255)
                   cGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGt                                                                                 <  2:663479/100‑1 (MQ=255)
                     aaCCAAGATAGCATTTTCGCTGACAGCGGCGGCAAATAAATGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTAc                                                                               <  1:1094687/100‑1 (MQ=255)
                       ccAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACaa                                                                             >  1:732508/1‑100 (MQ=255)
                        cAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAAc                                                                            <  1:64907/100‑1 (MQ=255)
                        cAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATCCTCATGGCTGATTCATCTGTACAAc                                                                            <  1:237523/100‑1 (MQ=255)
                         aaGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACACTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAAcc                                                                           >  2:1304824/1‑100 (MQ=255)
                             tAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAAccagcc                                                                       >  2:265308/1‑100 (MQ=255)
                             tAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAAccagcc                                                                       <  2:902942/100‑1 (MQ=255)
                             tAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAAccagcc                                                                       <  2:1608165/100‑1 (MQ=255)
                               gCATGTTCGCTGACAGCGGCGGCAAATAAAGGAACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAAccagccag                                                                     >  2:472230/1‑100 (MQ=255)
                                  tGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGa                                                                  <  1:235171/100‑1 (MQ=255)
                                  tGTTCGCTGACAGCAACAGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGa                                                                  >  2:722759/1‑100 (MQ=255)
                                     tCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCCATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAg                                                               <  2:793664/100‑1 (MQ=255)
                                       gCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATCCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGAt                                                             <  2:1311734/100‑1 (MQ=255)
                                          gACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTAc                                                          >  2:1678389/1‑100 (MQ=255)
                                            cAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCg                                                        >  2:946358/1‑100 (MQ=255)
                                               cggcggCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCa                                                     >  2:2077213/1‑100 (MQ=255)
                                                ggcggcAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCAc                                                    <  1:1152840/100‑1 (MQ=255)
                                                ggcggcAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCAc                                                    <  1:1269304/100‑1 (MQ=255)
                                                 gcggcAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATCCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACg                                                   >  1:306651/1‑100 (MQ=255)
                                                    gcAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATCCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGata                                                <  1:1245296/100‑1 (MQ=255)
                                                        aTAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAAc                                            <  1:1631041/100‑1 (MQ=255)
                                                           aaGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACgtt                                         >  1:2243051/1‑100 (MQ=255)
                                                               tACCTTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACgttgt                                       <  1:248662/98‑1 (MQ=255)
                                                                aCCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACgttgttgt                                    <  1:261653/100‑1 (MQ=255)
                                                                     gggAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGttcttc                               <  2:1414366/100‑1 (MQ=255)
                                                                        aaTGCTATTTTGTCGTGTTCCGCAATTCTCATGTCTGATTCAGCTGTACAACCAGCCAGTGATATATTACCGCCACGATATAACGTTGTTGTTCTTCGCt                            <  2:539340/100‑1 (MQ=255)
                                                                          tGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCt                            <  1:1403385/98‑1 (MQ=255)
                                                                           gCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCTTCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTaa                         <  2:1301500/100‑1 (MQ=255)
                                                                           gCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATACCGTTGTTGTTCTTCGCTTaa                         >  1:1608904/1‑100 (MQ=255)
                                                                           gCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTaa                         >  2:1721278/1‑100 (MQ=255)
                                                                                tttGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGaaa                    >  2:692527/1‑100 (MQ=255)
                                                                                  tGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAaca                  >  1:1171686/1‑100 (MQ=255)
                                                                                  tGTCGTTTTCCGCAATCCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAaca                  >  1:2133455/1‑100 (MQ=255)
                                                                                   gTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAacac                 >  1:701010/1‑100 (MQ=255)
                                                                                    tCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACAcc                >  1:1701346/1‑100 (MQ=255)
                                                                                      gTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATATATTACCGCCACGATATAACGTAGTTGGTCTTCGCTTAACGATacac                 >  1:350785/1‑97 (MQ=255)
                                                                                      gTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCt              >  2:349619/1‑100 (MQ=255)
                                                                                       ttttCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTg             <  1:2077235/100‑1 (MQ=255)
                                                                                          tCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCgg          <  1:1226552/100‑1 (MQ=255)
                                                                                             gCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAAt       >  2:1828900/1‑100 (MQ=255)
                                                                                              cAATCCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATg      <  1:2244333/100‑1 (MQ=255)
                                                                                               aaTTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGt     <  1:1398809/100‑1 (MQ=255)
                                                                                               aaTCCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATg      <  1:137341/99‑1 (MQ=255)
                                                                                                aTTCTCATGGCTGATTCATCTGTACAACAAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGt     <  2:1380484/99‑1 (MQ=255)
                                                                                                  tctcTTGGATGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTAt  <  1:946369/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
CCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTAT  >  NC_000913/1611866‑1612063

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