Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 5,485,083 C→A 13.3% I409I (ATC→ATA)  rcsC_18 → Sensor histidine kinase RcsC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1435,485,0830CA13.3% 89.7 / 5.4 30I409I (ATC→ATA) rcsC_18Sensor histidine kinase RcsC
Reads supporting (aligned to +/- strand):  ref base C (13/13);  new base A (2/2);  total (15/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 = 6.05e-01

TCCTCAACCTGGCGGTCAGTGAACGGCTGCGGCGCATAAGCGAGGAGCGTTTCAGCGCCAGCCGTGCACTGGCCGCCAGCGATGCGGAAATCAACGCCAAAGCCGAGTTTCTGGCCAAGATCAGCCACGAAATCCGCACCCCGATGAACGGCGTGCTCGGCATGACTGAACTGCTGCTGGGCACACCGCTGTCGGTGAAGCAACGTGATTATGTGCAGACCATCCACAGCGCGGGCAACGAGCTGCTGACGCTGATCA  >  ACB143/5484950‑5485207
                                                                                                                                     |                                                                                                                            
tcctcAACCTGGCGGTCAGTGAACGGCTGCGGCGCATAAGCGAGGAGCGTTTCAGCGCCAGCCGTGCACTGGCCGCCAGCGATGCGGAAATCAACGCCAAAGCCGAGTTTCTGGCCAAGATCAGCCACGAAATCCg                                                                                                                            >  2:192564/1‑136 (MQ=255)
     aaCCTGGCGGTCAGTGAACGGCTGCGGCGCATAAGCGAGGAGCGTTTCAGCGCCAGCCGTGCACTGGCCGCCAGCGATGCGGAAATCAACGCCAAAGCCGAGTTTCTGGCCAAGATCAGCCACGAAATCCGCAccc                                                                                                                       >  7:102313/1‑136 (MQ=255)
            cGGTCAGTGAACGGCTGCGGCGCATAAGCGAGGAGCGTTTCAGCGCCAGCCGTGCACTGGCCGCCAGCGATGCGGAAATCAACGCCAAAGCCGAGTTTCTGGCCAAGATCAGCCACGAAATCCGCACCCCGATGaa                                                                                                                >  8:212781/1‑136 (MQ=255)
             ggTCAGTGAACGGCTGCGGCGCATAAGCGAGGAGCGTTTCAGCGCCAGCCGTGCACTGGCCGCCAGCGATGCGGAAATCAACGCCAAAGCCGAGTTTCTGGCCAAGATCAGCCACGAAATCCGCACCCCGATGAAc                                                                                                               >  6:241274/1‑136 (MQ=255)
                    gAACGGCTGCGGCGCATAAGCGAGGAGCGTTTCAGCGCCAGCCGTGCACTGGCCGCCAGCGATGCGGAAATCAACGCCAAAGCCGAGTTTCTGGCCAAGATCAGCCACGAAATCCGCACCCCGATGAACGGCGTGc                                                                                                        >  3:144751/1‑136 (MQ=255)
                     aaCGGCGGGGGCGAATATGCGAGGAGGGTTTCAGCGCCCGCGGGCCCCGGGCCCCCAGCGATGCGGAAATCAACGCCAAAGCCGAGTTTCTGGCCAAGATCAGCCACGAAATCCGCACCCCGATGAACGGCGTGCt                                                                                                       <  1:69871/136‑1 (MQ=255)
                      ccgggggcggccCATAAGCGTGGCGCGTTGCCGGGCATGCGGGGCAGGGGCCGCCAGCGATGCGGAAATCAACGCCAAAGCCGAGTTTCTGGCCAAGATCAGCCACGAAATCCGCACCCCGATGAACGGCGTGCt                                                                                                       <  8:230792/129‑1 (MQ=255)
                         gCTGCGGCGCATAAGCGAGGAGCGTTTCAGCGCCAGCCGTGCACTGGCCGCCAGCGATGCGGAAATCAACGCCAAAGCCGAGTTTCTGGCCAAGATCAGCCACGAAATCCGCACCCCGATGAACGGCGTGCTCgg                                                                                                    <  7:7589/135‑1 (MQ=255)
                                             agCGTTTCAGCGCCAGCCGTGCACTGGCCGCCAGCGATGCGGAAATCAACGCCAAAGCCGAGTTTCTGGCCAAGATCAGCCACGAAATCCGCACCCCGATGAACGGCGTGCTCGGCATGACTGAACTGCTGCTggg                                                                               >  1:25424/1‑136 (MQ=255)
                                             agCGTTTCACCGCCAGCCGTCCACGGGCCGCCAGCGATGCGGAAATCAACGCCAAAGCCGAGTTTCTGGCCAAGATCAGCCACGAAATCCGCACCCCGATGAACGGCGTGCTCGGCATGACTGAACTGCTGCTggg                                                                               <  5:79011/136‑1 (MQ=255)
                                                  ttCAGCGCCAGCCGTGCACTGGCCGCCAGCGATGCGGAAATCAACGCCAAAGCCGAGTTTCTGGCCAAGATCAGCCACGAAATCCGCACCCCGATGAACGGCGTGCTCGGCATGACTGAACTGCTGCTGGGCacac                                                                          >  3:205617/1‑136 (MQ=255)
                                                    cAGCGCCAGCCGTGCACTGGCCGCCAGCGATGCGGAAATCAACGCCAAAGCCGAGTTTCTGGCCAAGATCAGCCACGAAATCCGCACCCCGATGAACGGCGTGCTCGGCATGACTGAACTGCTGCTGGGCACACCg                                                                        >  4:115287/1‑136 (MQ=255)
                                                                    cTGGCCGCCAGCGATGCGGAAATCAACGCCAAAGCCGAGTTTCTGGCCAAGATCAGCCACGAAATCCGCACCCCGATGAACGGCGTGCTCGGCATGACTGAACTGCTGCTGGGCACACCGCTGTCGGTGAAGCAAc                                                        >  3:278370/1‑136 (MQ=255)
                                                                      ggCCGCCAGCGATGCGGAAATCAACGCCAAAGCCGAGTTTCTGGCCAAGATCAGCCACGAAATCCGCACCCCGATGAACGGCGTGCTCGGCATGACTGAACTGCTGCTGGGCACACCGCTGTCGGTGAAGCAACGt                                                      <  3:98618/136‑1 (MQ=255)
                                                                          gccAGCGATGCGGAAATCAACGCCAAAGCCGAGTTTCTGGCCAAGATCAGCCACGAAATCCGCACCCCGATGAACGGCGTGCTCGGCATGACTGAACTGCTGCTGGGCACACCGCTGTCGGTGAAGCAACGTGAt                                                   >  3:69974/1‑135 (MQ=255)
                                                                          gccAGCGATGCGGAAATCAACGCCAAAGCCGAGTTTCTGGCCAAGATCAGCCACGAAATCCGCACCCCGATGAACGGCGTGCTCGGCATGACTGAACTGCTGCTGGGCACACCGCTGTCGGTGAAGCAACGTGAt                                                   >  4:213553/1‑135 (MQ=255)
                                                                          gccAGCGATGCGGAAATCAACGCCAAAGCCGAGTTTCTGGCCAAGATCAGCCACGAAATCCGCACCCCGATGAACGGCGTGCTCGGCATGACTGAACTGCTGCTGGGCACACCGCTGTCGGTGAAGCAACGTGAt                                                   >  7:77624/1‑135 (MQ=255)
                                                                           ccAGCGATGCGGAAATCAACGCCAAAGCCGAGTTTCTGGCCAAGATCAGCCACGAAATCCGCACCCCGATGAACGGCGTGCTCGGCATGACTGAACTGCTGCTGGGCACACCGCTGTCGGTGAAGCAACGTGATTa                                                 <  8:129154/136‑1 (MQ=255)
                                                                             aGCGATGCGGTAATCAACGCCAAAGCCGAGTTTCTGGCCAAGATCAGCCACGAAATCCGCACCCCGATGAACGGCGTGCTCGGCATGACTGAACTGCTGCTGGGCACACCGCTGTCGGTGAAGCAACGTGATTAtg                                               <  7:277395/136‑1 (MQ=255)
                                                                             aGCGATGCGGAAATCAACTCCAAAGCCGAGTTTCTGGCCAAGATCAGCCACGAAATCCGCACCCCGATGAACGGCGTGCTCGGCATGACTGAACTGCTGCTGGGCACACCGCTGTCGGTGAAGCAACGTGATTAtg                                               <  8:34412/136‑1 (MQ=255)
                                                                             aGCGATGCGGAAATCAACGCCAAAGCCGAGTTTCTGGCCAAGATCAGCCACGAAATACGCACCCCGATGAACGGCGTGCTCGGCATGACTGAACTGCTGCTGGGCACACCGCTGTCGGTGAAGCAACGTg                                                     <  2:282057/130‑1 (MQ=255)
                                                                             aGCGATGCGGAAATCAACGCCAAAGCCGAGTTTCTGGCCAAGATCAGCCACGAAATACGCACCCCGATGAACGGCGTGCTCGGCATGACTGAACTGCTGCTGGGCACACCGCTGTCGGTGAAGCAACGTg                                                     <  2:277546/130‑1 (MQ=255)
                                                                             aGCGATGCGGAAATCAACGCCAAAGCCGAGTTTCTGGCCAAGATCAGCCACGAAATACGCACCCCGATGAACGGCGTGCTCGGCATGACTGAACTGCTGCTGGGCACACCGCTGTCGGTGAAGCAACGTg                                                     >  1:282057/1‑130 (MQ=255)
                                                                             aGCGATGCGGAAATCAACGCCAAAGCCGAGTTTCTGGCCAAGATCAGCCACGAAATACGCACCCCGATGAACGGCGTGCTCGGCATGACTGAACTGCTGCTGGGCACACCGCTGTCGGTGAAGCAACGTg                                                     >  1:277546/1‑130 (MQ=255)
                                                                                gATGCGGAAATCAACGCCAAAGCCGAGTTTCTGGCCAAGATCAGCCACGAAATCCGCACCCCGATGAACGGCGTGCTCGGCATGACTGAACTGCTGCTGGGCACACCGCTGTCGGTGAAGCAACGTGATTATGTGc                                            <  6:144161/136‑1 (MQ=255)
                                                                                                  aaaGCCGAGGTTCTGGCCAAGATCAGCCACGAAATCCGCACCCCGATGAACGGCGTGCTCGGCATGACTGAACTGCTGCTGGGCACACCGCTGTCGGTGAAGCAACGTGATTATGTGCAGACCATCCACAGCGCgg                          <  1:130424/136‑1 (MQ=255)
                                                                                                      ccGAGTTTCTGGCCAAGATCAGCCACGAAATCCGCACCCCGATGAACGGCGTGCTCGGCATGACTGAACTGCTGCTGGGCACACCGCTGTCGGTGAAGCAACGTGATTATGTGCAGACCATCCACAGCGCGGGCaa                      >  4:82265/1‑136 (MQ=255)
                                                                                                      ccGAGTTTCTGGCCAAGATCAGCCACGAAATCCGCACCCCGATGAACGGCGTGCTCGGCATGACTGAACTGCTGCGGGGCACACCGCTGTCGGTGAAGCAACGTGATTATGTGCAGACCATCCACAGCGCGGGCaa                      <  8:55981/136‑1 (MQ=255)
                                                                                                                     aGATCAGCCACGAAATCCGCACCCCGATGAACGGCGTGCTCGGCATGACTGAACTGCTGCTGGGCACACCGCTGTCGGTGAAGCAACGTGATTATGTGCAGACCATCCACAGCGCGGGCAACGAGCTGCTGACGCt       <  5:7561/136‑1 (MQ=255)
                                                                                                                          aGCCACGAAATCCGCACCCCGATGAACGGCGTGCTCGGCATGACTGAACTGCTGCTGGGCACACCGCTGTCGGTGAAGCAACGTGATTATGTGCAGACCATCCACAGCGCGGGCAACGAGCTGCTGACGCTGATCa  <  7:184723/136‑1 (MQ=255)
                                                                                                                                     |                                                                                                                            
TCCTCAACCTGGCGGTCAGTGAACGGCTGCGGCGCATAAGCGAGGAGCGTTTCAGCGCCAGCCGTGCACTGGCCGCCAGCGATGCGGAAATCAACGCCAAAGCCGAGTTTCTGGCCAAGATCAGCCACGAAATCCGCACCCCGATGAACGGCGTGCTCGGCATGACTGAACTGCTGCTGGGCACACCGCTGTCGGTGAAGCAACGTGATTATGTGCAGACCATCCACAGCGCGGGCAACGAGCTGCTGACGCTGATCA  >  ACB143/5484950‑5485207

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: ‑