Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,792,240 Δ1 bp 17.1% intergenic (‑220/+27) ydiV ← / ← nlpC anti‑FlhD4C2 factor, inactive EAL family phosphodiesterase/putative C40 clan peptidase lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,792,2400T.17.1% 230.3 / 44.4 82intergenic (‑220/+27)ydiV/nlpCanti‑FlhD4C2 factor, inactive EAL family phosphodiesterase/putative C40 clan peptidase lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (30/37);  new base . (8/6);  total (38/44)
Fisher's exact test for biased strand distribution p-value = 5.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.25e-01

TCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTG  >  NC_000913/1792142‑1792339
                                                                                                  |                                                                                                    
tCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTa                                                                                                     <  1:1759806/100‑1 (MQ=255)
     tttttAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTatcatc                                                                                               >  2:1727789/1‑100 (MQ=255)
     tttttAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTatcatc                                                                                               >  2:548826/1‑100 (MQ=255)
      ttttAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTGTGTATGCAAAATACTACGGGAGTGAGGTGTATATGGACGGCGATGTATCATCg                                                                                              <  2:1365356/100‑1 (MQ=255)
       tttAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGc                                                                                             >  1:264292/1‑100 (MQ=255)
       tttAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGc                                                                                             >  1:1695813/1‑100 (MQ=255)
        ttAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGcc                                                                                            <  1:1074427/100‑1 (MQ=255)
        ttAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGcc                                                                                            >  1:2255049/1‑100 (MQ=255)
         tAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCg                                                                                           <  2:1932698/100‑1 (MQ=255)
         tAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCg                                                                                           >  2:502121/1‑100 (MQ=255)
          aGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGATACATCGCCGt                                                                                          <  2:330929/100‑1 (MQ=39)
          aGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGATACATCGCCGt                                                                                          <  1:1757939/100‑1 (MQ=39)
           gATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGATACATCGCCGTc                                                                                         <  1:31147/100‑1 (MQ=255)
           gATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGATACATCGCCGTc                                                                                         <  1:1585326/100‑1 (MQ=39)
            aTTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGATACATCGCCGTcc                                                                                        <  2:1863235/100‑1 (MQ=39)
            aTTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGATACATCGCCGTcc                                                                                        <  1:1101720/100‑1 (MQ=39)
             ttATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGATACATCACCGTCCa                                                                                       <  1:1715673/100‑1 (MQ=38)
              tATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCAt                                                                                      <  2:1408100/100‑1 (MQ=255)
              tATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCAt                                                                                      >  2:1724742/1‑100 (MQ=255)
              tATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCAt                                                                                      <  1:295541/100‑1 (MQ=255)
                  cTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCAttttt                                                                                  <  2:1571455/100‑1 (MQ=255)
                  cTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCAttttt                                                                                  >  2:2165128/1‑100 (MQ=255)
                   tGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTc                                                                                 <  1:1913327/100‑1 (MQ=255)
                      ttATAAACGAATAGTCCATTTACGGAATTTTTGTCTGCAAAAGACTACTGTATTCAGTAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGAc                                                                              <  2:2294958/100‑1 (MQ=255)
                        ataAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACAt                                                                            >  1:1732980/1‑100 (MQ=255)
                         taAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATa                                                                           >  2:2379422/1‑100 (MQ=255)
                          aaaCGAATAATCCATTTACGGACTTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGGCGGCGATGTATCATCCCCGTCCATTTTTCGACATAc                                                                          >  2:617149/1‑100 (MQ=255)
                          aaaCGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATAc                                                                          <  2:2134762/100‑1 (MQ=255)
                           aaCGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACt                                                                         >  1:1821036/1‑100 (MQ=255)
                           aaCGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGAAGATGTATTATCGACGTCCATTTTTCGACATACt                                                                         >  2:2181058/1‑100 (MQ=255)
                            aCGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTa                                                                        <  2:2319108/100‑1 (MQ=255)
                               aataatCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGAt                                                                     <  1:25534/100‑1 (MQ=255)
                                      cATTTACGGAATTTTTGTTTGTAAAATACTACTGTATTCAGGAGTAAATGTACGTCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGtctc                                                              >  2:2280470/1‑100 (MQ=255)
                                         ttACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGcc                                                           <  2:2001188/100‑1 (MQ=255)
                                          tACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCt                                                          >  2:112233/1‑100 (MQ=255)
                                             ggAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGcc                                                       <  2:1353109/100‑1 (MQ=255)
                                                    ttGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAAt                                                 <  2:1257640/99‑1 (MQ=255)
                                                      gTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAAtttt                                              <  1:1817040/100‑1 (MQ=255)
                                                       tCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAAttttt                                             <  2:1236438/100‑1 (MQ=255)
                                                        cTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTc                                            <  2:1360811/100‑1 (MQ=255)
                                                          gCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTg                                          <  2:658675/100‑1 (MQ=255)
                                                          gCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTg                                          <  1:2269723/100‑1 (MQ=255)
                                                               aTACTCCTGTATTCAGGAGTAAATGGCAGGCTATGTATCATCGTCGTCCATTTTTCGAGCTACGACATTCGTCCCGCTTGCCAGAAAATTTTCTGCCa                                       <  2:2372344/98‑1 (MQ=255)
                                                                  ctactGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGa                                  >  1:2074422/1‑100 (MQ=255)
                                                                    actGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACa                                >  2:399689/1‑100 (MQ=255)
                                                                    actGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACa                                >  2:1221818/1‑100 (MQ=255)
                                                                     ctGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACAt                               >  1:382813/1‑100 (MQ=255)
                                                                     ctGTATCCAGGAGTAAATGGATGGCGCTGTATCAGCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCGGAAATTTTTCTGCCAATAg                                   <  2:1463531/96‑1 (MQ=255)
                                                                      tGTATTCAGGAGTAAATGGACAGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACAtt                              <  2:1719990/100‑1 (MQ=255)
                                                                       gTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTa                             >  1:900602/1‑100 (MQ=255)
                                                                        tATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTAt                            <  2:13334/100‑1 (MQ=255)
                                                                        tATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTAt                            <  1:502123/100‑1 (MQ=255)
                                                                         aTTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATc                           <  1:2036171/100‑1 (MQ=255)
                                                                         aTTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATc                           >  1:447685/1‑100 (MQ=255)
                                                                          ttCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCa                          >  2:354302/1‑100 (MQ=255)
                                                                          ttCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCa                          >  1:339211/1‑100 (MQ=255)
                                                                            cAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAg                        >  1:810957/1‑100 (MQ=255)
                                                                            cAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAg                        >  2:1710112/1‑100 (MQ=255)
                                                                             aGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGt                       <  1:2058393/100‑1 (MQ=255)
                                                                              ggAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTg                      <  1:1033928/100‑1 (MQ=255)
                                                                              ggAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTg                      <  1:730974/100‑1 (MQ=255)
                                                                                aGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAg                    >  2:1397472/1‑100 (MQ=255)
                                                                                 gTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAgg                   <  2:1085900/100‑1 (MQ=255)
                                                                                  tAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGa                  <  2:311976/100‑1 (MQ=255)
                                                                                  tAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGa                  <  2:197534/100‑1 (MQ=255)
                                                                                    aaTGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAAc                <  1:1068702/100‑1 (MQ=255)
                                                                                    aaTGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAAc                >  2:1218487/1‑100 (MQ=255)
                                                                                      tGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGc              <  2:852339/100‑1 (MQ=255)
                                                                                       ggACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCa             <  1:1201931/100‑1 (MQ=255)
                                                                                       ggACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCa             <  1:194395/100‑1 (MQ=255)
                                                                                       ggACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATCTTTCTGCCAATAGACATTATCAAGTGAGGAACGCa             >  1:177080/1‑100 (MQ=255)
                                                                                       ggACGGCGATGATACATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCa             >  1:382075/1‑100 (MQ=39)
                                                                                        gACGGCGATGTATCATCGCCGTCCATTTTTTTATATATTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCAt            >  2:5389/1‑100 (MQ=255)
                                                                                        gACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCAt            >  1:2372957/1‑100 (MQ=255)
                                                                                        gACGGCGATGATACATCGCCGTCCATTTTTCTACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCAt            >  1:715899/1‑100 (MQ=38)
                                                                                        gACGGCGATGATACATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCAt            >  1:2236596/1‑100 (MQ=39)
                                                                                        gACGGCGATGATACATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCAt            >  1:1222050/1‑100 (MQ=39)
                                                                                        gACGGCGATGATACATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCAt            >  1:1848807/1‑100 (MQ=39)
                                                                                        gACGGCGATGATACATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCAt            >  1:2000857/1‑100 (MQ=39)
                                                                                        gACGGCGATGATACATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCAAGAAATTTTTTTGAGAATAGACATTATCAAGTGAGGAACGAAt            >  2:1468068/1‑100 (MQ=37)
                                                                                         aCGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATc           <  2:2255029/100‑1 (MQ=255)
                                                                                         aCGGCGATGATACATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATc           >  1:33839/1‑100 (MQ=39)
                                                                                                  tATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTg  >  2:2268486/1‑100 (MQ=255)
                                                                                                  |                                                                                                    
TCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTG  >  NC_000913/1792142‑1792339

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: