Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,289,430 T→G 9.1% K1053Q (AAG→CAG)  PP_3761 ← 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_0029474,289,4300TG 9.1% 116.7 / 2.6 42K1053Q (AAG→CAG) PP_3761two‑component system sensor histidine kinase/response regulator
Reads supporting (aligned to +/- strand):  ref base T (23/15);  new base G (2/2);  total (25/17)
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.89e-01

GGCGGATCAGGCGGGTGGCCTCGAAGCCGTCCATCTCCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGCGTTGACGTTCGTTGGACAGCTGCGATTCGACCTTGT  >  NC_002947/4289292‑4289576
                                                                                                                                           |                                                                                                                                                  
ggCGGATCAGGCGGGTGGCCTCGAAGCCGTCCATCTCCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTtgt                                                                                                                                                  <  5:246652/141‑1 (MQ=255)
ggCGGATCAGGCGGGTGGCCTCGAAGCCGTCCATCTCCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTtgt                                                                                                                                                  >  6:246652/1‑141 (MQ=255)
 gCGGATCAGGCGGGTGGCCTCGAAGCCGTCCATCTCCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAgg                                                                                                                                          <  4:20507/148‑1 (MQ=255)
   ggATCAGGCGGGTGGCCTCGAAGCCGTCCATCTCCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGgcg                                                                                                                                       >  8:153373/1‑149 (MQ=255)
        aGGCGGGTGGCCTCGAAGCCGTCCATCTCCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTGGAGCACAAGGgcgccatt                                                                                                                                  <  6:320019/149‑5 (MQ=255)
             ggTGGCCTCGAAGCCGTCCATCTCCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAggg                                                                                                                             <  7:153373/149‑1 (MQ=255)
                   cTCGAAGCCGTCCATCTCCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCtt                                                                                                                                                    <  8:169756/120‑1 (MQ=255)
                   cTCGAAGCCGTCCATCTCCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCtt                                                                                                                                                    >  7:169756/1‑120 (MQ=255)
                   cTCGAAGCCGTCCATCTCCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCa                                                                                                                           >  3:235629/1‑145 (MQ=255)
                    tCGAAGCCGTCCATCTCCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGa                                                                                                                       >  5:123111/1‑148 (MQ=255)
                           cGTCCATCTCCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACg                                                                                                                >  2:164408/1‑148 (MQ=255)
                                aTCTCCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACAtc                                                                                                          >  1:302178/1‑149 (MQ=255)
                                     cGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACAtcgtc                                                                                                       <  4:235629/147‑1 (MQ=255)
                                     cGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGAACATcgtcgt                                                                                                     >  8:260274/1‑149 (MQ=255)
                                                             cAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGGATGCGGCCGATTTCGACAATGCCGCCCTTGTGCTCCAGGGCCCTGGTCAGGCCAAAGATATAACGCACATCGTCGTCCACCACCAGCACATTGCGGcc                                                                               >  1:81426/1‑148 (MQ=255)
                                                                   gATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGgcgc                                                                                                                                      <  3:87906/86‑1 (MQ=255)
                                                                   gATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGgcgc                                                                                                                                      >  4:87906/1‑86 (MQ=255)
                                                                            gtGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACGTTGCGGCCCTCGAACACCTTGTcgc                                                              <  7:46042/149‑1 (MQ=255)
                                                                               ctgctCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTa                                                                                                                  <  2:246111/94‑1 (MQ=255)
                                                                               ctgctCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTa                                                                                                                  >  1:246111/1‑94 (MQ=255)
                                                                               ctgctCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGc                                                           >  7:179565/1‑149 (MQ=255)
                                                                               ctgctCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGc                                                           >  7:174843/1‑149 (MQ=255)
                                                                                       ggCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGcgcgcg                                                   <  8:174843/149‑1 (MQ=255)
                                                                                       ggCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGcgcgcg                                                   <  8:179565/149‑1 (MQ=255)
                                                                                         cgctcgctGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCCCATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGt                                                 <  2:81426/149‑1 (MQ=255)
                                                                                           cTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCt                                               <  5:253079/149‑1 (MQ=255)
                                                                                              gATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCtt                                              >  1:216525/1‑147 (MQ=255)
                                                                                                      cgcgCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGc                                      <  1:164408/147‑1 (MQ=255)
                                                                                                      cgcgCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGc                                      <  7:260274/147‑1 (MQ=255)
                                                                                                        cgcCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGCGtt                                   >  3:283038/1‑148 (MQ=255)
                                                                                                                   ggCCGATTTCGACAATGGCGCCCTTGTGCTCCTGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGCGTTGACGTTCGTTg                        <  6:261322/149‑1 (MQ=255)
                                                                                                                   ggCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAacac                                                                      >  4:205816/1‑103 (MQ=255)
                                                                                                                   ggCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAacac                                                                      <  3:205816/103‑1 (MQ=255)
                                                                                                                   ggCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGCGTTGACGTTCGTTg                        >  4:294749/1‑149 (MQ=255)
                                                                                                                   ggCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGCGTTGACGTTCGTTg                        >  7:81553/1‑149 (MQ=255)
                                                                                                                       tccttCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGCGTTGACGTTCGTTGGa                      >  2:240773/4‑147 (MQ=255)
                                                                                                                             gACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGCGTTGACGTTCGTTGGACAGCTGc               >  5:148335/1‑148 (MQ=255)
                                                                                                                             gACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGCGTTGACGTTCGTTGGACAGCTGc               >  8:60310/1‑148 (MQ=255)
                                                                                                                               ctccagcgcCCTGGAGCACAAGGTCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGCGTTGACGTTCGTTGGACAGCTGCGa             <  5:145193/143‑1 (MQ=255)
                                                                                                                                 aTGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGCGTTGACGTTCGTTGGACAGCTGCGATTc          <  2:216525/149‑1 (MQ=255)
                                                                                                                                  cagcgcCCTGGAGCACAAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGCGTTGACGTTCGTTGGACAGCTGCgg             >  6:145193/3‑144 (MQ=255)
                                                                                                                                  cagcgcCCTGGAGCACAAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGCGTTGACGTTCGTTGGACAGCTGCGATTc          >  1:263466/3‑148 (MQ=255)
                                                                                                                                   ggCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGCGTTGACGTTCGTTGGACAGCTGCGATTCGa        >  7:155715/1‑149 (MQ=255)
                                                                                                                                   ggCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGCGTTGACGTTCGTTGGACAGCTGCGATTCGa        >  4:83532/1‑149 (MQ=255)
                                                                                                                                   ggCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGACCACCAGCAGCACCTTGCCGCCCTCGAACACCTTGTCGCGGCTGCGCCCGGTCTTGAGCATGCGTTGACGATCGTTGGACAGCTGCGATACGa        >  3:294875/1‑149 (MQ=255)
                                                                                                                                         cTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGCGTTGACGTTCGTTGGACAGCTGCGATTCGACCTtgt  >  4:7699/1‑149 (MQ=255)
                                                                                                                                           |                                                                                                                                                  
GGCGGATCAGGCGGGTGGCCTCGAAGCCGTCCATCTCCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGCGTTGACGTTCGTTGGACAGCTGCGATTCGACCTTGT  >  NC_002947/4289292‑4289576

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