Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,641,490:1 +G 23.7% coding (113/1425 nt) creC → sensory histidine kinase in two‑component regulatory system with CreB or PhoB, regulator of the CreBC regulon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,641,4901.G23.7% 57.3 / 43.2 38coding (113/1425 nt)creCsensory histidine kinase in two‑component regulatory system with CreB or PhoB, regulator of the CreBC regulon
Reads supporting (aligned to +/- strand):  ref base . (13/16);  new base G (3/6);  total (16/22)
Fisher's exact test for biased strand distribution p-value = 7.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGG  >  W3110S.gb/4641423‑4641557
                                                                    |                                                                   
gcggtggcAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACg                                                                    >  1:3014212/1‑69 (MQ=255)
      gcAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACgg                                                                   >  1:2476141/1‑64 (MQ=255)
          ccTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGGAGGGGACGtt                                                         <  1:2764283/71‑1 (MQ=255)
           cTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGGAGGGGACGTTg                                                        <  1:702707/71‑1 (MQ=255)
                tCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGa                                                   >  1:1654089/1‑70 (MQ=255)
                 cGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGac                                                  <  1:3495013/70‑1 (MQ=255)
                  gTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACgg                                                                   <  1:96093/52‑1 (MQ=255)
                     cTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACggga                                                                 >  1:3113445/1‑49 (MQ=255)
                       ggCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGaca                                                 >  1:4236610/1‑65 (MQ=255)
                           aTTTTTGTCAAAGAAGTTAACCCGGGCGTGCGAAGAGCAACGGGAGGGGACGTTGATc                                                     <  1:1166116/58‑1 (MQ=255)
                           aTTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGtt                                       <  1:119610/71‑1 (MQ=255)
                            tttttGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACg                                                           <  1:4442759/50‑1 (MQ=255)
                             ttttGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATc                                                     <  1:2703494/55‑1 (MQ=255)
                             ttttGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGc                                     >  1:3373216/1‑71 (MQ=255)
                              tttGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAAc                                                                      >  1:262731/1‑38 (MQ=255)
                              tttGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGGAGGGGACGTTGAt                                                      >  1:3995661/1‑54 (MQ=255)
                              tttGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGa                                                                  >  1:4109331/1‑41 (MQ=255)
                              tttGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAAc                                          >  1:2855445/1‑65 (MQ=255)
                              tttGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCt                                    >  1:2405645/1‑71 (MQ=255)
                              tttGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCt                                    >  1:2847228/1‑71 (MQ=255)
                              tttGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCt                                    >  1:4344218/1‑71 (MQ=255)
                               ttGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGGAGGGGACGTTGAt                                                      >  1:1441167/1‑53 (MQ=255)
                                    aaaGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGa                                                   >  1:4317005/1‑50 (MQ=255)
                                     aagaagTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTgg                                  <  1:1535266/66‑1 (MQ=255)
                                       gaagTTAAACCGGGCGTGCGAAGAGCAACGGGAGGGGACGTTGATCGa                                                   <  1:3529456/48‑1 (MQ=255)
                                        aagTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTg                          <  1:526347/71‑1 (MQ=255)
                                          gTTAAACCGGGCGTGCGAAGAGCAACGGGAGGGGACGTTGATCGACACCGCAACGTTGCTgg                                  <  1:1582320/62‑1 (MQ=255)
                                          gTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGgc                        <  1:3209038/71‑1 (MQ=255)
                                           ttAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCgg                               <  1:389701/63‑1 (MQ=255)
                                            tAAACCGGGCGTGCGAAGAGCAACGGGAGGGGACGTTGATCGacac                                                >  1:3327565/1‑46 (MQ=255)
                                             aaaCCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGt                                        >  1:973285/1‑52 (MQ=255)
                                             aaaCCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGcgcg                     >  1:4533478/1‑71 (MQ=255)
                                              aaCCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGAAACCGCAACGTTGCTGGCGGAGCTGGcgcg                     <  1:1991658/70‑1 (MQ=255)
                                                 cGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACCCCGCAACGTTGCTGGCGGAGCTGGCGCGt                    >  1:1336774/1‑68 (MQ=255)
                                                   ggCGTGCGAAGAGCAACGGAGTGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGa               <  1:2126260/71‑1 (MQ=255)
                                                      gTGCGAAGAGCAACGGGAGGGGACGTTGATCGACACCGCAACGTTGCTgg                                  <  1:3814844/50‑1 (MQ=255)
                                                       tGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTg           <  1:3415274/71‑1 (MQ=255)
                                                         cGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGct         <  1:2191440/71‑1 (MQ=255)
                                                             gagCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGctctct     <  1:777470/71‑1 (MQ=255)
                                                              agCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCt                                    <  1:254093/39‑1 (MQ=255)
                                                               gCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTg    <  1:2661454/70‑1 (MQ=255)
                                                                cAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTggg  <  1:947520/71‑1 (MQ=255)
                                                                cAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTggg  >  1:1279666/1‑71 (MQ=255)
                                                                    |                                                                   
GCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGG  >  W3110S.gb/4641423‑4641557

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