Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,856,887 C→T 70.1% intergenic (‑81/‑206) ygbA ← / → mutS uncharacterized protein/methyl‑directed mismatch repair protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,856,8870CT70.1% 101.7 / 55.4 83intergenic (‑81/‑206)ygbA/mutSuncharacterized protein/methyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base C (13/11);  new base T (23/34);  total (38/45)
Fisher's exact test for biased strand distribution p-value = 3.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACTCGCCGTTC  >  NC_000913/2856791‑2856982
                                                                                                |                                                                                               
tACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCaa                                                                                              >  2:907119/1‑100 (MQ=255)
  cGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAAcc                                                                                            <  2:2103709/100‑1 (MQ=255)
    cTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGa                                                                                          >  1:1466049/1‑100 (MQ=255)
       gCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATAc                                                                                       <  1:1349920/100‑1 (MQ=255)
            aCATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGTGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAAttt                                                                                  <  2:217277/100‑1 (MQ=255)
            aCATGCTGCCTACACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAAtt                                                                                   <  2:2072848/99‑1 (MQ=255)
             cATGCTGCCTCCACCTCATTAAGATGTATATATATTACATCTTAATCTTAAAAGGCACTATGACTCCAAAGAAGAAAGGTTAGTCAACCGATACAAtttt                                                                                 >  2:424100/1‑100 (MQ=255)
              aTGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTg                                                                                <  2:195522/100‑1 (MQ=255)
              aTGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTg                                                                                <  2:737368/100‑1 (MQ=255)
               tgctgcCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGc                                                                               >  1:870835/1‑100 (MQ=255)
                 ctgcCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGTCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGt                                                                             <  2:1291008/100‑1 (MQ=255)
                  tgcCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTa                                                                            <  1:2164778/100‑1 (MQ=255)
                    ccTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACt                                                                          <  2:758571/100‑1 (MQ=255)
                    ccTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACt                                                                          >  1:1609888/1‑100 (MQ=255)
                     cTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACtt                                                                         >  1:494480/1‑100 (MQ=255)
                       ccACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGc                                                                       <  1:159102/100‑1 (MQ=255)
                        cACCTCATTAAGATGTATTTATATTACATCTTAATCTTAACGGGCACTATGACTCTAAAGAAGAAGGGTTCGACAACCGCTATCATTTTGCGTACTTGCt                                                                      >  2:13245/1‑100 (MQ=255)
                         aCCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCtt                                                                     >  1:1462185/1‑100 (MQ=255)
                         aCCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCtt                                                                     <  2:2153297/100‑1 (MQ=255)
                          ccTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTc                                                                    >  2:723912/1‑100 (MQ=255)
                           cTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCtt                                                                     <  2:955939/98‑1 (MQ=255)
                            tCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGTAGTGTTATCCAAGCGATACAATTTTGCGTACTTGCTTCa                                                                   <  2:1333230/99‑1 (MQ=255)
                             cATTAAGATGTATTTATATTACATCTTAATCTTACAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTCCTTTCTTCATa                                                                 <  2:1176454/100‑1 (MQ=255)
                             cATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACAATGAATAAAAAAAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATa                                                                 <  2:435120/100‑1 (MQ=255)
                              aTTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATaa                                                                <  2:1884573/100‑1 (MQ=255)
                               ttAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAg                                                               >  2:2073028/1‑100 (MQ=255)
                               ttAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAg                                                               >  2:1867766/1‑100 (MQ=255)
                               ttAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGACAACCGATACAATTTTGCGTACTTGCTTCATAAg                                                               >  2:1731995/1‑100 (MQ=255)
                               ttAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATAACTCCAAAGAAGAAAGGTTAGCAAACCGATAAACTTTAGCGTACATACTTCATaa                                                                >  2:1359171/1‑99 (MQ=255)
                                tAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGc                                                              <  2:1622655/100‑1 (MQ=255)
                                   gATGTATTTATAGTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATc                                                           >  1:1473410/1‑100 (MQ=255)
                                    aTGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCa                                                          <  2:997590/100‑1 (MQ=255)
                                     tGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCAc                                                         >  2:1399419/1‑100 (MQ=255)
                                     tGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCAc                                                         >  1:529942/1‑100 (MQ=255)
                                       tATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGc                                                       >  1:1429736/1‑100 (MQ=255)
                                       tATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGc                                                       <  1:2188498/100‑1 (MQ=255)
                                        aTTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCa                                                      >  1:724008/1‑100 (MQ=255)
                                        aTTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGATTAACAAACAGATACAATTTTGCGTACTTGGTTTATAAGAATCACgaa                                                      >  2:2124249/1‑98 (MQ=255)
                                           tataTTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCaaa                                                    <  2:1857690/99‑1 (MQ=255)
                                            atatTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAg                                                  <  1:1540834/100‑1 (MQ=255)
                                             tatTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGc                                                 <  1:354590/100‑1 (MQ=255)
                                              atTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGATGGGTTAGCTAACCGATACAATTTTGTGTACTTGTTTCATAAGCATCACGCAAAAGCt                                                <  2:912274/100‑1 (MQ=255)
                                              atTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCt                                                >  2:2227493/1‑100 (MQ=255)
                                               ttACATCTTAATCTTAAATGACACTATGACTACAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCAAGCAAAAGCTg                                               >  2:948442/1‑100 (MQ=255)
                                               ttACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTg                                               >  2:448203/1‑100 (MQ=255)
                                                  cATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCa                                             >  1:1881193/1‑99 (MQ=255)
                                                   atcttaatcttaAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCaaa                                           >  1:6878/1‑100 (MQ=255)
                                                   atcttaatcttaAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCaaa                                           <  1:513560/100‑1 (MQ=255)
                                                    tcttaatcttaaAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCaaaa                                          >  1:1383967/1‑100 (MQ=255)
                                                         atcgtaaAGGGCACTATTACTCCAAAGAAGAGGGGTTAGTCAACCTATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCa                                     <  2:1943123/100‑1 (MQ=255)
                                                          tcttaaAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCAt                                    >  2:1554712/1‑100 (MQ=255)
                                                           cttaaAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATc                                   <  1:435758/100‑1 (MQ=255)
                                                           cttaaAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATc                                   <  1:478796/100‑1 (MQ=255)
                                                            ttaaAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCt                                  <  1:492051/100‑1 (MQ=255)
                                                            ttaaAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCa                                     >  1:775757/1‑97 (MQ=255)
                                                            ttaaAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCt                                  >  2:1728254/1‑100 (MQ=255)
                                                             taaAGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCtt                                 >  1:73243/1‑100 (MQ=255)
                                                                aGGGCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTc                               <  2:1866962/99‑1 (MQ=255)
                                                                 gggCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTccc                             <  2:1734137/100‑1 (MQ=255)
                                                                 gggCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTccc                             <  1:1338487/100‑1 (MQ=255)
                                                                  ggCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCg                            <  2:1609877/100‑1 (MQ=255)
                                                                  ggCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCg                            >  1:254115/1‑100 (MQ=255)
                                                                  ggCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCg                            <  1:809198/100‑1 (MQ=255)
                                                                  ggCACTATGACTCCAAAGAAGAAGGCTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCg                            <  2:475536/100‑1 (MQ=255)
                                                                   gCACTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCgg                           >  2:696980/1‑100 (MQ=255)
                                                                      cTATGACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAAc                        >  1:1895336/1‑100 (MQ=255)
                                                                      cTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAAc                        >  2:1606097/1‑100 (MQ=255)
                                                                          gACTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGc                    <  1:1099917/100‑1 (MQ=255)
                                                                           aCTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAATCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCa                   <  2:449517/100‑1 (MQ=255)
                                                                            cTCCAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTACTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCAt                  >  2:368911/1‑100 (MQ=255)
                                                                              ccAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTGCCGGAACCAGCATCa                >  2:1553307/1‑100 (MQ=255)
                                                                              ccAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCa                <  1:244928/100‑1 (MQ=255)
                                                                              ccAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACCTGCTTCATAAGCATCACGCAAAAGCTGCAACACAGCATCTTTCCCGGACCCAGCATc                 >  2:365981/1‑99 (MQ=255)
                                                                               cAAAGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCaa               <  2:671933/100‑1 (MQ=255)
                                                                               cAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCaa               >  1:1922173/1‑100 (MQ=255)
                                                                                aaaGAAGAAGGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAg              <  1:1473098/100‑1 (MQ=255)
                                                                                   gaagaagGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAAc           <  2:619753/100‑1 (MQ=255)
                                                                                      gaagGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACTCg        <  2:1895324/100‑1 (MQ=255)
                                                                                       aagGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCCTCAAGAACTCGc       <  1:1469913/100‑1 (MQ=255)
                                                                                       aagGGTTAGTCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACTCGc       <  2:1922507/100‑1 (MQ=255)
                                                                                         gggTTAGTCACCCGATACAATTTTGCGGACTTGCCGCATGAGCATCCCGCAAAAGCTGCAAGCCAGCAGGTTTCCCGGAACCAGCATCAAGAACTCGcc      <  2:795587/99‑1 (MQ=255)
                                                                                           gTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACTCGCCGtt   <  2:149041/100‑1 (MQ=255)
                                                                                             tAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACTCGCCGTTc  >  1:109095/1‑99 (MQ=255)
                                                                                                |                                                                                               
TACGCTTGCCGGACATGCTGCCTCCACCTCATTAAGATGTATTTATATTACATCTTAATCTTAAAGGGCACTATGACTCCAAAGAAGAAGGGTTAGCCAACCGATACAATTTTGCGTACTTGCTTCATAAGCATCACGCAAAAGCTGCAAAACAGCATCTTTCCCGGAACCAGCATCAAGAACTCGCCGTTC  >  NC_000913/2856791‑2856982

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