Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,276,134 C→T 6.9% R222C (CGC→TGC)  yjcC → putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,276,1340CT 6.9% 156.1 / 2.6 58R222C (CGC→TGC) yjcCputative membrane‑anchored cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base C (39/15);  new base T (2/2);  total (41/17)
Fisher's exact test for biased strand distribution p-value = 5.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

AATGGGGTGTGTATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCT  >  NC_000913/4275988‑4276280
                                                                                                                                                  |                                                                                                                                                  
aaTGGGGTGTGTATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGCCAAAATGc                                                                                                                                                  >  1:799304/1‑149 (MQ=255)
aaTGGGGTGTGTATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAATGc                                                                                                                                                  >  1:358264/1‑149 (MQ=255)
      gtgtgtATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAAtt                                                                                                                                            >  1:116193/1‑149 (MQ=255)
                 aaggTGACGAAAACCTTTTTCACGTAAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCAATTGAAGGATTTAACCGTACAAACAAATGGCTATTAATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCAatgt                                                                                                                                 <  2:739008/147‑1 (MQ=255)
                   ggTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGTATTAAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCAttgttg                                                                                                                               <  1:627049/149‑1 (MQ=255)
                    gTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGc                                                                                                                              <  2:570263/149‑1 (MQ=255)
                          aaaaCCTTTTTCTCGTTAATCAAATAGGCCTCGGAAGCAACGTTTTCTCCGCTGATTCATTTGATGGAGTTAACCGTACAATGAAATGTCTTTTTATATGCGACAGTTTATTCGACAAAACGCCCAATTTCAGCCATTGTTGCGTCTTc                                                                                                                        >  2:510456/1‑149 (MQ=255)
                                    tctcGTTAAGCAAAGAGGCCTCGGCAGCAACTTCTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGt                                                                                                              <  2:604428/149‑1 (MQ=255)
                                               aaaGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCa                                                                                                   >  2:483256/1‑149 (MQ=255)
                                                  tagGCCTCGGCAGCAACGTTTTCTCCGCTTATTCATTTGAAGGATTTAACCGTACAAAGAAATCGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATGGTTGCGACTTCATATCAACGTCTTATAACCCAttt                                                                                                <  2:714495/148‑1 (MQ=255)
                                                     gCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTta                                                                                             >  2:529730/1‑149 (MQ=255)
                                                      ccTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTtat                                                                                            >  1:499683/1‑149 (MQ=255)
                                                      ccTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTtaa                                                                                            >  2:444270/1‑148 (MQ=255)
                                                                aaCGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTa                                                                                  >  2:421646/1‑149 (MQ=255)
                                                                       tctcCGCTGATTCATTTTAAGTATTTAAGCGTACAAAGAAATGGCTTTTTATATGAGAAAGTTTATTCGACAAAACGCCCTATTGCAGCAATTGTTGCGACTTCATATCAACGTCTTATAACCCATTATTATAATCATCTTATTTTTGc                                                                           <  1:853312/149‑1 (MQ=255)
                                                                        ctcCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCg                                                                          <  2:992442/149‑1 (MQ=255)
                                                                        ctcCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGATGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCg                                                                          >  2:425288/1‑149 (MQ=255)
                                                                         tcCGCTGATTCATTTGAAGGATTTAACCGTACTAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGt                                                                         >  1:949538/1‑149 (MQ=255)
                                                                                aTTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTTTATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAAACCATTTTTATAATCTACTTATTTTTGCGTTGCccga                                                                  >  2:1002070/1‑148 (MQ=255)
                                                                                aTTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAATCGCCCTATTGCAGCCATTGTTGCGAGTTCTTATCAACGGCATATACCCCATATTTAGAACAATCATATTTTTGCGTTCCccgc                                                                  >  1:448835/1‑149 (MQ=255)
                                                                                 ttCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCccgcc                                                                 >  1:650243/1‑149 (MQ=255)
                                                                                  tCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCccgccg                                                                >  1:448264/1‑149 (MQ=255)
                                                                                        gAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTAttt                                                          >  2:382227/1‑149 (MQ=255)
                                                                                         aaGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAGCCCAGTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTAtttt                                                         >  2:505449/1‑149 (MQ=255)
                                                                                            gtattAACCGTACAAAGAAATGGCTATTTATAAGCGACATTTTATTCGACAAAACGCCCAATTTCACCCATTGTTGCGTCTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTggg                                                      <  1:444270/146‑1 (MQ=255)
                                                                                             aTTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTgggg                                                     >  2:753091/1‑149 (MQ=255)
                                                                                              tttAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATCACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGa                                                    >  1:775533/1‑149 (MQ=255)
                                                                                               ttAAACGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAg                                                   >  2:539490/1‑149 (MQ=255)
                                                                                                 aaCCGTACAAAGAAATTTCTATTTATATGCGACAGTTTATTCGACAAAATGCCGAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATATTTGCGTTGCCCGCCGGTATTTTGGGGAggc                                                 <  2:799304/149‑3 (MQ=255)
                                                                                                 aaCCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAATGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTc                                                 <  2:358264/149‑1 (MQ=255)
                                                                                                 aaCCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATGTTATTTTTGCGTTGCCCGCGGGTATTTTGGGGAGtg                                                 >  2:164604/1‑148 (MQ=255)
                                                                                                 aaCCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACACATTTTTATAATCATCTTATTTTTGCTTTGCCCGCCGGTATTTTGGGGAGTc                                                 >  1:34306/1‑149 (MQ=255)
                                                                                                 aaCCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCACAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTCATAACCCATGTTTATAATCATCTTTTTTTTGCGTTGCACGCCGGGAGTTTGGGGAggc                                                 >  2:656433/1‑147 (MQ=255)
                                                                                                 aaCCCTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTc                                                 <  1:51024/149‑1 (MQ=255)
                                                                                                    cGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTg                                              >  1:174851/1‑149 (MQ=255)
                                                                                                     gTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGt                                             >  1:725058/1‑149 (MQ=255)
                                                                                                         aaaGAAATGGCTATTTATATGCGCCAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTAGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTAGTTTAGCGTTGCACGCCGGTATTTTGGGGAGTCTTGTTCTg                                         >  2:37114/1‑149 (MQ=255)
                                                                                                         aaaGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTg                                         >  1:278298/1‑149 (MQ=255)
                                                                                                            gAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTCTATCCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCtt                                      <  2:266677/149‑2 (MQ=255)
                                                                                                             aaaTGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTAt                                     >  1:289465/1‑149 (MQ=255)
                                                                                                                  gCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCGATTTTTATAATGATCTTATTTTTGCGTTGCCCGCCGTTATTTTGGGGAGTCTTGTTATGCTATTaata                                >  2:586240/1‑146 (MQ=255)
                                                                                                                    tATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTTCGTTGCCAGCCGTTATTTTGGGGAGTCTTGTTCTTCTATTACTCTg                              >  1:578068/1‑149 (MQ=255)
                                                                                                                         atatGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGGAATGTGGGGAGTCTTGTTCTGCTATTACTCTGGCTAc                         <  1:347525/149‑1 (MQ=255)
                                                                                                                          tatGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACg                        >  2:625901/1‑149 (MQ=255)
                                                                                                                          tatGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGAGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTAATTTTGCGTTGCCAGCCGGGATTTTGGGGAGTCTTGGTCTGCTATTACTCTGGCTACg                        >  2:734540/1‑149 (MQ=255)
                                                                                                                             gCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTAt                     >  1:293093/1‑149 (MQ=255)
                                                                                                                               gACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTc                   >  1:1026079/1‑149 (MQ=255)
                                                                                                                                aCAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCg                  >  2:288247/1‑149 (MQ=255)
                                                                                                                                    tttATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACaa              >  1:466219/1‑149 (MQ=255)
                                                                                                                                      tATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACaaaa            <  1:529730/149‑1 (MQ=255)
                                                                                                                                       aatCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGGCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAAc           <  2:116193/147‑1 (MQ=255)
                                                                                                                                        ttCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTAGTTCTCTGGCTACTGATTCGACACAACt          >  2:958577/1‑149 (MQ=255)
                                                                                                                                           gACAAACCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGGTATTACTCTGGCTACGTATTCGACAAAACTAtt       <  2:949538/149‑1 (MQ=255)
                                                                                                                                           gACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACTAAACTAtt       <  2:650243/149‑1 (MQ=255)
                                                                                                                                            aCAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTAttt      >  2:854228/1‑149 (MQ=255)
                                                                                                                                              aaaaCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTAt    >  2:858446/1‑149 (MQ=255)
                                                                                                                                              aaaaCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTAGTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACGATTTAt    >  2:743157/1‑149 (MQ=255)
                                                                                                                                                aaCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTAtct  >  1:316461/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
AATGGGGTGTGTATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCT  >  NC_000913/4275988‑4276280

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