Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 3,638,024 (TTG)2→3 coding (675/2163 nt) uvrD ← DNA‑dependent ATPase I and helicase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,638,0161.T76.5% 37.1 / 10.1 17coding (683/2163 nt)uvrDDNA‑dependent ATPase I and helicase II
Reads supporting (aligned to +/- strand):  ref base . (1/3);  new base T (13/0);  total (14/3)
Fisher's exact test for biased strand distribution p-value = 5.88e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*W3110S.gb3,638,0162.G76.5% 43.0 / 8.6 17coding (683/2163 nt)uvrDDNA‑dependent ATPase I and helicase II
Reads supporting (aligned to +/- strand):  ref base . (1/3);  new base G (13/0);  total (14/3)
Fisher's exact test for biased strand distribution p-value = 5.88e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*W3110S.gb3,638,0163.T76.5% 37.1 / 10.1 17coding (683/2163 nt)uvrDDNA‑dependent ATPase I and helicase II
Reads supporting (aligned to +/- strand):  ref base . (1/3);  new base T (13/0);  total (14/3)
Fisher's exact test for biased strand distribution p-value = 5.88e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

CCGACGATCATCACTTTGCCGGTGTCGCCCGCCAGCAGGCGGATCCACGCGTACTGAA‑‑‑TGTTGTTGGTATCCTGGAATTCGTCCACCAGGATATTGGTAAAACGTTCGCGGTAGTG  >  W3110S.gb/3637959‑3638074
                                                          |||                                                          
ccGACGATCATCACTTTGCCGGTGTCGCCCGCCAGCAGGCGGATCCACGCGTACTGAA‑‑‑tgtt                                                        >  1:116064/1‑62 (MQ=255)
ccGACGATCATCACTTTGCCGGTGTCGCCCGCCAGCAGGCGGATCCACGCGTACTGAA‑‑‑tgtt                                                        >  1:190329/1‑62 (MQ=255)
ccGACGATCATCACTTTGCCGGTGTCGCCCGCCAGCAGGCGGATCCACGCGTACTGAA‑‑‑tgtt                                                        >  1:2136580/1‑62 (MQ=255)
ccGACGATCATCACTTTGCCGGTGTCGCCCGCCAGCAGGCGGATCCACGCGTACTGAA‑‑‑tgtt                                                        >  1:2271070/1‑62 (MQ=255)
ccGACGATCATCACTTTGCCGGTGTCGCCCGCCAGCAGGCGGATCCACGCGTACTGAA‑‑‑tgtt                                                        >  1:360810/1‑62 (MQ=255)
ccGACGATCATCACTTTGCCGGTGTCGCCCGCCAGCAGGCGGATCCACGCGTACTGAA‑‑‑tgtt                                                        >  1:320214/1‑62 (MQ=255)
   acgaTCATCACTTTGCCGGTGTCGCCCGCCAGCAGGCGGATCCACGCGTACTGAA‑‑‑tgttgt                                                      >  1:2793546/1‑61 (MQ=255)
   acgaTCATCACTTTGCCGGTGTCGCCCGCCAGCAGGCGGATCCACGCGTACTGAA‑‑‑Tgttgtt                                                     >  1:1149360/1‑62 (MQ=255)
                   cGGTGTCGCCCGCCAGCAGGCGGATCCACGCGTACTGAA‑‑‑TGTTGTTGGTATCCTGGAATTc                                      <  1:38277/61‑1 (MQ=255)
                      tgtCGCCCGCCAGCAGGCGGATCCACGCGTACTGAATGTTGTTGTTGGTATCCTGGAATTc                                      >  1:662460/1‑61 (MQ=255)
                      tgtCGCCCGCCAGCAGGCGGATCCACGCGTACTGAATGTTGTTGTTGGTATCCTGGAATTCg                                     >  1:977433/1‑62 (MQ=255)
                      tgtCGCCCGCCAGCAGGCGGATCCACGCGTACTGAATGTTGTTGTTGGTATCCTGGAATTCg                                     >  1:1114678/1‑62 (MQ=255)
                      tgtCGCCCGCCAGCAGGCGGATCCACGCGTACTGAATGTTGTTGTTGGTATCCTGGAATTCg                                     >  1:8856/1‑62 (MQ=255)
                      tgtCGCCCGCCAGCAGGCGGATCCACGCGTACTGAATGTTGTTGTTGGTATCCTGGAATTCg                                     >  1:639246/1‑62 (MQ=255)
                      tgtCGCCCGCCAGCAGGCGGATCCACGCGTACTGAATGTTGTTGTTGGTATCCTGGAATTCg                                     >  1:489321/1‑62 (MQ=255)
                      tgtCGCCCGCCAGCAGGCGGATCCACGCGTACTGAATGTTGTTGTTGGTATCCTGGAATTCg                                     >  1:3084200/1‑62 (MQ=255)
                      tgtCGCCCGCCAGCAGGCGGATCCACGCGTACTGAATGTTGTTGTTGGTATCCTGGAATTCg                                     >  1:2790188/1‑62 (MQ=255)
                      tgtCGCCCGCCAGCAGGCGGATCCACGCGTACTGAATGTTGTTGTTGGTATCCTGGAATTCg                                     >  1:2176959/1‑62 (MQ=255)
                      tgtCGCCCGCCAGCAGGCGGATCCACGCGTACTGAATGTTGTTGTTGGTATCCTGGAATTCg                                     >  1:2140785/1‑62 (MQ=255)
                      tgtCGCCCGCCAGCAGGCGGATCCACGCGTACTGAATGTTGTTGTTGGTATCCTGGAATTCg                                     >  1:1980365/1‑62 (MQ=255)
                      tgtCGCCCGCCAGCAGGCGGATCCACGCGTACTGAATGTTGTTGTTGGTATCCTGGAATTCg                                     >  1:1938074/1‑62 (MQ=255)
                      tgtCGCCCGCCAGCAGGCGGATCCACGCGTACTGAATGTTGTTGTTGGTATCCTGGAATTCg                                     >  1:1557220/1‑62 (MQ=255)
                                                       gAA‑‑‑TGTTGTTGGTATCCTGGAATTCGTCCACCAGGATATTGGTAAAACGTTCGCGGTAgtg  <  1:3060657/61‑1 (MQ=255)
                                                       gAA‑‑‑TGTTGTTGGTATCCTGGAATTCGTCCACCAGGATATTGGTAAAACGTTCGCGGTAgtg  <  1:1755507/61‑1 (MQ=255)
                                                          |||                                                          
CCGACGATCATCACTTTGCCGGTGTCGCCCGCCAGCAGGCGGATCCACGCGTACTGAA‑‑‑TGTTGTTGGTATCCTGGAATTCGTCCACCAGGATATTGGTAAAACGTTCGCGGTAGTG  >  W3110S.gb/3637959‑3638074

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