Predicted mutation
evidence seq id position mutation annotation gene description
RA ACB122 3,065,508 G→A V53M (GTG→ATG)  rcsC_6 → Sensor histidine kinase RcsC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1223,065,5080GA95.0% 57.1 / NA 20V53M (GTG→ATG) rcsC_6Sensor histidine kinase RcsC
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (9/10);  minor base C (1/0);  total (10/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGCGGGCGTTGGCCGGGCTGCTGGGGCTGGGTGACCACTCGGCGCGTAAAAGCCATTACCCGGAGCTGTCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCGTGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCGGCTATGACGACCCGCAAGCGGTGCTGTTTTCCCTGACTGACCTGG  >  ACB122/3065376‑3065625
                                                                                                                                    |                                                                                                                     
ccgcGGGCGTTGGCCGGGCTGCTGGGGCGGGGTGCCCACTCGGCGCGGAAAAGCCAGTACCCGGCGCTGTCCGCCCGCCTGGCCGAACTGGCAGCCGAACGCCACCGTGACAAATGGCTGGTCGCAAACGCCCTg                                                                                                                     >  6:198845/3‑135 (MQ=255)
    gggCGGTGGCGGGGCGGCTGGGGCTGGGTGACCACTCGGCGCGTAAAAGCCATTCCCCGGAGCTGGCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACg                                                                                                                 <  5:198845/135‑1 (MQ=255)
              cGGGCTGCTGGGGCTGGGTGACCACTCGGCGCGTAAAAGCCATTACCCGGAGCTGTCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCa                                                                                                       <  5:184957/135‑1 (MQ=255)
              cGGGCTGCTGGGGCTGGGTGACCACTCGGCGCGTAAAAGCCATTACCCGGAGCTGTCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCa                                                                                                       <  8:40766/135‑1 (MQ=255)
                ggCTGCTGGGGCTGGGTGACCACTCGGCGCGTAAAAGCCATTACCCGGAGCTGTCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAgg                                                                                                     <  2:95777/135‑1 (MQ=255)
                     ctggggctgggTGACCACTCGGCGCGTAAAAGCCATTACCCGGAGCTGTCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGc                                                                                                <  5:50118/135‑1 (MQ=255)
                         ggctgggTGACCACTCGGCGCGTAAAAGCCATTACCCGGAGCTGTCCGCCCGCCTGGCCGAACTGGCAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCCCGGTATTTTCCAGGCCCGCCt                                                                                              >  6:124086/1‑133 (MQ=255)
                              ggTGACCACTCGGCGCGTAAAAGCCATTACCCGGAGCTGTCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTGACAAATGGCTGTTCGAAAACGCCATGCAAGGTATTTTCCAGGCCAGCCTTCAGGAc                                                                                       >  2:111707/1‑135 (MQ=255)
                                      cTCGGCGCGTAAAAGCCATTACCCGGAGCTGTCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACTCAATGCACGCTATTTTCCAGGCCAGCCTTCAGGACGGCATGCg                                                                               >  7:175158/1‑135 (MQ=255)
                                           cgcgTAAAAGCCATTACCCGGAGCTGTCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTgccg                                                                          >  8:180277/1‑135 (MQ=255)
                                                aaaaGCCATTACCCGGAGCTGTCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAgg                                                                                                     <  5:140834/103‑1 (MQ=255)
                                                aaaaGCCATTACCCGGAGCTGTCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAgg                                                                                                     >  6:140834/1‑103 (MQ=255)
                                                        ttACCCGGAGCTGTCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCt                                                             >  1:64372/1‑135 (MQ=255)
                                                               gAGCTGTCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGgcgcgc                                                      <  4:203762/135‑1 (MQ=255)
                                                                 gCTGTCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCAt                                                    <  4:179622/135‑1 (MQ=255)
                                                                        gcccgccTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCgg                                              >  3:3925/1‑134 (MQ=255)
                                                                        gcccgccTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCgg                                              >  1:129396/1‑134 (MQ=255)
                                                                                     aaCTGGAACCCGAACGCAACCGTTACAAATGGCTGTTCGAAACCGCCATGCACGGTATTTTCCAGGCCAGCCTTCTGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCGCGGCTATGACGACCCGc                                <  1:205010/135‑1 (MQ=255)
                                                                                                               aaaTGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCGGCTATGACGACCCGCAAGCGGTGCTGTTTTCCctgactga       <  7:182747/134‑1 (MQ=255)
                                                                                                                   ggCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCGGCTATGACGACCCGCAAGCGGTGCTGTTTTCCCTGACTGACCTgg  >  4:156442/1‑135 (MQ=255)
                                                                                                                                    |                                                                                                                     
CAGCGGGCGTTGGCCGGGCTGCTGGGGCTGGGTGACCACTCGGCGCGTAAAAGCCATTACCCGGAGCTGTCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCGTGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCGGCTATGACGACCCGCAAGCGGTGCTGTTTTCCCTGACTGACCTGG  >  ACB122/3065376‑3065625

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