Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 4,269,081 2 bp→CC coding (1176‑1177/1416 nt) dnaB → replicative DNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,269,0810TC77.8% 31.5 / 4.3 18R392R (CGT→CGCdnaBreplicative DNA helicase
Reads supporting (aligned to +/- strand):  ref base T (0/4);  new base C (0/14);  total (0/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*W3110S.gb4,269,0820GC77.8% 33.7 / 4.3 18A393P (GCC→CCC) dnaBreplicative DNA helicase
Reads supporting (aligned to +/- strand):  ref base G (0/4);  new base C (0/14);  total (0/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAAGAACTGAACGTGCCGGTGGTGGCGCTGTCCCAGTTGAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTAT  >  W3110S.gb/4269021‑4269128
                                                            ||                                              
gAAAGAACTGAACGTGCCGGTGGTGGCGCTGTCCCAGTTGAACCGTTCTCTGGAACAACGt                                                 <  1:2879612/61‑1 (MQ=255)
      aCTGAACGTGCCGGTGGTGGCGCTGTCCCAGTTGAACCGTTCTCTGGAACAACGt                                                 <  1:2252961/55‑1 (MQ=255)
          aaCGTGCCGGTGGTGGCGCTGTCCCAGTTGAACCGTTCTCTGGAACAACGt                                                 <  1:949895/51‑1 (MQ=255)
          aaCGTGCCGGTGGTGGCGCTGTCCCAGTTGAACCGTTCTCTGGAACAACGt                                                 <  1:136822/51‑1 (MQ=255)
          aaCGTGCCGGTGGTGGCGCTGTCCCAGTTGAACCGTTCTCTGGAACAACGt                                                 <  1:137396/51‑1 (MQ=255)
          aaCGTGCCGGTGGTGGCGCTGTCCCAGTTGAACCGTTCTCTGGAACAACGt                                                 <  1:2842087/51‑1 (MQ=255)
          aaCGTGCCGGTGGTGGCGCTGTCCCAGTTGAACCGTTCTCTGGAACAACGt                                                 <  1:1666900/51‑1 (MQ=255)
          aaCGTGCCGGTGGTGGCGCTGTCCCAGTTGAACCGTTCTCTGGAACAACGt                                                 <  1:1794319/51‑1 (MQ=255)
          aaCGTGCCGGTGGTGGCGCTGTCCCAGTTGAACCGTTCTCTGGAACAACGt                                                 <  1:2754677/51‑1 (MQ=255)
          aaCGTGCCGGTGGTGGCGCTGTCCCAGTTGAACCGTTCTCTGGAACAACGt                                                 <  1:2503407/51‑1 (MQ=255)
          aaCGTGCCGGTGGTGGCGCTGTCCCAGTTGAACCGTTCTCTGGAACAACGt                                                 <  1:1891975/51‑1 (MQ=255)
          aaCGTGCCGGTGGTGGCGCTGTCCCAGTTGAACCGTTCTCTGGAACAACGt                                                 <  1:2471017/51‑1 (MQ=255)
          aaCGTGCCGGTGGTGGCGCTGTCCCAGTTGAACCGTTCTCTGGAACAACGt                                                 <  1:2400806/51‑1 (MQ=255)
           aCGTGCCGGTGGTGGCGCTGTCCCAGTTGAACCGTTCTCTGGAACAACGt                                                 <  1:2907715/50‑1 (MQ=255)
                        ggCGCTGTCCCAGTTGAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTcc                         <  1:1232897/61‑1 (MQ=255)
                        ggCGCTGTCCCAGTTGAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTcc                         <  1:218905/61‑1 (MQ=255)
                        ggCGCTGTCCCAGTTGAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTcc                         <  1:1798107/61‑1 (MQ=255)
                        ggCGCTGGCCCAGTTGAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTcc                         <  1:2718774/61‑1 (MQ=255)
                                              tcgctGGAACAACGCCCCAACAAACGTCCCGTGATGTCCGACCTGCGTGAATCCGGCGCTAt  <  1:967102/59‑1 (MQ=255)
                                              tcgctGGAACAACGCCCCAACAAACGTCCCGTGATGTCCGACCTGCGTGAATCCGGCGCTAt  <  1:1156263/59‑1 (MQ=255)
                                              tcgctGGAACAACGCCCCAACAAACGTCCCGTGATGTCCGACCTGCGTGAATCCGGCGCTAt  <  1:940666/59‑1 (MQ=255)
                                              tcgctGGAACAACGCCCCAACAAACGTCCCGTGATGTCCGACCTGCGTGAATCCGGCGCTAt  <  1:1339546/59‑1 (MQ=255)
                                              tcgctGGAACAACGCCCCAACAAACGTCCCGTGATGTCCGACCTGCGTGAATCCGGCGCTAt  <  1:29036/59‑1 (MQ=255)
                                              tcgctGGAACAACGCCCCAACAAACGTCCCGTGATGTCCGACCTGCGTGAATCCGGCGCTAt  <  1:1452226/59‑1 (MQ=255)
                                              tcgctGGAACAACGCCCCAACAAACGTCCCGTGATGTCCGACCTGCGTGAATCCGGCGCTAt  <  1:1537191/59‑1 (MQ=255)
                                              tcgctGGAACAACGCCCCAACAAACGTCCCGTGATGTCCGACCTGCGTGAATCCGGCGCTAt  <  1:2702613/59‑1 (MQ=255)
                                              tcgctGGAACAACGCCCCAACAAACGTCCCGTGATGTCCGACCTGCGTGAATCCGGCGCTAt  <  1:2657334/59‑1 (MQ=255)
                                              tcgctGGAACAACGCCCCAACAAACGTCCCGTGATGTCCGACCTGCGTGAATCCGGCGCTAt  <  1:1857772/59‑1 (MQ=255)
                                              tcgctGGAACAACGCCCCAACAAACGTCCCGTGATGTCCGACCTGCGTGAATCCGGCGCTAt  <  1:2481426/59‑1 (MQ=255)
                                              tcgctGGAACAACGCCCCAACAAACGTCCCGTGATGTCCGACCTGCGTGAATCCGGCGCTAt  <  1:2057201/59‑1 (MQ=255)
                                              tcgctGGAACAACGCCCCAACAAACGTCCCGTGATGTCCGACCTGCGTGAATCCGGCGCTAt  <  1:2059971/59‑1 (MQ=255)
                                              tcgctGGAACAACGCCCCAACAAACGTCCCGTGATGTCCGACCTGCGTGAATCCGGCGCTAt  <  1:2072536/59‑1 (MQ=255)
                                                            ||                                              
GAAAGAACTGAACGTGCCGGTGGTGGCGCTGTCCCAGTTGAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTAT  >  W3110S.gb/4269021‑4269128

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: