Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | REL606 | 1,333,604 | G→T | 33.6% | pseudogene (73/129 nt) | yciX → | predicted protein; b4522(b4574) |
Read alignment evidence... | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|
seq id | position | change | freq | score | reads | annotation | genes | product | ||
* | REL606 | 1,333,604 | 0 | G→T | 33.6% | 4.9 ‑4.0 0.0 | 22 | pseudogene (73/129 nt) | yciX | predicted protein; b4522(b4574) |
Reads supporting (aligned to +/- strand): new base (8/0): ref base (6/8): total (14/8) | ||||||||||
Fisher's exact test strand distribution p-value = 1.77e-02 | ||||||||||
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00 |
TTAAAACGTAAAGCCTGGCTGGCGGTTTTTCTTGGTTCTGCACTTTTCTGGGTGGTTGTCGCACTGCTG > REL606/1333571‑1333639
|
ttAAAACGTACAGCCTGGCCGGCGGTTTTTCTTGGt < 1:985513/36‑1 (MQ=37)
tAAAACGTAAAGCCTGGCTGGCGGTTTTTTTTTGtt > 1:590996/1‑36 (MQ=37)
tAAAACGTAAAGCCTGGCTGGCGGGTTTTTTTTGtt > 1:802377/1‑36 (MQ=21)
aaCGTAAAGCCTGGCTGGCGGTTTTTCTTGGTTCTg < 1:3535100/36‑1 (MQ=255)
aCGTAAAGCCTTGCTGGCGGTTTTTCTTGGTTCTGc < 1:651358/36‑1 (MQ=255)
aaaGCCTGGCTGGCGGTTTTTTTTTTTTCTGTACtt > 1:1762075/1‑36 (MQ=16)
aaGCCTGGCTGGCGGTTTTTTTTTGTTTTTCACttt > 1:4605790/1‑36 (MQ=16)
gCCTGGCTGGCGGTTTTTCTTGGTTCTGCACttttt > 1:1181244/1‑35 (MQ=255)
ccTGGCTGGCGGTTTTTTTTTGTTCTTCACtttttt > 1:4027652/1‑34 (MQ=21)
ctggctggCGGTTTTTTTTTTTTCTTCACTTTTCtt > 1:5364358/1‑35 (MQ=16)
ctggctggCGGTTTTTCTTGGTTCTGCACTTTTCTg < 1:4090500/36‑1 (MQ=255)
ctggctggCGGGTTTTTTTTGTTGTGCCCttttttt > 1:1612488/1‑33 (MQ=9)
tggctggcGGTTTTTCTTGGTTCTGCACTTTTCTgg < 1:889600/36‑1 (MQ=255)
ctggcGGTTTTTTTTTGTTTTTCACTTTTCTTggtt > 1:5845521/1‑35 (MQ=9)
gcGGTTTTTCTTGGTTCTGCACTTTTCTGGGTGGtt < 1:4535365/36‑1 (MQ=255)
ggTTTTTCTTGGTTCTGCACTTTTCTTGGTGGTTGt > 1:892875/1‑36 (MQ=255)
ggTTTTTCTTGGTTCTGCACTTTTCTGGGTGGTTGt > 1:3026124/1‑36 (MQ=255)
gTTTTTCTTGGTTCTGCACTTTTCTGGGTGGTTGTc > 1:5374164/1‑36 (MQ=255)
ttttCTTGGTTCTGCACTTTTCTGGGTGGTTGTCGc > 1:400940/1‑36 (MQ=255)
ttttCTTGGTTCTGCACTTTTCTGGGTGGTTGTCGc < 1:178216/36‑1 (MQ=255)
tCTTGGTTCTGCACTTTTCTGGGTGGTTGTCGCAct < 1:2503137/36‑1 (MQ=255)
ttGGTTCTGCACTTTTCTGGGTGTTTGTCGCActgc > 1:2083825/1‑36 (MQ=255)
ggTTCTGCACTTTTCTGGGTGGTTGTCGCActgctg > 1:601300/1‑36 (MQ=255)
|
TTAAAACGTAAAGCCTGGCTGGCGGTTTTTCTTGGTTCTGCACTTTTCTGGGTGGTTGTCGCACTGCTG > REL606/1333571‑1333639
Base quality scores: 
ATCG
< 3 ≤
ATCG
< 5 ≤
ATCG
< 10 ≤
ATCG
< 26 ≤
ATCG
< 37 ≤
ATCG
|