Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 3,065,508 G→A 89.1% 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,5080GA89.1% 127.3 / 7.6 46V53M (GTG→ATG) rcsC_6Sensor histidine kinase RcsC
Reads supporting (aligned to +/- strand):  ref base G (2/3);  new base A (22/19);  total (24/22)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

CAGCGGGCGTTGGCCGGGCTGCTGGGGCTGGGTGACCACTCGGCGCGTAAAAGCCATTACCCGGAGCTGTCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCGTGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCGGCTATGACGACCCGCAAGCGGTGCTGTTTTCCCTGACTGACCTGGCCGCCAAC  >  ACB122/3065376‑3065633
                                                                                                                                    |                                                                                                                             
cagcGGGCGTTGGCCGGGCTGCTGGGGCTGGGTGACCACTCGGCGCGTAAAAGCCATTACCCGGAGCTGTCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATg                                                                                                                             >  8:59742/1‑135 (MQ=255)
cagcGGGCGTTGGCCGGGCTGCTGGGGCTGGGTGACCACTCGGCGCGTAAAAGCCATTACCCGGAGCTGTCCGCCCCCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATg                                                                                                                             >  1:201192/1‑135 (MQ=255)
 agcGGGCGTTGGCCGGGCTGCTGGGGCTGGGTGACCACTCGGCGCGTAAAAGCCATTACCCGGAGCTGTCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGc                                                                                                                            <  4:30681/135‑1 (MQ=255)
   cGGGCGTTGGCCGGGCTGCTGGGGCTGGGTGACCACTCGGCGCGTAAAAGCCATTACCCGGAGCTGTCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACg                                                                                                                         <  8:152004/136‑1 (MQ=255)
        gTTGGCCGGGCTGCTGGGGCTGGGTGACCACTCGGCGCGTAAAAGCCATTACCCGGAGCTGTCCGCCCGCCTGGACGAACTGGAAGCCGAACGCACCCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTag                                                                                                                     >  5:318401/1‑134 (MQ=255)
              cGGGCTGCTGGGGCTGGGTGACCACTCGGCGCGTAAAAGCCATTACCCGGAGCTGTCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCa                                                                                                               <  2:311732/135‑1 (MQ=255)
                ggcggcGGGGGCTGGGGGACCACTCGGCGCGTAAAAGCCATTACCCGGAGCTGTCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAgg                                                                                                             <  5:218592/135‑1 (MQ=255)
                ggCTGCTGGGGCTGGGTGACCACTCGGCGCGTAAAAGCCATTACCCGGAGCTGTCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAgg                                                                                                             <  2:315651/135‑1 (MQ=255)
                      tggggctgggTGACCACTCGGCGCGTAAAAGCCATTACCCGGAGCTGTCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGcc                                                                                                       <  6:172436/135‑1 (MQ=255)
                        gggctgggTGACCACTCGGCGCGTAAAAGCCATTACCCGGAGCTGTCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCtt                                                                                                     <  3:206064/135‑1 (MQ=255)
                         ggctgggTGACCACTCGGCGCGTAAAAGCCATTACCCGGAGCTGTCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCa                                                                                                   <  8:49062/136‑1 (MQ=255)
                              ggTGACCACTCGGCGCGTAAAAGCCATTACCCGGAGCTGTCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGc                                                                                                        >  4:10301/1‑126 (MQ=255)
                              ggTGACCACTCGGCGCGTAAAAGCCATTACCCGGAGCTGTCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGc                                                                                                        <  3:10301/126‑1 (MQ=255)
                                tGACCACTCGGCGCGTAAAAGCCATTACCCGGAGCTGTCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACg                                                                                              >  3:177569/1‑134 (MQ=255)
                                  aCCACTCGGCGCGTAAAAGCCATTACCCGGAGCTGTCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCa                                                                                           >  4:193253/1‑135 (MQ=255)
                                  aCCACTCGGCGCGTAAAAGCCATTACCCGGAGCTGTCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCa                                                                                           >  4:246335/1‑135 (MQ=255)
                                      cTCGGCGCGTAAAAGCCATTACCCGGAGCTGTCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCg                                                                                       <  3:39180/135‑1 (MQ=255)
                                            gcgTAAAAGCCATTACCCGGAGCTGTCCGCCCGCCTGGACGAACAGGAAGCCGACCGCAACCGTTACAAATGGGTCTTCGAACACGCCATCCACGGTATTATCCAGGCCAGCCTTTACGACCGCATGCggcccg                                                                                  >  6:88480/1‑129 (MQ=255)
                                                    gCCATTACCCGGAGCTGTCAGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTAGAAAACGCCGTGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCAGgc                                                                        >  8:135468/1‑136 (MQ=255)
                                                           cccGGAGCTGTCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGgcg                                                                 >  7:235795/1‑136 (MQ=255)
                                                             cGGAGCTGTCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGgcgcg                                                               >  7:368607/1‑136 (MQ=255)
                                                                aGCTGTCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGGTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGACATGCGTGCCGCCAACCCGGCGCTGGCGCGCAt                                                            <  7:306635/136‑1 (MQ=255)
                                                                 gCTGTCCGCCCGCCTGGACGAACTGGAAGCCGACCGCAACCGTTACAACTGGCTGTTCGAAAACGCCGTGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCAt                                                            <  1:331337/135‑1 (MQ=255)
                                                                 gCTGTCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCAt                                                            <  4:342730/135‑1 (MQ=255)
                                                                 gCTGTCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCAt                                                            <  2:70238/135‑1 (MQ=255)
                                                                 gCTGTCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCAt                                                            <  3:390123/135‑1 (MQ=255)
                                                                 gCTGGCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCGTGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCAt                                                            <  5:24302/135‑1 (MQ=255)
                                                                  cTGTCCGCCCGCCTGGACGAACTGGAAGCCGACCGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATg                                                           >  1:396859/1‑135 (MQ=255)
                                                                           cgccTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCGTGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCGGCTAt                                                  <  4:263166/135‑1 (MQ=255)
                                                                             ccTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCGGCTATgac                                               >  8:191605/1‑136 (MQ=255)
                                                                             ccTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCGGCTATga                                                >  1:59839/1‑135 (MQ=255)
                                                                               tGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCGCGGCTATgacg                                              >  6:238391/1‑135 (MQ=255)
                                                                                    gAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCGGCTATGACGACCCGc                                        <  7:59742/136‑1 (MQ=255)
                                                                                     aaCTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCGGCTATGACGACCCGc                                        <  4:372893/135‑1 (MQ=255)
                                                                                         ggAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAgg                                                                                                             >  5:221858/1‑62 (MQ=255)
                                                                                         ggAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCGGCTATGACGACCCGCAAGc                                    >  4:357479/1‑135 (MQ=255)
                                                                                         ggAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCGGCTATGACGACCCGCAAGCg                                   >  7:6888/1‑136 (MQ=255)
                                                                                         ggAAGCCGAACGCAACAGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAgg                                                                                                             <  6:221858/62‑1 (MQ=255)
                                                                                            aGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCGGCTATGACGACCCGCAAGCGGt                                 <  4:177569/135‑1 (MQ=255)
                                                                                                    gCAACCGTTACAAATGGCTGTTCGAAAACGCCGTGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCGGCTATGACGACCCGCAAGCGGTGCTGtttt                         >  1:372143/1‑135 (MQ=255)
                                                                                                                        ttCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCGGCTATGACGACCCGCAAGCGGTGCTGTTTTCCCTGACTGACCTGgccgcc     >  3:367236/1‑135 (MQ=255)
                                                                                                                         tCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCGGCTATGACGACCCGCAAGCGGTGCTGTTTTCCCTGACTGACCTGGCCGCCa    >  1:57526/1‑135 (MQ=255)
                                                                                                                         tCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCGCGGCTATGACGACCCGCAAGCGGTGCTGTTTTCCCTGACTGACCTGgccgcc     <  6:304009/134‑1 (MQ=255)
                                                                                                                           gAAAACGCCATGCACGGTATTTTCCAGGCTAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCGGCTATGACGACCCGCAAGCGGTGCTGTTTTCCCTGACTGACCTGGCCGCCAAc  >  1:234760/1‑135 (MQ=255)
                                                                                                                           gAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCGGCTATGACGACCCGCAAGCGGTGCTGTTTTCCCTGACTGACCTGGCCGCCAAc  >  1:329719/1‑135 (MQ=255)
                                                                                                                           gAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGAGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCGGCTATGACGACCCGCAAGCGGTGCTGTTTTCCCTGACTGACCTGGCCGCCAAc  >  1:257340/1‑135 (MQ=255)
                                                                                                                                    |                                                                                                                             
CAGCGGGCGTTGGCCGGGCTGCTGGGGCTGGGTGACCACTCGGCGCGTAAAAGCCATTACCCGGAGCTGTCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCGTGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCGGCTATGACGACCCGCAAGCGGTGCTGTTTTCCCTGACTGACCTGGCCGCCAAC  >  ACB122/3065376‑3065633

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