Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,641,490:1 +G 40.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.G40.7% ‑3.6 / 51.9 27coding (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/3);  new base G (3/8);  total (16/11)
Fisher's exact test for biased strand distribution p-value = 1.47e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTG  >  W3110S.gb/4641428‑4641555
                                                               |                                                                 
ggcAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGa                                                           >  1:1533217/1‑71 (MQ=255)
ggcAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGa                                                           >  1:3377045/1‑71 (MQ=255)
      cTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGa                                                                >  1:3588057/1‑60 (MQ=255)
         gTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATc                                                   <  1:128026/70‑1 (MQ=255)
           tCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGa                                                 >  1:3819083/1‑70 (MQ=255)
                cTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGGAGGGGACGTTGATCGaca                                               >  1:4074718/1‑68 (MQ=255)
                cTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATc                                                   >  1:3237603/1‑63 (MQ=255)
                   gCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGacac                                              >  1:1311115/1‑65 (MQ=255)
                   gCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCaa                                         <  1:3591925/70‑1 (MQ=255)
                    ccATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGGATGGGACGTTGATCGACACCGCAAc                                        <  1:3820502/71‑1 (MQ=255)
                     cATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGacac                                              >  1:1274482/1‑63 (MQ=255)
                     cATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGt                                      <  1:2003258/71‑1 (MQ=255)
                      aTTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGGAGGGGACGTTGATCGACACCGCAACGt                                      <  1:2631186/71‑1 (MQ=255)
                      aTTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGGAGGGGACGTTGATCGACACCGCAACGt                                      <  1:3958926/71‑1 (MQ=255)
                       tttttGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGGAGGGGACGTTGATCGACACCGCAACGtt                                     <  1:3250910/71‑1 (MQ=255)
                       tttttGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGGAGGGGACGTTGATCGACACCGCAACGtt                                     <  1:3538638/71‑1 (MQ=255)
                         tttGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAg                                                               >  1:1206718/1‑42 (MQ=255)
                         tttGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCt                                  >  1:2697373/1‑71 (MQ=255)
                         tttGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCt                                  >  1:3154250/1‑71 (MQ=255)
                                       tAAACCGGGCGTGCGAAGAGCAACGGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGgcg                     >  1:2626561/1‑71 (MQ=255)
                                       tAAACCGGGCGTGCGAAGAGCAACGGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGgcg                     >  1:1504940/1‑71 (MQ=255)
                                                cGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGAtt           >  1:2586410/1‑71 (MQ=255)
                                                cGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGAtt           >  1:2084811/1‑71 (MQ=255)
                                                      aaGAGCAACGGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGct       <  1:3896003/70‑1 (MQ=255)
                                                        gagCAACGGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGctctc    <  1:3346890/71‑1 (MQ=255)
                                                         agCAACGGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGct       <  1:3185380/67‑1 (MQ=255)
                                                         agCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTg  >  1:406042/1‑71 (MQ=255)
                                                               |                                                                 
GGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTG  >  W3110S.gb/4641428‑4641555

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