Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 593,467 T→G D435A (GAC→GCC)  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_000913593,4670TG100.0% 99.6 / NA 27D435A (GAC→GCC) cusScopper‑sensing histidine kinase in two‑component regulatory system with CusR
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (10/17);  total (10/17)

GCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGG  >  NC_000913/593332‑593593
                                                                                                                                       |                                                                                                                              
gCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCaa                                                                                                                             <  1:342426/139‑1 (MQ=255)
    gTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAAcgcg                                                                                                                         >  1:55105/1‑139 (MQ=255)
           gATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGaaa                                                                                                                  <  1:464463/139‑1 (MQ=255)
              aaCAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACgg                                                                                                               >  2:234689/1‑139 (MQ=255)
                         tCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATcg                                                                                                    <  2:182651/139‑1 (MQ=255)
                                    tCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCt                                                                                         <  2:55105/139‑1 (MQ=255)
                                        aCGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCggg                                                                                     <  1:243373/139‑1 (MQ=255)
                                                 cAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATggg                                                                            <  1:748620/139‑1 (MQ=255)
                                                                         caaTCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCgggg                                                                  <  1:551644/125‑1 (MQ=255)
                                                                         caaTCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCgggg                                                                  >  2:551644/1‑125 (MQ=255)
                                                                         caaTCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAAc                                                    <  1:419995/139‑1 (MQ=255)
                                                                          aaTCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACt                                                   <  2:698391/139‑1 (MQ=255)
                                                                          aaTCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACt                                                   <  2:567811/139‑1 (MQ=255)
                                                                            tCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATgg                                                                             >  1:472040/1‑111 (MQ=255)
                                                                            tCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATgg                                                                             <  2:472040/111‑1 (MQ=255)
                                                                                        tCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATc                                     <  1:234689/139‑1 (MQ=255)
                                                                                                    tACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACa                         <  1:753724/139‑1 (MQ=255)
                                                                                                         cTACCTTCACCTTTTCGCTGGGGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCTGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACa                         <  1:568472/134‑1 (MQ=255)
                                                                                                         cTACCTTCACCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACa                         >  2:568472/1‑134 (MQ=255)
                                                                                                                 aCCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAAc                                                    >  1:622040/1‑99 (MQ=255)
                                                                                                                 aCCTTTTCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAAc                                                    <  2:622040/99‑1 (MQ=255)
                                                                                                                      ttCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCt                                                                                         <  1:644342/57‑1 (MQ=255)
                                                                                                                      ttCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCt                                                                                         >  2:644342/1‑57 (MQ=255)
                                                                                                                      ttCGCTGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGtctc       >  2:219235/1‑139 (MQ=255)
                                                                                                                           tGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGt                                                                                                <  1:591345/45‑1 (MQ=255)
                                                                                                                           tGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGt                                                                                                >  2:591345/1‑45 (MQ=255)
                                                                                                                           tGGCGGGAGGGGGCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCAATACAATTGTATCTCCgg  >  1:400091/1‑139 (MQ=255)
                                                                                                                                       |                                                                                                                              
GCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACGGTCAAACAATCGCGGTAAGTGCTCGGGCGCAATGGGCGTACCGGGGTTTTCGACGATAACTTGCACCAGGTGATCGACCGTCTGACAGCGCACTACAATTGTCTCTCCGG  >  NC_000913/593332‑593593

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