Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,289,435 T→A 8.7% E1051V (GAG→GTG)  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,4350TA 8.7% 128.0 / 3.2 46E1051V (GAG→GTG) PP_3761two‑component system sensor histidine kinase/response regulator
Reads supporting (aligned to +/- strand):  ref base T (25/17);  new base A (2/2);  total (27/19)
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.02e-01

GCGGATCAGGCGGGTGGCCTCGAAGCCGTCCATCTCCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCC‑GTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGCGTTGACGTTCGTTGGACAGCTGCGATTCGACCTTGTGCAGGAA  >  NC_002947/4289293‑4289583
                                                                                                                                               |                                                                                                                                                    
gCGGATCAGGCGGGTGGCCTCGAAGCCGTCCATCTCCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCC‑GTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAgg                                                                                                                                                 <  4:20507/148‑1 (MQ=255)
  ggATCAGGCGGGTGGCCTCGAAGCCGTCCATCTCCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCC‑GTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGgcg                                                                                                                                              >  8:153373/1‑149 (MQ=255)
       aGGCGGGTGGCCTCGAAGCCGTCCATCTCCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCC‑GTTGCGGCCGATTTCGACAATGGCGCCCTGGAGCACAAGGgcgccatt                                                                                                                                         <  6:320019/149‑5 (MQ=255)
            ggTGGCCTCGAAGCCGTCCATCTCCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCC‑GTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAggg                                                                                                                                    <  7:153373/149‑1 (MQ=255)
                  cTCGAAGCCGTCCATCTCCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCC‑GTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCa                                                                                                                                  >  3:235629/1‑145 (MQ=255)
                   tCGAAGCCGTCCATCTCCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCC‑GTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGa                                                                                                                              >  5:123111/1‑148 (MQ=255)
                          cGTCCATCTCCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCC‑GTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACg                                                                                                                       >  2:164408/1‑148 (MQ=255)
                               aTCTCCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCC‑GTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACAtc                                                                                                                 >  1:302178/1‑149 (MQ=255)
                                    cGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCC‑GTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACAtcgtc                                                                                                              <  4:235629/147‑1 (MQ=255)
                                    cGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCC‑GTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGAACATcgtcgt                                                                                                            >  8:260274/1‑149 (MQ=255)
                                                            cAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGGATGCGGCCGATTTCGACAATGCCGCCCTTGTGCTCCAGGGCCCTGGTCAGGCCAAAGATATAACGCACATCGTCGTCCACCACCAGCACATTGCGGcc                                                                                      >  1:81426/1‑148 (MQ=255)
                                                                  gATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCC‑GTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGgcgc                                                                                                                                             >  4:87906/1‑86 (MQ=255)
                                                                  gATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCC‑GTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGgcgc                                                                                                                                             <  3:87906/86‑1 (MQ=255)
                                                                           gtGCTGCTCCAGGCGCTCGATGGCTTCGCGCCC‑GTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACGTTGCGGCCCTCGAACACCTTGTcgc                                                                     <  7:46042/149‑1 (MQ=255)
                                                                              ctgctCCAGGCGCTCGATGGCTTCGCGCCC‑GTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTa                                                                                                                         <  2:246111/94‑1 (MQ=255)
                                                                              ctgctCCAGGCGCTCGATGGCTTCGCGCCC‑GTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTa                                                                                                                         >  1:246111/1‑94 (MQ=255)
                                                                              ctgctCCAGGCGCTCGATGGCTTCGCGCCC‑GTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGc                                                                  >  7:179565/1‑149 (MQ=255)
                                                                              ctgctCCAGGCGCTCGATGGCTTCGCGCCC‑GTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGc                                                                  >  7:174843/1‑149 (MQ=255)
                                                                                      ggCGCTCGATGGCTTCGCGCCC‑GTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGcgcgcg                                                          <  8:174843/149‑1 (MQ=255)
                                                                                      ggCGCTCGATGGCTTCGCGCCC‑GTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGcgcgcg                                                          <  8:179565/149‑1 (MQ=255)
                                                                                        cgctcgctGGCTTCGCGCCC‑GTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCCCATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGt                                                        <  2:81426/149‑1 (MQ=255)
                                                                                          cTCGATGGCTTCGCGCCC‑GTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCt                                                      <  5:253079/149‑1 (MQ=255)
                                                                                             gATGGCTTCGCGCCC‑GTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCtt                                                     >  1:216525/1‑147 (MQ=255)
                                                                                                     cgcgCCC‑GTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGc                                             <  7:260274/147‑1 (MQ=255)
                                                                                                     cgcgCCC‑GTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGc                                             <  1:164408/147‑1 (MQ=255)
                                                                                                       cgcCC‑GTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGCGtt                                          >  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)
                                                                                                                                           gtgCTCCAGGGCGCTGGTCAGTGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGCGTTGACGTTCGTTGGACAGCTGCGATTCGACCTTGTGc       <  2:302178/148‑1 (MQ=255)
                                                                                                                                           gtgCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGCGTTGACGTTCGTTGGACAGCTGCGATTCGACCTtgtg        >  5:197462/1‑147 (MQ=255)
                                                                                                                                           gtgCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGCGTTGACGTTCGTTGGACAGCTGCGATTCGACCTtgtg        <  8:81553/147‑1 (MQ=255)
                                                                                                                                              cTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGCGTTGACGTTCGTTGGACAGCTGCGATTCGACCTTGTGCAgg    >  7:189940/1‑148 (MQ=255)
                                                                                                                                               tCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTGTTGAGCATGCGTTGACGTTCGTTGGACAGCTGCGATTCGACCTTGTGCAGGaa  >  6:187566/1‑149 (MQ=255)
                                                                                                                                               tCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGCGTTGACGTTCGTTGGACAGCTGCGATTCGACCTTGTGCAGGaa  >  3:280582/1‑149 (MQ=255)
                                                                                                                                               tCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGCGTTGACGTTCGTTCGACAGCTGCGAtt                  <  6:307419/133‑1 (MQ=255)
                                                                                                                                               tCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGCGTTGACGTTCGTTCGACAGCTGCGAtt                  >  5:307419/1‑133 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
GCGGATCAGGCGGGTGGCCTCGAAGCCGTCCATCTCCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCC‑GTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGCGTTGACGTTCGTTGGACAGCTGCGATTCGACCTTGTGCAGGAA  >  NC_002947/4289293‑4289583

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