Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,636,827 C→T D44D (GAC→GAT)  creC → sensory histidine kinase CreC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,636,8270CT96.7% 89.5 / ‑4.4 31D44D (GAC→GAT) creCsensory histidine kinase CreC
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (16/13);  total (18/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.30e-01

CACCGTCGATACCCACATTAAAACGCTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACAC‑CGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCT  >  NC_000913/4636585‑4637069
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                   
cACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATAC‑CGCaa                                                                                                                                                                                                                                             >  2:144197/1‑250 (MQ=255)
    gttgATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGGGGGGACGTTAATCGAAAC‑CGCAACGTTg                                                                                                                                                                                                                                        >  2:169736/4‑251 (MQ=255)
               caTTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGAAGTTAATCGATAC‑CGCAACGTTGATGGCGGAGCt                                                                                                                                                                                                                             >  1:82798/1‑251 (MQ=255)
                          tGCGTGCCAAGCGGCGGGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAATTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATAC‑CGCAACGTTGCTGGCGGAGCTGGCGCGTccc                                                                                                                                                                                                                   <  2:216935/250‑1 (MQ=255)
                               gCCAAGCTGCGCGCCATCAACCCCTATCTTTCACCGATTAATACTCATCCCGGCATGGGATATAGCCGGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGTCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATAC‑CGCAACGTTGCTGGCGGATCTGGCGCGGCCCGAttt                                                                                                                                                                                                              <  1:63392/250‑1 (MQ=255)
                                            ccATCAGCCCCGATATTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGAGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGCAGTTAAACCGGGTGTGCGAAGGGCACAGGAGGGGGCGTTAATCG‑CTCACGCACCGTTGCTGGCGGTGCCGGCGCGCACCGGTTTGCTCTCTGGGGAcc                                                                                                                                                                                                >  2:15057/1‑251 (MQ=255)
                                                             tCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATAC‑CGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTg                                                                                                                                                                               >  1:31497/1‑251 (MQ=255)
                                                              cACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATAC‑CGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTgg                                                                                                                                                                              <  1:65066/251‑1 (MQ=255)
                                                                       aTACTCATCGCGGCATGGGATATAGCCTGAGGGGACTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTAATGGTGGCGGTGGCGGCCTGGCTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGCCGTTAATCGATCC‑CGAAACGTTGCTGGCGGAGCTGGAGCGTCACGATTTGATCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGt                                                                                                                                                                     >  1:267396/1‑251 (MQ=255)
                                                                                       gggATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATAC‑CGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTacaac                                                                                                                                                      <  2:31497/250‑1 (MQ=255)
                                                                                                        gCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATAC‑CGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCa                                                                                                                                    <  1:169736/251‑1 (MQ=255)
                                                                                                                               ggTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGCAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATAC‑CGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAAATGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAAAAAAGTg                                                                                                             >  1:164339/1‑251 (MQ=255)
                                                                                                                                   gctgCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATAC‑CGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTgcgc                                                                                                          <  1:126389/248‑1 (MQ=255)
                                                                                                                                      gCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATAC‑CGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAAt                                                                                                       <  1:246974/250‑1 (MQ=255)
                                                                                                                                           gCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATAC‑CGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAatat                                                                                                 <  2:189116/251‑1 (MQ=255)
                                                                                                                                                   ttgcTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTGTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATAC‑CGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCt                                                                                          <  2:82798/247‑1 (MQ=255)
                                                                                                                                                     aCTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATAA‑CGCAACGTTGCTGGCGGCGCTGGCGCGTACCGGTTTGCTCTCTGGGGAACCAACGCATGGGCAACTGGCGCTAGCGTTTAATCAGCTACAACATCACCCGTTTCTCGCCAATATCGGTGGCCTTAACAAAGTGCGCAATGAATAACAAGTCtat                                                                                        >  1:136695/1‑250 (MQ=255)
                                                                                                                                                           ggcggtggcGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATAC‑CGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCCGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCg                                                                                 >  1:194477/1‑251 (MQ=255)
                                                                                                                                                           ggcggtggcGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATAC‑CGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCg                                                                                 <  1:144197/251‑1 (MQ=255)
                                                                                                                                                              ggtggcGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATAC‑CGCAACGTTGCTGGCGGCGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGAtt                                                                              >  2:24358/1‑250 (MQ=255)
                                                                                                                                                              ggcggcggCCTTGTTCGTACTGGATCTTTTTTTCAGAGTAGGTAAACCGGGCGTGCGCCGAGCACCGGAGGGGCCGTTACTCGATAC‑CGCAACGTTGCTGGCGGCGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACCACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACACAGTGCGCCCTGAATCTCATGTCTATATGAcccatt                                                                              >  2:29484/8‑247 (MQ=255)
                                                                                                                                                                    ggCCTGGTTCGTACTGGCTCTTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATAC‑CGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTTGCATTCACCAAGTGCGCAATGAATATCCTGTCTATATGACCGCTTCGCAgg                                                                        >  2:227799/2‑251 (MQ=255)
                                                                                                                                                                        tGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATAGATAC‑CGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAgggta                                                                     >  1:108669/1‑248 (MQ=255)
                                                                                                                                                                             cGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATAC‑CGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCAGGTCTATATGACCGATTCGCAGGGTAAAGTAt                                                               <  2:44931/251‑1 (MQ=255)
                                                                                                                                                                                                  agaagTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATAC‑CGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAg                                          <  1:24358/251‑1 (MQ=255)
                                                                                                                                                                                                      gTTAAACCGGGAGTGCGAAGAGCAACGGAGGGGACGTTAAGCGATAC‑CGCAACGTTGCTGGCGGAGCTGGCGTGTCCCGATTTGCTCTCTGGGTACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGt                                      <  2:77076/251‑1 (MQ=255)
                                                                                                                                                                                                                gCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATAC‑CGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCCTTAACAAAGGGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGat                            >  1:149215/1‑251 (MQ=255)
                                                                                                                                                                                                                  gTGCGAAGAGCAACGGAGGGGACGTTAATCGATAC‑CGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGatt                           >  1:207078/1‑250 (MQ=255)
                                                                                                                                                                                                                           gCAACGGAGGGGACGTTAATCGATAC‑CGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCCTGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGCTTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTg                 >  2:28596/1‑251 (MQ=255)
                                                                                                                                                                                                                             aaCGGAGGGGACGTTAATCGATAC‑CGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTgg                >  1:73579/1‑250 (MQ=255)
                                                                                                                                                                                                                                          ttaaTCGATAC‑CGCAACGATGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTAGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGAACGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCt  >  1:52499/4‑251 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                   
CACCGTCGATACCCACATTAAAACGCTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACAC‑CGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCT  >  NC_000913/4636585‑4637069

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