Read alignment evidence... | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|
seq id | position | change | freq | score | reads | annotation | genes | product | ||
* | REL606 | 3,660,230 | 0 | T→G | 100.0% | 7.3 | 15 | T332P (ACC→CCC) | xylA | xylose isomerase |
Rejected: E-value exceeds prediction threshold. | ||||||||||
Reads supporting (aligned to +/- strand): new base (9/0): ref base (4/2): total (13/2) | ||||||||||
Fisher's exact test strand distribution p-value = 1.43e-01 | ||||||||||
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00 |
GTACTTTGGCATCGAAGTTCAGACCACCGGTGGTGAAACCGCCTGCTTTGAGAATTTCATACATCACC > REL606/3660197‑3660264 | gTACTTTGGCATCGAAGTTCAGCCCACCGgtggtg < 1:7351476/35‑1 (MQ=255) gTACTTTGGCATCGAAGTTCAGCCCACCGgtggtg < 1:6816980/35‑1 (MQ=255) cTTTGGCATCGAAGTTCAGACCACCGGTGGGGaaa > 1:16914230/1‑35 (MQ=255) cTTTGGCATCGAAGTTCAGACCACCGGTGGGGaaa > 1:5058150/1‑35 (MQ=255) cTTGGGCATCGAAGTTCAGCCCACCGGTGGTGaaa < 1:1824826/35‑1 (MQ=38) tttGGCATCGAAGTTCAGCCCACCGGTGGTGAAAc < 1:12262739/35‑1 (MQ=255) tttGGCATCGAAGTTCAGACCACCGGTGGGGaaaa > 1:760180/1‑34 (MQ=255) tttGGCATCGAAGTTCAGACCACCGGTGGGGGAAc > 1:12202052/1‑35 (MQ=37) tttGGCATCGAAGTTCAGACCACCGGTGGGGAAAc > 1:9083929/1‑35 (MQ=255) ttGGCATCGAAGTTCAGACCACCGGTGGGGGAAcc > 1:12074276/1‑35 (MQ=37) ttGGCATCGAAGTTCAGACCACCGGTGGGGGAAcc > 1:7155803/1‑35 (MQ=37) ttGGCATCGAAGTTCAGACCACCGGTGGGGAAAca > 1:6007984/1‑34 (MQ=255) ggCATCGAAGTTCAGACCACCGGTGGGGAAACCGc > 1:18006598/1‑35 (MQ=255) ggCATCGAAGTTCAGACCACCGGTGGAGGAAACGc > 1:15147792/1‑35 (MQ=25) gCATCGAAGTTCAGACCACCGGTGGGGGAACCGcc > 1:5493418/1‑35 (MQ=37) gCATCGAAGTTCAGACCACCGGTGGGGAAAACGcc > 1:2126261/1‑35 (MQ=37) gCATCGAAGTTCAGACCACCGGTGGGGAAAACGcc > 1:14463657/1‑35 (MQ=37) gCATCGAAGTTCAGACCACCGGTGGGGAAAACGcc > 1:10448651/1‑35 (MQ=37) cATCGAAGTTCAGACCACCGGTGGGGGAACCGccc > 1:17954832/1‑34 (MQ=37) cATCGAAGTTCAGACCACCGGTGGGGAAACCGCCt > 1:17441631/1‑35 (MQ=255) cATCGAAGTTCAGACCACCGGTGGGGAAACCCccc > 1:16103988/1‑34 (MQ=37) cATCGAAGTTCAGACCACCGGTGGGGAAAACGCCt > 1:6301507/1‑35 (MQ=37) aTCGAAGTTCAGACCACCGGTGGGGGAACCGCCTg > 1:14793622/1‑35 (MQ=37) tCGAAGTTCAGACCACCGGTGGGGGAAACGCCTGc > 1:14123183/1‑35 (MQ=25) tCGAAGTTCAGACCACCGGTGGGGAAAACGCCTGc > 1:3977700/1‑35 (MQ=37) cGAAGTTCAGACCACCGGTGGGGAAACCGCCTGCt > 1:5864991/1‑35 (MQ=255) cGAAGTTCAGACCACCGGTGGGGAAACCGCCTGCt > 1:20075575/1‑35 (MQ=255) cGAAGTTCAGACCACCGGTGGGGAAACCGCCTGCt > 1:6806547/1‑35 (MQ=255) gAAGTTCAGACCACCGGTGGGGGAACCGCCTGCtt > 1:831635/1‑35 (MQ=37) gAAGTTCAGACCACCGGTGGGGGAACCGCCTGCtt > 1:11602797/1‑35 (MQ=37) aaGTTCAGACCACCGGTGGGGAAACCGCCTGCttt > 1:3922162/1‑35 (MQ=255) aGTTCAGACCACCTGTGGTGAAACCGCCTGCTTTg > 1:18018905/1‑35 (MQ=255) aGTTCAGACCACCGGTGGGGGAACCGCCTGCtttt > 1:9099438/1‑34 (MQ=37) aGTTCAGACCACCGGTGGGGGAAACCCCTGCtttt > 1:6731719/1‑34 (MQ=16) gTTCAGACCACCGGTGGGGAAACCCCCGCCtttag > 1:6484091/1‑33 (MQ=16) ttCAGACCACCGGTGGGGGAACCGCCCGCTTTtag > 1:7873503/1‑35 (MQ=19) aGACCACCGGTGGGGAAACCGCCTGCTTTGAGAAt > 1:3678373/1‑35 (MQ=255) gACCACCGGTGTTGAAACCGCCTGCTTTGAGAAtt > 1:10980168/1‑35 (MQ=255) gACCACCGGTGGGGAAACCGCCTGCTTTTAGAAtt > 1:14039436/1‑35 (MQ=37) ccGGTGGTGAAACCGCCTGCTTTGAGAATTTCATa > 1:1206900/1‑35 (MQ=255) ccGGTGGTGAAACCGCCTGCTTTGAGAATTTCATa > 1:8431158/1‑35 (MQ=255) ggtggtGAAACCGCCTGCTTTGAGAATTTCATACa > 1:4295019/1‑35 (MQ=255) ggtggaGAACCCGCCTGCTTTGAGAATTTCATACa < 1:1775487/35‑1 (MQ=37) gggggtaAACCCGCCTGCTTTGAGAATTTCATACa < 1:9967918/32‑1 (MQ=25) gtggtgAAACCGCCTGCTTTGAGAATTTCATACAt > 1:9843137/1‑35 (MQ=255) tggtgAAACCGCCTGCTTTGAGAATTTCATACATc > 1:2138006/1‑35 (MQ=255) tggtgAAACCGCCTGCTTTGAGAATTTCATACATc > 1:1216119/1‑35 (MQ=255) tggagAACCCGCCTGCTTTGAGAATTTCATACATc < 1:18063526/35‑1 (MQ=37) ggtgAACCCGCCTGCTTTGAGAATTTCATACATCa < 1:20206518/35‑1 (MQ=255) ggtgAAACCGCCTGCTTTGAGAATTTCATACATCa > 1:18241310/1‑35 (MQ=255) ggtgAAACCGCCTGCTTTGAGAATTTCATACATCa > 1:16615970/1‑35 (MQ=255) ggtgAAACCGCCTGCTTTGAGAATTTCATACATCa > 1:1044075/1‑35 (MQ=255) gtgAAACCGCCTGCTTTGAGAATTTCATACATCAc > 1:19305241/1‑35 (MQ=255) gtgAAACCGCCTGCTTTGAGAATTTCATACATCAc > 1:19266666/1‑35 (MQ=255) gtgAAACCGCCTGCTTTGAGAATTTCATACATCAc > 1:8240841/1‑35 (MQ=255) gtgAAACCGCCTGCTTTGAGAATTTCATACATCAc > 1:15103237/1‑35 (MQ=255) gtgAAACCGCCTGCTTTGAGAATTTCATACATCAc > 1:12399001/1‑35 (MQ=255) gtgAAACCGCATGCTTTGAGAATTTCATACATCAc < 1:4939667/35‑1 (MQ=255) tgAAACCGCCTGCTTTGAGAATTTCATACATCAcc > 1:581769/1‑35 (MQ=255) tgAAACCGCCTGCTTTGAGAATTTCATACATCAcc > 1:14790025/1‑35 (MQ=255) | GTACTTTGGCATCGAAGTTCAGACCACCGGTGGTGAAACCGCCTGCTTTGAGAATTTCATACATCACC > REL606/3660197‑3660264 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ 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. |