Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | REL606 | 3,157,725 | C→A | 13.0% | P759Q (CCA→CAA) | ebgA → | cryptic beta‑D‑galactosidase, alpha subunit |
Read alignment evidence... | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|
seq id | position | change | freq | score | reads | annotation | genes | product | ||
* | REL606 | 3,157,725 | 0 | C→A | 13.0% | 4.8 ‑4.5 ‑0.0 | 51 | 2276 | ebgA | cryptic beta‑D‑galactosidase, alpha subunit |
Reads supporting (aligned to +/- strand): new base (0/7): ref base (24/20): total (24/27) | ||||||||||
Fisher's exact test strand distribution p-value = 1.07e-02 | ||||||||||
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.75e-01 |
GGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACC > REL606/3157689‑3157760
|
ggCAGGTGAATGGCG.ATTGCTTCTTACTCGCGGGcc > 1:1605177/1‑36 (MQ=9)
cAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCa < 1:2481884/36‑1 (MQ=255)
cAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCa < 1:2795640/36‑1 (MQ=255)
aGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCaa > 1:4496561/1‑36 (MQ=255)
ggTGAATGGCGAATCGCTGCTGACTCGCGAGCCaaa < 1:637020/36‑1 (MQ=255)
ggTGAATGGCGAATCGCTGCTGACTCGCGAGCCaaa < 1:4821933/36‑1 (MQ=255)
ggTGAATGGCGAATCGCTGCTGACTCGCGAGCCaaa > 1:1939930/1‑36 (MQ=255)
ggTGAATGGCGAATCGCTGCTGACTCGCGAGCCaaa < 1:7615663/36‑1 (MQ=255)
gTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAg > 1:5027352/1‑36 (MQ=255)
gTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAg > 1:3883333/1‑36 (MQ=255)
tGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGa < 1:6802429/36‑1 (MQ=255)
tGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGa > 1:3264558/1‑36 (MQ=255)
tGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGa > 1:2138718/1‑36 (MQ=255)
tGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGa < 1:3095192/36‑1 (MQ=255)
aaTGGCGAATCGCTGCTGACTCGCGAGCCAAAGagc > 1:8414669/1‑34 (MQ=255)
aaTGGCGAATCGCTGCTGACTCGCGAGCCAAAGATc > 1:7827904/1‑36 (MQ=255)
gCGAATCGCTGCTGACTCGCGAGCCAAAGATCAAct < 1:7402892/36‑1 (MQ=255)
cGAATCGCTGCTGACTCGCGAGCCAAAGATCAActt < 1:5789617/36‑1 (MQ=255)
gAATCGCTGCTGACTCGCGAGCCAAAGATCAActtc < 1:6837683/36‑1 (MQ=255)
gAATCGCTGCTGACTCGCGAGCCAAAGATCAActtc < 1:4171610/36‑1 (MQ=255)
gAATCGCTGCTGACTCGCGAGCCAAAGATCAActtc < 1:1682896/36‑1 (MQ=255)
aaTCGCTGCTGACTCGCGAGCCAAAGATCAActtct > 1:593872/1‑36 (MQ=255)
aaTCGCTGCTGACTCGCGAGCCAAAGATCAActtct < 1:3974825/36‑1 (MQ=255)
aaTCGCTGCTGACTCGCGAGCCAAAGATCAActtct > 1:833775/1‑36 (MQ=255)
aTCGCTGCTGACTCGCGAGCCAAAGATCAActtctt < 1:7993250/36‑1 (MQ=255)
cGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCa < 1:1713536/36‑1 (MQ=255)
gctgctGACTCGCGAGCCAAAGATCAACTTCTTCaa > 1:5345431/1‑36 (MQ=255)
ctgctgACTCGCGAGCCAAAGATCAACTTCTTCAAg > 1:8214121/1‑36 (MQ=255)
ctgctgACTCGCGAGCCAAAGATCAACTTCTTCAAg < 1:6716415/36‑1 (MQ=255)
ctgctgACTCGCGAGCCAAAGATCAACTTCTTCAAg > 1:6787966/1‑36 (MQ=255)
ctgctgACTCGCGAGCCAAAGATCAACTTCTTCAAg > 1:809999/1‑36 (MQ=255)
ctgctgACTCGCGAGCCAAAGATCAACTTCTTCAAg < 1:8059951/36‑1 (MQ=255)
ctgctgACTCGCGAGCCAAAGATCAACTTCTTCAAg < 1:3020884/36‑1 (MQ=255)
gctgACTCGCGAGCCAAAGATCAACTTCTTCAAGcc > 1:1277112/1‑36 (MQ=255)
ctgACTCGCGAGCAAAAGATCAACTTCTTCAAGCCg < 1:1304428/36‑1 (MQ=255)
tgACTCGCGAGCCAAAGATCAACTTCTTCAAGCCga > 1:5055852/1‑36 (MQ=255)
tgACTCGCGAGCCAAAGATCAACTTCTTCAAGCCga > 1:4779987/1‑36 (MQ=255)
tgAATCGCGAGCCAAAGATCAACTTCTTCAAGccta < 1:968614/36‑3 (MQ=38)
aCTCGCGAGCCAAAGATCAACTTCTTCAAGCAgatg > 1:3986893/1‑36 (MQ=255)
aCTCGCGAGAAAAAGATCAACTTCTTCAAGCCgatg < 1:3339472/36‑1 (MQ=37)
tCGCGAGCCAAAGATCAACTTCTTCAAGCCgatgat > 1:2321322/1‑36 (MQ=255)
cgcgAGCCAAACATAAACTTCTTCAAGCCGATGATc < 1:7269077/36‑1 (MQ=37)
cgcgAGCAAAAGATCAACTTCTTCAAGCCGATGATc < 1:1942250/36‑1 (MQ=255)
cgcgAGAAAAAGATCAACTTCTTCAAGCCGATGATc < 1:8264184/36‑1 (MQ=37)
cgcgAGAAAAAGATAAACTTCTTCAAGCCGATGATc < 1:3833658/36‑1 (MQ=25)
cgAGCCAAAGATCAACTTCTTCAAGCCGATGATCGa > 1:7221201/1‑36 (MQ=255)
cgAGCAAAAGATCAACTTCTTCAAGCCGATGATCGa < 1:721836/36‑1 (MQ=255)
gAGCCAAAGATCAACTTCTTCAATCCGATGATAGAc > 1:8183398/1‑36 (MQ=37)
gAGCCAAAGATCAACTTCTTCAAGCCGATGATCGAc > 1:3487193/1‑36 (MQ=255)
aGCCAAAGATCAACTTCTTCAAGCCGATGATCGACa > 1:740766/1‑36 (MQ=255)
aGCCAAAGATCAACTTCTTCAAGCCGATGATCGACa > 1:5958385/1‑36 (MQ=255)
aGCAAAAGATCAACTTCTACAAGCCGATGATCGACa < 1:3156827/36‑1 (MQ=37)
gacAAAGATCAACTTCTTAAAGCCGATGATCGCCaa < 1:6864121/34‑1 (MQ=37)
cAAAGATCAACTTCTTCAAGCCGATGATCGACAAcc > 1:1823399/1‑36 (MQ=255)
cAAAGATCAACTTCTTCAAGCCGATGATCGACAAcc > 1:2577235/1‑36 (MQ=255)
cAAAGATCAACTTCTTCAAGCCGATGATCGACAAcc > 1:1653345/1‑36 (MQ=255)
|
GGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACC > REL606/3157689‑3157760
Base quality scores: 
ATCG
< 3 ≤
ATCG
< 10 ≤
ATCG
< 21 ≤
ATCG
< 30 ≤
ATCG
< 37 ≤
ATCG
|