Read alignment evidence... | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|
seq id | position | change | freq | score | reads | annotation | genes | product | ||
* | REL606 | 4,231,619 | 0 | T→C | 8.4% | 0.7 ‑5.3 0.0 | 76 | 182 | ubiC | chorismate pyruvate lyase |
Rejected: E-value exceeds prediction threshold. | ||||||||||
Reads supporting (aligned to +/- strand): new base (0/8): ref base (37/30): total (37/39) | ||||||||||
Fisher's exact test strand distribution p-value = 5.19e-03 | ||||||||||
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00 |
GATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCC.GAAGAACTGCCGCTGCTGCCGAAAGAGTCTC > REL606/4231584‑4231654
|
gatgatCCGCGAAGGGTTTGTCGAGCAGAATGAAAt < 1:1296134/36‑1 (MQ=255)
atgatCCGCGAAGGGTTTGTCGAGCAGAATGAAATc > 1:1820223/1‑36 (MQ=255)
atgatCCGCGAAGGGTTTATCGAGCAGAATGAAATc < 1:1923218/36‑1 (MQ=255)
tgatCCGCGAAGGGTTTGTCGAGCAGAATGAAATcc < 1:5184081/36‑1 (MQ=255)
tgatCCGCGAAGGGTTTGTCGAGCAGAATGAAATcc < 1:1617024/36‑1 (MQ=255)
tgatCCGCGAAGGGTTTGTCGAGCAGAATGAAATcc > 1:7075270/1‑36 (MQ=255)
tgatCCGCGAAGGGTTTGTCGAGCAGAATGAAATcc > 1:6011475/1‑36 (MQ=255)
tgatCCGCGAAGGGCTTGTCGAGCAGAATGAAATcc > 1:3070492/1‑36 (MQ=255)
atCCGCGAAGGGTTTGTCGAGCAGAATGAAATcccc > 1:438127/1‑36 (MQ=255)
tCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCC.g < 1:3146230/36‑1 (MQ=255)
tCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCC.g < 1:62223/36‑1 (MQ=255)
ccGCGAAGGGTTTGTCGAGCAGAATGAAATCCCC.ga < 1:3300668/36‑1 (MQ=255)
cgcgAAGGGTTTGTCGAGCAGAATGAAATCCCC.gaa > 1:2207599/1‑36 (MQ=255)
gcgAAGGGTTTGTCGAGCAGAATGAAATCCCC.gaag < 1:5701730/36‑1 (MQ=255)
gcgAAGGGTTTGTCGAGCAGAATGAAATCCCC.gaag < 1:3003671/36‑1 (MQ=255)
gcgAAGGGTTTGTCGAGCAGAATGAAATCCCC.gaag > 1:7588897/1‑36 (MQ=255)
cgAAGGGTTTGTCGAGCAGAATGAAATCCCC.gaaga > 1:4050961/1‑36 (MQ=255)
cgAAGGGTTTGTCGAGCAGAATGAAATCCCC.gaaga < 1:5372313/36‑1 (MQ=255)
cgAAGGGTTTGTCGAGCAGAATGAAATCCCC.gaaga < 1:3791239/36‑1 (MQ=255)
cgAAGGGTTTGTCGAGCAGAATGAAATCCCC.gaaga > 1:7086415/1‑36 (MQ=255)
cgAAGGGTTTGTCGAGCAGAATGAAATCCCC.gaaga < 1:6614411/36‑1 (MQ=255)
cgAAGGGTTTGTAGAGCAGAATGAAATCCCC.gaaga < 1:1729607/36‑1 (MQ=255)
gAAGGGTTTGTCGAGCAGAATGAAATCCCC.gaagaa < 1:4448258/36‑1 (MQ=255)
gAAGGGTTTGTCGAGCAGAATGAAATCCCC.gaagaa > 1:5500895/1‑36 (MQ=255)
gAAGGGTTTGTCGAGCAGAATGAAATCCCC.gaagaa < 1:3687686/36‑1 (MQ=255)
gAAGGGTATGTCGAGCAGAATGAAATCCCC.gaagaa < 1:5766309/36‑1 (MQ=255)
aaGGGTTTGTCGAGCAGAATGAAATCCCC.GAAGAAc < 1:5821870/36‑1 (MQ=255)
aGGGTTTGTCGAGCAGAATGAAATCCCC.GAAGAACt > 1:2235948/1‑36 (MQ=255)
ggTTTGTCGAGCAGAATGAAATCCCC.GAAGAACTGc > 1:2553364/1‑36 (MQ=255)
ggTTTGTCGAGCAGAATGAAATCCCC.GAAGAACTGc < 1:2958566/36‑1 (MQ=255)
gTTTGTCGAGCAGAATGAAATCCCC.GAAGAACTGcc < 1:387084/36‑1 (MQ=255)
gTTTGTCGAGCAGAATGAAATCCCC.GAAGAACTGcc > 1:1756389/1‑36 (MQ=255)
gTTTGTCGAGCAGAATGAAATCCCC.GAAGAACTGcc < 1:4765371/36‑1 (MQ=255)
tttGTCGAGCAGAATGAAATCCCC.GAAGAACTGCCg > 1:487098/1‑36 (MQ=255)
tttGTCGAGCAGAATGAAATCCCC.GAAGAACTGCCg > 1:1241407/1‑36 (MQ=255)
tttGTCGAGAAGAATGAAATCCCC.GAAGAACTGCCg < 1:1956424/36‑1 (MQ=255)
ttGTCGAGCAGAATGAAATCCCC.GAAGAACTGCCgc > 1:1627714/1‑36 (MQ=255)
ttGTCGAGCAGAATGAAATCCCC.GAAGAACTGCCgc > 1:6345918/1‑36 (MQ=255)
ttGTCGAGCAGAATGAAATCCCC.GAAGAACTGCCgc < 1:3042340/36‑1 (MQ=255)
tGTCGAGCAGAATGAAATCCCC.GAAGAACTGCCgct > 1:3308157/1‑36 (MQ=255)
gTCGAGCAGAATGAAATCCCC.GAAGAACTGCCgctg > 1:4774873/1‑36 (MQ=255)
gTCGAGCAGAATGAAACCCCC.GAAGAACTGCCgctg < 1:7000058/36‑1 (MQ=255)
tcaaGCACAAAGAAACCCCC.GAAGATCTGCCgctgc < 1:5432075/33‑1 (MQ=12)
cGAGCAGAATGAAATCCCC.GAAGAACTGCCgctgct > 1:3643663/1‑36 (MQ=255)
cGAGAAGAATGAAATCCCC.GAAGAACTGCCgctgct < 1:1354217/36‑1 (MQ=255)
taGCAGAATGAAATCCCC.TAAGAACTGCCGctgctg < 1:1028800/35‑1 (MQ=38)
gAGCAGAATGAAATCCCC.GAAGAACTGCCGctgctg > 1:928545/1‑36 (MQ=255)
gAGAAGAATGAAATCCCC.GAAGAACTGCCGctgctg < 1:2289630/36‑1 (MQ=255)
aGCAGAATGAAATCCCC.GAAGAACTGCCGCtgctgc > 1:6947472/1‑36 (MQ=255)
aGCAGAATGAAATCCCC.GAAGAACTGCCGCtgctgc < 1:39330/36‑1 (MQ=255)
aGCAGAATGAAATCCCC.GAAGAACTGCCGCtgctgc > 1:7107741/1‑36 (MQ=255)
aGCAGAAGGAAACCCCC.GAGGATCTGCCGCtgctgc < 1:4732461/36‑1 (MQ=16)
gCAGAATGAAATCCCC.GAAGAACTGCCGCTGCTGcc < 1:3850201/36‑1 (MQ=255)
gCAGAATGAAATCCCC.GAAGAACTGCCGCTGCTGcc > 1:5517815/1‑36 (MQ=255)
gCAGAATGAAATCCCC.GAAGAACTGCCGCTGCTGcc > 1:7817380/1‑36 (MQ=255)
gCAGAATGAAATCCCC.GAAGAACTGCCGCTGCTGcc > 1:5175558/1‑36 (MQ=255)
gCAAAATAAAACCCCC.GAAGAACTGCCGCTGCTGcc < 1:928737/36‑1 (MQ=25)
aaaaTGAAACCCCC.GAAAAGCTGCCGCTGCTGCCGa < 1:4969901/34‑1 (MQ=21)
aGAATGAAATCCCC.GAAGAACTGCCGCTGCTGCCGa > 1:3490657/1‑36 (MQ=255)
aGAATGAAATCCCC.GAAGAACTGCCGCTGCTGCCGa > 1:2057780/1‑36 (MQ=255)
aGAATGAAATCCCC.GAAGAACTGCCGCTGCTGCCGa > 1:1826831/1‑36 (MQ=255)
gAATGAAATCCCC.GAAGAACTGCCGCTGCTGCCGaa > 1:6252613/1‑36 (MQ=255)
aaaTGAAACCCCC.GAAGAACTGCCGCTGCTGCCGaa < 1:7624821/35‑1 (MQ=255)
aaTGAAATCCCC.GAAGAACTGCCGCTGCTGCCGaaa > 1:3660818/1‑36 (MQ=255)
aaTAAAATCCCC.AAAGAACTGCCGCTGCTGCCGaaa < 1:7293996/36‑1 (MQ=37)
cataaaTCCCCAGAATAACTGCCGCTGCTGCCGaaa < 1:426905/33‑1 (MQ=19)
aTGAAATCCCC.GAAGAACTGCCGCTGCTGCCGAAag > 1:148389/1‑36 (MQ=255)
taaaaCCCCC.AAAGAACTGCCGCTGCTGCCGAAaga < 1:6289103/34‑1 (MQ=25)
tGAAATCCCC.GAAGAACTGCCGCTGCTGCTGAAaga > 1:7914348/1‑36 (MQ=255)
cgAAAACCCC.GAAAAACTGCCGCTGCTGCCGAAaga < 1:1651026/35‑1 (MQ=37)
aaaaTCCCC.GAAGAACTGCCGCTGCTGCCGAAagag < 1:714315/35‑1 (MQ=255)
aaaCCCCC.AAAGAACTGCCGCTGCTGCCGAAAGAGt < 1:6624950/36‑1 (MQ=37)
aaTCCCC.GAAGAACTGCCGCTGCTGCCGAAAGAGtc > 1:2215552/1‑36 (MQ=255)
aaTCCCC.AAAAAACTGCCGCTGCTGCCGAAAGAGtc < 1:992855/36‑1 (MQ=37)
aTCCCC.GAAGAACTGCCGCTGCTGCCGAAAGAGtct > 1:5729973/1‑36 (MQ=255)
aTCCCC.GAAGAACTGCCGCTGCTGCCGAAAGAGtct > 1:3992030/1‑36 (MQ=255)
aTCCCC.GAAGAACTGCCGCTGCTGCCGAAAGAGtct > 1:2535551/1‑36 (MQ=255)
tCCCC.GAAGAACTGCCGCTGCTGCCGAAAGAGtctc > 1:1871226/1‑36 (MQ=255)
|
GATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCC.GAAGAACTGCCGCTGCTGCCGAAAGAGTCTC > REL606/4231584‑4231654
Base quality scores: 
ATCG
< 3 ≤
ATCG
< 10 ≤
ATCG
< 22 ≤
ATCG
< 32 ≤
ATCG
< 37 ≤
ATCG
|