Read alignment evidence... | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|
seq id | position | change | freq | score | reads | annotation | genes | product | ||
* | REL606 | 2,757,213 | 0 | A→T | 4.9% | 4.1 ‑17.3 ‑0.0 | 241 | I7F (ATC→TTC) | ygbK | hypothetical protein |
Rejected: Polymorphism prediction not supported by minimum number of reads on both strands. | ||||||||||
Rejected: Prediction has biased read strand distribution supporting polymorphism. | ||||||||||
Rejected: Polymorphic base substitution creates a homopolymer stretch. | ||||||||||
Reads supporting (aligned to +/- strand): new base (0/13): ref base (173/55): total (173/68) | ||||||||||
Fisher's exact test strand distribution p-value = 2.94e-08 | ||||||||||
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00 |
Only 100 of 271 total aligned reads displayed.
TACCGGGAGCGAAATCATGATCAAGATTGGCGTTATCGCCGATGATTTTACCGGCGCGACGGATATCGC > REL606/2757179‑2757247 | tACCGGGAGCGAAATCATGATCAAGATTGGCGTTa < 1:2368582/35‑1 (MQ=255) aCCGGGAGCGAAATCATGATCAAGATTGGCGTTAt > 1:7029158/1‑35 (MQ=255) aCCGGGAGCGAAATCATGATCAAGATTGGCGTTAt < 1:1257009/35‑1 (MQ=255) aCCGGGAGCGAAATCATGATCAAGATTGGCGTTAt > 1:2696640/1‑35 (MQ=255) aCCGGGAGCGAAATCATGATCAAGATTGGCGTTAt < 1:10470521/35‑1 (MQ=255) ccGGGAGCGGAATCATGATCAAGATTGGCGTTATc < 1:15330558/35‑1 (MQ=255) ccGGGAGCGAAATCATGATCAAGATTGGCGTTATc > 1:16777833/1‑35 (MQ=255) ccGGGAGCGAAATCATGATCAAGATTGGCGTTATc < 1:4049962/35‑1 (MQ=255) ccGGGAGCGAAATCATGATCAAGATTGGCGTTATc > 1:3300480/1‑35 (MQ=255) cGGGAGCGAAATCATGATCAAGATTGGCGTTATCg < 1:11457460/35‑1 (MQ=255) cGGGAGCGAAATCATGATCAAGATTGGCGTTATCg > 1:1334919/1‑35 (MQ=255) gggAGCGAAATCATGATCAAGATTGGCGTTATCGc < 1:14641673/35‑1 (MQ=255) gggAGCGAAATCATGATCAAGATTGGCGTTATCGc > 1:14724872/1‑35 (MQ=255) gggAGCGAAATCATGATCAAGATTGGCGTTATCGc < 1:11794328/35‑1 (MQ=255) gggAGCGAAATCATGATCAAGATTGGCGTTATCGc < 1:13862240/35‑1 (MQ=255) gggAGCGAAATCATGATCAAGATTGGCGTTATCGc > 1:16128305/1‑35 (MQ=255) ggAGCGAAATCATGATCAAGATTGGCGTTATCGcc > 1:1589769/1‑35 (MQ=255) gAGCGAAATCATGATCAAGATTGGCGTTATCGCCg > 1:870311/1‑35 (MQ=255) gAGCGAAATCATGATCAAGATTGGCGTTATCGCCg > 1:6949836/1‑35 (MQ=255) gAGCGAAATCATGATCAAGATTGGCGTTATCGCCg < 1:19384309/35‑1 (MQ=255) aGCGAAATCATGATCAAGATTGGCGTTATCGCCga > 1:12109353/1‑35 (MQ=255) aGCGAAATCATGATCAAGATTGGCGTTATCGCCga < 1:11014811/35‑1 (MQ=255) aGCGAAATCATGATCAAGATTGGCGTTATCGCCga > 1:5851663/1‑35 (MQ=255) gCGAAATCATGATCAAGATTGGCGTTATCGCCgat > 1:5662145/1‑35 (MQ=255) gCGAAATCATGATCAAGATTGGCGTTATCGCCgat > 1:20160110/1‑35 (MQ=255) gCGAAATCATGATCAAGATTGGCGTTATCGCCgat > 1:9678387/1‑35 (MQ=255) gCGAAATCATGATCAAGATTGGCGTTATCGCCgat > 1:218511/1‑35 (MQ=255) cGAAATCATGATCAAGATTGGCGTTATCGCCgatg > 1:18286697/1‑35 (MQ=255) gAAATCATGATCAAGATTGGCGTTATCGCCgatga > 1:2584322/1‑35 (MQ=255) gAAATCATGATCAAGATTGGCGTTATCGCCgatga > 1:10866409/1‑35 (MQ=255) gAAATCATGATCAAGATTGGCGTTATCGCCgatga > 1:1036683/1‑35 (MQ=255) aaaTCATGATCAAGATTGGCGTTATCGCCgatgat > 1:20380232/1‑35 (MQ=255) aaaTCATGATCAAGATTGGCGTTATCGCCgatgat > 1:137844/1‑35 (MQ=255) aaaTCATGATCAAGATTGGCGTTATCGCCgatgat > 1:2832132/1‑35 (MQ=255) aaTCATGATCAAGATTGGCGTTATCGCCGATGAtt < 1:1672349/35‑1 (MQ=255) aTCATGATCAAGATTGGCGTTATCGCCGATGAttt > 1:17875081/1‑35 (MQ=255) tCATGATCAAGATTGGCGTTATCGCCGATGAtttt > 1:5358952/1‑35 (MQ=255) tCATGATCAAGATTGGCGTTATCGCCGATGAtttt > 1:1485042/1‑35 (MQ=255) tCATGATCAAGATTGGCGTTATCGCCGATGAtttt > 1:10214416/1‑35 (MQ=255) cATGATCAAGATTGGCGTTATCGCCGATGATTTTa > 1:1715190/1‑35 (MQ=255) aTGATCAAGATTGGCGTTATCGCCGATGATTTTAc > 1:9011892/1‑35 (MQ=255) tGATCAAGATTGGCGTTATCGCCGATGATTTTAcc > 1:5093565/1‑35 (MQ=255) tGATCAAGATTGGCGTTATCGCCGATGATTTTAcc > 1:3199695/1‑35 (MQ=255) gATCAAGATTGGCTTTTTCGCCGATGATTTTACCg < 1:1491561/35‑1 (MQ=37) gATCAAGATTGGCGTTATCGCCGATGATTTTACCg > 1:13656541/1‑35 (MQ=255) gATCAAGATTGGCGTTATCGCCGATGATTTTACCg < 1:12308987/35‑1 (MQ=255) gATCAAGATTGGCGTTATCGCCGATGATTTTACCg < 1:19547427/35‑1 (MQ=255) ataaaGATTGGCTTTTTCGCCGATGATTTTACCgg < 1:14464233/32‑1 (MQ=25) agcAAGATTGGCTTTTTCGCCGATGATTTTACCgg < 1:8182469/33‑1 (MQ=25) aTCAAGATTGGCGTTATCGCCGATGATTTTACCgg > 1:8210529/1‑35 (MQ=255) tCAAGATTGGCGTTATCGCCGATGATTTCACCGgc > 1:11842855/1‑35 (MQ=255) cAAGATTGGCGTTATCGCCGATGATTTTACCGgcg > 1:5581738/1‑35 (MQ=255) cAAGATTGGCGTTACCGCCGATGATTTTACCGgcg > 1:7612663/1‑35 (MQ=255) aaGATTGGGTTTATCGCCGATGATTTTACCGgcgc < 1:17262105/35‑1 (MQ=25) aaGATTGGCGTTATCGCCGATGATTTTACCGgcgc > 1:3401703/1‑35 (MQ=255) aaGATTGGCGTTATCGCCGATGATTTTACCGgcgc > 1:5182908/1‑35 (MQ=255) aaGATGGGCTTTTTCGCCGATGATTTTACCGgcgc < 1:10723828/35‑1 (MQ=25) gATTGGCGTTATCGCCGATGATTTTACCGGCGCGa > 1:7741093/1‑35 (MQ=255) gATTGGCGTTATCGCCGATGATTTTACCGGCGCGa > 1:8473316/1‑35 (MQ=255) gATTGGCGTTATCGCCGATGATTTTACCGGCGCGa > 1:951077/1‑35 (MQ=255) gATGGGCTTTTTCGCCGATGATTTTACCGGCGCGa < 1:3584314/35‑1 (MQ=25) aTTGTCGTTATCGCCGATGATTTTACCGGCGCGAc > 1:6886694/1‑35 (MQ=255) aTTGGCGTTATCGCCGATGATTTTACCGGCGCGAc > 1:5170128/1‑35 (MQ=255) aTTGGCGTTATCGCCGATGATTTTACCGGCGCGAc > 1:4811432/1‑35 (MQ=255) aTTGGCGTTATCGCCGATGATTTTACCGGCGCGAc > 1:4852836/1‑35 (MQ=255) aTTGGCGTTATCGCCGATGATTTTACCGGCGCGAc > 1:9427089/1‑35 (MQ=255) ttGGCTTTATCGCCGATGATTTTACCGGCGCGACg < 1:15991060/35‑1 (MQ=255) ttGGCGTTATCGCCGATGATTTTACCGGCGCGACg > 1:20035782/1‑35 (MQ=255) tgggCTTTTTCGCCGATGATTTTACCGGCGCGACg < 1:2416009/33‑1 (MQ=25) tgggCGTTATCGCCGATGATTTTACCGGCGCGACg < 1:4592624/33‑1 (MQ=255) tGGCGTTATCGCCGATGATTTTACCGGCGCGACgg > 1:11605629/1‑35 (MQ=255) gggCGTTATCGCCGATGATTTTACCGGCGCGACgg < 1:15094285/34‑1 (MQ=255) ggCTTTTTCGCCGATGATTTTACCGGCGCGACGGa < 1:10023762/35‑1 (MQ=37) ggCGTTTTCGCCGATGATTTTACCGGCGCGACGGa < 1:8633056/35‑1 (MQ=255) ggCGTTATCGCCGATGATTTTACCGGCGCGACGGa > 1:19871086/1‑35 (MQ=255) gctttATCGCCGATGATTTTACCGGCGCGACGGat < 1:9897551/32‑1 (MQ=255) gCGTTATCGCCGATGATTTTACCGGCGCGACGGat > 1:6222833/1‑35 (MQ=255) gCGTTATCGCCGATGATTTTACCGGCGCGACGGat > 1:4406800/1‑35 (MQ=255) ctttATCGCCGATGATTTTACCGGCGCGACGGata < 1:8586011/33‑1 (MQ=255) ctttATCGCCGATGATTTTACCGGCGCGACGGata < 1:17522879/33‑1 (MQ=255) cGTTATCGCCGATGATTTTACCGGCGCGACGGata > 1:15783579/1‑35 (MQ=255) cGTTATCGCCGATGATTTTACCGGCGCGACGGata > 1:796661/1‑35 (MQ=255) cGTTATCGCCGATGATTTTACCGGCGCGACGGata > 1:3962649/1‑35 (MQ=255) gTTTTCGCCGATGATTTTACCGGCGCGACGGatat < 1:13972007/35‑1 (MQ=255) gTTTTCGCCGATGATTTTACCGGCGCGACGGatat < 1:14325449/35‑1 (MQ=255) gTTATCGCCGATGATTTTACCGGCGCGACGGatat > 1:10153589/1‑35 (MQ=255) gTTATCGCCGATGATTTTACCGGCGCGACGGatat > 1:7301625/1‑35 (MQ=255) gTTATCGCCGATGATTTTACCGGCGCGACGGatat > 1:20038793/1‑35 (MQ=255) ttATCGCCGATGATTTTACCGGCGCGACGGATATc > 1:19161305/1‑35 (MQ=255) ttATCGCCGATGATTTTACCGGCGCGACGGATATc > 1:12658618/1‑35 (MQ=255) ttATCGCCGATGATTTTACCGGCGCGACGGATATc > 1:12962321/1‑35 (MQ=255) ttATCGCCGATGATTTTACCGGCGCGACGGATATc > 1:20296852/1‑35 (MQ=255) tATCGCCGATGATTTTACCGGCGCGACGGATATCg > 1:13542792/1‑35 (MQ=255) tATCGCCGATGATTTTACCGGCGCGACGGATATCg > 1:14979972/1‑35 (MQ=255) tATCGCCGATGATTTTACCGGCGCGACGGATATCg > 1:16917316/1‑35 (MQ=255) tATCGCCGATGATTTTACCGGCGCGACGGATATCg > 1:18057810/1‑35 (MQ=255) tATCGCCGATGATTTTACCGGCGCGACGAATATCg < 1:1890960/35‑1 (MQ=255) aTCGCCGATGATTTTACCGGCGCGACGGATATCGc > 1:16442457/1‑35 (MQ=255) aTCGCCGATGATTTTACCGGCGCGACGGATATCGc < 1:305223/35‑1 (MQ=255) aTCGCCGATGATTTTACCGGCGCGACGGATATCGc > 1:20577907/1‑35 (MQ=255) | TACCGGGAGCGAAATCATGATCAAGATTGGCGTTATCGCCGATGATTTTACCGGCGCGACGGATATCGC > REL606/2757179‑2757247 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |
Reads not counted as support for junction |
read_name Not counted due to insufficient overlap past the breakpoint. |
read_name Not counted due to not crossing MOB target site duplication. |