Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 593,357 T→C T472A (ACA→GCA)  cusS ← sensory histidine kinase CusS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913593,3570TC94.3% 104.9 / ‑2.1 35T472A (ACA→GCA) cusSsensory histidine kinase CusS
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (19/14);  total (20/15)
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.69e-01

AGGTTTTTTATGAATTGTGACAGCTTTTTACCATTAATAGGTATGACTATTGCGGCACGTTATTTTTACACTGGTTATAAAAGTTGCCGTTTGCTGAAGGATTAAGCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACG  >  NC_000913/593227‑593483
                                                                                                                                  |                                                                                                                              
aGGTTTTTTATGAATTGTGACAGCTTTTTACCATTAATAGGTATGACTATTGCGGCACGTTATTTTTACACTGGTTATAAAAGTTGCCGTTTGCTGAAGGATTAAGCGGGTAATGTGATAACAAACCTTGccc                                                                                                                              >  2:459709/1‑133 (MQ=255)
     ttttATGAATTGTGACAGCTTTTTACCATTAATAGGTATGACTATTGCGGCACGTTATTTTTACACTGGTTATAAAAGTTGCCGTTTGCTGAAGGATTAAGCGGGTAATGTGATAACAAACCTTGCCCCCcgc                                                                                                                         >  2:295988/1‑133 (MQ=255)
      tttATGAATTGTGACAGCTTTTTACCATTAATAGGTATGACTATTGCGGCACGTTATTTTTACACTGGTTATAAAAGTTGCCGTTTGCTGAAGGATTAAGCGGGTGATGTGATAACAAACCTTGCCCCCcgcg                                                                                                                        >  1:264854/1‑133 (MQ=255)
              ttGTGACAGCCTTTTACCATTAATTGGTTTGACTATTGCGGTACGTTATTTTTACACTGGTTATAAAAGTTGCCGTTTGCTGAAGGATTAAGCGGGTAATGTGATAACAAACCTTGCCCCCCGCGCATCCGAc                                                                                                                <  2:181657/133‑1 (MQ=255)
                         ttttACCATTAATAGGTATGACTATTGCGGCACGTTATTTTTACACTGGTTATAAAAGTTGCCGTTTGCTGAAGGATTAAGCGGGTAATGTGATAACAAACCTTGCCCCCCGCGCATCCGACGTTACCGCAAc                                                                                                     >  1:337156/1‑133 (MQ=255)
                                         tATGACTATTGCGGCACGTTATTTTTACACTGGTTATAAAAGTTGCCGTTTGCTGAAGGATTAAGCGGGTAATGTGATAACAAACCTTGCCCCCCGCGCAGCCGGCGTCACCGCAACTGTGCCTTTATGcgcg                                                                                     >  2:192784/1‑133 (MQ=255)
                                           tGACTATTGCGGCACGTTATTTTTACACTGGTTATAAAAGTTGCCGTTTGCTGAAGGATTAAGCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGac                                                                                   >  1:26529/1‑133 (MQ=255)
                                             aCTATTGCAGCACGTTTTATTTTCCCTGGTTATAAAAGTTGCCGTTTGCTGAAGGATTAAGCGGGTAATGTGATAACAAACCTTGCCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGacaa                                                                                 <  1:209835/133‑1 (MQ=255)
                                                 ttGCGGCACGTTATTTTTACACTGGTTATAAAAGTTGCCGTTTGCTGAAGGATTAAGCGGGTAATGTGATAACAAACCTTGCCCCCCGCGCGTCCGACGTAACCACAAACGTGCCTTGATGCACGACAACAAt                                                                             >  2:148382/1‑133 (MQ=255)
                                                  tGCGGCACGTTATTTTTACACTGGTTATAAAAGTTGCCGTTTGCTGAAGGATTAAGCGGGTAATGTGATAACAAACCTTGCCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATc                                                                            >  1:417072/1‑133 (MQ=255)
                                                       cACGTTATTTTTACACTGGTTATAAAAGTTGCCGTTTGCTGAAGGATTAAGCGGGTAATGTGATAACAAACCTTGCCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGAttt                                                                       >  1:383329/1‑133 (MQ=255)
                                                         cGTTATTTTTACACTGGTTATAAAAGTTGCCGTTTGCTGAAGGATTAAGCGGGTAATGTGATAACAAACCTTGCCCCCCGCGCATCCGACGCTACCGCAACCGTGCTTATATGCGCGGCAACTATCGATTTCa                                                                     >  2:227192/1‑133 (MQ=255)
                                                                tttACACTTGTTATACAATTTGGCGTTTGCTGAAGGATTAAGCGGGTAATGTGATAACAAACCTTGCCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGc                                                              <  1:192784/133‑1 (MQ=255)
                                                                tttACACTTGTTATACAACTTGCCGTTTGCTGAAGGATTAAGCGGGTAATGTGATAACAAACCTTGCCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGc                                                              <  1:148382/133‑1 (MQ=255)
                                                                tttACACTGGTTATAAAAGTTGCCGTTTGCTGAAGGATTAAGCGGGTAATGTGATAACAACCCTTGCCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGc                                                              <  1:36381/133‑1 (MQ=255)
                                                                tttACACTGGTTATAAAAGTTGCCGTTTGCTGAAGGATTAAGCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGc                                                              <  2:26529/133‑1 (MQ=255)
                                                                tttACACTGGTTATAAAAGTTGCCGTTTGCTGAAGGATTAAGCGGGTAATGTGATAACAAACCTTGCCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGc                                                              <  2:417072/133‑1 (MQ=255)
                                                                tttACACTGGTTATAAAAGTTGCCGTTTGCTGAAGGATTAAGCGGGTAATGTGATAACAAACCTTGCCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGc                                                              <  1:26047/133‑1 (MQ=255)
                                                                   acacTGGTTATAAAAGTTGCCGTTTGCTGAAGGATTAAGCGGGTAATGTGATAACAAACCTTGCCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAg                                                           >  2:81099/1‑133 (MQ=255)
                                                                        ggTTATAAAAGTTGCCGTTTGCTGAAGGATTAAGCGGGTAATGTGATAACAAACCTTGCCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGc                                                              >  2:378172/1‑125 (MQ=255)
                                                                        ggTTATAAAAGTTGCCGTTTGCTGAAGGATTAAGCGGGTAATGTGATAACAAACCTTGCCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGc                                                              <  1:378172/125‑1 (MQ=255)
                                                                        ggTTATAAAAGTTGCCGTTTGCTGAAGGATTAAGCGGGTAATGTGATAACAAACCTTGCCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCaa                                                      <  2:383329/133‑1 (MQ=255)
                                                                        ggTTATAAAAGTTGCCGTTTGCTGAAGGATTAAGCGGGTAATGTGATAACAAACCTTGCCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCaa                                                      >  2:363688/1‑133 (MQ=255)
                                                                             taAAAGTTGCCGTTTGCTGAAGGATTAAGCGGGTAATGTGATAACAAACCTTGCCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCg                                                 <  1:81099/133‑1 (MQ=255)
                                                                                   ttGCCGTTTGCTGAAGGATTAAGCGGGTAATGTGATAACAAACCTTGCCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCt                                           >  2:4099/1‑133 (MQ=255)
                                                                                   ttGCCGTTTGCTGAAGGATTAAGCGGGTAATGTGATAACAAACCTTGCCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCt                                           >  1:344833/1‑133 (MQ=255)
                                                                                   ttGCCGTTTCTTGAAGGATTAAGTAGGTAATGCGATAACAAACCTTGCCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCt                                           <  2:244016/133‑1 (MQ=255)
                                                                                                       aaGCGGGTAATGTGATAACAAACCTTGCCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGa                       >  2:313207/1‑133 (MQ=255)
                                                                                                       aaGCGGGTAATGTGATAACAAACCTTGCCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGa                       >  2:401455/1‑133 (MQ=255)
                                                                                                               aaTGTGATAACAAACCTTGCCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCaa               <  1:401455/133‑1 (MQ=255)
                                                                                                               aaTGTGATAACAAACCTTGCCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCaa               <  1:363688/133‑1 (MQ=255)
                                                                                                               aaTGTGATAACAAACCTTGCCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCaa               <  1:295988/133‑1 (MQ=255)
                                                                                                                       aaCAAACCTTGCCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAg       >  1:96469/1‑133 (MQ=255)
                                                                                                                        aCAAACCTTGCCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGa      <  2:344833/133‑1 (MQ=255)
                                                                                                                        aCAAAACTTGCCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAgg      >  1:389127/1‑132 (MQ=255)
                                                                                                                           aaCCTTGCCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAAc   >  2:502694/1‑133 (MQ=255)
                                                                                                                            aCCTTGCCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACg  >  1:485181/1‑133 (MQ=255)
                                                                                                                                  |                                                                                                                              
AGGTTTTTTATGAATTGTGACAGCTTTTTACCATTAATAGGTATGACTATTGCGGCACGTTATTTTTACACTGGTTATAAAAGTTGCCGTTTGCTGAAGGATTAAGCGGGTAATGTGATAACAAACCTTGTCCCCCGCGCATCCGACGTTACCGCAACCGTGCCTTTATGCGCGACAACAATCGATTTCACTATCGCCAGGCCAATACCGCTACCTTCACCTTTTCGCTGGCGGGAGGGGTCAACGCGATAGAAACG  >  NC_000913/593227‑593483

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