Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,486,807 C→T 9.9% C812C (TGC→TGT)  evgS → hybrid sensory histidine kinase in two‑component regulatory system with EvgA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,486,8070CT 9.9% 107.1 / 2.9 41C812C (TGC→TGT) evgShybrid sensory histidine kinase in two‑component regulatory system with EvgA
Reads supporting (aligned to +/- strand):  ref base C (20/16);  new base T (2/2);  total (22/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.51e-01

AATTGGTGAAATCCTTGATGTCGACAAAATTGAATCGGGTAACTATCAACTTCAACCACAATGGGTCGATATCCCTACTTTAGTCCAGAACACTTGTCACTCTTTCGGTGCGATTGCTGCAAGCAAATCGATCGCATTAAGTTGCAGCAGTACGTTTCCTGAACATTACCTGGTTAAGATCGACCCTCAGGCGTTTAAGCAGGTCTTATCAAATTTACTGAGTAATGCTCTCAAATTTACCACCGAGGGGGCAGTAAAAATTACGACCTCCCTGGGTCACATTGATG  >  NC_000913/2486663‑2486949
                                                                                                                                                |                                                                                                                                              
aaTTGGTGAAATCCTTGATGTCGACAAAATTGAATCGGGTAACTATCAACTTCAACCACAATGGGTCGATATCCCTACTTTAGTCCAGAACACTTGTCACTCTTTCGGTGCGATTGCTGCAAGCAAATCGATCGCATTAAGTTgcagca                                                                                                                                            >  2:733366/1‑149 (MQ=255)
   tGGTGAAATCCTTGATGTCGACAAAATTGAATCGGGTAACTATCAACTTCAACCACAATGGGTCGATATCCCTACTTTAGTCCAGAACACTTGTCACTCTTTCGGTGCGATTGCTGCAAGCAAATCGATCGCATTAAGTTGTAGCAGTa                                                                                                                                         >  1:36790/1‑149 (MQ=255)
   tGGTGAAATCCTTGATGTCGACAAAATTGAATCGGGTAACTATCAACTTCAACCACAATGGGTCGATATCCCTACTTTAGTCCAGAACACTTGTCACTCTTTCGGTGCGATTGCTGCAAGCAAATCGATCGCATTAAGTTGTAGCAGTa                                                                                                                                         >  1:36794/1‑149 (MQ=255)
    ggTGAAATCCTTGATGTCGACAAAATTGAATCGGGTAACTATCAACTTCAACCACAATGGGTCGACATCCCTACTTTAGTCCAGAACACTTGTCACTCTTTCGGTGCGATTGCTGCAAGCAAATCGATCGCATTAAGTTGCAGCAGTAc                                                                                                                                        <  1:51307/149‑1 (MQ=255)
    ggTGAAATCCTTGACGTCGACAAAAGTGAATCGTGTAACTCTCAACTTCACCCACAATGGGTCGCTCTCCCTACTTGAGTCCAGAACACTTGTCACTCTTTCGGTGCTATTGCTGCAAGCAAATCGATCGCATTAAGTTGCAGCAGTAc                                                                                                                                        <  2:613466/149‑1 (MQ=255)
         aaTCCTTGATGTCGACAAAATTGAATCGGGTAACTATCAACTTCAACCACAATGGGTCGATATCCCTACTTTAGTCCAGAACACTTGTCACTCTTTCGGTGCGATTGCTGCAAGCAAATCGATCGCATTAAGTTGCAGCAGTACGTTTc                                                                                                                                   <  1:488560/149‑1 (MQ=255)
           tCCTTGATGTCGACAAAATTGAATCGGGCAACTATCAACTTCAACCACACTGGGTCGATATCCCTACTTTAGGCCAGAACACTTGTCACTCTTTCGGTGCGATTGCTGCAAGCAAATCGATCGCATTAAGTTGCAGCAGTACGTTTCCt                                                                                                                                 <  1:7190/149‑1 (MQ=255)
               tGATGTCGACAAAATTGAATCGGGTAACTATCAACTTCAACCACAATGGGTCGATATGCCTACTTTAGTCCAGAACACTTGTCACTCTTTCGGTGCGATTGCTGGAAGCAAATCGATCGCATTCAGTTGCCGCAGTACGTTTCCtacac                                                                                                                             >  2:698781/1‑145 (MQ=255)
               tGATGTCGACAAAATTGAATCGGGTAACTATCAACTTCAACCACAATGGGTCGATATCCCTACTTTAGTCCAGAACACTTGTCACTCTTTCGGTGCGATTGCTGCAAGCAAATCGATCGCATAAAGTTGCAGCAGGACGTTTCCTGAAc                                                                                                                             >  2:455751/1‑149 (MQ=255)
                   gTCGACCAAATTGAATCGGGTAACTATCAACTTCAACCACAATGGGTCGATATCCCTACTTTAGTCCCGAACCCTTGTCACTCTTTCGGTGCGATTGCTGCAGGCAAATAGATCGCATTAAGTTGCAGCCGTGCGTTTACCGAACAtac                                                                                                                         >  2:276756/1‑147 (MQ=255)
                    tCGACAAAATTGAATCGGGTAACTATCAACTTCAACCACAATGGGTCGATATCCCTACTCTAGTCCAGAACACTTGTCACTCTTTCGGTGCGATTGCTGCAAGCAAATCGATCGCATTAAGTTGCAGCAGTACGTTTCCTGAACATTAc                                                                                                                        >  1:532027/1‑149 (MQ=255)
                             ttGAATCGGGTAACTATCAACTTCAACCACAATGGGTCGATATCCCTACTTTAGTCCAGAACACTTGTCACTCTTTCGGTGCGATTGCTGCAAGCAAATCGATCGCATTAAGTTGCAGCAGTACGTTTCCTGAACATTACCTGGTTAAg                                                                                                               >  1:60676/1‑149 (MQ=255)
                               gAATCGGGTAACTATCAACTTCAACCACAATGGGTCGATATCCCTACTTTAGTCCAGAACACTTGTCACTCTTTCGGTGCGATTGCTGCAAGCAAATCGATCGCATTAAGTTGCAGCAGTACGTTTCCTGAACATTACCTGGTTAAGAt                                                                                                             <  1:35922/149‑1 (MQ=255)
                                           tATCAACTTCAACCACAATGGGTCGATATCCCTACTTTAGTCCAGAACACTTGTCACTCTTTCGGTGCGATTGCTGCAAGCAAATCGATCGCAGTCACTTAAAACGGTACGTTTCCTGAACATTACCTGGTTAAGATCGACCCTCAGGc                                                                                                 <  1:790553/149‑1 (MQ=255)
                                              cAACTTCAACCACAATGGGTCGATATCCCTACTTTAGTCCAGAACACTTGTCACTCTTTCGGTGCGATTGCTGCAAGCAAATCGATCGCATTAAGTTGCAGCAGTACGTTTCCTGAACATTACCTGGTTAAGATCGACCCTCAGGCGtt                                                                                              <  1:644058/149‑1 (MQ=255)
                                                  ttCAACCACAATGGGTCGATATCCCTACTTTAGTCCAGAACACTTGTCACTCTTTCGGTGCGATTGCTGCAAGCAAATCGATCGCATTAAGTTGCAGCAGTACGTTTCCTGAACATTACCTGGTTAAGATCGACCCTCAGGCGTTTAAg                                                                                          >  1:36072/1‑149 (MQ=255)
                                                  ttCAACCACAATGGGTCGATATCCCTACTTTAGTCCAGAACACTTGTCACTCTTTCGGTGCGATTGCTGCAAGCAAATCGATCGCATTAAGTTGCAGCAGTACGTTTCCTGAACATTACCTGGTTAAGATCGACCCTCAGGCGTTTAAg                                                                                          >  2:546040/1‑149 (MQ=255)
                                                      aCCACAATGGGTCGATATCCCTACTTTAGTCCAGAACACTTGTCACTCTTTCGGTGCGATTGCTGCAAGCAAATCGATCGCATTAAGTTGCAGCAGTACGTTTCCTGAACATTACCTGGTTAAGATCGACCCTCAGGCGTTTAAGCAgg                                                                                      <  2:587193/149‑1 (MQ=255)
                                                       ccACAATGGGTCGATATCCCTACTTTAGTCCAGAACACTTGTCACTCTTTCGGTGCGATTGCTGCAAGCAAATCGATCGCATTAAGTTGCAGCAGTACGTTTCCTGAACATTACCTGGTTAAGATCGACCCTCAGGCGTTTAAGCAGGt                                                                                     >  2:67780/1‑149 (MQ=255)
                                                        cacaATGGGTCGATATCCCTACTTTAGTCCAGAACACTTGTCACTCTTTCGGTGCGATTGCTGCAAGCAAATCGATCGCATTAAGTTGCAGCAGTACGTTTCCTGAACATTACCTGGTTAAGATCGACCCACAGGCGTTTAAGCAGGTc                                                                                    >  2:20723/1‑149 (MQ=255)
                                                                 tCGATATCCCTACTTTAGTCCAGAACACTTGTCACTCTTTCGGTGCGATTGCTGCAAGCAAATCGATCGCATTAAGTTGCAGCAGTACGTTTCCTGAACATTACCTGGTTAAGATCGACCCTCAGGCGTTTAAGCAGGTCTTATCAAAt                                                                           <  1:230068/149‑1 (MQ=255)
                                                                  cGATATCCCTACTTTAGTCCAGAACACTTGTCACTCTTTCGGTGCGATTGCTGCAAGCAAATCGATCGCATTAAGTTGCAGCAGTACGTTTCCTGAACATTACCTGGTTAAGATCGACCCTCAGGCGTTTAAGCAGGTCTTATCAAAtt                                                                          <  1:165884/149‑1 (MQ=255)
                                                                      atCCCTACTTTAGTCCAGAACACTTGTCACTCTTTCGGTGCGATTGCTGCAAGCAAATCGATCGCATTAAGTTGCAGCAGTACGTTTCCTGAACATTACCTGGTTAAGATCGACCCTCAGGCGTTTAAGCAGGTCTTATCAAATTTACt                                                                      <  2:743878/149‑1 (MQ=255)
                                                                      atCCCTACTTTAGTCCAGAACACTTGTCACTCTTTCGGTGCGATTGCTGCAAGCAAATCGATCGCATTAAGTTGCAGCAGTACGTTTCCTGAACATTACCTGGTTAAGATCGACCCTCAGGCGTTCAAGCAGGTCTTATCAAATTTACt                                                                      >  2:247310/1‑149 (MQ=255)
                                                                      atCCCTACTTTAGTCCAGAACACTTGTCACTCTTTCGGTGCGATTGCTGCAAGCAAATCGATCGCATTAAGTTGCAGCAGTACGTTTCCTGAACATTACCTGGTTAAGATCGACCCCCAGGCGTTTAAGCAGGTCTTATCAAATTTGCt                                                                      >  2:47423/1‑149 (MQ=255)
                                                                        cccTACTTTAGTCCAGAACACTTGTCACTCTTTCGGTGCGATTGCTGCAAGCAAATCGATCGCATTAAGTTGCAGCAGTACGTTTCCTGAACATTACCTGGTTAAGATCGACCCTCAGGCGTTTAAGCAGGCCTTATCAAATTTACTGa                                                                    >  2:13666/1‑149 (MQ=255)
                                                                             cTTTAGTCCAGAACACTTGTCACTCTTTCGGTGCGATTGCTGCAAGCAAATCGATCGCATTAAGTTGCAGCAGTACGTTTCCTGAACAATACCTGGTTAAGATCGACCCTCAGGCGTTTAAACAGGGCTTATCAAATTTACTGAGTAAt                                                               >  2:709829/1‑149 (MQ=255)
                                                                                tAGTCCAGAACACTTGTCACTCTTTCGGTGCGATTGCTGCAAGCAAATCGATCGCATTAAGTTGCAGCAGTACGTTTCCTGAACATTACCTGGTTAAGATCGACCCTCAGGCGTTTAAGCAGGTCTTATCAAATTTACTGAGTAATGct                                                            >  2:31792/1‑149 (MQ=255)
                                                                                   tCCAGAACATTTGTCACTTTTTCGATGCGATTGCTGCAAGCAAATTGATCGCATTAAGTTGCAGCAGTACGTTTCCTGAACATTACCTGGTTAAGATCGACCCTCAGGCGTTTAAGCAGGTCTTATCAAATTTACTGAGTAATGCtctc                                                         <  2:593178/149‑1 (MQ=255)
                                                                                       gAACACTTGTCACTCTTTCGGTGCGATTGCTGCAAGCAAATCGATCGCATTAAGTTGCAGCAGTACGTTTCCTGAACATTACCTGGTTAAGATCGACCCTCAGGCGTTTAAGCAGGTCTTATCAAATTTACTGAGTAATGCTCTCAAAt                                                     <  2:398803/149‑1 (MQ=255)
                                                                                       gAACACTTGTCACTCTTTCGGTGCGATTGCTGCAAGCAAATCGATCGCATTAAGTTGCAGCAGTACGTTTCCTGAACATTACCTGGTTAAGATCGACCCTCAGGCGTTTAAGCAGGTCTTATCAAATTTACTGAGTAATGCTCTCAAAt                                                     <  2:574834/149‑1 (MQ=255)
                                                                                                tCACTCTTTCGGTGCGATTGCTGCAAGCAAATCGATCGCATTAAGTTGCAGCAGTACGTTTCCTGAACATTACCTGGTTAAGATCGACCCCCAGGCGTTTAAGCAGGTCTTATCAAATTTACTGAGTAATGCCCTCAAATTTACCACCg                                            >  2:460611/1‑149 (MQ=255)
                                                                                                      tttCGGTGCGATTGCTGCAAGCAAATCGATCGCATTAAGTTGTAGCAGTACGTTTCCTGAACATTACCTGGTTAAGATCGACCCTCAGGCGTTTAAGCAGGTCTTATCAAATTTACTGAGTAATGCTCTCAAATTTACCACCGAggggg                                      <  2:36790/149‑1 (MQ=255)
                                                                                                      tttCGGTGCGATTGCTGCAAGCAAATCGATCGCATTAAGTTGTAGCAGTACGTTTCCTGAACATTACCTGGTTAAGATCGACCCTCAGGCGTTTAAGCAGGTCTTATCAAATTTACTGAGTAATGCTCTCAAATTTACCACCGAggggg                                      <  2:36794/149‑1 (MQ=255)
                                                                                                        tCGGTGCGATTGCTGCAAGCAAATGGATCGCATTAAGTTGCAGCAGTACGTTTCCTGAACATTACCTGGTTAAGATCGACCCTCAGGCGTTTAAGCAGGTCTTATCAAATTTACTGAGTAATGCTCTCAAATTTACCACCGAGGGGGCa                                    <  2:17695/149‑1 (MQ=255)
                                                                                                                 ttGCTGCAAGCAAATCGATCGCATTAAGTTGCAGCAGTACGTTTCCTGAACATTACCTGGTTAAGATCGACCCTCAGGCGTTTAAGCAGGTCTTATCAAATTTACTGAGTAATGCTCTCAAATTAACCACCGAGGGGGCAGTAAAAAtt                           >  2:716264/1‑149 (MQ=255)
                                                                                                                  tgctgcAAGCAAATCGATCGCATTAAGTTGCAGCAGTACGTTTCCTGAACATTACCTGGTTAAGATCGACCCTCAGGCGTTTAAGCAGGTCTTATCAAATTTACTGAGTAATGCTCTCAAATTTACCACCGAGGGGGCAGTAAAAATTa                          <  2:277026/149‑1 (MQ=255)
                                                                                                                     tgcAAGCAAATCGATCGCATTAAGTTGCAGCAGTACGTTTCCTGAACATTACCTGGTTAAGATCGACCCTCAGGCGTTTAAGCAGGTCTTATCAAATTTACTGAGTAATGCTCTCAAATTTACCACCGAGGGGGCAGTAAAAATTACGa                       >  1:797869/1‑149 (MQ=255)
                                                                                                                       caagcaaATCGATCGCATTAAGTTGCAGCAGTACGTTTCCTGAACATTACCTGGTTAAGATCGACCCTCAGGCGTTTAAGCAGGTCTTATCAAATTTACTGAGTAATGCTCTCAAATTTACCACCGAGGGGGCAGTAAAAATTACGAcc                     >  1:701808/1‑149 (MQ=255)
                                                                                                                                 gatcgCATTAAGTTGCAGCAGTACGTTTCCTGAACATTACCTGGTTAAGATCGACCCTCAGGCGTTTAAGCAGGTCTTATCAAATTTACTGAGTAATGCTCTCAAATTTACCACCGACGGGGCCGGAAAAATTACGACCTCCCTGGGTc           >  2:129795/1‑149 (MQ=255)
                                                                                                                                         tAAGTTGCAGCAGTACGTTTCCTGAACATTACCTGGTTAAGATCGACCCTCAGGCGTTTAAGCAGGTCTTATCAAATTTACTGAGTAATGCTCTCAAATTTACCACCGAGGGGGCCGTAAAAATTACGACCTCCCTGGGTCACATtgat   >  1:19349/1‑149 (MQ=255)
                                                                                                                                          aaGTTGCAGCAGTACGTTTCCTGAACATTACCTGGCTAAGATCGACCCTCAGGCGTTTAAGCAGGTCTTATCAAATTTACTGAGTAATGCTCTCAAATTTACCACCGAGGGGGCAGTAAAAATTACGACCTCCCTGGGTCACATtgatg  <  1:457593/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                              
AATTGGTGAAATCCTTGATGTCGACAAAATTGAATCGGGTAACTATCAACTTCAACCACAATGGGTCGATATCCCTACTTTAGTCCAGAACACTTGTCACTCTTTCGGTGCGATTGCTGCAAGCAAATCGATCGCATTAAGTTGCAGCAGTACGTTTCCTGAACATTACCTGGTTAAGATCGACCCTCAGGCGTTTAAGCAGGTCTTATCAAATTTACTGAGTAATGCTCTCAAATTTACCACCGAGGGGGCAGTAAAAATTACGACCTCCCTGGGTCACATTGATG  >  NC_000913/2486663‑2486949

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