Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 594,154 A→T 7.9% V206E (GTA→GAA)  cusS ← copper‑sensing histidine kinase in two‑component regulatory system with CusR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913594,1540AT 7.9% 143.7 / 3.1 51V206E (GTA→GAA) cusScopper‑sensing histidine kinase in two‑component regulatory system with CusR
Reads supporting (aligned to +/- strand):  ref base A (22/25);  new base T (2/2);  total (24/27)
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.08e-01

TACGCTCGATCATATGGTTGAACGACAGTACCAGCTGTTCCAGCTCAATGGGCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATATTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTACGATAAAGACGATCAGGATGCTGATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTTATGTAATGAAGATGAAAATCGATCGAAAGCGCGATGTAGAGCGTATAAATCGGTTTGCCGTCCACCAACG  >  NC_000913/594008‑594301
                                                                                                                                                  |                                                                                                                                                   
tACGCTCCATCATATGGTTGATCGACAGTACCAGCTGTTCCTGCTCATTGGGCACGGTCTGCGGGTCGAGGCGAACCTCGAGATCTTTCGAGGTAATATTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTAcc                                                                                                                                                   <  2:4806/149‑1 (MQ=255)
        aTCATATGGTTGAACGACAGTACCAGCTGTTCCAGCTCAATGGGCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATATTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAg                                                                                                                                           <  7:85066/149‑1 (MQ=255)
                 ttGAACGACAGTACCAGCTGTTCCAGCTCAATGGGCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATATTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTACGATa                                                                                                                                    >  7:189727/1‑147 (MQ=255)
                   gAACGACAGTACCAGCTGTTCCAGCTCAATGGGCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATATTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTCCGATaatgc                                                                                                                                >  5:142565/1‑146 (MQ=255)
                      cGACAGTACCAGCTGTTCCAGCTCAATGGGCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATATTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTACGATAAAGACGa                                                                                                                             >  5:184163/1‑149 (MQ=255)
                       gACAGTACCAGCTGTTCCAGCTCAATGGGCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATATTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTACGATAAAGACGAt                                                                                                                            >  4:161233/1‑149 (MQ=255)
                           gTACCAGCTGTTCCAGCTCAATGGGCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATATTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTACGATAAAGACGATCAgg                                                                                                                        <  8:48566/149‑1 (MQ=255)
                              ccAGCTGTTCCAGCTCAATGGGCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATATTATGGATTTGACGGCTGACGCTGCCGATAGGCGCGTGACCTTAATGTACCGCCAACAGTAAGAAAAAGACGATCAGGATg                                                                                                                     >  3:216190/1‑149 (MQ=255)
                                    gTTCCAGCTCAATGGGCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATATTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTACGATAAAGACGATCAGGATGCTGAt                                                                                                                >  7:74853/1‑148 (MQ=255)
                                        cAGCTCAATGGGCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATATTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCCACAGTACGATAAAGACGATCAGGATGCTGATTACCg                                                                                                           >  3:158287/1‑149 (MQ=255)
                                            tCAATGGGCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATATTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTATGATAAAGACGATCAGGATGCTGATTACCGATGc                                                                                                       >  2:11452/1‑149 (MQ=255)
                                                             cGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATATTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTACGATAAAGACGATCAGGATGCTGATTACCGATGCGGTCATAATAAGttt                                                                                        <  7:111990/147‑1 (MQ=255)
                                                                 tCGAGGCGAACGTCGAGATCTTTCGAGGTAATATTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTACGATAAAGACGATCAGGATGCTGATTACCGATGCGGTCATAATAAGTTTATtcat                                                                                  <  5:8540/149‑1 (MQ=255)
                                                                      gCGAACGTCGAGATCTTTCGAGGTAATATTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTACGATAAAGACGATCAGGATGCTGATTACCGATGCGGTCATAATAAGTTTATTCATCAAAt                                                                             >  2:98971/1‑149 (MQ=255)
                                                                      gCGAACGTCGAGATCTTTCGAGGTAATATTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTACGATAAAGACGATCAGGATGCTGATTACCGATGCGGTCATAATAAGTTTATTCATCAAAt                                                                             <  3:161233/149‑1 (MQ=255)
                                                                      gCGAACGTCGAGATCTTTCGAGGTAATATTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTACGATAAAGACGATCAGGATGCTGATTACCGATGCGGTCATAATAAGTTTATTCATCAAAt                                                                             >  4:13462/1‑149 (MQ=255)
                                                                           cGTCGAGATCTTTCGAGGTAATATTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTACGATAAAGACGATCAGGATGCTGATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCAtt                                                                         >  3:152105/1‑148 (MQ=255)
                                                                                 gaTCTTTCGAGGTAATATTCTGGATTTGACGGCTGACGCTGCGGATCTGCGCGTGACCTTTATGTTCCGCCAACAGTACGATAAAGACGATCAGGATGCTGATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTTATTTaa                                                                  >  2:90143/1‑149 (MQ=255)
                                                                                 gaTCTTTCGAGGTAATATTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTTCCGCCAACAGTACGATAAAGACGATCAGGATGCTGATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTTATGTaa                                                                  >  2:93550/1‑149 (MQ=255)
                                                                                 gaTCTTTCGAGGTAATATTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTACGATAAAGACGATCAGGATGCTGa                                                                                                                 >  4:214530/1‑102 (MQ=255)
                                                                                 gaTCTTTCGAGGTAATATTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTACGATAAAGACGATCAGGATGCTGa                                                                                                                 <  3:214530/102‑1 (MQ=255)
                                                                                      ttcgttgTAATATTCTGGATTTGAAGGCTGACGCTGCGGATGGGCGCGTGACCTTTATGTACCGCCAAGAGTACGATAAAGACGATCAGGATGCTGATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTTATGTAATGaa                                                              <  6:232422/142‑1 (MQ=255)
                                                                                       tCGATGTAATATTCTTTATTTTACGGCTGACGCTGCGGATATGCGCGTGACCTTTATGTACCGCCAACAGTACGATAAAGACGATCAGGATGCTGATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTTATGTAATGa                                                               <  6:207036/146‑1 (MQ=255)
                                                                                       tCGAGGTAATATTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTACGATAAAGACGATCAGGATGCTGATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTTATGTAATGa                                                               <  7:79671/146‑1 (MQ=255)
                                                                                           ggTAATATTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTACGATAAAGACGATCAGGATGCTGATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTTATGTAATGAAGATGa                                                         >  8:76564/1‑148 (MQ=255)
                                                                                            gTAATATTCTGGATTTGACTGCTGACGCTGCTGATCGGAGCGTGACCTTTATGTACCGCCAACAGTACGATAAAGACGATCAGGATTCTGATTACTGATGCGGTGATAATAAGTTTATTCATCAAATCATTTATGTAATGAAGATGaa                                                        <  4:216190/148‑1 (MQ=255)
                                                                                               atatTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTACGATAAAGACGATCAGGATGCTGATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTTATGTAATGAAGATGAAAat                                                     <  8:74853/148‑1 (MQ=255)
                                                                                                tatTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTAGCGCCAACAGTACGATAAAGACGATCAGGATGCTGATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTTATGTAATGAAGATGAAAatc                                                    >  2:229698/1‑148 (MQ=255)
                                                                                                  ttCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTACGATAAAGACGATCAGGATGCTGATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTTATGTAATGAAGATGAAAatcga                                                  <  1:98971/148‑1 (MQ=255)
                                                                                                     tGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTACGATAAAGACGATCAGGATGCTGATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTTATGTAATGAAGATGAAAatcgatcg                                               >  3:227402/1‑148 (MQ=255)
                                                                                                       gATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTACGATAAAGACGATCAGGATGCTGATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTTATGTAATGAAGATGAAAATCGATCGaaa                                            <  4:152105/149‑1 (MQ=255)
                                                                                                       gATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTACGATAAAGACGATCAGGATGCTGATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTTATGTAATGAAGATGAAAATCGATCGaaa                                            <  7:92648/149‑1 (MQ=255)
                                                                                                             aCGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTACGATAAAGACGATCAGGATGCTGATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTTATGTAATGAAGATGAAAATCGATCGAAAGCGCGa                                      >  2:96524/1‑149 (MQ=255)
                                                                                                                gCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTACGATAAAGACGATCAGGATGCTGATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTTATGTAATGAAGATGAAAATCGATCGAAAGCGCGATGt                                   <  1:96524/149‑1 (MQ=255)
                                                                                                                 cggACGCTGCGTAGGGGCGCGTGAACTTTATGTACCTCCAACAGTACGATAAAGACGATCAGGAATCTGATTACCTATTCGGTCATAATAAGTTTATTCATCAAATCATTTATGTAATTAAGATGAAAATCGATCGAAAGCGCTATGTa                                  <  6:59915/147‑1 (MQ=255)
                                                                                                                 cTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTACGATAAAGACGATCAGGATGCTGATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTTATGTAATGAAGATGAAAATCGATCGAAAGCGCGATg                                    >  3:201049/1‑147 (MQ=255)
                                                                                                                  tGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTACGATAAAGACGATCAGGATGCTGATTTCCGATGCGGTCATAATAAGTTTATTCATCAAATCATTTATGTAATGAAGATGAAAATCGATCGAAAGCGCGATg                                    >  7:7500/1‑146 (MQ=255)
                                                                                                                  tGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTACGATAAAGACGATCAGGATGCTGATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTTATGTAATGAAGATGAAAATCGATCGAAAGCGCGATg                                    >  3:180414/1‑146 (MQ=255)
                                                                                                                          cGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTACGATAAAGACGATCAGGATGCTGATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTTATGTAATGAAGATGAAAATCGATCGAAAGCGCGATGTAGAGCGtat                          <  6:184163/148‑1 (MQ=255)
                                                                                                                           ggATCGGCGCGTGACCTTTATGTACCGCCAACAGTACGATAAAGACGATCAGCATGCTGATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTTATGTAATGAAGATGAAAATCGATCGAAAGCGCGATGTAGAGCGTATaa                        <  2:218567/149‑1 (MQ=255)
                                                                                                                              tCGGCGCGTGACCTTTATGTACCGCCAACAGTACGATAAAGACGATCAGGATGCTGATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTTATGTAATGAAGATGAAAATCGATCGAAAGCGCGATGTAGAGCGTATAAATc                     <  8:189727/149‑1 (MQ=255)
                                                                                                                                 gcgcGTGACCTTTATGTACCGCCAACAGTACGATAAAGACGATCAGGATGCTGATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTTATGTAATGAAGATGAAAATCGATCGAAAGCGCGATGTAGAGCGTATAAATCGGt                  <  4:201049/149‑1 (MQ=255)
                                                                                                                                          cTTTATGTACCGCCAACAGTACGATAAAGACGATCAGGATGCTGATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTTATGTAATGAAGATGAAAATCGATCGAAAGCGCGATGTAGAGCGTATAAATCGGTTTGCCGTcc         >  5:97820/1‑149 (MQ=255)
                                                                                                                                          cTTTATGTACCGCCAACAGTACGATAAAGACGATCAGGATGCTGATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTTATGTAATGAAGATGAAAATCGATCGAAAGCGCGATGTAGAGCGTATAAATCGGTTTGCCGTc          >  7:175300/1‑148 (MQ=255)
                                                                                                                                          cTTTATGTACCGCCAACAGTACGATAAAGACGATCAGGATGCTGATTACCGATGCGGTCATAATAAGTTTATTCATAAAATCATTTATGTAATGAAGATGAAAATCGATCGAAAGCGCGATGTAGAGCGTATAAATCGGTTTGCCGTc          >  5:47001/1‑148 (MQ=255)
                                                                                                                                             tATGTTCCGCCAACAGTACGATAAAGACGATCAGGATGCTGATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTTATGTAATGAAGATGAAAATCGATCGAAAGCGCGATGTAGAGCGTATAAACCGGTTTGCCGTccac       <  1:90143/148‑1 (MQ=255)
                                                                                                                                             tATGTTCCGCCAACAGTACGATAAAGACGATCAGGATGCTGATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTTATGTAATGAAGATGAAAATCGATCGAAAGCGCGATGTAGAGCGTATAAACCGGTTTGCCGTccac       <  1:93550/148‑1 (MQ=255)
                                                                                                                                             tATGTACCGCCAACAGTATGATAAAGACGATCAGGATGCTGATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTTATGTAATGAAGATGAAAATCGATCGAAAGCGCGATGTAGAGCGTATAAATCGGTTTGCCGTccac       <  1:11452/148‑1 (MQ=255)
                                                                                                                                             tATGTACCGCCAACAGTACGATAAAGACGATCAGGATGCTGATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTTATGTAATGAAGATGAAAATCGATCGAAAGCGCGATGTAGAGCGTATAAATCGGTTTGCCGTcca        <  2:35594/147‑1 (MQ=255)
                                                                                                                                                 tACCGCCAACAGTACGATAAAGACGATCAGGATGCTGATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTTATGTAATGAAGATGAAAATCGATCGAAAGCGCGATGTAGAGCGTATAAATCGGTTTGCCGTCCACCAACg  <  6:142565/149‑1 (MQ=255)
                                                                                                                                                 tACCGCCAACAGTACGATAAAGACGATCAGGATGCTGATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTTATGTAATGAAGATGAAAATCGATCGAAAGCGCGATGTAGAGCGTATAAATCGGTTTGCCGTCCACCAACg  <  4:227402/149‑1 (MQ=255)
                                                                                                                                                  aCCGCCAACAGTACGATAAAGACGATCAGGATGCTGATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTTATGTaa                                                                  <  4:90320/84‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TACGCTCGATCATATGGTTGAACGACAGTACCAGCTGTTCCAGCTCAATGGGCACGGTCTGCGGGTCGAGGCGAACGTCGAGATCTTTCGAGGTAATATTCTGGATTTGACGGCTGACGCTGCGGATCGGCGCGTGACCTTTATGTACCGCCAACAGTACGATAAAGACGATCAGGATGCTGATTACCGATGCGGTCATAATAAGTTTATTCATCAAATCATTTATGTAATGAAGATGAAAATCGATCGAAAGCGCGATGTAGAGCGTATAAATCGGTTTGCCGTCCACCAACG  >  NC_000913/594008‑594301

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: