Read alignment evidence... | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|
seq id | position | change | freq | score | reads | annotation | genes | product | ||
* | REL606 | 3,947,797 | 0 | A→G | 5.4% | 0.9 0.0 ‑0.2 | 55 | 10 | aslA | acrylsulfatase‑like enzyme |
Rejected: E-value exceeds prediction threshold. | ||||||||||
Reads supporting (aligned to +/- strand): new base (2/1): ref base (35/17): total (37/18) | ||||||||||
Fisher's exact test strand distribution p-value = 1.00e+00 | ||||||||||
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.01e-01 |
GGCGACAGCTACAACAAGAAGTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATGTCGTC > REL606/3947762‑3947832
|
ggcgACAGCTACAACAAGAAGTTTGGGTGAAAACGa > 1:3415449/1‑36 (MQ=255)
ggcgACAGCTACAACAAGAAGTTTGGGTGAAAACGa > 1:2117166/1‑36 (MQ=255)
ggcgACAGCTACAACAAGAAGTTTGGGTGAAAACGa < 1:5362845/36‑1 (MQ=255)
gcgACAGCTACAACAAGAAGTTTGGGTGAAAACGaa > 1:4578155/1‑36 (MQ=255)
gcgACAGCTACAACAAGAAGTTTGGGTGAAAACGaa > 1:5164657/1‑36 (MQ=255)
gcgACAGCTACAACAAGAAGTTTGGGTGAAAACGaa < 1:1598173/36‑1 (MQ=255)
cgACAGCTACAACAAGAAGTTTGGGTGAAAACGaaa < 1:5069976/36‑1 (MQ=255)
gACAGCTACAACAAGAAGTTTGGGTGAAAACGaaaa > 1:384351/1‑36 (MQ=255)
gACAGCTACAACAAGAAGTTTGGGTGAAAACGaaaa < 1:2474653/36‑1 (MQ=255)
gACAGCTACAACAAGAAGTTTGGGTGAAAACGaaaa < 1:5725955/36‑1 (MQ=255)
aCAGCTACAACAAGAAGTTTGGGTGAAAACGAAAAt < 1:4351662/36‑1 (MQ=255)
cAGCTACAACAAGAAGTTTTGTTGAAAACGAAAAtt < 1:2693170/36‑1 (MQ=37)
cAGCTACAACAAGAAGTTTGGGTGAAAACGAAAAtt < 1:3124154/36‑1 (MQ=255)
aGCTACAACAAGAAGTTTGGGTGAAAACGAAAATTc > 1:3652649/1‑36 (MQ=255)
aGCTACAACAAGAAGTTTGGGTGAAAACGAAAATTc > 1:3750030/1‑36 (MQ=255)
aGCTACAACAAGAAGTTTGGGCGAAAACGAAAATTc > 1:3037292/1‑36 (MQ=255)
gCTACAACAAGAAGTTTGGGTGAAAACGGAAATTcc < 1:1860854/36‑1 (MQ=255)
gCTACAACAAGAAGTTTGGGTGAAAACGAAAATTcc > 1:5929357/1‑36 (MQ=255)
cTACAACAAGAAGTTTGGGTGAAAACGAAAATTCCa < 1:1657905/36‑1 (MQ=255)
tACAACAAGAAGTTTGGGTGAAAACGAAAATTCCAt > 1:4303667/1‑36 (MQ=255)
tACAACAAGAAGTTTGGGTGAAAACGAAAATTCCAt < 1:4469288/36‑1 (MQ=255)
caacaaGAAGTTTGGGTGAAAACGAAAATTCCATGc < 1:1355878/36‑1 (MQ=255)
caacaaGAAGTTTGGGTGAAAACGAAAATTCCATGc < 1:4222208/36‑1 (MQ=255)
aacaaGAAGTTTGGGTGAAAACGAAAATTCCATGCa > 1:2199152/1‑36 (MQ=255)
aacaaGAAGTTTGGGTGAAAACGAAAATTCCATGCa > 1:4655319/1‑36 (MQ=255)
acaaGAAGTTTGGGTGAAAACGAAAATTCCATGCaa > 1:2585242/1‑36 (MQ=255)
caaGAAGTTTGGGTGAAAACGAAAATTCCATGCaaa > 1:3007136/1‑36 (MQ=255)
aagaagTTTGGGTGAAAACGAAAATTCCATGCaaaa > 1:2373749/1‑36 (MQ=255)
aagaagTTTGGGTGAAAACGAAAATTCCATGCaaaa > 1:5906335/1‑36 (MQ=255)
aagaagTTTGGGTGAAAACGAAAATTCCATGCaaaa < 1:3590860/36‑1 (MQ=255)
agaagTTTGGGTGAAAACGAAAATTCCATGCAAAAt < 1:1827356/36‑1 (MQ=255)
agaagTTTGGGTGAAAACGAAAATTCCATGCAAAAt < 1:1089842/36‑1 (MQ=255)
gaagTTTGGGTGAAAACGGAAATTCCATGCAAAATg > 1:1314963/1‑36 (MQ=255)
gaagTTTGGGTGAAAACGAAAATTCCATGCAAAATg > 1:3030158/1‑36 (MQ=255)
gaagTTTGGGTGAAAACGAAAATTCCATGCAAAATg > 1:1013028/1‑36 (MQ=255)
agTTTGGGTGAAAACGAAAATTCCATGCAAAATGCt > 1:979710/1‑36 (MQ=255)
agTTTGGGTGAAAACGAAAATTCCATGCAAAATGCt < 1:2659770/36‑1 (MQ=255)
gTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTc > 1:951514/1‑36 (MQ=255)
gTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTc > 1:5479067/1‑36 (MQ=255)
gTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTc > 1:5587428/1‑36 (MQ=255)
gTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTc > 1:1680273/1‑36 (MQ=255)
tttGGGTGAAAACGAAAATTCCATGCAAAATGCTcc > 1:5505659/1‑36 (MQ=255)
tttGGGTGAAAACGAAAATTCCATGCAAAATGCTcc > 1:1476258/1‑36 (MQ=255)
tttGGGTGAAAACGAAAATTCCATGCAAAATGCTcc > 1:3811959/1‑36 (MQ=255)
tGGGTGAAAACGAAAATTCCATGCAAAATGCTCCgg > 1:2732589/1‑36 (MQ=255)
gggTGAAAACGAAAATTCCATGAAAAATGCTCAGGt > 1:3568690/1‑36 (MQ=37)
ggTGAAAACGAAAATTCCATGCAAAATGCTCCGGtt < 1:5292355/36‑1 (MQ=255)
ggTGAAAACGAAAATTCCATGCAAAATGCTCCGGtt > 1:665528/1‑36 (MQ=255)
gTGAAAACGAAAATTCCATGCAAAATGCTCCGGttt < 1:3494267/36‑1 (MQ=255)
gTGAAAACGAAAATTCCATGCAAAATGCTCCGGttt > 1:4709387/1‑36 (MQ=255)
tGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTc > 1:4625037/1‑36 (MQ=255)
gAAAACGGAAATTCCATGCAAAATGCTCCGGTTTCa > 1:3986076/1‑36 (MQ=255)
gAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCa > 1:833949/1‑36 (MQ=255)
gAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCa > 1:5261728/1‑36 (MQ=255)
gAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCa > 1:4307663/1‑36 (MQ=255)
aaaaCGAAAATTCCATGCAAAATGCTCCGGTTTCAt > 1:2579791/1‑36 (MQ=255)
aaaaCGAAAATTCCATGCAAAATGCTCCGGTTTCAt < 1:5108744/36‑1 (MQ=255)
aaaaCGAAAATTCCATGCAAAATGCTCCGGTTTCAt > 1:3153451/1‑36 (MQ=255)
aaaaCGAAAATTCCATGCAAAATGCTCCGGTTTCAt > 1:3253187/1‑36 (MQ=255)
aaaCGAAAATTCCATGCAAAATGCTCCGGTTTCATg < 1:4445335/36‑1 (MQ=255)
aaaCGAAAATTCCATGCAAAATGCTCCGGTTTCATg < 1:5454124/36‑1 (MQ=255)
aaCGAAAATTCCATGCAAAATGCTCCGGTTTCATgt > 1:2598278/1‑36 (MQ=255)
aaCGAAAATTCCATGCAAAATGCTCCGGTTTCATgt > 1:3807703/1‑36 (MQ=255)
aaCGAAAATTCCATGCAAAATGCTCCGGTTTCATgt > 1:3693603/1‑36 (MQ=255)
aCGAAAATTCCATGCAAAATGCTCCGGTTTCATgtc < 1:2968643/36‑1 (MQ=255)
aaaaTTCCATGCAAAATGCTCGGGTTTCATgtcgtc < 1:1811992/36‑1 (MQ=255)
aaaaTTCCATGCAAAATGCTCCGGTTTCATgtcgtc < 1:2049824/36‑1 (MQ=255)
aaaaTTCCATGCAAAATGCTCCGGTTTCATgtcgtc > 1:3867700/1‑36 (MQ=255)
aaaaTTCCATGCAAAATGCTCCGGTTTCATgtcgtc > 1:1268200/1‑36 (MQ=255)
aaaaTTCCACGCAAAATGCTCCGGTTTCATgtcgtc < 1:806183/36‑1 (MQ=255)
|
GGCGACAGCTACAACAAGAAGTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATGTCGTC > REL606/3947762‑3947832
Base quality scores: 
ATCG
< 3 ≤
ATCG
< 15 ≤
ATCG
< 24 ≤
ATCG
< 32 ≤
ATCG
< 37 ≤
ATCG
|