Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 2,267,040 +A intergenic (‑211/‑157) fruB ← / → setB fused fructose‑specific PTS enzyme IIA component and HPr component/lactose/glucose efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,267,0361.A90.0% 124.1 / 8.2 40intergenic (‑207/‑161)fruB/setBfused fructose‑specific PTS enzyme IIA component and HPr component/lactose/glucose efflux system
Reads supporting (aligned to +/- strand):  ref base . (0/4);  new base A (0/36);  total (0/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.01e-01

TGCTGATGGGAAAATCCGTTATCAGCAATTTCATTTCAGCAGCGTCAGCAGCCAATCAGCAACA  >  W3110S.gb/2267027‑2267090
          |                                                      
tGCTGATGGGAAAAATCCGTTATCAGCAATTTCATTTCAGCAGCGTCAGCAGCCAATCAGCAACa  <  1:631288/65‑1 (MQ=255)
tGCTGATGGGAAAAATCCGTTATCAGCAATTTCATTTCAGCAGCGTCAGCAGCCAATCAGCAACa  <  1:2257651/65‑1 (MQ=255)
tGCTGATGGGAAAAATCCGTTATCAGCAATTTCATTTCAGCAGCGTCAGCAGCCAATCAGCAACa  <  1:1079933/65‑1 (MQ=255)
tGCTGATGGGAAAAATCCGTTATCAGCAATTTCATTTCAGCAGCGTCAGCAGCCAATCAGCAACa  <  1:2359429/65‑1 (MQ=255)
tGCTGATGGGAAAAATCCGTTATCAGCAATTTCATTTCAGCAGCGTCAGCAGCCAATCAGCAACa  <  1:2496162/65‑1 (MQ=255)
tGCTGATGGGAAAAATCCGTTATCAGCAATTTCATTTCAGCAGCGTCAGCAGCCAATCAGCAACa  <  1:2530944/65‑1 (MQ=255)
tGCTGATGGGAAAAATCCGTTATCAGCAATTTCATTTCAGCAGCGTCAGCAGCCAATCAGCAACa  <  1:302723/65‑1 (MQ=255)
tGCTGATGGGAAAAATCCGTTATCAGCAATTTCATTTCAGCAGCGTCAGCAGCCAATCAGCAACa  <  1:370353/65‑1 (MQ=255)
tGCTGATGGGAAAAATCCGTTATCAGCAATTTCATTTCAGCAGCGTCAGCAGCCAATCAGCAACa  <  1:517176/65‑1 (MQ=255)
tGCTGATGGGAAAAATCCGTTATCAGCAATTTCATTTCAGCAGCGTCAGCAGCCAATCAGCAACa  <  1:2356184/65‑1 (MQ=255)
tGCTGATGGGAAAAATCCGTTATCAGCAATTTCATTTCAGCAGCGTCAGCAGCCAATCAGCAACa  <  1:65512/65‑1 (MQ=255)
tGCTGATGGGAAAAATCCGTTATCAGCAATTTCATTTCAGCAGCGTCAGCAGCCAATCAGCAACa  <  1:680973/65‑1 (MQ=255)
tGCTGATGGGAAAAATCCGTTATCAGCAATTTCATTTCAGCAGCGTCAGCAGCCAATCAGCAACa  <  1:681836/65‑1 (MQ=255)
tGCTGATGGGAAAAATCCGTTATCAGCAATTTCATTTCAGCAGCGTCAGCAGCCAATCAGCAACa  <  1:790499/65‑1 (MQ=255)
tGCTGATGGGAAAAATCCGTTATCAGCAATTTCATTTCAGCAGCGTCAGCAGCCAATCAGCAACa  <  1:799337/65‑1 (MQ=255)
tGCTGATGGGAAAAATCCGTTATCAGCAATTTCATTTCAGCAGCGTCAGCAGCCAATCAGCAACa  <  1:860174/65‑1 (MQ=255)
tGCTGATGGGAAAAATCCGTTATCAGCAATTTCATTTCAGCAGCGTCAGCAGCCAATCAGCAACa  <  1:912925/65‑1 (MQ=255)
tGCTGATGGGAAAAATCCGTTATCAGCAATTTCATTTCAGCAGCGTCAGCAGCCAATCAGCAACa  <  1:931786/65‑1 (MQ=255)
tGCTGATGGGAAAAATCCGTTATCAGCAATTTCATTTCAGCAGCGTCAGCAGCCAATCAGCAACa  <  1:1822559/65‑1 (MQ=255)
tGCTGATGGGAAAAATCCGTTATCAGCAATTTCATTTCAGCAGCGTCAGCAGCCAATCAGCAACa  <  1:1145098/65‑1 (MQ=255)
tGCTGATGGGAAAAATCCGTTATCAGCAATTTCATTTCAGCAGCGTCAGCAGCCAATCAGCAACa  <  1:13456/65‑1 (MQ=255)
tGCTGATGGGAAAAATCCGTTATCAGCAATTTCATTTCAGCAGCGTCAGCAGCCAATCAGCAACa  <  1:1598593/65‑1 (MQ=255)
tGCTGATGGGAAAAATCCGTTATCAGCAATTTCATTTCAGCAGCGTCAGCAGCCAATCAGCAACa  <  1:1637727/65‑1 (MQ=255)
tGCTGATGGGAAAAATCCGTTATCAGCAATTTCATTTCAGCAGCGTCAGCAGCCAATCAGCAACa  <  1:1738318/65‑1 (MQ=255)
tGCTGATGGGAAAAATCCGTTATCAGCAATTTCATTTCAGCAGCGTCAGCAGCCAATCAGCAACa  <  1:1752252/65‑1 (MQ=255)
tGCTGATGGGAAAAATCCGTTATCAGCAATTTCATTTCAGCAGCGTCAGCAGCCAATCAGCAACa  <  1:2204397/65‑1 (MQ=255)
tGCTGATGGGAAAAATCCGTTATCAGCAATTTCATTTCAGCAGCGTCAGCAGCCAATCAGCAACa  <  1:1834379/65‑1 (MQ=255)
tGCTGATGGGAAAAATCCGTTATCAGCAATTTCATTTCAGCAGCGTCAGCAGCCAATCAGCAACa  <  1:1867115/65‑1 (MQ=255)
tGCTGATGGGAAAAATCCGTTATCAGCAATTTCATTTCAGCAGCGTCAGCAGCCAATCAGCAACa  <  1:1894247/65‑1 (MQ=255)
tGCTGATGGGAAAAATCCGTTATCAGCAATTTCATTTCAGCAGCGTCAGCAGCCAATCAGCAACa  <  1:2083674/65‑1 (MQ=255)
tGCTGATGGGAAAAATCCGTTATCAGCAATTTCATTTCAGCAGCGTCAGCAGCCAATCAGCAACa  <  1:2097318/65‑1 (MQ=255)
tGCTGATGGGAAAAATCCGTTATCAGCAATTTCATTTCAGCAGCGTCAGCAGCCAATCAGCAACa  <  1:2195773/65‑1 (MQ=255)
tGCTGATGGGAAAAATCCGTTATCAGCAATTTCATTTCAGCAGCGTCAGCAGCAAATCAGCAACa  <  1:747265/65‑1 (MQ=255)
tGCTGATGGGAAAATCCGTTATCAGCAATTTCATTTCAGCAGCGTCAGCAGCCAATCAGCAACa  <  1:279160/64‑1 (MQ=255)
tGCTGATGGGAAAATCCGTTATCAGCAATTTCATTTCAGCAGCGTCAGCAGCCAATCAGCAACa  <  1:1401359/64‑1 (MQ=255)
tGCTGATGGGAAAATCCGTTATCAGCAATTTCATTTCAGCAGCGTCAGCAGCCAATCAGCAACa  <  1:1333603/64‑1 (MQ=255)
tGCTGATGGGAAAATCCGTTATCAGCAATTTCATTTCAGCAGCGTCAGCAGCCAATCAGCAACa  <  1:2275960/64‑1 (MQ=255)
 gCTGATGGGAAAAATCCGTTATCAGCAATTTCATTTCAGCAGCGTCAGCAGCCTATCAGCAACa  <  1:467411/64‑1 (MQ=255)
 gCTGATGGGAAAAATCCGTTATCAGCAATTTCATTTCAGCAGCGTCAGCAGCCAATCAGCAACa  <  1:2379354/64‑1 (MQ=255)
    gATGGGAAAAATCCGTTATCAGCAATTTCATTTCAGCAGCGTCAGCAGCCAATCAGCAACa  <  1:1381059/61‑1 (MQ=255)
          |                                                      
TGCTGATGGGAAAATCCGTTATCAGCAATTTCATTTCAGCAGCGTCAGCAGCCAATCAGCAACA  >  W3110S.gb/2267027‑2267090

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