Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,313,302 G→A 9.3% F18F (TTC→TTT)  rcsC ← hybrid sensory kinase in two‑component regulatory system with RcsB and YojN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,313,3020GA 9.3% 126.1 / 4.3 43F18F (TTC→TTT) rcsChybrid sensory kinase in two‑component regulatory system with RcsB and YojN
Reads supporting (aligned to +/- strand):  ref base G (19/20);  new base A (2/2);  total (21/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.92e-01

CCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGATCGGAACTCAGATTAAATTCCTGACGAATTTCCGATTCTCGCTGATGTAACGCATTAACGATGTAAAAAACGGATGAAAAAGCAATCAACAGCCAGAGCACTAACGCCAATGCTCTGAACATGTAGCGCGAGGCTTTCAGGGTTGTACGAAAAGAAGCAAGGTATTTCAAAGGGGCGAAGCTCCGCCTCAGGTGACCGATGGAGTGTGGTTAAGGTAGCGGTAAAAGCGTGTTACCGCAATGTTCTCTCTTCTCTGGA  >  CP009273/2313067‑2313443
                                                                                                                                                                                                                                           |                                                                                                                                             
ccTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGATCGGAACTCAGATTAAATTCCTGACGAATTTCCGATTCTCGCTGATGTAACGCATTAACGATGTAAAAAACGGATGAAAAAGCAATCAACAGCCAGAGCACTAACGCCAATGCTCTGaa                                                                                                                                             >  2:226903/1‑238 (MQ=255)
    cGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGATCGGAACTCAGATTAAATTCCTGACGAATTTCCGATTCTCGCTGATGTAACGCATTAACGATGTAAAAAACGGATGAAAAAGCAATCAACAGCCAGAGCACTAACGCCAATGCTCTGAACATg                                                                                                                                         <  1:226907/238‑1 (MQ=255)
     tcgtgacGTCAACGCGGGGAAATCACACCGTTTTCTGCCGATAAGCTATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCATTGAATAAAGCGCTGAGCCTGATCGGAACTCAGATTAAATTCCTGACGAATTTCCGATTCTCGCTGATGTAACTCATTAACGATGTAAAAAACGGATGAAAAAGCAATCCACAGCCAGAGCACTAACGCCAATGCTCTGAACATGTAGCGCGTGGCTTTCCGGGTTGTACGAAAAGAAGCa                                                                                                     <  2:109129/267‑1 (MQ=255)
                                        tGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCATAGGTGCGTTGAATAAAGCGCTGAGCCTGATCGGAACTCAGATTAAATTCCTGACGAATTTCCGATTCTCGCTGATGTAACGCATTAACGATGTAAAAAACGGATGAAAAAGCAATCAACAGCCAGAGCACTAACGCCAATGCTCTGAACATGTAGCGCGAGGCTTTCAGGGTTGTACGAAAAGAAGCAAGGTAt                                                                                               >  1:421420/1‑244 (MQ=255)
                                            gATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCATAGGTGCGTTGAATAAAGCGCTGAGCCTGATCGGAACTCAGATTAAATTCCTGACGAATTTCCGATTCTCGCTGATGTAACGCATTAACGATGTAAAAAACGGATGAAAAAGCAATCAACAGCCAGAGCACTAACGCCAATGCTCTGAACATGTAGCGCGAGGCTTTCAGGGTTGTACGAAAAGAAGCAAGGTATTTCa                                                                                           <  2:421415/244‑1 (MQ=255)
                                                                      cAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGATCGGAACTCAGATTAAATTCCTGACGAATTTCCGATTCTCGCTGATGTAACGCATTAAAGATGTAAAAAACGGATGAAAAAGCAATCAACAGCCAGAGCACTAACGCCAATGCTCTGAACATGTAGCGCGAGGCTTTCAGGGTTGTACGAAAAGAAGCAAGGTATTTCAAAGGGGCGAAGCTCCGCCTCAGGTGACCGATGGAGGGGGGTTactgtcgataatac                                   >  1:87521/1‑261 (MQ=255)
                                                                           ctTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGATCGGAACTCAGATTAAATTCCTGACGAATTTACGATTCTCGCTGATGTAACGCATTAACGATGTAAAAAACGGATGAAAAAGCAATCAACAGCCAGAGCACTAACGCCAATGCTCTGAACATGTAgcg                                                                                                                                    >  1:341933/1‑172 (MQ=255)
                                                                               catcacatcaGAGGTGCGTTGAATAAAGCGCTGAGCCTGATCGGAACTCAGATTAAATTCCTGACGAATTTACGATTCTCGCTGATGTAACGCATTAACGATGTAAAAAACGGATGAAAAAGCAATCAACAGCCAGAGCACTAACGCCAATGCTCTGAACATGTAGCGCGAg                                                                                                                                <  2:341928/172‑1 (MQ=255)
                                                                                                gTTGAATAAAGCGCTGAGCCTGATCGGAACTCAGATTAAATTCCTGACGAATTTCCGATTCTCGCTGATGTAACGCATTAACGATGTAAAAAACGGATGAAAAAGCAATCAACAGCCAGAGCACTAACGCCAATGCTATGa                                                                                                                                              <  1:80508/141‑1 (MQ=255)
                                                                                                 ctGAATAAAGCGCTGAGCCTGATCGGAACTCAGATTAAATTCCTGACGAATTTCCGATTCTCGCTGGTGTAACGCATTAACGATGTAAAAAACGGATGAAAAAGCAATCAACAGACAGAGCACTAACGCCAATGCTCTGAAAATGTAGCGCAAGGCTTTCAGGGTTGAACGAAAAGAAGCAAGGTATTTCAAAGGGGCGAAACTCAGCCTAAGGAGACCGATGAAGTGTGGTTAAGGTAGGGGGAAAAGCGTGTTACTGCAATGTtcgctgtccg       >  2:179301/2‑270 (MQ=255)
                                                                                                     aTAAAGCGCTGAGCCTGATCGGAACTCAGATTAAATTCCTGACGAATTTCCGATTCTCGCTGATGTAACGCATTAACGATGTAAAAAACGGATGAAAAAGCAATCAACAGCCAGAGCACTAACGCCAATGCTCTGAACATGTAGCGCGAGGCTTTCAGGGTTGTACGAAAAGAAGCAAGGTATTTCAAAGGGGCGAAGCTCCGCCTCAGGTGACCGATGGAGTGTGGTTAAGGTAGCGGTAAAAGCGTGTTACCGCa                     <  1:179305/257‑1 (MQ=255)
                                                                                                              tGAGCCTGATCGGAACTCAGATTAAATTCCTGACGAATTTCCGATTCTCGCTGATGTAACGCATTAACGATGTAAAAAACGGATGAAAAAGCAATCAACAGCCAGAGCACTAACGCCAATGCTCTGAACATGTAg                                                                                                                                      >  1:236499/1‑135 (MQ=255)
                                                                                                                  ccTGATCGGAACTCAGATTAAATTCCTGACGAATTTCCGATTCTCGCTGATGTAACGCATTAACGATGTAAAAAACGGATGAAAAAGCAATCAACAGCCAGAGCACTAACGCCAATGCTCTGAACATGTAGcgcg                                                                                                                                  <  2:236495/135‑1 (MQ=255)
                                                                                                                     gATCGGAACTCAGATTAAATTCCTGACGAATTTCCGATTCTCGCTGATGTAACGCATTAACGATGTAAAAAACGGATGAAAAAGCAATCAACAGCCAGCGCACTAACGCCAATGCTCTGAACATGTAGCGCGAGGCTTTCAGGGTTGTACGAAACGAAGCAAGGTATTTCAAAGGAGCTAAGCTCCGCCTCAGGTGACCGa                                                             >  2:614667/1‑201 (MQ=255)
                                                                                                                         ggAACTCAGATTAAATTCCTGACGAATTTCCGATTCTGGCTGNTGTAACGCATGTACTATGTAAAAAACGGATGAAAAAGCAATTAACAGCCAGAGCACTAACGCCAATGCTCTGAACATGTAGCGCGAGGCTTTCAGGGTTGTACGAAAAGAAGCAAGGTATTTCAAAGGGGCGAAGCTCCGCCTCAGGTGACCGATGGa                                                         <  1:614678/201‑1 (MQ=255)
                                                                                                                                        ttCCTGACGAATTTCCGATTCTCGCTGATGTAACGCATTAACGATGTAAAAAACGGATGAAAAAGCAATCAACAGCCAGAGCACTAACGCCAATGCTCTAAACATGTAGCGCGAGGCTTTCAGGGTTGTACGAAAAGAAGCAAGGTATTTCAAAGGGGCGAAGCTCCGCCTCAGGTGACCGATGGAGTGTGGTTAAGGTAGCGGTAAAAg                                 >  2:565549/1‑210 (MQ=255)
                                                                                                                                            tGACGAATTTCCGATTCTCGCTGATGTAACGCATTAACGATGTAAAAAACGGATGAAAAAGCAATCAACAGCCAGAGCACTAACGCCAATGCTCTAAACATGTAGCGCGAGGCTTTCAGGGTTGTACGAAAAGAAGCAAGGTATTTCAAAGGGGCGAAGCTCCGCCTCAGGTGACCGATGGAGTGTGGTTAAGGTAGCGGTAAAAGCgtg                             <  1:565560/210‑1 (MQ=255)
                                                                                                                                                                                 cGATGTAAAAAACGGATGAAAAAGCAATCAACAGCCAGAGCACTAACGCCAATGCTCTGAACATGTAGCGCGAGGCTTTCAGGGTTGTACGAAAAGAAGCAAGGTATTTCAAAGGGGCGAAGCTCCGCCTCAGGTGACCGATGGAGTGTGGTTAAGGTAGCGGTAAAAGCGTGTTACCGCAATGTTCTCTCTTCTCTGGa  <  1:24499/200‑1 (MQ=255)
                                                                                                                                                                                   aTGTAAAAAACGGATGAAAAAGCAATCAACAGCCAGAGCACTAACGCCAATGCTCTAAACATGTAGCGCGAGGCTTTCAGGGTTGTACGAAAAGAAGCAAGGTATTTCAAAGGGGCGAAGCTCCGCCTCAGGTGa                                                                 >  2:23093/1‑135 (MQ=255)
                                                                                                                                                                                    tGTAAAAAACGGATGAAAAAGCAATCAACAGCCAGAGCACTAACGCCAATGCTCTGAACATGTAGCGCGAGGCTTTCAGGGt                                                                                                                     >  2:174620/1‑82 (MQ=255)
                                                                                                                                                                                     gTAAAAAACGGATGAAAAAGCAATCAACAGCCAGAGCACTAACGCCAATGCTCTGAACATGTa                                                                                                                                       >  1:96159/1‑63 (MQ=255)
                                                                                                                                                                                      tAAAAAACGGATGAAAAAGCAATCAACAGCCAGAGCACTAACGCCAATGCTCTGAACATGTAGCGCGAGGCTTTCAGGGTTGTACGAAAAGAAGc                                                                                                      >  1:303398/1‑95 (MQ=255)
                                                                                                                                                                                       aaaaaaCGGATGAAAAAGCAATCAACAGCCAGAGCACTAACGCCAATGCTCTAAACATGTAGCGCGAGGCTTTCAGGGTTGTACGAAAAGAAGCAAGGTATTTCAAAGGGGCGAAGCTCCGCCTCAGGTGACCGa                                                             <  1:23093/135‑1 (MQ=255)
                                                                                                                                                                                        aaaaaCGGATGAAAAAGCAATCAACAGCCAGAGCACTAACGCCAATGCTCTGAACATGTAGCGCGAGGCTTTCAGGGTTGTa                                                                                                                 <  1:174624/82‑1 (MQ=255)
                                                                                                                                                                                         aaaaCGGATGAAAAAGCAATCAACAGCCAGAGCACTAACGCCAATGCTCTGAACATGTAgcgc                                                                                                                                   <  2:96159/63‑1 (MQ=255)
                                                                                                                                                                                          aaaCGGATGAAAAAGCAATCAACAGCCAGAGCACTAACGCCAATGCTCTGAACATGTAGCGCGAGGCTTTCAGGGTTGTACGAAAAGAAGCAAgg                                                                                                  <  2:303393/95‑1 (MQ=255)
                                                                                                                                                                                                       aGCAATCAACAGCCAGAGCACTAACGCCAATGCTCTGAACATGTAGCGCGAGGCTTTCAGGGTTGTACGAAaaga                                                                                                         >  2:322049/1‑75 (MQ=255)
                                                                                                                                                                                                          aaTCAACAGCCAGAGCACTAACGCCAATGCTCTGAACATGTAGCGCGAGGCTTTCAGGGTTGTACGAAaaga                                                                                                         >  2:642605/1‑72 (MQ=255)
                                                                                                                                                                                                           aTCAACAGCCAGAGCACTAACGCCAATGCTCTGAACATGTAGCGCGAGGCTTTCAGGGTTGTACGAAAAGAAGCa                                                                                                     <  1:322054/75‑1 (MQ=255)
                                                                                                                                                                                                              aaCAGCCAGAGCACTAACGCCAATGCTCTGAACATGTAGCGCGAGGCTTTCAGGGTTGTACGAAAAGAAGCa                                                                                                     <  1:642616/72‑1 (MQ=255)
                                                                                                                                                                                                                   ccAGAGCACTAACGCCAATGCTCTGAACATGTAGCGCGAgg                                                                                                                               >  2:474787/1‑41 (MQ=255)
                                                                                                                                                                                                                    cAGAGCACTAACGCCAATGCTCTGAACATGTAGCGCGAGGCTTTCAGGGTTGTACGAAAAGAAGc                                                                                                      >  2:322639/1‑65 (MQ=255)
                                                                                                                                                                                                                      gagCACTAACGCCAATGCTCTGAACATGTAGCGCGAGGCTTTCAGGGTTGTACGAAAAGAAGCAAGGTATTTCAAAGGGGCGAAGCTCCGCCTCAGGTGACCGATGGAGTGTGGTTAAGGTAGCGGTAAAAGCGTGTTAc                         >  1:139422/1‑140 (MQ=255)
                                                                                                                                                                                                                        gCACTAACGCCAATGCTCTGAACATGTAGCGCGAGGCTTTCAGGGTTGTACGAAAAGAAGCAAgg                                                                                                  <  1:322644/65‑1 (MQ=255)
                                                                                                                                                                                                                          aCTAACGCCAATGCTCTGAACATGTAGCGCGAGGCTTTCAGGGTTGTACGAAAAGAAGCAAGGTATTTCAAAGGGGCGAAGCTCCGCCTCAGGTGACCGATGGAGTGTGGTTAAGGTAGCGGTAAAAGCGTGTTACCGCa                     <  2:139421/140‑1 (MQ=255)
                                                                                                                                                                                                                              aCGCCAATGCTCTGAACATGTAGCGCGAGGCTTTCAGGGTTGTAc                                                                                                                >  1:215183/1‑45 (MQ=255)
                                                                                                                                                                                                                                 ccAATGCTCTGAACATGTAGCGCGAGGCTTTCAGGGTTGTACGAAAAGAAGCAAGGTATTTCAAAGGGGCGAAGCTCCGCCTc                                                                       >  2:533062/1‑83 (MQ=255)
                                                                                                                                                                                                                                  cAATGCTCTGAACATGTAGCGCGAGGCTTTCAGGGTTGTACGaaa                                                                                                            <  2:215179/45‑1 (MQ=255)
                                                                                                                                                                                                                                   aaTGCTCTGAACATGTAGCGCGAGGCTTTCAGGGTTGTACGAAAAGAAGc                                                                                                      >  2:372644/1‑50 (MQ=255)
                                                                                                                                                                                                                                   aaTGCTCTGAACATGTAGCGCGAGGCTTTCAGGGTTGTACGAAAAGAAGCAAGGTATTTCAAAGGGGCGAAGCt                                                                              >  2:264838/1‑74 (MQ=255)
                                                                                                                                                                                                                                     tGCTCTGAACATGTAGCGCGAGGCTTTCAGGGTTGTACGAAAAGAAGCAAGGTATTTCAAAGGGGCGAAGCTCCGCCTCAGGt                                                                   <  1:533073/83‑1 (MQ=255)
                                                                                                                                                                                                                                       ctctGAACATGTAGCGCGAGGCTTTCAGGGTTGTACGAAAAGAAGCAAgg                                                                                                  <  1:372649/50‑1 (MQ=255)
                                                                                                                                                                                                                                       ctctGAACATGTAGCGCGAGGCTTTCAGGGTTGTACGAAAAGAAGCAAGGTATTTCAAAGGGGCGAAGCTCCGc                                                                          <  1:264842/74‑1 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                             
CCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGATCGGAACTCAGATTAAATTCCTGACGAATTTCCGATTCTCGCTGATGTAACGCATTAACGATGTAAAAAACGGATGAAAAAGCAATCAACAGCCAGAGCACTAACGCCAATGCTCTGAACATGTAGCGCGAGGCTTTCAGGGTTGTACGAAAAGAAGCAAGGTATTTCAAAGGGGCGAAGCTCCGCCTCAGGTGACCGATGGAGTGTGGTTAAGGTAGCGGTAAAAGCGTGTTACCGCAATGTTCTCTCTTCTCTGGA  >  CP009273/2313067‑2313443

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