Read alignment evidence... | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|
seq id | position | change | freq | score | reads | annotation | genes | product | ||
* | REL606 | 2,864,066 | 0 | C→T | 12.8% | 4.3 ‑7.1 0.0 | 61 | intergenic (‑469/‑216) | ygdP/mutH | dinucleoside polyphosphate hydrolase/DNA mismatch repair protein |
Rejected: Prediction has biased read strand distribution supporting polymorphism. | ||||||||||
Reads supporting (aligned to +/- strand): new base (9/0): ref base (21/31): total (30/31) | ||||||||||
Fisher's exact test strand distribution p-value = 8.25e-04 | ||||||||||
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00 |
CAAGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTCTGTGGATAACATGGTGT > REL606/2864031‑2864101
|
cAAGGCTGTAAATGGAAACAGTAGTGGCGGTTTTTc < 1:605187/36‑1 (MQ=255)
cAAGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTc < 1:5796273/36‑1 (MQ=255)
cAAGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTc < 1:3500671/36‑1 (MQ=255)
cAAGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTc < 1:4219496/36‑1 (MQ=255)
aGGTTGTAAATGGAAACAGTAGTGTAGGTTTTTcac < 1:4353628/36‑1 (MQ=38)
aGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTtac > 1:5847607/1‑36 (MQ=255)
aGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTcac < 1:1490564/36‑1 (MQ=255)
gCTGTAAATGGAAACAGTAGTGGAGGTTTTTCACAg < 1:3398548/36‑1 (MQ=255)
gCTGTAAATGGAAACAGTAGTGGAGGTTTTTCACAg < 1:3845054/36‑1 (MQ=255)
tGTAAATGGAAACAGTAGTGGAGGTTTTTCACAGtt < 1:549586/36‑1 (MQ=255)
gTAAATGGAAACAGTAGTGGAGGTTTTTTACCGttt > 1:4758980/1‑35 (MQ=37)
gTAAATGGAAACAGTAGTGGAGGTTTTTTACAGTTa > 1:3811422/1‑36 (MQ=255)
gTAAATGGAAACAGTAGTGGAGGTTTTTTACAAttt > 1:3911040/1‑35 (MQ=37)
gTAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTa < 1:4194575/36‑1 (MQ=255)
tAAATGGAAACCGTAGTGGAGGTTTTTCACAGTTAt < 1:2518081/36‑1 (MQ=255)
tAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTAt < 1:2554100/36‑1 (MQ=255)
aaaTGGAAACAGTAGTGGAGGTTTTTCACAGTTATc > 1:341081/1‑36 (MQ=255)
aaaTGGAAACAGTAGTGGAGGTTTTTCACAGTTATc < 1:3546576/36‑1 (MQ=255)
aaaTGGAAACAGTAGTGGAGGTTTTTCACAGTTATc < 1:605545/36‑1 (MQ=255)
aaaTGGAAACAGTAGTGGAGGGTTTTCACAGTTATc > 1:274957/1‑36 (MQ=255)
aaaTGGAAACAGGAGTGGAGGGTGTTTACAGTTATc > 1:75140/1‑36 (MQ=19)
aaTGGAAACAGTAGTGGAGGTTTTTCACAGTTATcc < 1:5769673/36‑1 (MQ=255)
aTGGAAACAGTAGTGGAGGTTTTTTACAGTTTTccc > 1:2790795/1‑36 (MQ=38)
aTGGAAACAGTAGTGGAGGTTTTTCACAGTTATccc < 1:3488309/36‑1 (MQ=255)
aTGGAAACAGTAGTGCAGGTTTTTCACATTTATccc < 1:2309977/36‑1 (MQ=37)
tGGAAAAAGTAGTGGAGGATTTTTGCCGTTATGCCa > 1:1400957/1‑36 (MQ=9)
aaCAGTAGTGGAGGTTTTTCACAGTTATCCCAGCtt < 1:1224433/36‑1 (MQ=255)
aCAGTCGTGGAGGTTTTTCACAGTTATCCCAGCttt < 1:2885867/36‑1 (MQ=255)
aCAGTAGTGGAGGTTTTTTACAGTTATCCCAGCttt > 1:458513/1‑36 (MQ=255)
aCAGTAGTGGAGGTTTTTCACAGTTATNNCAGCttt > 1:54296/1‑36 (MQ=255)
aCAGTAGTGGAGGGTTTTCACAGTTATTCCAACttt > 1:4408555/1‑36 (MQ=37)
cAGTAGTGGATGTTTTTCACAGTTATCCCAGCTTTc < 1:5770711/36‑1 (MQ=255)
cAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTc < 1:5861069/36‑1 (MQ=255)
agtagtGGAGGTTTTTCGCAGTTATCCCAGCTTTCt < 1:4875608/36‑1 (MQ=255)
tagtGGAGGTTTTTTAAAGTTATTCCAGGTTTCttt > 1:2067689/1‑34 (MQ=16)
tagtGGAGGTTTTTCACAGTTATCCCAGCTTTCtgt < 1:5708893/36‑1 (MQ=255)
agtGGAGGTTTTTCACAGTTATCCCAGCTTTCtgtg < 1:5189661/36‑1 (MQ=255)
agtGGAGGTTTTTCACAGTTATCCCAGCTTTCtgtg < 1:4274463/36‑1 (MQ=255)
gtGGAGGTTTTTTACAGTTATCCCAGCTTTTTGTgg > 1:2750582/1‑36 (MQ=38)
gtGGAGGTTTTTCACAGTTATCCCAGCTTTCTGTgg > 1:2058973/1‑36 (MQ=255)
gtGGAGGTTTTTCACAGTTATCCCAGCTTTCTGTgg < 1:4739219/36‑1 (MQ=255)
gtGGAGGTTTTTCACAGTTATCCCAGCTTTCTGTgg > 1:890885/1‑36 (MQ=255)
tGGAGGTTTTTCACAGTTATCCCAGCTTTCTGTGGa > 1:1831605/1‑36 (MQ=255)
tGGAGGTTTTTCACAGTTATCCCAGCTTTCTGTGGa > 1:2205270/1‑36 (MQ=255)
tGGAGGTTTTTCACAGTTATCCCAGCTTTCTGTGGa < 1:754704/36‑1 (MQ=255)
tGGAGGTTTTTCACAGTTATCCCAGCCTTCTGTGGa < 1:947404/36‑1 (MQ=255)
ggAGGTTTTTCACAGTTATCCCAGCTTTCTGTGGAt > 1:1872060/1‑36 (MQ=255)
ggAGGTTTTTCACAGTTATCCCAGCTTTCTGTGGAt > 1:4048851/1‑36 (MQ=255)
ggAGGTTTTTCACAGTTATCCCAGCTTTCTGTGGAt < 1:4537687/36‑1 (MQ=255)
ggAGGTTTTTCACAGTTATCCCAGCTTTCTGTGGAt < 1:2699644/36‑1 (MQ=255)
gAGGTTTTTCACAGTTATCCCAGCTTTCTGTGGATa < 1:1590889/36‑1 (MQ=255)
gAGGTTTTTCACAGTTATCCCAGCTTTCTGTGGATa > 1:3981989/1‑36 (MQ=255)
aGGTTTTTCACAGTTATCCCAGCTTTCTGTGGATaa < 1:1109124/36‑1 (MQ=255)
aGGTTTTTCACAGTTATCCCAGCTTTCTGTGGATaa < 1:740937/36‑1 (MQ=255)
aGGTTTTTCACAGTTATCCCAGCTTTCTGTGGATaa < 1:3512216/36‑1 (MQ=255)
aGGTTTTTCACAGTTATCCCAGCTTTCTGTGGATaa > 1:3639489/1‑36 (MQ=255)
ggTTTTTCACAGTTATCCCAGCTTTCTGTGGATAAc > 1:1551128/1‑36 (MQ=255)
ggTTTTTCACAGTTATCCCAGCTTTCTGTGGATAAc > 1:3541176/1‑36 (MQ=255)
ggTTTTTCACAGTTATCCCAGCTTTCTGTGGATAAc > 1:5265422/1‑36 (MQ=255)
gTTTTTCACAGTTATCCCAGCTTTCTGTGGATAACa > 1:5545923/1‑36 (MQ=255)
gTTTTTCACAGTTATCCCAGCTTTCTGTGGATAACa > 1:5544078/1‑36 (MQ=255)
gTTTTTCACAGTTATCCCAGCTTTCTGTGGATAACa > 1:840865/1‑36 (MQ=255)
tttttCACAGTTATCCCTGCTTTCTGTGGATAACAt < 1:5944689/36‑1 (MQ=255)
ttttCACAGTTATCCCAGCTTTCTGTGGATAACATg < 1:2221794/36‑1 (MQ=255)
tttCACAGTTATCCCAGCTTTCTGTGGATAACATgg > 1:3341548/1‑36 (MQ=255)
tttCACAGTTATCCCAGCTTTCTGTGGATAACATgg < 1:4261158/36‑1 (MQ=255)
tCACAGTTATCCCAGCTTTCTGTGGATAACATGgtg > 1:3649961/1‑36 (MQ=255)
tCACAGTTATCCCAGCTTTCTGTGGATAACATGgtg > 1:3387614/1‑36 (MQ=255)
cacaGTTATCCCAGCTTTCTGTGGATAACATGgtgt < 1:3749862/36‑1 (MQ=255)
cacaGTTATCCCAGCTTTCTGTGGATAACATGgtgt > 1:657607/1‑36 (MQ=255)
|
CAAGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTCTGTGGATAACATGGTGT > REL606/2864031‑2864101
Base quality scores: 
ATCG
< 3 ≤
ATCG
< 10 ≤
ATCG
< 20 ≤
ATCG
< 30 ≤
ATCG
< 37 ≤
ATCG
|