Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,640,981 C→T 31.0% A98V (GCT→GTT)  creB → DNA‑binding response regulator in two‑component regulatory system with CreC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,640,9810CT31.0% 37.2 / 21.0 29A98V (GCT→GTT) creBDNA‑binding response regulator in two‑component regulatory system with CreC
Reads supporting (aligned to +/- strand):  ref base C (14/6);  new base T (2/7);  total (16/13)
Fisher's exact test for biased strand distribution p-value = 4.05e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.83e-01

TACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGC  >  W3110S.gb/4640915‑4641046
                                                                  |                                                                 
tACCTGTACTGTTCCTGACGTCCCTAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGTTg                                                                 <  1:2316448/69‑1 (MQ=37)
 aCCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCt                                                                  >  1:3295007/1‑67 (MQ=255)
     gTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGTTgacgac                                                            >  1:3591256/1‑69 (MQ=255)
      tACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACt                                                           >  1:3029448/1‑69 (MQ=255)
             ccTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGTTGACGACTACGTgg                                                     >  1:2964336/1‑68 (MQ=255)
                gACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGTTGACGACTACGTGGCTa                                                  <  1:929493/68‑1 (MQ=255)
                gACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTaa                                                 <  1:1472542/69‑1 (MQ=255)
                         aaGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGTTGACGACTACGTGGCTAAACCGTTTTc                                        <  1:350369/69‑1 (MQ=255)
                           gTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACt                                                           >  1:399331/1‑48 (MQ=255)
                           gTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGtttt                                         >  1:3432153/1‑66 (MQ=255)
                           gTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGtttt                                         >  1:1054346/1‑66 (MQ=255)
                           gTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGtttt                                         >  1:882896/1‑66 (MQ=255)
                           gTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGtttt                                         >  1:738610/1‑66 (MQ=255)
                           gTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGtttt                                         >  1:3537710/1‑66 (MQ=255)
                              aaGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGtttt                                         >  1:179617/1‑63 (MQ=255)
                                   gTCGATCGCCTGCTTGGGCTGGAAATTGGTGTTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGaa                              <  1:2958831/69‑1 (MQ=255)
                                     cGATCGCCTGCTTGGGCTGGAAATTGGTGTTGACGACTACGTGGCTAAAcc                                              <  1:660893/51‑1 (MQ=255)
                                            cTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGcc                     <  1:2270928/69‑1 (MQ=255)
                                            cTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGcc                     <  1:446118/69‑1 (MQ=255)
                                             tGCTTGGGCTGGAAATTGGTGTTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCa                    <  1:3542345/69‑1 (MQ=255)
                                              gCTTGGGCTGGAAATTGGTGTTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCa                    <  1:2992785/68‑1 (MQ=255)
                                               cTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAgg                  >  1:2185718/1‑69 (MQ=255)
                                                  gggCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGcc                     >  1:3281123/1‑63 (MQ=255)
                                                  gggCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGcc                     >  1:2103832/1‑63 (MQ=255)
                                                  gggCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGcc                     >  1:2693791/1‑63 (MQ=255)
                                                    gCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTgcg             >  1:3719970/1‑69 (MQ=255)
                                                    gCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTg               <  1:1493770/67‑1 (MQ=255)
                                                      tGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCa           <  1:2428101/69‑1 (MQ=255)
                                                                tGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGc  <  1:2562566/68‑1 (MQ=255)
                                                                  |                                                                 
TACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGC  >  W3110S.gb/4640915‑4641046

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