Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | REL606 | 4,204,731 | A→G | 11.4% | E658G (GAA→GGA) | metH → | B12‑dependent methionine synthase |
Read alignment evidence... | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|
seq id | position | change | freq | score | reads | annotation | genes | product | ||
* | REL606 | 4,204,731 | 0 | A→G | 11.4% | 2.6 ‑3.7 0.0 | 49 | 1973 | metH | B12‑dependent methionine synthase |
Reads supporting (aligned to +/- strand): new base (6/0): ref base (20/22): total (27/22) | ||||||||||
Fisher's exact test strand distribution p-value = 2.48e-02 | ||||||||||
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00 |
CGCTAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAG > REL606/4204696‑4204766
|
cGCTAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGa < 1:3686070/36‑1 (MQ=255)
cGCTAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGa > 1:5329785/1‑36 (MQ=255)
gCTAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGaa < 1:2664099/36‑1 (MQ=255)
cTAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGGAg > 1:5836209/1‑36 (MQ=255)
cTAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAg > 1:4111063/1‑36 (MQ=255)
cTAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAg < 1:6970442/36‑1 (MQ=255)
cTAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAg < 1:1820106/36‑1 (MQ=255)
tAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGt < 1:8418451/36‑1 (MQ=255)
aaCGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTg > 1:8161583/1‑36 (MQ=255)
aCGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGa > 1:2189794/1‑36 (MQ=255)
aCGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGa > 1:3977591/1‑36 (MQ=255)
aCGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGa < 1:5575736/36‑1 (MQ=255)
aCGCCCAGCAGGCGGAGTGGCGCTCGCGGGTAGTGa > 1:3896326/1‑36 (MQ=37)
cGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGaa > 1:2059399/1‑36 (MQ=255)
cccAGCAGGCGGAGTGGCGCTCGTGGGGAGGGGATa > 1:3613135/1‑36 (MQ=25)
ccAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATaa > 1:7725842/1‑36 (MQ=255)
ccAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATaa < 1:7953578/36‑1 (MQ=255)
cagcagGCGGAGTGGCGCTCGTGGGAAGTGAATaaa < 1:1507739/36‑1 (MQ=255)
cagcagGCGGAGTGGCGCTCGTGGGAAGTGAATaaa < 1:8299003/36‑1 (MQ=255)
cagcagGCGGAGTGGCGCTCGGGGGAAGGGAATaaa > 1:2849332/1‑36 (MQ=38)
agcagGCGGAGTGGCGCTCGTGGGGAGTGAATAAAc > 1:5605443/1‑36 (MQ=255)
agcagGCGGAGTGGCGCTCGTGGGAAGTGAATAAAc > 1:336569/1‑36 (MQ=255)
agcagGCGGAGTGGCGCTCGTGGGAAGTGAATAAAc < 1:2109251/36‑1 (MQ=255)
agcagGCGGAGTGGCGCTCGTGGGAAGTGAATAAAc > 1:2551918/1‑36 (MQ=255)
gcagGCGGAGTGGCGCTCGTGGGAAGTGAATAAACg > 1:4227274/1‑36 (MQ=255)
gcagGCGGAGTGGCGCTCGTGGGAAGTGAATAAACg > 1:875324/1‑36 (MQ=255)
gcagGCGGAGTGGCGCTCGGGGGAAGTGAATAAACg > 1:2726109/1‑36 (MQ=255)
cagGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGt < 1:7646436/36‑1 (MQ=255)
agGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTc > 1:2378656/1‑36 (MQ=255)
gCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTg > 1:4960757/1‑36 (MQ=255)
cGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTgg > 1:7634123/1‑36 (MQ=255)
cGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTgg < 1:3016492/36‑1 (MQ=255)
ggAGTGGCGCTCGTGGGGAGTGAATAAACGTCTGGa > 1:434354/1‑36 (MQ=255)
ggAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGa < 1:139840/36‑1 (MQ=255)
ggAGTGGCGCTCGGGGGGAGTGGATAAACGGCTggg > 1:6936008/1‑35 (MQ=16)
aGTGGCGCTCGTGGGGAGTGAATAAAAGTCTGGGAt > 1:4665053/1‑36 (MQ=37)
gTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATa > 1:2967036/1‑36 (MQ=255)
gTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATa < 1:6690199/36‑1 (MQ=255)
tGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATAc < 1:5990015/36‑1 (MQ=255)
gcgcTCGTGGGAAGTGAATAAACGTCTGGAATACTc < 1:5998488/36‑1 (MQ=255)
gcgcTCGTGGGAAGTGAATAAACGTCTGGAATACTc < 1:3652081/36‑1 (MQ=255)
gcTCGTGGGAAGTGAATAAACGTCTGGAATACTCGc > 1:738577/1‑36 (MQ=255)
cTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCt < 1:4092165/36‑1 (MQ=255)
cGTGGGAAGTGAATAAGCGTCTGGAATACTCGCTgg < 1:6108231/36‑1 (MQ=255)
gTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGt > 1:6078887/1‑36 (MQ=255)
gTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGt > 1:2972960/1‑36 (MQ=255)
tGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTc < 1:5632577/36‑1 (MQ=255)
tggAAGTGAATAAACGTCTGGAATACTCGCTGGTCa < 1:6542361/35‑1 (MQ=255)
gggAAGTGAATAAACGTCTGGAATACTCGCTGGTCa > 1:2625739/1‑36 (MQ=255)
ggAAGTGAATAAACGTCTGGAATACTCGCTGGTCaa < 1:643397/36‑1 (MQ=255)
gAAGTGAATAAACGTCTGGAATACTCGCTGGTCaaa < 1:796757/36‑1 (MQ=255)
gAAGTGAATAAACGTCTGGAATACTCGCTGGTCaaa < 1:2339996/36‑1 (MQ=255)
aaGTGAATAAACGTCTGGAATACTCGCTGGTCAAAg > 1:7251172/1‑36 (MQ=255)
aaGTGAATAAACGTCTGGAATACTCGCTGGTCAAAg < 1:5194051/36‑1 (MQ=255)
aaGTGAATAAACGTCTGGAATACTCGCTGGTCAAAg < 1:1401665/36‑1 (MQ=255)
|
CGCTAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAG > REL606/4204696‑4204766
Base quality scores: 
ATCG
< 3 ≤
ATCG
< 14 ≤
ATCG
< 24 ≤
ATCG
< 30 ≤
ATCG
< 37 ≤
ATCG
|