Read alignment evidence... | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|
seq id | position | change | freq | score | reads | annotation | genes | product | ||
* | REL606 | 2,735,994 | 0 | C→G | 13.2% | 5.8 ‑13.4 0.0 | 73 | 145 | ascB | cryptic 6‑phospho‑beta‑glucosidase |
Rejected: Polymorphism prediction not supported by minimum number of reads on both strands. | ||||||||||
Reads supporting (aligned to +/- strand): new base (11/0): ref base (14/47): total (26/47) | ||||||||||
Fisher's exact test strand distribution p-value = 1.47e-06 | ||||||||||
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00 |
CACATGGCGAGCATCGAATGGC.GGTGA.AACTGGGGCTGGAAAAACGTTTTCAGTTGCGCGATGACGAGTTT > REL606/2735959‑2736029
|
cacaTGGCGAGCATCGAATGGC.GGGGA.AACCAGGGc > 1:1045539/1‑36 (MQ=25)
acaTGGCGAGCATCGAATGGC.GGTGA.AACTGGGGCt < 1:3809306/36‑1 (MQ=255)
aTGGCGAGCATCGAATGGC.GGTGA.AACTGGGGCTgg < 1:5190139/36‑1 (MQ=255)
aTGGCGAGCATCGAATGGC.GGGGG.AAATGGGGGTgg > 1:2663627/1‑36 (MQ=16)
aTGGCGAGCATCGAATGGC.GGGGA.AAATGGGGGTgg > 1:3153511/1‑36 (MQ=25)
tGGCGAGCATCGAATGGC.GGTGA.AACTGGGGCTGGa < 1:4615166/36‑1 (MQ=255)
tGGCGAGCATCGAATGGC.GGTGA.AACTGGGGCTGGa < 1:309059/36‑1 (MQ=255)
tGGCGAGCATCGAATGGC.GGGGA.AAAAGGGGGTggg > 1:2156053/1‑35 (MQ=16)
ggggAGCATCGAATGAC.GGTGA.AACTGGGGCTGGaa < 1:96130/33‑1 (MQ=37)
ggCGAGCATCGAATGGC.GGTGA.AACTGGGGCTGGaa < 1:2218393/36‑1 (MQ=255)
ggCGAGCATCGAATGGC.GGTGA.AACTGGGGCTGGaa < 1:2414174/36‑1 (MQ=255)
gCGAGCATCGAATGGG.GGTGA.AACTGGGGCTGGaaa < 1:5840954/36‑1 (MQ=255)
gCGAGCATCGAATGGC.GGTGG.AAACGGGGGTgggga > 1:2574793/1‑33 (MQ=12)
cGAGCATCGAATGGG.GGTGA.AACTGGGGCTGGaaaa < 1:5186746/36‑1 (MQ=255)
cGAGCATCGAATGGC.GGTGA.AACTGGGGCTGGaaaa < 1:1186745/36‑1 (MQ=255)
gAGCATCGAATGGC.GGTGG.AACTGGGGGTGGgaaaa > 1:1869216/1‑36 (MQ=25)
gAGCATCGAATGGC.GGTGG.AAATGGGGGTGGgaaaa > 1:2974262/1‑36 (MQ=19)
gAGCATCGAATGGC.GGTGG.AAAAGGGGGTGGgaaaa > 1:894135/1‑36 (MQ=16)
gAGCATCGAATGGC.GGTAA.AACTGGGGGGGGgaaaa > 1:626920/1‑36 (MQ=19)
aGCATCGAATGGC.GGTGAGAATGGGGGGTGGaaaaa > 1:5269231/1‑36 (MQ=16)
aGCATCGAATGGC.GGTGA.AACTGGGGCTGGAAAAAc < 1:5097304/36‑1 (MQ=255)
gCATCGAATGGC.GGGGC.AAAAGGGGCTGGaaaaaaa > 1:1503270/1‑34 (MQ=16)
cATCGAATGGC.GGCGA.AACTGGGGGGGGGAAAACGt > 1:4858330/1‑36 (MQ=16)
aTCGAATGGCGGGGGA.AACCGGGGTGGGAAAAAGGt > 1:3321385/1‑36 (MQ=2)
aTCGAATGGC.GGTGA.AACTGGGGCTGGAAAAACGtt < 1:1583146/36‑1 (MQ=255)
aTCGAATGGC.GGTGA.AACTGGGGCTGGAAAAACGtt < 1:4426043/36‑1 (MQ=255)
tCGAATGGC.GGTGG.AACTGGGGGTGGGAAAACGgtt > 1:1841974/1‑36 (MQ=16)
cGAATGGC.GGTGA.AACTGGGGCTGGAAAAACGtttt < 1:5331580/36‑1 (MQ=255)
cGAATGGC.GGTGA.AACTGGGGCTGGAAAAACGtttt < 1:3743109/36‑1 (MQ=255)
aaTGGC.GGTGA.AACTGGGGCTGGAAAAACGTTTTCa < 1:1477905/36‑1 (MQ=255)
aaTGGC.GGTGA.AACTGGGGCTGGAAAAACGTTTTCa < 1:3897914/36‑1 (MQ=255)
tGGC.GGTGA.AACTGGGGCTGGAAAAACGTTTTCAGt < 1:5598876/36‑1 (MQ=255)
tGGC.GGTGA.AACTGGGGCTGGAAAAACGTTTTCAGt < 1:1769221/36‑1 (MQ=255)
tGGC.GGTGA.AACTGGGGCTGGAAAAACGTTTTCAGt < 1:501513/36‑1 (MQ=255)
ggC.GGTGA.AACTGGGGCTGGAAAAACGTTTTCAGtt < 1:281118/36‑1 (MQ=255)
ggC.GGTGA.AACTGGGGCTGGAAAAACGTTTTCAGtt < 1:5273897/36‑1 (MQ=255)
gC.GGTGA.AACTGGGGCTGGAAAAACGTTTTCAGTTg < 1:4441975/36‑1 (MQ=255)
gC.GGTGA.AACTGGGGCTGGAAAAACGTTTTCAGTTg < 1:1790599/36‑1 (MQ=255)
gC.GGTGA.AACTGGGGCTGGAAAAACGTTTTCAGGTg > 1:2843954/1‑36 (MQ=255)
c.GGTGA.AACTGGGGCTGGAAAAACGTTTTCAGTTgc < 1:4810524/36‑1 (MQ=255)
c.GGTGA.AACTGGGGCTGGAAAAACGTTTTCAGTTgc < 1:5906418/36‑1 (MQ=255)
ccgTGA.AACTGGGGCTGGAAAAGCGTTTTCAGTTgc < 1:79645/34‑1 (MQ=37)
ggTGA.AACTGGGGCTGGAAAAACGTTTTCAGTTgcg > 1:5043584/1‑36 (MQ=255)
ggTGA.AACTGGGGCTGGAAAAACGTTTTCAGTTgcg > 1:1693302/1‑36 (MQ=255)
gTGA.AACTTGGGCTGGAAAAACGTTTTCAGTTgcgc < 1:1142138/36‑1 (MQ=255)
gTGA.AACTGGGGCTGGAAAAACGTTTTCAGTTgcgc > 1:1271737/1‑36 (MQ=255)
gTGA.AACTGGGGCTGGAAAAACGTTTTCAGTTgcgc > 1:2888334/1‑36 (MQ=255)
tGA.AACTGTGGCTGGAAAAGCGTTTTCAGTTGcgcg < 1:1660651/36‑1 (MQ=25)
tGA.AACTGGGGCTGGAAAAACGTTTTCAGTTGcgcg < 1:4266001/36‑1 (MQ=255)
a.aaCTGGGTCTGGAGAAACGTTTTCAGTTGCGCGAt < 1:37462/36‑1 (MQ=37)
a.aaCTGGGGCTGGAAAAACGTTTTCAGTTGCGCGAt > 1:638271/1‑36 (MQ=255)
a.aaCTGGGGCTGGAAAAACGTTTTCAGTTGCGCGAt < 1:2384924/36‑1 (MQ=255)
a.aaCTGGGGCTGGAAAAACGTTTTCAGTTGCGCGAt > 1:2532993/1‑36 (MQ=255)
a.aaCTGGGGCTGGAAAAACGTTTTCAGNNGCGCGAt > 1:2437183/1‑36 (MQ=255)
aaCTGGGGCTGGAAAAACGTTTTCAGTTGCGCGATg > 1:3569296/1‑36 (MQ=255)
aaCTGGGGCTGGAAAAACGTTTTCAGTTGCGCGATg > 1:287523/1‑36 (MQ=255)
aCTGGGGCTGGAAAAACGTTTTCAGTTGCGCGATGa > 1:2122137/1‑36 (MQ=255)
aCTGGGGCTGGAAAAACGTTTTCAGTTGCGCGATGa < 1:4791490/36‑1 (MQ=255)
aCTGGGGCTGGAAAAACGTTTTCAGATGCGCGATGa < 1:3646449/36‑1 (MQ=255)
cTGGGGCTGGAAAAACGTTTTCAGTTGCGCGATGAc < 1:4224646/36‑1 (MQ=255)
cTGGGGCTGGAAAAACGTTTTCAGTTGCGCGATGAc > 1:3207572/1‑36 (MQ=255)
tGGGGCTGGAAAAACGTTTTCAGTTGCGCGATGACg < 1:3780398/36‑1 (MQ=255)
gggggCTGGAAAAACGTTTTCAGTTGCGCGATGACg < 1:3344639/35‑1 (MQ=255)
gggCTGGAAAAACGTTTTCAGTTGCGCGATGACGAg < 1:3269157/36‑1 (MQ=255)
ggCTGGAAAACCGTTTTCATTTGCGCGATGACGAGt < 1:2728480/36‑1 (MQ=37)
ggCTGGAAAAACGTTTTCAGTTGCGCGATGACGAGt < 1:5227466/36‑1 (MQ=255)
ggCTGGAAAAACGTTTTCAGTTGCGCGATGACGAGt < 1:3182888/36‑1 (MQ=255)
ggCTGGAAAAACGTTTTCAGTTGCGCGATGACGAGt > 1:1207338/1‑36 (MQ=255)
ggCTGGAAAAACGTTTTCAGTTGCGCGATGACGAGt < 1:2263815/36‑1 (MQ=255)
gCTTGAAAAACGTTTTCAGTTGCGCGATGACGAGtt < 1:2987946/36‑1 (MQ=255)
gCTGGAAAACCGTTTTCAGTGGCGCGATGACGAGtt < 1:2156379/36‑1 (MQ=38)
gCTGGAAAAACGTTTTCAGTTGCGCGATGCCGAGtt < 1:5567448/36‑1 (MQ=255)
gCTGGAAAAACGTTTTCAGTTGCGCGATGACGAGtt < 1:5567741/36‑1 (MQ=255)
gCTGGAAAAACGTTTTCAGTTGCGCGATGACGAGtt < 1:853519/36‑1 (MQ=255)
cTGGAAAAACGTTTTCAGTTGCGCGATGACGAGttt > 1:4607416/1‑36 (MQ=255)
|
CACATGGCGAGCATCGAATGGC.GGTGA.AACTGGGGCTGGAAAAACGTTTTCAGTTGCGCGATGACGAGTTT > REL606/2735959‑2736029
Base quality scores: 
ATCG
< 3 ≤
ATCG
< 7 ≤
ATCG
< 16 ≤
ATCG
< 29 ≤
ATCG
< 37 ≤
ATCG
|