Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | REL606 | 2,574,203 | A→T | 32.2% | intergenic (+45/+773) | sseA → / ← sseB | 3‑mercaptopyruvate sulfurtransferase/rhodanase‑like enzyme, sulfur transfer from thiosulfate |
Read alignment evidence... | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|
seq id | position | change | freq | score | reads | annotation | genes | product | ||
* | REL606 | 2,574,203 | 0 | A→T | 32.2% | 6.6 ‑2.1 ‑0.1 | 17 | intergenic (+45/+773) | sseA/sseB | 3‑mercaptopyruvate sulfurtransferase/rhodanase‑like enzyme, sulfur transfer from thiosulfate |
Reads supporting (aligned to +/- strand): new base (5/0): ref base (6/5): total (12/5) | ||||||||||
Fisher's exact test strand distribution p-value = 1.19e-01 | ||||||||||
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.97e-01 |
CAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTG > REL606/2574171‑2574238
|
cAATAAAAAACCGCCGAATTTGGCGGTTTTTTTttt > 1:3773867/1‑35 (MQ=255)
aaTAAAAAACCGCCAAATTCGGCGGCTTTTTTtttt > 1:1754050/1‑34 (MQ=17)
aTAAAAAACCGCCGAATTTGGCGGGTTTTTGttttt > 1:4936253/1‑33 (MQ=14)
aaaaaCCGCCGAATTTGGCGGCTTTTTTTTTTTAGt > 1:5105746/1‑36 (MQ=12)
gccgCATTTGGCGGTTTTTTATTGCTAGTCTGGTTc < 1:2703443/36‑1 (MQ=255)
gccgAATTTGGCGGTTTTTTTTTTTTTGTTTGGTTc > 1:1358190/1‑36 (MQ=9)
gccgAATTTGGCGGTTTTTTTTTTCTATTTTTGTTc > 1:972590/1‑36 (MQ=12)
gAATTTGGCGGTTTTTTATTTCTAGTCTGGTTcgcg > 1:216633/1‑36 (MQ=255)
ggTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCca < 1:1483135/36‑1 (MQ=255)
ttttttATTGCTAGTCTGGTTCGCGGCCTTTCcagc > 1:529511/1‑36 (MQ=255)
ttttATTGCTAGTCTGGTTCGCGGCCTTTCcagcag < 1:3466720/36‑1 (MQ=255)
ttttATTGCTAGTCTGGTTCGCGGCCTTTCcagcag < 1:1421073/36‑1 (MQ=255)
ttttATTGCTAGTCTGGTTCGCGACCTTTCcagcag > 1:2342438/1‑36 (MQ=255)
tttATTGCTAGTCTGGTTCGCGGCATTTCCAGCAgg > 1:1969427/1‑36 (MQ=255)
ttATTGCTAGTCTGTTTCTCGGCCTTTCCAGCAGGt < 1:3124950/36‑1 (MQ=18)
ttATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGt > 1:1472300/1‑36 (MQ=255)
tATTGCTAGTCTGGTTCGCGGCCTTTCGAGCAGGtt > 1:5892136/1‑36 (MQ=255)
aTTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTg > 1:3748747/1‑36 (MQ=255)
aTTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTg > 1:5010379/1‑36 (MQ=255)
|
CAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTG > REL606/2574171‑2574238
Base quality scores: 
ATCG
< 3 ≤
ATCG
< 7 ≤
ATCG
< 15 ≤
ATCG
< 27 ≤
ATCG
< 37 ≤
ATCG
|