Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 3,035,419 C→G R448R (CGG→CGC recJ ← ssDNA exonuclease, 5' ‑‑> 3'‑specific

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,035,4190CG95.6% 79.0 / ‑3.6 45R448R (CGG→CGCrecJssDNA exonuclease, 5' ‑‑> 3'‑specific
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base G (0/43);  total (0/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.13e-01

AACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTATCCTCTTCCAGCGACAAACCCGCCGC  >  W3110S.gb/3035404‑3035469
               |                                                  
aaCCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTATCCTCTTCCAGCGACAAACccg      <  1:45682/62‑1 (MQ=255)
 aCCAGTTCGCCAAACCGTTGTTGAAAGAGGTTGAATTTATCCTCTTCCAGCGACAAACccg      <  1:45488/61‑1 (MQ=255)
    aGTTCGCCAAAGCTTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:71122/62‑1 (MQ=255)
    aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:64338/62‑1 (MQ=255)
    aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:47136/62‑1 (MQ=255)
    aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:49461/62‑1 (MQ=255)
    aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:52705/62‑1 (MQ=255)
    aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:53906/62‑1 (MQ=255)
    aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:54901/62‑1 (MQ=255)
    aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:55338/62‑1 (MQ=255)
    aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:59781/62‑1 (MQ=255)
    aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:60056/62‑1 (MQ=255)
    aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:61847/62‑1 (MQ=255)
    aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:46102/62‑1 (MQ=255)
    aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:66779/62‑1 (MQ=255)
    aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:67629/62‑1 (MQ=255)
    aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:68773/62‑1 (MQ=255)
    aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:75144/62‑1 (MQ=255)
    aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:79044/62‑1 (MQ=255)
    aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:79103/62‑1 (MQ=255)
    aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:79776/62‑1 (MQ=255)
    aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:81561/62‑1 (MQ=255)
    aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:81876/62‑1 (MQ=255)
    aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:30144/62‑1 (MQ=255)
    aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:11911/62‑1 (MQ=255)
    aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:13097/62‑1 (MQ=255)
    aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:18141/62‑1 (MQ=255)
    aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:18331/62‑1 (MQ=255)
    aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:19119/62‑1 (MQ=255)
    aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:19779/62‑1 (MQ=255)
    aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:2086/62‑1 (MQ=255)
    aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:22305/62‑1 (MQ=255)
    aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:27628/62‑1 (MQ=255)
    aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:46519/62‑1 (MQ=255)
    aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:35288/62‑1 (MQ=255)
    aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:37777/62‑1 (MQ=255)
    aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:39286/62‑1 (MQ=255)
    aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:39482/62‑1 (MQ=255)
    aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:41195/62‑1 (MQ=255)
    aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:41347/62‑1 (MQ=255)
    aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:45858/62‑1 (MQ=255)
    aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:10180/62‑1 (MQ=255)
     gTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:35421/61‑1 (MQ=255)
     gTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:79214/61‑1 (MQ=255)
     gTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:4149/61‑1 (MQ=255)
               |                                                  
AACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTATCCTCTTCCAGCGACAAACCCGCCGC  >  W3110S.gb/3035404‑3035469

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: