Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,890,197 A→G 8.8% W40R (TGG→CGG)  PP_1695 ← sodium‑solute symporter/two‑component system sensor histidine kinase/response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,890,1970AG 8.8% 123.9 / 2.2 45W40R (TGG→CGG) PP_1695sodium‑solute symporter/two‑component system sensor histidine kinase/response regulator
Reads supporting (aligned to +/- strand):  ref base A (18/23);  new base G (2/2);  total (20/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

CTTCTGCAGCACCCAAGGGGCGAAGATCAGCAGCAGTACCGGGCCCAGGTAGATGGGCAGGAACGCCCAAAGCTGTTCGGCGGCCTGGCCGAC‑CGCGCCGAAGAAGGTCCAGCTGGTGCAGTACACGGCCAGCGACAGGCTGTACACCCAGGCACGCAGCTTTGGCGGCAATGGCGTGCTGCGGCGGTCGCCATAGAAGGCGATGGCGAACATGATGGCCATATAGGCCAGGGCGACCACGGCGATCAGCCCGCTGGACAACGACATGAAGACTCCGAAGCAAAAAGGTAACGCG  >  NC_002947/1890048‑1890342
                                                                                                                                                      |                                                                                                                                                 
cttcgacagcaCCCACGGGGCGAAGATCAGCAGCAGTACCGGGCCCAGGTAGATGGGCAGGAACGCCCAAAGCTGTTCGGCGGCCTAGACTAC‑CGCGTCGAAGAAGGACCAGCTGGTGCAGTACACGGTCAGCGACAGGCTGTACACCCa                                                                                                                                                   <  1:1077780/144‑1 (MQ=255)
      cagcaCCCAAGGGGCGAAGATCAGCAGCAGTACCGGGCCCAGGTAGATGGGCAGGAACGCCCAAAGCTGTTCGGCGGCCTGGCCTAC‑CGCGCCGAAGAAGGTCCAGCTGGTGCAGTACACGGCCAGCGACAGGCTGTACACCCAGGCACg                                                                                                                                             <  2:266685/150‑1 (MQ=255)
       agcaCCCAAGGGGCGAAGATCAGCAGCAGTACCGGGCCCAGGTAGATGGGCAGGAACGCCCAAAGCTGTTCGGCGGCCTGGCCGAC‑CGCGCCGAAGAAGGTCCAGCTGGTGCAGTACACGGCCAGCGACAGGCTGTACACCCAGGCACGc                                                                                                                                            >  2:196980/1‑150 (MQ=255)
            ccAAGGGGCGAAGATCAGCAGCAGTACCGGGCCCAGGTAGATGGGCAGGAACGCCCAAAGCTGTTCGGCGGCCTGGCCGAC‑CGCGCCGAAGAAGGTCCAGCTGGTGCAGTACACGGCCAGCGACAGGCTGTACACCCAGGCACGCAGCtt                                                                                                                                       >  1:968178/1‑150 (MQ=255)
               aGGGCCGAAGATCAGCAGCAGTAACGGGCCCAGGTAGATGGGCAGGAACGCCCAAAGCTGTTCGGCGGCCTGGCCGAC‑CGCGCCGAAGAAGGTCCAGCTGGTGCAGTACACGGCCAGCGACAGGCTGTACACCCAGGCACGCAGCTTTgg                                                                                                                                    <  1:1003971/150‑1 (MQ=255)
                     gAAGATCAGCAGCAGTACCGGGCCCAGGTAGATGGGCAGGAACGCCCAAAGCTGTTCGGCGGCCTGGCCGAC‑CGCGCCGAAGAAGGTCCAGCTGGTGCAGTACACGGCCAGCGACAGGCTGTACACCCAGGCACGCAGCTTTGGCGGCaa                                                                                                                              >  2:843361/1‑150 (MQ=255)
                      aaGATCAGCAGCAGTACCGGGCCCAGGTAGATGGGCAGGAACGCCCAAAGCTGTTCGGCGGCCTGGCCGAC‑CGCGCCGAAGAAGGTCCAGCTGGTGCAGTACACGGCCAGCGACAGGCTGTACACCCAGGCACGCAGCTTTGGCGGCAAt                                                                                                                             >  2:408765/1‑150 (MQ=255)
                       aGATCAGCAGCAGTACCGGGCCCAGGTAGATGGGCAGGAACGCCCAAAGCTGTTCGGCGGCCTGGCCGAC‑CGCGCCGAAGAAGGTCCAGGTGGTGCAGTACACGGCCAGCGACAGGCTGTACACCCAGGCACGCAGCTTTGGCGGCAATg                                                                                                                            <  2:538989/150‑1 (MQ=255)
                       aGATCAGCAGCAGTACCGGGCCCAGGTAGATGGGCAGGAACGCCCAAAGCTGTTCGGCGGCCTGGCCGAC‑CGCGCCGAAGAAGGTCCAGCTGGTGCAGTACACGGCCAGCGACAGGCTGTACACCCAGGCACGCAGCTTTGGCGGCAATg                                                                                                                            <  1:49058/150‑1 (MQ=255)
                        gATCAGCAGCAGTACCGGGCCCAGGTAGATGGGCAGGAACGCCCAAAGCTGTTCGGCGGCCTGGCCGAC‑CGCGCCGAAGAAGGTCCAGCTGGTGCAGTACACGGCCAGCGACAGGCTGTACACCCAGGCACGCAGCTTTGGCGGCAATgg                                                                                                                           >  2:316161/1‑150 (MQ=255)
                           cagcagcagTACCGGGCCCAGGTAGATGGGCAGGACCGCCCAAAGCTGTTCGGCGGCCTGGCCGTC‑CGGGCCGAAGAGGGTCCAGCTGGTGCAGTACACGGCCAGCGACAGGCTGTACACCCAGGCACGCAGCTTTGGCGGCAATGGCGt                                                                                                                        <  2:124930/150‑1 (MQ=255)
                               agcagTACCGGGCCCAGGTAGATGGGCAGGAACGCCCAAAGCTGTTCGGCGGCCTGGCCGAC‑CGCGCCGAAGAAGGTCCAGCTGGTGCAGTACACGGCCAGCGACAGGCTGTACACACAGGCACGCAGCTTTGGCGGCAATGGCGtgctg                                                                                                                    >  1:4545/1‑150 (MQ=255)
                                   gTACCGGGCCCAGGTAGATGGGCAGGAACGCCCAAAGCTGTTCGGCGGCCTGGCCGAC‑CGCGCCGAAGAAGGTCCAGCTGGTGCAGTACACGGCCAGCGACAGGCTGTACACCCAGGCACGCAGCTTTGGCGGCAATGGCGTGCTgcggc                                                                                                                <  1:319002/150‑1 (MQ=255)
                                         ggCCCAGGTAGATGGGCAGGAACGCCCAAAGCTGTTCGGCGGCCTGGCCGAC‑CGCGCCGAAGAAGGTCCAGCTGGTGCAGTACACGGCCAGCGACAGGCTGTACACCCAGGCACGCAGCTTTGGCGGCAATGGCGTGCTGCGGCGGTCGc                                                                                                          <  2:4545/150‑1 (MQ=255)
                                             cAGGTAGATGGGCAGGAACGCCCAAAGCTGTTCGGCGGCCTGGCCGAC‑CGCGCCGAAGAAGGTCCAGCTGGTGCAGTACACGGCCAGCGACAGGCTGTACACCCAGGCACGCAGCTTTGGCGGCAATGGCGTGCTGCGGCGGTCGCCATa                                                                                                      <  1:196980/150‑1 (MQ=255)
                                                gTAGATGGGCAGGAACGCCCAAAGCTGTTCGGCGGCCTGGCCGACGTGCGAC‑TTTAACTGCTAGCTGGTGCAGTACACGGCCAGCGACAGGCTGTACACCCAGGCACGCAGCTTTGGCGGCAATGGCGTGCTGCGGCGGTCGCCATAGaa                                                                                                   <  2:870055/150‑1 (MQ=255)
                                                       ggCAGGAACGCCCAAAGCTGTTCGGCGGCCTGGCCGAC‑CGCGCCGAAGAAGGTCCAGCTGGTGCAGTACACGGCCAGCGACAGGCTGTACACCCAGGCACGCAGCTTTGGCGGCAATGGCGTGCTGCGGCGGTCGCCATAGAAGGCGATg                                                                                            <  1:324239/150‑1 (MQ=255)
                                                       ggCAGGAACGCCCAAAGCTGTTCGGCGGCCTGGCCGAC‑CGCGCCGAAGAAGGTCCAGCTGGTGCAGTACACGGCCAGCGACAGGCTGTACACCCAGGCACGCAGCTTTGGCGGCAATGGCGTGCTGCGGCGGTCGCCATAGAAGGCGATg                                                                                            >  2:739233/1‑150 (MQ=255)
                                                       ggCAGGAACGCCCAAAGCTGTTCGGCGGCCTGGCCGAC‑CGCGCCGAAGAAGGTCCAGCTGGTGCAGTACACGGCCAGCGACAGGCTGTACACCCAGGCACGCAGCTTTGGCGGCAATGGCGTGCTGCGGCGGTCGCCATAGAAGGCGATg                                                                                            <  1:739233/150‑1 (MQ=255)
                                                          aGGAACGCCCAAAGCTGTTCGGCGGCCTGGCCGAC‑CGCGCCGAAGAAGGTCCAGCTGGTGCAGTACACGGCCAGCGACAGGCTGTACACCCAGGCACGCAGCTTTGGCGGCAATGGCGTGCTGCGGCGGTCGCCATAGAAGGCGATGGCg                                                                                         >  2:901567/1‑150 (MQ=255)
                                                            gAACGCCCAAAGCTGTTCGGCGGCCTGGCCGAC‑CGCGCCGAAGAAGGTCCAGCTGGTGCAGTACACGGCCAGCGACAGGCTGTACACCCAGGCACGCAGCTTTGGCGGCAATGGCGTGCTGCGGCGGTCGCCATAGAAGGCGATGGCGaa                                                                                       >  2:985375/1‑150 (MQ=255)
                                                                 cccAAGGCTGTTCGGCGGCCTGGCCGAC‑CGCGCCGAAGAAGGTCCAGCTGGTGCAGTACACGGCCAGCGACAGGCTGTACACCCAGGCACGCAGCTTTGGCGGCAATGGCg                                                                                                                         <  1:616142/111‑1 (MQ=255)
                                                                 cccAAGGCTGTTCGGCGGCCTGGCCGAC‑CGCGCCGAAGAAGGTCCAGCTGGTGCAGTACACGGCCAGCGACAGGCTGTACACCCAGGCACGCAGCTTTGGCGGCAATGGCg                                                                                                                         >  2:616142/1‑111 (MQ=255)
                                                                   cAAAGCTGTTCGGCGGCCTGGCCGAC‑CGCGCCGAAGAAGGTCCAGCTGGTGCAGTACACGGCCAGCGACAGGCTGTACACCCAGGCACGCAGCTTTGGCGGCAATGGCGTGCTGCGGCGGTCGCCATAGAAGGCGATGGCGAACatgatg                                                                                <  1:901567/150‑1 (MQ=255)
                                                                    aaaGCTGTTCGGCGGCCTGGCCGAC‑CGCGCCGAAGAAGGTCCAGCTGGTGCAGTACACGGCCAGCGACAGGCTGTACACCCAGGCACGCAGCTTTGGCGGCAATGGCGTGCTGCGGCGGTCGCCATAGAAGGCGATGGCGAACATGATgg                                                                               <  1:703462/150‑1 (MQ=255)
                                                                      aGCTGTTCGGCGGCCTGGCCGAC‑CGCGCCGAAGAAGGTCCAGCTGGTGCAGTACACGGCCAGCGACAGGCTGTACACCCAGGCACGCAGCTTTGGCGGCAATGGCGTGCTGCGGCGGTCGCCATAGAAGGCGATGGCGAACATGATGGcc                                                                             <  2:1113723/150‑1 (MQ=255)
                                                                          gTTCGGCGGCCTGGCCGAC‑CGCGCCGAAGAAGGTCCAGCTGGTGCAGTACACGGCCAGCGACAGGCTGTACACCCGGGCACGCAGCTTTGGCGGCAATGGCGTGCTGCGGCGGt                                                                                                             >  2:629700/1‑114 (MQ=255)
                                                                          gTTCGGCGGCCTGGCCGAC‑CGCGCCGAAGAAGGTCCAGCTGGTGCAGTACACGGCCAGCGACAGGCTGTACACCCGGGCACGCAGCTTTGGCGGCAATGGCGTGCTGCGGCGGt                                                                                                             >  2:629483/1‑114 (MQ=255)
                                                                          gTTCGGCGGCCTGGCCGAC‑CGCGCCGAAGAAGGTCCAGCTGGTGAAGTACACGGCCAGCGACAGGCTGTACACCCGGGCACGCAGCTTTGGCGGCAATGGCGTGCTGCGGCGGt                                                                                                             <  1:629483/114‑1 (MQ=255)
                                                                          gTTCGGCGGCCTGGCCGAC‑CGCGCCGAAGAAGGTCCAGCTGGTGAAGTACACGGCCAGCGACAGGCTGTACACCCGGGCACGCAGCTTTGGCGGCAATGGCGTGCTGCGGCGGt                                                                                                             <  1:629700/114‑1 (MQ=255)
                                                                                     tGGCCGAC‑CGCGCCGAAGAAGGTCCAGCTGGTGCAGTACACGGCCAGCGACAGGCTGTACACCCAGGCACGCAGCTTTGGCGGCAATGGCGTGCTGCGGCGGTCGCCATAGAAGGCGATGGCGAACATGATGGCCATATAGGCCAGGGCg                                                              >  1:585581/1‑150 (MQ=255)
                                                                                         cGAC‑CGCGCCGAAGAAGGTCCAGCTGGTGCAGTACACGGCCAGCGACAGGCTGTACACCCAGGCACGCAGCTTTGGCGGCAATGGCGTGCTGCGGCGGTCGCCATAGAAGGCGATGGCGAACATGATGGCCATATAGGCCAGGGCGACCa                                                          <  1:589036/150‑1 (MQ=255)
                                                                                            c‑cGCGCCGAAGAAGGTCCAGCTGGTGCAGTACACGGCCAGCGACAGGCTGTACACCCAGGCACGCAGCTTTGGCGGCAATGGCGTGCTGCGGCGGTCGCCATAGAAGGCGATGGCGAACATGATGGCCATATAGGCCAGGGCGACCACgg                                                       <  1:815960/150‑1 (MQ=255)
                                                                                                   cGAAGAAGGTTCAGCTGGTGCAGTACACGTCCAGCGACAGGCTGTACACCCAGGGACGCAGCTTTGACGGTAATTGCGGCCTGTGGCGATCGCCATAGAAGGCGATGGCGAATATGATGGCCATATATGCCAGGGCGACCACGATGATCa                                                 <  2:443129/150‑1 (MQ=255)
                                                                                                            tCCAGCTGGTGCAGTACACGGCCAGCGACAGGCTGTACACCCAGGCACGCAGCTTTGGCGGCAATGGCGTGCTGCGGCGGTCGCCATAGAAGGCGATGGCGAACATGATGGCCATATAGGCCAGGGCGACCACGGCGATCAGCCCGCTgg                                        >  1:310636/1‑150 (MQ=255)
                                                                                                                gCTGGTGCAGTACACGGCCAGCGACAGGCTGTACACCCAGGCACGCAGCTTTGGCGGCAATGGCGTGCTGCGGCGGTCGCCATAGAAGGCGATGGCGAACATGATGGCCATATAGGCCAGGGCGACCACGGCGATCAGCCCGCTg                                         <  1:1127648/145‑1 (MQ=255)
                                                                                                                gCTGGTGCAGTACACGGCCAGCGACAGGCTGTACACCCAGGCACGCAGCTTTGGCGGCAATGGCGTGCTGCGGCGGTCGCCATAGAAGGCGATGGCGAACATGATGGCCATATAGGCCAGGGCGACCACGGCGATCAGCCCGCTg                                         >  2:1127648/1‑145 (MQ=255)
                                                                                                                  gatggCAGTACACGGCCAGCGACAGGCTGTACACCCAGGCACGCAGCTTTGGCGGCAATGGCGTGCTGCGGCGGTCGCCATAGAAGGCGATGGCGAACATGATGGCCATATAGGCCAGGGCGACCACGGCGATCAGCCCGCTGGACAACg                                  >  1:879518/5‑150 (MQ=255)
                                                                                                                         gTACACGGCCAGCGACAGGCTGTACACCCAGGCACGCAGCTTTGGCGGCAATGGCGTGCTGCGGCGGTCGCCATAGAAGGCGATGGCGAACATGATGGCCATATAGGCCAGGGCGACCACGGCGATCAGCCCGCTGGACAACGACATGaa                           >  1:1128024/1‑150 (MQ=255)
                                                                                                                           acacGGCCAGCGACAGGCTGTACACCCAGGCACGCAGCTTTGGCGGCGATGGCGTGCTGCGGCGGTCGCCATAGAAGGCGATGGCGAACATGATGGCCATATAGGCCAGGGCGACCACGGCGATCAGCCCGCTGGACAACGACATGAAGa                         <  2:919259/150‑1 (MQ=255)
                                                                                                                              cGGCCAGCGACAGGCTGTACACCCAGGCACGCAGCTTTGGCGGCAATGGCGTGCTGCGGCGGTCGCCATAGAAGGCGATGGCGAACATGATGGCCATATAGGCCAGGGCGACCACGGCGATCAGCCCGCTGGACAACGACATGAAGACTc                      <  2:968178/150‑1 (MQ=255)
                                                                                                                                 ccAGCAACAGGCTGTACACCCAGGCACGCAGCTTTGGCGGCAATGGCCTGCTGCGGCGGTCGCCATAGAAGGCGATGGCGAACATGATGGCCATATAGGCCAGGGCGACCACGGCGATCAGCCCGCTGGACAACGACATGAAGACTCCGa                   <  1:799521/150‑1 (MQ=255)
                                                                                                                                           gCTGTACACCCAGGCACGCAGCTTTGGCGGCAATGGCGTGCTGCGGCGGTCGCCATAGAAGGCGATGGCGAACATGATGGCCATATAGGCCAGGGCGACCACGGCGATCAGCCCGCTGGACAACGACATGAAGACTCCGAAGCAAAAAgg         >  1:697049/1‑150 (MQ=255)
                                                                                                                                             tGTACACCCAGGCACGCAGCTTTGGCGGCAATGGCGTGCTGCGGCGGTCGCCATAGAAGGCGATGGCGAACATGATGGCCATATAGGCCAGGGCGACCACGGCGATCAGCCCGCTGGACAACGACATGAAGACTCCGAAGCAAAAAGGTa       >  2:210320/1‑150 (MQ=255)
                                                                                                                                               tACACCCAGGCACGCAGCTTTGGCGGCAATGGCGTGCTGCGGCGGTCGCCATAGAAGGCGATGGCGAACATGATGGCCATATAGGCCAGGGCGACCACGGCGATCAGCCCGCTGGACAACGACATGAAGACTCCGAAGCAAAAAGGTAAc     >  2:965049/1‑150 (MQ=255)
                                                                                                                                                  acCCAGGCACGCAGCTTTGGCGGCAATGGCGTGCTGCGGCGGTCGCCATAGAAGGCGATGGCGAACATGATGGCCATATAGGCCAGGGCGACCACGGCGATCAGCCCGCTGGACAACGACATGAAGACTCCGAAGCAAAAAGGTAAcgcg  <  1:316161/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                 
CTTCTGCAGCACCCAAGGGGCGAAGATCAGCAGCAGTACCGGGCCCAGGTAGATGGGCAGGAACGCCCAAAGCTGTTCGGCGGCCTGGCCGAC‑CGCGCCGAAGAAGGTCCAGCTGGTGCAGTACACGGCCAGCGACAGGCTGTACACCCAGGCACGCAGCTTTGGCGGCAATGGCGTGCTGCGGCGGTCGCCATAGAAGGCGATGGCGAACATGATGGCCATATAGGCCAGGGCGACCACGGCGATCAGCCCGCTGGACAACGACATGAAGACTCCGAAGCAAAAAGGTAACGCG  >  NC_002947/1890048‑1890342

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