Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,640,981 C→T 30.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,9810CT30.0% 43.2 / 19.3 30A98V (GCT→GTT) creBDNA‑binding response regulator in two‑component regulatory system with CreC
Reads supporting (aligned to +/- strand):  ref base C (14/7);  new base T (5/4);  total (19/11)
Fisher's exact test for biased strand distribution p-value = 6.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCG  >  W3110S.gb/4640916‑4641047
                                                                 |                                                                  
aCCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTgac                                                                >  1:1870899/1‑70 (MQ=255)
aCCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTgac                                                                >  1:3426818/1‑70 (MQ=255)
aCCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTgac                                                                >  1:3284961/1‑70 (MQ=255)
          ttacTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGTTGACGACTACGTGGc                                                     >  1:2994674/4‑71 (MQ=255)
          ttCTTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGTTGACGACTACGTGGc                                                     >  1:3219820/1‑71 (MQ=255)
          ttCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGTTGACGACTACGTGGc                                                     >  1:1096937/1‑71 (MQ=255)
          ttCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGTTGACGACTACGTGGc                                                     >  1:3542796/1‑71 (MQ=255)
          ttCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGTTGACGACTACGTGGc                                                     >  1:2788785/1‑71 (MQ=255)
                 cGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCg                                              <  1:1229226/71‑1 (MQ=255)
                             aaGAGGGCGATCGCCTGCTTGGGCTGGAAATTGGGGTTGACGACTACGTGGCTaa                                                  <  1:369378/55‑1 (MQ=37)
                               gagGTCGATCGCCTGCTTGGGCTGGAAATTGGTGTTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGa                                <  1:1700125/71‑1 (MQ=255)
                                    cGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAgt                             >  1:721428/1‑69 (MQ=255)
                                    cGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAgt                             >  1:1714610/1‑69 (MQ=255)
                                    cGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAgt                             >  1:1942158/1‑69 (MQ=255)
                                        cGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTgcgc                       >  1:1111815/1‑71 (MQ=255)
                                        cGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTgcgc                       >  1:3197208/1‑71 (MQ=255)
                                        cGCCTGCTTGGGCTGGAAATTGATGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTgcgc                       >  1:3456434/1‑71 (MQ=255)
                                          ccTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGcc                      <  1:2378424/70‑1 (MQ=255)
                                          ccTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGcc                      <  1:756971/70‑1 (MQ=255)
                                          ccTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGcc                      <  1:939773/70‑1 (MQ=255)
                                                  ggCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCg                                              <  1:2418152/38‑1 (MQ=255)
                                                  ggCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCg                                              <  1:1709306/38‑1 (MQ=255)
                                                  ggCGGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGGgcgc             <  1:3696425/71‑1 (MQ=255)
                                                   gCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAgtg                            >  1:3590656/1‑55 (MQ=255)
                                                   gCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCGAGGGTGCGCa            >  1:1159185/1‑71 (MQ=255)
                                                       gAAATTGGTGTTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCt         <  1:3098025/70‑1 (MQ=255)
                                                       gAAATTGGTGTTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCt         <  1:3716444/70‑1 (MQ=255)
                                                             ggTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTAc      >  1:3540045/1‑67 (MQ=255)
                                                             ggTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGc   >  1:1547753/1‑70 (MQ=255)
                                                             ggTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCg  >  1:573920/1‑71 (MQ=255)
                                                                 |                                                                  
ACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCG  >  W3110S.gb/4640916‑4641047

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: