Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,163,369 A→C 100% Q165P (CAA→CCA)  baeS → sensory histidine kinase BaeS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,163,3690AC97.2% 285.3 / ‑0.0 72Q165P (CAA→CCA) baeSsensory histidine kinase BaeS
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base C (33/37);  total (35/37)
Fisher's exact test for biased strand distribution p-value = 2.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.16e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TCTGGTCAACGGTGCGGAAGTTGGCGCGGTGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCAAACCAGCTGGTTGATTGTCGCCCTGGCAACGTTACTCGCGGCACTTGCCACTTTTCTGCTGGCGCGCGGTTTACTGGCACCGGTAAAACGAC  >  NC_000913/2163274‑2163461
                                                                                               |                                                                                            
tCTGGTCAACGGTGCGGAAGTTGGCGCGGTGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCCAAc                                                                                           <  1:539737/99‑1 (MQ=255)
  tGGTCAACGGTGCGGAAGTTGGCGCGGTGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCCAACCa                                                                                         >  2:899165/1‑99 (MQ=255)
   ggTCAACGGTGCGGAAGTTGGCGCGGTGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCCAACCAg                                                                                        >  2:245335/1‑99 (MQ=255)
        aCGGTGCGGAAGTTGGCGCGGTGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCCAACCAGCTGGt                                                                                   <  2:1639407/99‑1 (MQ=255)
           gTGCGGAAGTTGGCGCGGTGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCCAACCAGCTGGTTGa                                                                                >  2:1980154/1‑99 (MQ=255)
           gTGCGGAAGTTGGCGCGGTGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCCAACCAGCTGGTTGa                                                                                >  1:1968519/1‑99 (MQ=255)
           gTGCGGAAGTTGGCGCGGTGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCCAACCAGCTGGTTGa                                                                                <  2:455437/99‑1 (MQ=255)
            tGCGGAAGTTGGCGCGGTGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCCAACCAGCTGGTTGAt                                                                               <  2:2087983/99‑1 (MQ=255)
             gCGGAAGTTGGCGCGGTGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCCAACCAGCTGGTTGAtt                                                                              <  1:1075694/99‑1 (MQ=255)
             gCGGAAGTTGGCGCGGTGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCCAACCAGCTGGTTGAtt                                                                              <  1:1032526/99‑1 (MQ=255)
              cGGAAGTTGGCGCGGTGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCCAACCAGCTGGTTGATTg                                                                             <  1:1202997/99‑1 (MQ=255)
               ggAAGTTGGCGCGGTGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCCAACCAGCTGGTTGATTGt                                                                            >  2:333411/1‑99 (MQ=255)
                gAAGTTGGCGCGGTGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCCAACCAGCTGGTTGATTGTc                                                                           >  2:626612/1‑99 (MQ=255)
                  aGTTGGCGCGGTGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCCAACCAGCTGGTTGATTGTCGc                                                                         >  1:79838/1‑99 (MQ=255)
                       gcgcGGTGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCCAACCAGCTGGTTGATTGTCGCCCTgg                                                                    <  1:1443607/99‑1 (MQ=255)
                         gcgGTGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCCAACCAGCTGGTTGATTGTCGCCCTGGCa                                                                  >  2:287559/1‑99 (MQ=255)
                           ggTGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCCAACCAGCTGGTTGATTGTCGCCCTGGCAAc                                                                <  2:697300/99‑1 (MQ=255)
                           ggTGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCCAACCAGCTGGTTGATTGTCGCCCTGGCAAc                                                                <  2:1196197/99‑1 (MQ=255)
                           ggTGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCCAACCAGCTGGTTGATTGTCGCCCTGGCAAc                                                                <  1:818381/99‑1 (MQ=255)
                            gTGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCCAACCAGCTGGTTGATTGTCGCCCTGGCAACg                                                               <  1:1134532/99‑1 (MQ=255)
                             tGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCCAACCAGCTGGTTGATTGTCGCCCTGGCAACGt                                                              >  2:2244800/1‑99 (MQ=255)
                             tGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCCAACCAGCTGGTTGATTGTCGCCCTGGCAACGt                                                              <  1:713655/99‑1 (MQ=255)
                             tGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCCAACCAGCTGGTTGATTGTCGCCCTGGCAACGt                                                              >  2:382289/1‑99 (MQ=255)
                                tCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCCAACCAGCTGGTTGATTGTCGCCCTGGCAACGTTAc                                                           <  2:2146839/99‑1 (MQ=255)
                                 cGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCCAACCAGCTGGTTGATTGTCGCCCTGGCAACGTTACt                                                          <  2:1648256/99‑1 (MQ=255)
                                      cccccGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCCAACCAGCTGGTTGATTGTCGCCCTGGCAACGTTACTCGCgg                                                     <  2:1606406/99‑1 (MQ=255)
                                       ccccGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCCAACCAGCTGGTTGATTGTCGCCCTGGCAACGTTACTCGCGGc                                                    >  2:261651/1‑99 (MQ=255)
                                       ccccGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCCAACCAGCTGGTTGATTGTCGCCCTGGCAACGTTACTCGCGGc                                                    >  2:963768/1‑99 (MQ=255)
                                        cccGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCCAACCAGCTGGTTGATTGTCGCCCTGGCAACGTTACTCGCGGCa                                                   >  1:1654992/1‑99 (MQ=255)
                                         ccGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCCAACCAGCTGGTTGATTGTCGCCCTGGCAACGTTACTCGCGGCAc                                                  <  1:756441/99‑1 (MQ=255)
                                          cGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCCAACCAGCTGGTTGATTGTCGCCCTGGCAACGTTACTCGCGGCACt                                                 >  1:1473083/1‑99 (MQ=255)
                                          cGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCCAACCAGCTGGTTGATTGTCGCCCTGGCAACGTTACTCGCGGCACt                                                 <  1:2114468/99‑1 (MQ=255)
                                          cGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCCAACCAGCTGGTTGATTGTCGCCCTGGCAACGTTACTCGCGGCACt                                                 >  2:2052624/1‑99 (MQ=255)
                                           gTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCCAACCAGCTGGTTGATTGTCGCCCTGGCAACGTTACTCGCGGCACtt                                                <  2:205711/99‑1 (MQ=255)
                                           gTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCCAACCAGCTGGTTGATTGTCGCCCTGGCAACGTTACTCGCGGCACtt                                                <  2:1889116/99‑1 (MQ=255)
                                                gCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCCAACCAGCTGGTTGATTGTCGCCCTGGCAACGTTACTCGCGGCACTTGCCAc                                           <  1:256484/99‑1 (MQ=255)
                                                   gTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCCAACCAGCTGGTTGATTGTCGCCCTGGCAACGTTACTCGCGGCACTTGCCACttt                                        >  1:827193/1‑99 (MQ=255)
                                                   gTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCCAACCAGCTGGTTGATTGTCGCCCTGGCAACGTTACTCGCGGCACTTGCCACttt                                        <  1:1407315/99‑1 (MQ=255)
                                                   gTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCCAACCAGCTGGTTGATTGTCGCCCTGGCAACGTTACTCGCGGCACTTGCCACttt                                        >  1:2167161/1‑99 (MQ=255)
                                                     tAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCCAACCAGCTGGTTGATTGTCGCCCTGGCAACGTTACTCGCGGCACTTGCCACTTTTc                                      <  2:855930/99‑1 (MQ=255)
                                                     tAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCCAACCAGCTGGTTGATTGTCGCCCTGGCAACGTTACTCGCGGCACTTGCCACTTTTc                                      <  2:2128768/99‑1 (MQ=255)
                                                     tAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCCAACCAGCTGGTTGATTGTCGCCCTGGCAACGTTACTCGCGGCACTTGCCACTTTTc                                      <  2:2116937/99‑1 (MQ=255)
                                                       aCGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCCAACCAGCTGGTTGATTGTCGCCCTGGCAACGTTACTCGCGGCACTTGCCACTTTTctg                                    >  1:1177091/1‑99 (MQ=255)
                                                           gcAATACTGATATCAATTTCGATAAACAACAGCGGCCAACCAGCTGGTTGATTGTCGCCCTGGCAACGTTACTCGCGGCACTTGCCACTTTTCTGCTgg                                <  2:855570/99‑1 (MQ=255)
                                                            cAATACTGATATCAATTTCGATAAACAACAGCGGCCAACCAGCTGGTTGATTGTCGCCCTGGCAACGTTACTCGCGGCACTTGCCACTTTTCTGCTGgc                               >  1:1897532/1‑99 (MQ=255)
                                                              aTACTGATATCAATTTCGATAAACAACAGCGGCAAACCAGCTGGTTGATTGTCGCCCTGGCAACGTTACTCGCGGCACTTGCCACTTTTCTGCTGgcgc                             >  2:1523339/1‑99 (MQ=255)
                                                               tACTGATATCAATTTCGATAAACAACAGCGGCCAACCAGCTGGTTGATTGTCGCCCTGGCAACGTTACTCGCGGCACTTGCCACTTTTCTGCTGgcgcg                            >  2:2203955/1‑99 (MQ=255)
                                                                aCTGATATCAATTTCGATAAACAACAGCGGCCAACCAGCTGGTTGATTGTCGCCCTGGCAACGTTACTCGCGGCACTTGCCACTTTTCTGCTGgcgcgc                           >  2:335192/1‑99 (MQ=255)
                                                                 cTGATATCAATTTCGATAAACAACAGCGGCCAACCAGCTGGTTGATTGTCGCCCTGGCAACGTTACTCGCGGCACTTGCCACTTTTCTGCTGGcgcgcg                          >  2:1859955/1‑99 (MQ=255)
                                                                 cTGATATCAATTTCGATAAACAACAGCGGCCAACCAGCTGGTTGATTGTCGCCCTGGCAACGTTACTCGCGGCACTTGCCACTTTTCTGCTGGcgcgcg                          >  2:1228239/1‑99 (MQ=255)
                                                                 cTGATATCAATTTCGATAAACAACAGCGGCCAACCAGCTGGTTGATTGACGCCCTGGCAACGTTACTCGCGGCACTTGCCACTTTTCTGCTGGcgcgcg                          >  2:373291/1‑99 (MQ=255)
                                                                  tGATATCAATTTCGATAAACAACAGCGGCCAACCAGCTGGTTGATTGTCGCCCTGGCAACGTTACTCGCGGCACTTGCCACTTTTCTGCTGGCGCGCgg                         <  2:1789941/99‑1 (MQ=255)
                                                                  tGATAGCAATTTCGATAAACAACAGCGGCCATCCAGCTGGTTGATTGTCGCCCTGGCAACGTTACTCGCGGCACTTGCCACTTTTCTGCTGGCGCGCgg                         <  2:2132631/99‑1 (MQ=255)
                                                                   gATATCAATTTCGATAAACAACAGCGGCCAACCAGCTGGTTGATTGTCGCCCTGGCAACGTTACTCGCGGCACTTGCCACTTTTCTGCTGGCGCGCGGt                        >  2:1980169/1‑99 (MQ=255)
                                                                      atCAATTTCGATAAACAACAGCGGCCAACCAGCTGGTTGATTGTCGCCCTGGCAACGTTACTCGCGGCACTTGCCACTTTTCTGCTGGCGCGCGGTTTa                     >  1:1716993/1‑99 (MQ=255)
                                                                       tCAATTTCGTTTGACAACAGCGGACAACCCGCTGGTTGATTGTCGCCCTGGAAACGTTACCCGCGGCACTTGCGACTTTTCTGCTGGCGCGCAGTTTAc                    <  2:1260670/99‑1 (MQ=255)
                                                                       tCAATTTCGATAAACAACAGCGGCCAACCAGCTGGTTGATTGTCGCCCTGGCAACGTTACTCGCGGCACTTGCCACTTTTCTGCTGGCGCGCGGTTTAc                    <  1:860207/99‑1 (MQ=255)
                                                                       tCAATTTCGATAAACAACAGCGGCCAACCAGCTGGTTGATTGTCGCCCTGGCAACGTTACTCGCGGCACTTGCCACTTTTCTGCTGGCGCGCGGTTTAc                    <  2:2221682/99‑1 (MQ=255)
                                                                       tCAATTTCGATAAACAACAGCGGCCAACCAGCTGGTTGATTGTCGCCCTGGCAACGTTACTCGCGGCACTTGCCACTTTTCTGCTGGCGCGCGGTTTAc                    >  1:1303838/1‑99 (MQ=255)
                                                                        cAATTTCGATAAACAACAGCGGCCAACCAGCTGGTTGATTGTCGCCCTGGCAACGTTACTCGCGGCACTTGCCACTTTTCTGCTGGCGCGCGGTTTACt                   <  1:1711324/99‑1 (MQ=255)
                                                                         aaTTTCGATAAACAACAGCGGCCAACCAGCTGGTTGATTGTCGCCCTGGCAACGTTACTCGCGGCACTTGCCACTTTTCTGCTGGCGCGCGGTTTACTg                  <  1:1778850/99‑1 (MQ=255)
                                                                             tCGATAAACAACAGCGGCCAACCAGCTGGTTGATTGTCGCCCTGGCAACGTTACTCGCGGCACTTGCCACTTTTCTGCTGGCGCGCGGTTTACTGGCAc              >  1:1707134/1‑99 (MQ=255)
                                                                               gATAAACAACAGCGGCCAACCAGCTGGTTGATTGTCGCCCTGGCAACGTTACTCGCGGCACTTGCCACTTTTCTGCTGGCGCGCGGTTTACTGGCACCg            >  1:1505952/1‑99 (MQ=255)
                                                                                aTAAACAACAGCGGCCAACCAGCTGGTTGATTGTCGCCCTGGCAACGTTACTCGCGGCACTTGCCACTTTTCTGCTGGCGCGCGGTTTACTGGCACCgg           <  2:444256/99‑1 (MQ=255)
                                                                                 tAAACAACAGCGGCCAACCAGCTGGTTGATTGTCGCCCTGGCAACGTTACTCGCGGCACTTGCCACTTTTCTGCTGGCGCGCGGTTTACTGGCACCGGt          >  2:1702363/1‑99 (MQ=255)
                                                                                  aaaCAACAGCGGCAAACCAGCTGGTTGATTGTCGCCCTGGCAACGTTACTCGCGGCACTTGCCACTTTTCTGCTGGCGCGCGGTTTACTGGCACCGGTa         >  1:367576/1‑99 (MQ=255)
                                                                                    acaacaGCGGCCAACCAGCTGGTTGATTGTCGCCCTGGCAACGTTACTCGCGGCACTTGCCACTTTTCTGCTGGCGCGCGGTTTACTGGCACCGGTaaa       >  2:2000537/1‑99 (MQ=255)
                                                                                      aacaGCGGCCAACCAGCTGGTTGATTGTCGCCCTGGCAACGTTACTCGCGGCACTTGCCACTTTTCTGCTGGCGCGCGGTTTACTGGCACCGGTAAAAc     <  2:1465729/99‑1 (MQ=255)
                                                                                        caGCGGCCAACCAGCTGGTTGATTGTCGCCCTGGCAACGTTACTCGCGGCACTTGCCACTTTTCTGCTGGCGCGCGGTTTACTGGCACCGGTAAAACGa   <  1:1965811/99‑1 (MQ=255)
                                                                                        caGCGGCCAACCAGCTGGTTGATTGTCGCCCTGGCAACGTTACTCGCGGCACTTGCCACTTTTCTGCTGGCGCGCGGTTTACTGGCACCGGTAAAACGa   >  1:474108/1‑99 (MQ=255)
                                                                                        caGCGGCCAACCAGCTGGTTGATTGTCGCCCTGGCAACGTTACTCGCGGCACTTGCCACTTTTCTGCTGGCGCGCGGTTTACTGGCACCGGTAAAACGa   >  1:483484/1‑99 (MQ=255)
                                                                                         aGCGGCCAACCAGCTGGTTGATTGTCGCCCTGGCAACGTTACTCGCGGCACTTGCCACTTTTCTGCTGGCGCGCGGTTTACTGGCACCGGTAAAACGAc  <  1:447220/99‑1 (MQ=255)
                                                                                               |                                                                                            
TCTGGTCAACGGTGCGGAAGTTGGCGCGGTGATCGCCTCCCCCGTTGAGCGGTTAACGCGCAATACTGATATCAATTTCGATAAACAACAGCGGCAAACCAGCTGGTTGATTGTCGCCCTGGCAACGTTACTCGCGGCACTTGCCACTTTTCTGCTGGCGCGCGGTTTACTGGCACCGGTAAAACGAC  >  NC_000913/2163274‑2163461

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