Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,535,403 C→A 41.2% D156Y (GAT→TAT)  envZ ← sensory histidine kinase EnvZ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,535,4030CA41.2% 91.9 / 128.0 80D156Y (GAT→TAT) envZsensory histidine kinase EnvZ
Reads supporting (aligned to +/- strand):  ref base C (32/15);  new base A (17/16);  total (49/31)
Fisher's exact test for biased strand distribution p-value = 1.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.52e-01

TCGACCAACGGTCGGTTCTGGATACGAATAAACAGCCACGCCCCGCCTATCGCCAATAGCATAATCGCCAGCGTATAGCGGAACAGCGGAGAGAAATCGCCCTGATGAATTTCGGTCAGCGGCACGCGTACCCAGATATTGGGCGACAGCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAA  >  NC_000913/3535306‑3535498
                                                                                                 |                                                                                               
tCGACCAACGGTCGGTTCTGGATACGAATAAACAGCCACGCCCCGCCTATCGCCAATAGCATAATCGCCAGCGTATAGCGGAACAGCGGAGAGAAATCg                                                                                                <  2:1961325/99‑1 (MQ=255)
   aCCAACGGTCGGTTCTGGATACGAATAAACAGCCACGCCCCGCCTATCGCCAATAGCATAATCGCCAGCGTATAGCGGAACAGCGGAGAGAAATAGccc                                                                                             >  1:905822/1‑99 (MQ=255)
    ccAACAGTCGGTTCTGGATACGAATAAACAGCCACGCCCCGCCTATCGCCAATAGCATAATCGCCAGCGTATAGCGGAACAGCGGAGAGAAATCGCCCt                                                                                            <  1:2083796/99‑1 (MQ=255)
      aaCGGTCGGTTCTGGATACGAATAAACAGCCACGCCCCGCCTATCGCCAATAGCATAATCGCCAGCGTATAGCGGAACAGCGGAGAGAAATAGCCCtga                                                                                          <  1:947369/99‑1 (MQ=255)
      aaCGGTCGGTTCTGGATACGAATAAACAGCCACGCCCCGCCTATCGCCAATAGCATAATCGCCAGCGTATAGCGGAACAGCGGAGAGAAATAGCCCtga                                                                                          >  1:187639/1‑99 (MQ=255)
           tcggtTCTGGATACGAATAAACAGCCACGCCCCGCCTATCGCCAATAGCATAATCGCCAGCGTATAGCGGAACAGCGGAGAGAAATAGCCCTGATGAAt                                                                                     <  2:1729983/99‑1 (MQ=255)
            cggtTCTGGATACGAATAAACAGCCACGCCCCGCCTATCGCCAATAGCATAATCGCCAGCGTATAGCGGAACAGCGGAGAGAAATAGCCCTGATGAAtt                                                                                    <  2:2292236/99‑1 (MQ=255)
              gtTCTGGATACGAATAAACAGCCACGCCCCGCCTATCGCCAATAGCATAATCGCCAGCGTATAGCGGAACAGCGGAGAGAAATAGCCCTGATGAATTTc                                                                                  >  2:2212791/1‑99 (MQ=255)
               ttCTGGATACGAATAAACAGCCACGCCCCGCCTATCGCCAATAGCATAATCGCCAGCGTATAGCGGAACAGCGGAGAGAAATAGCCCTGATGAATTTCg                                                                                 >  2:15065/1‑99 (MQ=255)
                tCTGGATACGAATAAACAGCCACGCCCCGCCTATCGCCAATAGCATAATCGCCAGCGTATAGCGGAACAGCGGAGAGAAATCGCCCTGATGAATTTCgg                                                                                >  1:1861472/1‑99 (MQ=255)
                         gAATAAACAGCCACGCCCCGCCTATCGCCAATAGCATAATCGCCAGCGTATAGCGGAACAGCGGAGAGAAATAGCCCTGATGAATTTCGGTCAGCGGCa                                                                       <  2:1942633/99‑1 (MQ=255)
                          aaTAAACAGCCACGCCCCGCCTATCGCCAATAGCATAATCGCCAGCGTATAGCGGAACAGCGGAGAGAAATCGCCCTGATGAATTTCGGTCAGCGGCAc                                                                      <  1:1597093/99‑1 (MQ=255)
                           aTAAACAGCCACGCCCCGCCTATCGCCAATAGCATAATCGCCAGCGTATAGCGGAACAGCGGAGAGAAATCGCCCTGATGAATTTCGGTCAGCGGCAcg                                                                     >  1:1704906/1‑99 (MQ=255)
                             aaaCAGCCACGCCCCGCCTATCGCCAATAGCATAATCGCCAGCGTATAGCGGAACAGCGGAGAGAAATAGCCCTGATGAATTTCGGTCAGCGGCAcgcg                                                                   <  2:1978731/99‑1 (MQ=255)
                             aaaCAGCCACGCCCCGCATATCGCCAATAGCATAATCGCCAGCGTATAGCGGGACAGCGGAGAGAAATCGCCCGGATGAAGTTCGGTCAGCGGCAcgcg                                                                   <  2:2276066/99‑1 (MQ=255)
                                cAGCCACGCCCCGCCTATCGCCAATAGCATAATCGCCAGCGTATAGCGGAACAGCGGAGAGAAATCGCCCTGATGAATTTCGGTCAGCGGCACGCGTAc                                                                >  2:1428748/1‑99 (MQ=255)
                                cAGCCACGCCCCGCCTATCGCCAATAGCATAATCGCCAGCGTATAGCGGAACAGCGGAGAGAAATAGCCCTGATGAATTTCGGTCAGCGGCACGCGTAc                                                                >  2:362785/1‑99 (MQ=255)
                                cAGCCACGCCCCGCCTATCGCCAATAGCATAATCGCCAGCGTATAGCGGAACAGCGGAGAGAAATAGCCCTGATGAATTTCGGTCAGCGGCACGCGTAc                                                                <  2:2178812/99‑1 (MQ=255)
                                 aGCCACGCCCCGCCTATCGCCAATAGCATAATCGCCAGCGTATAGCGGAACAGCGGAGAGAAATAGCCCTGATGAATTTCGGTCAGCGGCACGCGTAcc                                                               <  1:601609/99‑1 (MQ=255)
                                  gCCACGCCCCGCCTATCGCCAATAGCATAATCGCCAGCGTATAGCGGAACAGCGGAGAGAAATAGCCCTGATGAATTTCGGTCAGCGGCACGCGTAccc                                                              >  2:2055473/1‑99 (MQ=255)
                                   ccACGCCCCGCCTATCGCCAATAGCATAATCGCCAGCGTATAGCGGAACAGCGGAGAGAAATCGCCCTGATGAATTTCGGTCAGCGGCACGCGTACCCa                                                             >  2:1826349/1‑99 (MQ=255)
                                     aCGCCCCGCCTATCGCCAATAGCATAATCGCCAGCGTATAGCGGAACAGCGGAGAGAAATCGCCCTGATGAATTTCGGTCAGCGGCACGCGTACCCAGa                                                           <  1:1434803/99‑1 (MQ=255)
                                      cGCCCCGCCTATCGCCAATAGCATAATCGCCAGCGTATAGCGGAACAGCGGAGAGAAATAGCCCTGATGAATTTCGGTCAGCGGCACGCGTACCCAGat                                                          >  1:1860420/1‑99 (MQ=255)
                                         cccGCCTATCGCCAATAGCATAATCGCCAGCGTATAGCGGAACAGCGGAGAGAAATCGCCCTGATGAATTTCGGTCAGCGGCACGCGTACCCAGATAtt                                                       >  1:1362145/1‑99 (MQ=255)
                                          ccGCCTATCGCCAATAGCATAATCGCCAGCGTATAGCGGAACAGCGGAGAGAAATCGCCCTGATGAATTTCGGTCAGCGGCACGCGTACCCAGATATTg                                                      >  2:335546/1‑99 (MQ=255)
                                          ccGCCTATCGCCAATAGCATAATCGCCAGCGTATAGCGGAACAGCGGAGAGAAATAGCCCTGATGAATTTCGGTCAGCGGCACGCGTACCCAGATATTg                                                      >  1:1825955/1‑99 (MQ=255)
                                           cGCCTATCGCCAATAGCATAATCGCCAGCGTATAGCGGAACAGCGGAGAGAAATCGCCCTGATGAATTTCGGTCAGCGGCACGCGTACCCAGATATTgg                                                     <  1:1498135/99‑1 (MQ=255)
                                             ccTATCGCCAATAGCATAATCGCCAGCGTATAGCGGAACAGCGGAGAGAAATCGCCCTGATGAATTTCGGTCAGCGGCACGCGTACCCAGATATTGGGc                                                   >  2:421652/1‑99 (MQ=255)
                                              cTATCGCCAATAGCATAATCGCCAGCGTATAGCGGAACAGCGGAGAGAAATCGCCCTGATGAATTTCGGTCAGCGGCACGCGTACCCAGAAATTGGGCg                                                  >  1:214797/1‑99 (MQ=255)
                                               tATCGCCAATAGCATAATCGCCAGCGTATAGCGGAACAGCGGAGAGAAATCGCCCTGATGAATTTCGGTCAGCGGCACGCGTACCCAGATATTGGGCGa                                                 >  2:2083390/1‑99 (MQ=255)
                                                aTCGCCAATAGCATAATCGCCAGCGTATAGCGGAACAGCGGAGAGAAATCGCCCTGATGAATTTCGGTCAGCGGCACGCGTACCCAGATATTGGGCGAc                                                >  2:1508730/1‑99 (MQ=255)
                                                 tCGCCAATAGCATAATCGCCAGCGTATAGCGGAACAGCGGAGAGAAATCGCCCTGATGAATTTCGGTCAGCGGCACGCGTACCCAGATATTGGGCGACa                                               <  2:2072998/99‑1 (MQ=255)
                                                 tCGCCAATAGCATAATCGCCAGCGTATAGCGGAACAGCGGAGAGAAATCGCCCTGATGAATTTCGGTCAGCGGCACGCGTACCCAGATATTGGGCGACa                                               >  2:1121477/1‑99 (MQ=255)
                                                 tCGCCAATAGCATAATCGCCAGCGTATAGCGGAACAGCGGAGAGAAATCGCCCTGATGAATTTCGGTCAGCGGCACGCGTACCCAGATATTGGGCGACa                                               >  1:928402/1‑99 (MQ=255)
                                                 tCGCCAATAGCATAATCGCCAGCGTATAGCGGAACAGCGGAGAGAAATCGCCCTGATGAATTTCGGTCAGCGGCACGCGTACCCAGATATTGGGCGACa                                               >  1:1950082/1‑99 (MQ=255)
                                                 tCGCCAATAGCATAATCGCCAGCGTATAGCGGAACAGCGGAGAGAAATAGCCCTGATGAATTTCGGTCAGCGGCACGCGTACCCAGATATTGGGCGACa                                               >  1:1181125/1‑99 (MQ=255)
                                                     cAATAGCATAATCGCCAGCGTATAGCGGAACAGCGGAGAGAAATAGCCCTGATGAATTTCGGTCAGCGGCACGCGTACCCAGATATTGGGCGACAGCCa                                           <  2:1233225/99‑1 (MQ=255)
                                                        tAGCATAATCGCCAGCGTATAGCGGAACAGCGGAGAGAAATCGCCCTGATGAATTTCGGTCAGCGGCACGCGTACCCAGATATTGGGCGACAGCCAGGt                                        >  1:1423889/1‑99 (MQ=255)
                                                          gCATAATCGCCAGCGTATAGCGGAACAGCGGAGAGAAATAGCCCTGATGAATTTCGGTCAGCGGCACGCGTACCCAGATATTGGGCGACAGCCAGGttt                                      >  2:1409793/1‑99 (MQ=255)
                                                          gCATAATCGCCAGCGTATAGCGGAACAGCGGAGAGAAATAGCCCTGATGAATTTCGGTCAGCGGCACGCGTACCCAGATATTGGGCGACAGCCAGGttt                                      <  1:1662002/99‑1 (MQ=255)
                                                            aTAATCGCCAGCGTATAGCGGAACAGCGGAGAGAAATCGCCCTGATGAATTTCGGTCAGCGGCACGCGTACCCAGATATTGGGCGACAGCCAGGTTTTc                                    >  2:2316276/1‑99 (MQ=255)
                                                             tAATCGCCAGCGTATAGCGGAACAGCGGAGAGAAATCGCCCTGATGAATTTCGGTCAGCGGCACGCGTACCCAGATATTGGGCGACAGCCAGGTTTTCa                                   >  1:1292249/1‑99 (MQ=255)
                                                             tAATCGCCAGCGTATAGCGGAACAGCGGAGAGAAATCGCCCTGATGAATTTCGGTCAGCGGCACGCGTACCCAGATATTGGGCGACAGCCAGGTTTTCa                                   >  2:755045/1‑99 (MQ=255)
                                                             tAATCGCCAGCGTATAGCGGAACAGCGGAGAGAAATCGCCCTGATGAATTTCGGTCAGCGGCACGCGTACCCAGATATTGGGCGACAGCCAGGTTTTCa                                   >  2:1908672/1‑99 (MQ=255)
                                                             tAATCGCCAGCGTATAGCGGAACAGCGGAGAGAAATAGCCCTGATGAATTTCGGTCAGCGGCACGCGTACCCAGATATTGGGCGACAGCCAGGTTTTCa                                   >  1:2189121/1‑99 (MQ=255)
                                                             tAATCGCCAGCGTATAGCGGAACAGCGGAGAGAAATAGCCCTGATGAATTTCGGTCAGCGGCACGCGTACCCAGATATTGGGCGACAGCCAGGTTTTCa                                   <  2:1196115/99‑1 (MQ=255)
                                                               aTCGCCAGCGTATAGCGGAACAGCGGAGAGAAATCGCCCTGATGAATTTCGGTCAGCGGCACGCGTACCCAGATATTGGGCGACAGCCAGGTTTTCAGc                                 <  2:395743/99‑1 (MQ=255)
                                                                tCGCCAGCGTATAGCGGAACAGCGGAGAGAAATCGCCCTGATGAATTTCGGTCAGCGGCACGCGTACCCAGATATTGGGCGACAGCCAGGTTTTCAGcc                                >  1:730174/1‑99 (MQ=255)
                                                                tCGCCAGCGTATAGCGGAACAGCGGAGAGAAATCGCCCTGATGAATTTCGGTCAGCGGCACGCGTACCCAGATATTGGGCGACAGCCAGGTTTTCAGcc                                >  1:651161/1‑99 (MQ=255)
                                                                  gCCAGCGTATAGCGGAACAGCGGAGAGAAATCGCCCTGATGAATTTCGGTCAGCGGCACGCGTACCCAGATATTGGGCGACAGCCAGGTTTTCAGCCAg                              >  2:454466/1‑99 (MQ=255)
                                                                   ccAGCGTATAGCGGAACAGCGGAGAGAAATCGCCCTGATGAATTTCGGTCAGCGGCACGCGTACCCAGATATTGGGCGACAGCCAGGTTTTCAGCCAga                             <  1:2127978/99‑1 (MQ=255)
                                                                      gCGTATAGCGGAACAGCGGAGAGAAATCGCCCTGATGAATTTCGGTCAGCGGCACGCGTACCCAGATATTGGGCGACAGCCAGGTTTTCAGCCAgacga                          >  2:860064/1‑99 (MQ=255)
                                                                       cGTATAGCGGAACAGCGGAGAGAAATAGCCCTGATGAATTTCGGTCAGCGGCACGCGTACCCAGATATTGGGCGACAGCCAGGTTTTCAGCCAgacgac                         <  2:1837377/99‑1 (MQ=255)
                                                                        gTATAGCGGAACAGCGGAGAGAAATCGCCCTGATGAATTTCGGTCAGCGGCACGCGTACCCAGATATTGGGCGACAGCCAGGTTTTCAGCCAGACGACa                        <  1:1830981/99‑1 (MQ=255)
                                                                        gTATAGCGGAACAGCGGAGAGAAATAGCCCTGATGAATTTCGGTCAGCGGCACGCGTACCCAGATATTGGGCGACAGCCAGGTTTTCAGCCAGACGACa                        >  1:1932204/1‑99 (MQ=255)
                                                                         tataGCGGAACAGCGGAGAGAAATCGCCCTGATGAATTTCGGTCAGCGGCACGCGTACCCAGATATTGGGCGACAGCCAGGTTTTCAGCCAGACGACAg                       >  1:217611/1‑99 (MQ=255)
                                                                         tataGCGGAACAGCGGAGAGAAATCGCCCTGATGAATTTCGGTCAGCGGCACGCGTACCCAGATATTGGGCGACAGCCAGGTTTTCAGCCAGACGACAg                       >  1:92692/1‑99 (MQ=255)
                                                                          ataGCGGAACAGCGGAGAGAAATCGCCCTGATGAATTTCGGTCAGCGGCACGCGTACCCAGATATTGGGCGACAGCCAGGTTTTCAGCCAGACGACAgg                      >  2:418287/1‑99 (MQ=255)
                                                                           taGCGGAACAGCGGAGAGAAATCGCCCTGATGAATTTCGGTCAGCGGCACGCGTACCCAGATATTGGGCGACAGCCAGGTTTTCAGCCAGACGACAGGc                     <  2:739239/99‑1 (MQ=255)
                                                                           taGCGGAACAGCGGAGAGAAATCGCCCTGATGAATTTCGGTCAGCGGCACGCGTACCCAGATATTGGGCGACAGCCAGGTTTTCAGCCAGACGACAGGc                     >  1:1206877/1‑99 (MQ=255)
                                                                           taGCGGAACAGCGGAGAGAAATAGCCCTGATGAATTTCGGTCAGCGGCACGCGTACCCAGATATTGGGCGACAGCCAGGTTTTCAGCCAGACGACAGGc                     <  2:1159052/99‑1 (MQ=255)
                                                                           taGCGGAACAGCGGAGAGAAATAGCCCTGATGAATTTCGGTCAGCGGCACGCGTACCCAGATATTGGGCGACAGCCAGGTTTTCAGCCAGACGACAGGc                     <  2:1894417/99‑1 (MQ=255)
                                                                               ggAACAGCGGAGAGAAATCGCCCTGATGAATTTCGGTCAGCGGCACGCGTACCCAGATATTGGGCGACAGCCAGGTTTTCAGCCAGACGACAGGCGAAc                 >  2:1055899/1‑99 (MQ=255)
                                                                                gAACAGCGGAGAGAAATAGCCCTGATGAATTTCGGTCAGCGGCACGCGTACCCAGATATTGGGCGACAGCCAGGTTTTCAGCCAGACGACAGGCGAACt                <  1:206625/99‑1 (MQ=255)
                                                                                 aaCAGCGGAGAGAAATCGCCCTGATGAATTTCGGTCAGCGGCACGCGTACCCAGATATTGGGCGACAGCCAGGTTTTCAGCCAGACGACAGGCGAACtt               >  1:1170138/1‑99 (MQ=255)
                                                                                  aCAGCGGAGAGAAATCGCCCTGATGAATTTCGGTCAGCGGCACGCGTACCCAGATATTGGGCGACAGCCAGGTTTTCAGCCAGACGACAGGCGAACttt              <  2:2162983/99‑1 (MQ=255)
                                                                                  aCAGCGGAGAGAAATCGCCCTGATGAATTTCGGTCAGCGGCACGCGTACCCAGATATTGGGCGACAGCCAGGTTTTCAGCCAGACGACAGGCGAACttt              <  1:1038071/99‑1 (MQ=255)
                                                                                  aCAGCGGAGAGAAATAGCCCTGATGAATTTCGGTCAGCGGCACGCGTACCCAGATATTGGGCGACAGCCAGGTTTTCAGCCAGACGACAGGCGAACttt              >  1:1366500/1‑99 (MQ=255)
                                                                                   cAGCGGAGAGAAATAGCCCTGATGAATTTCGGTCAGCGGCACGCGTACCCAGATATTGGGCGACAGCCAGGTTTTCAGCCAGACGACAGGCGAACtttt             >  1:2313729/1‑99 (MQ=255)
                                                                                     gCGGAGAGAAATCGCCCTGATGAATTTCGGTCAGCGGCACGCGTACCCAGATATTGGGCGACAGCCAGGTTTTCAGCCAGACGACAGGCGAACTTttgt           >  1:2328390/1‑99 (MQ=255)
                                                                                     gCGGAGAGAAATCGCCCTGATGAATTTCGGTCAGCGGCACGCGTACCCAGATATTGGGCGACAGCCAGGTTTTCAGCCAGACGACAGGCGAACTTttgt           >  1:1427011/1‑99 (MQ=255)
                                                                                       ggAGAGAAATCGCCCTGATGAATTTCGGTCAGCGGCACGCGTACCCAGATATTGGGCGACAGCCAGGTTTTCAGCCAGACGACAGGCGAACTTttgttg         >  2:329924/1‑99 (MQ=255)
                                                                                         agagaAATAGCCCTGATGAATTTCGGTCAGCGGCACGCGTACCCAGATATTGGGCGACAGCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGAc       >  2:1656508/1‑99 (MQ=255)
                                                                                          gagaAATCGCCCTGATGAATTTCGGTCAGCGGCACGCGTACCCAGATATTGGGCGACAGCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGAcc      >  1:1859255/1‑99 (MQ=255)
                                                                                          gagaAATAGCCCTGATGAATTTCGGTCAGCGGCACGCGTACCCAGATATTGGGCGACAGCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGAcc      >  2:1924468/1‑99 (MQ=255)
                                                                                           agaAATAGCCCTGATGAATTTCGGTCAGCGGCACGCGTACCCAGATATTGGGCGACAGCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCt     <  2:238162/99‑1 (MQ=255)
                                                                                            gaAATAGCCCTGATGAATTTCGGTCAGCGGCACGCGTACCCAGATATTGGGCGACAGCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTc    <  1:375858/99‑1 (MQ=255)
                                                                                             aaaTAGCCCTGATGAATTTCGGTCAGCGGCACGCGTACCCAGATATTGGGCGACAGCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCa   >  1:1951131/1‑99 (MQ=255)
                                                                                              aaTCGCCCTGATGAATTTCGGTCAGCGGCACGCGTACCCAGATATTGGGCGACAGCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCaa  <  2:1009026/99‑1 (MQ=255)
                                                                                              aaTCGCCCTGATGAATTTCGGTCAGCGGCACGCGTACCCAGATATTGGGCGACAGCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCaa  <  1:276039/99‑1 (MQ=255)
                                                                                                 |                                                                                               
TCGACCAACGGTCGGTTCTGGATACGAATAAACAGCCACGCCCCGCCTATCGCCAATAGCATAATCGCCAGCGTATAGCGGAACAGCGGAGAGAAATCGCCCTGATGAATTTCGGTCAGCGGCACGCGTACCCAGATATTGGGCGACAGCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAA  >  NC_000913/3535306‑3535498

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: