Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,448,304 C→T intergenic (+69/‑143) ytfP → / → chpS gamma‑glutamylamine cyclotransferase family protein YtfP/ChpS antitoxin of the ChpB‑ChpS toxin‑antitoxin system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,448,3040CT97.0% 87.4 / ‑3.7 34intergenic (+69/‑143)ytfP/chpSgamma‑glutamylamine cyclotransferase family protein YtfP/ChpS antitoxin of the ChpB‑ChpS toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (16/16);  total (17/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.51e-01

TGACAACGCCACGCTGGCCGAACTTGATGCCTTGCGCACCAGGGGCGGTGAATACGCGCGCCAGTTGATTCAGACGCCGTACGGGAGTGCATGGATGTACGTTTATCAACGACCCGTCGATGGATTAAAGCTAATTGAAAGCGGCGACTGGTTAGACAGGGATAAGTAACCATATGCATACGCCACCTTCGGGTGGCGTTGTTTTTTGCGAGACGACTCGCATTCTGTTTTGTAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCG  >  NC_000913/4448067‑4448540
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                            
tGACAACGCCACGCTGGCCGAACTTGATGCCTTGCGCACAAGGGGAGGTGAATACGCGCGCCAGTTGATTCAGACGCCGTACGGGAGTGCATGGATGTACGTTTATCAACGACCCGTCGATGGATTAAAGCTAATTGAAAGCGGCGACTGGTTAGACAGGGATAAGTAACCATATGCATACGCCACCTTCGGGTGGCGTTGGTTTTTGCGAGACGACTCGCCTTTTGTTTTGTAATTTCCTCACCTTTTGc                                                                                                                                                                                                                                 >  1:181109/1‑251 (MQ=255)
  acaacgcccactcggCCGAACTTGATGACCTTCGCACCAGGGGCGGTGAATTCGCGCGCCATTTGATTCAGACGCCGTACGGGAGTTCATGGAGGTACGTTTATCAACGACCCGTCGATGGATTAACGCTAATTGAAATCTGCGACTGGTTAGACTGGGATAAGTAACCAGATGCATACGCCCCCTTCGGGTTGCGGTGTTCTTTGCGAGACTACTCGCATTCTGTGTTGTAATTTCCTCACCTTTTGCtt                                                                                                                                                                                                                               <  2:645805/238‑1 (MQ=255)
           acgctcgccaaactgtaggctttGCGTACCAGGGGCGGTGAATACGCGCGCCATTTGATTCAGACGCCGTACGGGAGTGCATGGATGTTCGTTGATCAACGACCCGTCGATGGATTAAAGCTAATTTAAAGCGGCGACTGGTTAGACAGGGATAAGTACCCATATGCATACGCCATCTTCGGGTTGCGTTGTTTTTTGCGAGACGACTCGCATTCTGTTTTGTAATTTCCTCACCTTTTGCTTTtctctc                                                                                                                                                                                                                        <  2:666687/229‑1 (MQ=255)
            gCTGGCCGAACTTGATGCCTTGCGCACCAGGGGCGGTGAATACGCGCGCCAGTTGATTCAGACGCCGTACGGGAGTGCATGGATGTACGTTTATCAACGACCCGTCGATGGATTAAAGCTAATTGAAAGCGGCGACTGGTTAGACAGGGATAAGTAACCATATGCATACGCCACCTTCGGGTGGCGTTGTTTTTTGCGAGACGACTCGCATTCTGTTTTGTAATTTCCTCACCTTTTGCTTTTCTCTCCgg                                                                                                                                                                                                                     >  1:72632/1‑250 (MQ=255)
                  cGAACTTGATGCCTTGCGCACCAGGGGCGGTGAATACGCGCGCCAGTTGATTCAGACGCCGTACGGGAGTGCATGGATGTACGTTTATCAACGACCCGTCGATGGATTAAAGCTAATTGAAAGCGGCGACTGGTTAGACAGGGATAAGTAACCATATGCATACGCCACCTTCGGGTGGCGTTGTTTTTTGCGAGACGAATCGCATTCTGTTTTGTAATTTCCTCACCTTTTGCTTTTCTCTCCGAGCCGCt                                                                                                                                                                                                               >  1:91531/1‑251 (MQ=255)
                    aaCTTGATGCCTTGCGCACCAGGGGCGGTGAATACGCGCGCCAGTTGATTCAGACGCCGTACGGGAGTGCATGGATGTACGTTTATCAACGACCCGTCGATGGATTAAAGCTAATTGAAAGCGGCGACTGGTTAGACAGGGATAAGTAACCATATGCATACGCCACCTTCGGGTGGCGTTGTTTTTTGCGAGACGACTCGCATTCTGTTTTGTAATTTCCTCACCTTTTGCTTTTCTCTCCGAGCCGCttt                                                                                                                                                                                                             <  1:227449/251‑1 (MQ=255)
                      cTTGATGCCTTGCGCACCAGGGGCGGTGAATACGCGCGCCAGTTGATTCAGACGCCGTACGGGAGTGCATGGATGTACGTTTATCAACGACCCGTCGATGGATTAAAGCTAATTGAAAGCGGCGACTGGTTAGACAGGGAAAAGTAACCATATGCATACGCCACCTTCGGGTGGCGTTGGTTTTTGCGAGACGACTCGCATTCTGTTTTGTAATTTCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTcc                                                                                                                                                                                                           >  1:144947/1‑251 (MQ=255)
                          gggCCTTGCGCACCAGGGGCGGTGAATACTCGCGCCAGTTGATTCAGACGCCGTACGGGAGTGCAGGGATGGTCGTTTATCAACGACCCGTCGATGGATTAAAGCTAATTGAAAGCGGCGACTGGTTAGACAGGGATAAGTAACCATATGCATACGCCACCTTCGGGTGGCGTTGTTTTTTGCGAGACGACTCGCATTCTGTTTTTTAATTTCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCatat                                                                                                                                                                                                       <  2:157794/249‑1 (MQ=255)
                                   cgcACCAGGGGCGGTGAATACGCGCGCCAGTTGATTCAGACGCCGTACGGGAGTGCATGGATGTACGTTTATCAACGACCCGTCGATGGATTAAAGCTAATTGAAAGCGGCGACTGGTTAGACAGGGATAAGTAACCATATGCATACGCCACCTTCGGGTGGCGTTGTTTTTTGCGAGACGACTCGCATTCTGTTTTGTAATTTCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAAc                                                                                                                                                                                               <  1:515341/251‑1 (MQ=255)
                                                            gcCAGTTGATTCAGACGCCGTACGGGAGTGCATGGATGTACGTTTATCAACGACCCGTCGATGGATTAAAGCTAATTGAAAGCGGCGACTGGTTAGACAGGGATAAGTAACCATATGCATACGCCACCTTCGGGTGGCGTTGTTTTTTGCGAGACGACTCGCATTCTGTTTTGTAATTTCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTc                                                                                                                                                                       >  1:137399/1‑250 (MQ=255)
                                                                                  cGGGAGTGCATGGATGTACGTTTATCAACGACCCGTCGATGGATTAAAGCTAATTGAAAGCGGCGACTGGTTAGACAGGGATAAGTAACCATATGCATACGCCACCTTCGGGTGGCGTTGTTTTTTGCGAGACGACTCGCATTCTGTTTTGTAATTTCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGt                                                                                                                                                >  2:116407/1‑251 (MQ=255)
                                                                                     gAGTGCATGGATGTACGTTTATCAAAGACCCGTCGATGGATTAAAGCTAATTGAAAGCGGCGACTGGTTAGACAGGGATAAGTAACCATATGCATACGCCACCTTCGGGTGGCGTTGTTTTTTGCGAGACGACTCGCATTCTGTTTTGTAATTTCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATCTCTATTAACGCATAAAAAACTCTGCTGGCATTAACAACTGCGCAGGGGTAAAACGTTTc                                                                                                                                             >  1:449801/1‑251 (MQ=255)
                                                                                         gCATGGATGTACGTTTATCAACGACCCGTCGATGGATTAAAGCTAATTGAAAGCGGCGACTGGTTAGACAGGGATAAGTAACCATATGCATACGCCACCTTCGGGTGGCGTTGTTTTTTGCGAGACGACTCGCATTCTGTTTTGTAATTTCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGt                                                                                                                                         <  1:574615/251‑1 (MQ=255)
                                                                                                 gTACGTTTATCAACGACCCGTCGATGGATTAAAGCTAATTGAAAGCGGCGACTGGTTAGACAGGGATAAGTAACCATATGCATACGCCACCTTCGGGTGGCGTTGTTTTTTGCGAGACGACTCGCATTCTGTTTTGTAATTTCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGt                                                                                                                                 >  1:156208/1‑251 (MQ=255)
                                                                                                    cGTTTATCAACGACCCGTCGATGGATTAAAGCTAATTGAAAGCGGCGACTGGTTAGACAGGGATAAGTAACCATATGCATACGCCACCTTCGGGTGGCGTTGTTTTTTGCGAGACGACTCGCATTCTGTTTTGTAATTTCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAg                                                                                                                              <  1:249388/251‑1 (MQ=255)
                                                                                                         aTCAACTACCCGTCGATGTATTAATGCTAATTGAAAGCGGCGACTGGTTAGACAGGGATAAGTAACCATATGCATACGCCACCTTCGGGTGGCGTTGTTTTTTGCGAGACGACTCGCATTCTGTTTTGTAATTTCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTTTAGCACCGTGAGTtata                                                                                                                         <  1:196910/251‑1 (MQ=255)
                                                                                                         aTCAACGACCCGTCGATGGATTAAAGCTAATTGAAAGCGGCGACTGGTTAGACAGGGATAAGTAACCATATGCATACGCCACCTTCGGGTGGCGTTGTTTTTTGCGAGACGACTCGCATTCTGTTTTGTAATTTCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTtat                                                                                                                          >  1:535545/1‑250 (MQ=255)
                                                                                                           cAACGACCCGTCGATGGATTAAAGCTAAGTGAAAGCGGCGACTGGTTAGACAGGGATAAGTAACCATATGCATACGCCACCTTCGGGTGGCGTTGTTTTTTGCGAGACGACTCGCATTCTGTTTTGTAATTTCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACt                                                                                                                       <  1:496933/251‑1 (MQ=255)
                                                                                                                aCCCGTCGATGGATTAAAGCTAATTGAAAGCGGCGACTGGTTAGACAGGGATAAGTAACCATATGCATACGCCACCTTCGGGTGGCGTTGTTTTTTGCGAGACGACTCGCATTCTGTTTTGTAATTTCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAAc                                                                                                                                                  >  1:661285/1‑219 (MQ=255)
                                                                                                                      agATGGATAAAAGCTAATTGAAAGCGGCGACTGGTTAGACAGGGATAAGTAGCCATATGCATACGCCCTCTTCGGGTGGCGTTGTTTTTTGCGAGACGACTCGAATTCTGTTTTGTAATTTCCTCCCCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACtt                                                                                                            <  2:137399/250‑1 (MQ=255)
                                                                                                                               aaaGCTAATTGAAAGCGGCGACTGGTTAGACAGGGATAAGTAACCATATGCATACGCCACCTTCGGGTGGCGTTGTTTTTTGCGAGACGACTCGCATTCTGTTTTGTAATTTCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTg                                                                                                   >  2:647725/1‑251 (MQ=255)
                                                                                                                                                    cTGGTTAGACAGGGATAAGTAACCATATGCATACGCCACCTTCGGGTGGCGTTGTTTTTTGCGAGACGACTCGCATTCTGTTTTGTAATTTCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATaaa                                                                               >  1:639899/1‑250 (MQ=255)
                                                                                                                                                           gACAGGGATAAGTTACCATATGCATCCGCCGCCTTCGGGTGGCGTTGTTTTTTGCGAGACGACTCGCATTCTGTTTTGTAATTTCCTCTCCTTTTGCTTTTCTCGCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATggg                                                                       <  1:306199/251‑1 (MQ=255)
                                                                                                                                                                   tAAGTAACCATATGCATACGCCACCTTCGGGTGGCGTTGTTTTTTGCGAGACGACTCGCATTCTGTTTTGTAATTTCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTACAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTCTTACCATAAAAAGATGGGGGAACAGt                                                               >  2:646633/1‑251 (MQ=255)
                                                                                                                                                                     aGTAACCATATGCATACGCCACCTTCGGGTGGCGTTGTTTTTTGCGAGACGACCCGCATTCTGTTTTGTAATTTCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGc                                                             <  1:201358/251‑1 (MQ=255)
                                                                                                                                                                              atGCATACGCCACCTTCGGGTGGCGTTGTTTTTTGCGAGACGACTCGCATTCTGTTTTGTAATTTCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTACGGAGGAGCAGATGCGTATTCCCATAAAAAGATGGGGGAAAAGTGAAGGAATGGt                                                    >  2:268420/1‑251 (MQ=255)
                                                                                                                                                                                  aTACGCCACCTGCGGGTGGCGTTGTTTGTTGCGAGACGACTCGCATTCTTTTTTGTAATTTCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAATAAACTCTGCTGGCATTCACAAATGCGCAGGGGGAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCAtt                                                <  1:356241/251‑1 (MQ=255)
                                                                                                                                                                                    aCGCCACCTTCGGGTGGCGTTGTTTTTTGCGAGACGACTCGCATTCTGTTTTGTAATTCCCTAACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTcc                                              <  1:172222/251‑1 (MQ=255)
                                                                                                                                                                                          ccTTCGGATGGCGTTGTTTTTTGCGAGACGACTCGCATTCTGTTTTGTAATTTCCTCACCTTTTGCTTTTATATCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGAATTCACAAATGCGCAGGAGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTAGTATAACTTAAGGAAGTGCAGATCAGAATTACCATAAAAAGATGAGGAAACAGTGCAGGTATggaactttaaaatg                                         >  2:377167/1‑238 (MQ=255)
                                                                                                                                                                                               gggTGGCGTTGTTTTTTGCGAGACGACTCGCATTTTGTTTTGTAATTTCCTCCCCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACTTTTCCTGTAGCACCGTGAGTTATACCTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTaa                                   <  1:455631/251‑1 (MQ=255)
                                                                                                                                                                                                tgTGGCGTTGTTTTTTGCGACACGACTCGCATTCTGTTTTGTAGTTACCTCACCTTTTGCTTTTCTATCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACACATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGATATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAAAAGTGCAGTGATGGACATTCCCAATATCGTAAt                                  >  2:118369/2‑251 (MQ=255)
                                                                                                                                                                                                           ttttttcgattcgacTCGCATTCTGTTTTGTAATTTCCTCACCGTGTGCCTTTCTCTACGAGCCGCTTTCCCTATCTATTAACGCATAAAAAACTCTTCTGGCAGTCACAAATGCGCAGGGGTAAAACGTTTCCTGGAGCACCGTGAGTTATAATTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTACAAAGATTGTAATGAAAGaactt                        <  2:49467/239‑1 (MQ=255)
                                                                                                                                                                                                                                tCTGTTTTGTAATTTCCTCGCCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCGATTAACGCATAAAAAACGCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGc   <  1:8773/250‑1 (MQ=255)
                                                                                                                                                                                                                                tCTGTTTTGTAATTTCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCATGAGTTATACTTTGTATGACTTAAGGAGGTGCAGATGCGGATTACCATAAAAAGATGGGAGAACAGAGCAGGTATGGTCATTCCCAATCTCGTAATGAAAGCACTTAACTTACAGACGGGGCagaggg  >  2:94622/1‑249 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                            
TGACAACGCCACGCTGGCCGAACTTGATGCCTTGCGCACCAGGGGCGGTGAATACGCGCGCCAGTTGATTCAGACGCCGTACGGGAGTGCATGGATGTACGTTTATCAACGACCCGTCGATGGATTAAAGCTAATTGAAAGCGGCGACTGGTTAGACAGGGATAAGTAACCATATGCATACGCCACCTTCGGGTGGCGTTGTTTTTTGCGAGACGACTCGCATTCTGTTTTGTAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCG  >  NC_000913/4448067‑4448540

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