Read alignment evidence... | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|
seq id | position | change | freq | score | reads | annotation | genes | product | ||
* | REL606 | 2,724,423 | 0 | A→C | 8.5% | 8.7 ‑47.9 0.0 | 263 | T173P (ACC→CCC) | srlR | DNA‑bindng transcriptional repressor |
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/39): ref base (171/52): total (171/92) | ||||||||||
Fisher's exact test strand distribution p-value = 1.60e-21 | ||||||||||
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00 |
Only 100 of 299 total aligned reads displayed.
CGAGCATTTCACCTTCGATAAATTGTTTATGGGCACCGACGGCATCGATCTCAATGCGGGCGTAACCAC > REL606/2724389‑2724457 | cGAGCATTTCACCTTCGATAAATTGTTTATGGGCa < 1:1828396/35‑1 (MQ=255) cGAGCATTTCACCTTCGATAAATTGTTTATGGGCa > 1:15758580/1‑35 (MQ=255) cGAGCATTTCACCTTCGATAAATTGTTTATGGGCa > 1:9153987/1‑35 (MQ=255) cGAGCATTTCACCTTCGATAAATTGTTTATGGGCa < 1:2577720/35‑1 (MQ=255) cGAGCATTTCACCTTCGATAAATTGTTTATGGGCa > 1:4184760/1‑35 (MQ=255) cGAGCATTTCACCTTCGATAAATTGTTTATGGGCa > 1:6472668/1‑35 (MQ=255) cGAGCATTTCAACTTCGATAAATTGTTTATGGGCa < 1:13936453/35‑1 (MQ=255) gggCATTTCCCCTTCGATAAATTGTTTATGGGCAc < 1:23214/33‑1 (MQ=38) gAGCATTTCACCTTCGATAAATTGTTTATGGGCAc > 1:18672201/1‑35 (MQ=255) gAGCATTTCACCTTCGATAAATTGTTTATGGGCAc > 1:19922742/1‑35 (MQ=255) gAGCATTTCACCTTCGATAAATTGTTTATGGGCAc > 1:16542815/1‑35 (MQ=255) gAGCATTTCACCTTCGATAAATTGTTTATGGGCAc < 1:13202474/35‑1 (MQ=255) gAGCATTTCACCTTCGATAAATTGTTTATGGGCAc > 1:10047716/1‑35 (MQ=255) aGCATTTCCCCTTCGATAAATTGTTTATGGGCAcc < 1:909628/35‑1 (MQ=255) aGCATTTCCCCTTCGATAAATTGTTTATGGGCAcc < 1:1268961/35‑1 (MQ=255) aGCATTTCACCTTCGATAAATTGTTTATGGGCAcc < 1:6243484/35‑1 (MQ=255) aGCATTTCACCTTCGATAAATTGTTTATGGGCAcc < 1:14291428/35‑1 (MQ=255) gCATTTCCCCTTCGATAAATTGTTTATGGGCACCg < 1:20385350/35‑1 (MQ=255) gCATTTCCCCTTCGATAAATTGTTTATGGGCACCg < 1:16088396/35‑1 (MQ=255) gCATTTCACCTTCGATAAATTGTTTATGGTCACCg > 1:1669307/1‑35 (MQ=255) gCATTTCACCTTCGATAAATTGTTTATGGGCACCg > 1:15959378/1‑35 (MQ=255) gCATTTCACCTTCGATAAATTGTTTATGGGCACCg > 1:9340219/1‑35 (MQ=255) gCATTTCACCTTCGATAAATTGTTTATGGGCACCg > 1:15659702/1‑35 (MQ=255) gCATTTCACCTTCGATAAATTGTTTATGGGCACCg > 1:6764591/1‑35 (MQ=255) gCATTTCACCTTCGATAAATTGTTTATGGGCACCg < 1:12472625/35‑1 (MQ=255) cATTTCCCCTTCGATAAATTGTTTATGGGCACCGa < 1:19145836/35‑1 (MQ=255) cATTTCACCTTCGATAAATTGTTTATGGGCACCGa < 1:1201948/35‑1 (MQ=255) cATTTCACCTTCGATAAATTGTTTATGGGCACCGa > 1:16560350/1‑35 (MQ=255) cATTTCACCTTCGATAAATTGTTTATGGGCACCGa > 1:7218006/1‑35 (MQ=255) cATTTCAACTTCGATAAATTGTTTATGGGCACCGa < 1:14834217/35‑1 (MQ=255) aTTTCCCCTTCGATAAATTGTTTATGGGCACCGAc < 1:13580489/35‑1 (MQ=255) tttCACCTTCGATAAATTGTTTATGGGCACCGCCg > 1:10292230/1‑35 (MQ=255) ttCACCTTCGATAAATTGTTTATGGGCACCGACgg > 1:6375954/1‑35 (MQ=255) ttCACCTTCGATAAATTGTTTATGGGCACCGACgg > 1:5551480/1‑35 (MQ=255) ccccTTCGATAAATTGTTTATGGGCACCGACGGCa < 1:17397909/33‑1 (MQ=255) cACCTTCGATAAATTGTTTATGGGCACCGACGGCa > 1:13118329/1‑35 (MQ=255) ccTTCGATAAATTGTTTATGGGCACCGACGGCATc > 1:8715081/1‑35 (MQ=255) ccTTCGATAAATTGTTTATGGGCACCGACGGCATc > 1:10202168/1‑35 (MQ=255) cTTCGATAAATTGTTTATGGGCACCGACGGCATCg > 1:10775840/1‑35 (MQ=255) cTTCGATAAATTGTTTATGGGCACCGACGGCATCg > 1:4697738/1‑35 (MQ=255) cGATAAATTGTTTATGGGCACCGCCCGCATTGaga < 1:9755981/35‑3 (MQ=21) cGATAAATTGTTTATGGGCACCGACGGCATCGAtc > 1:4953406/1‑35 (MQ=255) cGATAAATTGTTTATGGGCACCGACGGCATCGAtc > 1:7352697/1‑35 (MQ=255) cGAAAAATTGTTTATGGGCACCGACGGCATCGAtc > 1:16811313/1‑35 (MQ=255) gATAAATTTTTTGTGGGCCCCGACGGCATCGAtct < 1:20073535/35‑1 (MQ=25) gATAAATTGTTTATGGGCCCCGACGGCATCGAtct < 1:5996820/35‑1 (MQ=255) gATAAATTGTTTATGGGCACCGACGGCATCGAtct > 1:5101506/1‑35 (MQ=255) gATAAATTGTTTATGGGCACCGACGGCATCGAtct > 1:16329218/1‑35 (MQ=255) gATAAATTGTTTATGGGCACCGACGGCATCGAtct > 1:17311218/1‑35 (MQ=255) gATAAATTGTTTATGGGCACCGACGGCATCGAtct > 1:18492599/1‑35 (MQ=255) gATAAATTGTTTATGGGCACCGACGGCATCGAtct > 1:18151980/1‑35 (MQ=255) aTAAATTGTTTGTGGGCCCCGACGGCATCGAtctc < 1:18977863/35‑1 (MQ=37) aTAAATTGTTTATGGGCACCGACGGCATCGAtctc > 1:9026694/1‑35 (MQ=255) aTAAATTGTTTATGGGCACCGACGGCATCGAtctc > 1:13793969/1‑35 (MQ=255) aTAAATTGTTTATGGGCACCGACGGCATCGAtctc > 1:12741056/1‑35 (MQ=255) aTAAATTGTTTATGGGCACCGACGGCATCGAtctc > 1:5890729/1‑35 (MQ=255) tAAATTGTTTATGGGCACCGACGGCATCGATCTCa > 1:17699489/1‑35 (MQ=255) tAAATTGTTTATGGCCCCCGACGGCATCGATCTCa < 1:11738498/35‑1 (MQ=37) aaaTTGTTTATGGGCACCGACGGCATCGATCTCaa > 1:15531428/1‑35 (MQ=255) aaTTGTTTATGGGCCCCGACGGCATCGATCTCAAt < 1:1232660/35‑1 (MQ=255) aaTTGTTTATGGGCCCCGACGGCATCGATCTAAAt < 1:9401354/35‑1 (MQ=37) ttttttATGGGCCCCGACGGCATCGATCTCAATGc < 1:4109452/32‑1 (MQ=37) ttGTTTATGGGCACCGACGGCATCGATCTCAATGc > 1:3509572/1‑35 (MQ=255) tGTTTATGGGCCCCGACGGCATCGATCTCAATGCg < 1:12155031/35‑1 (MQ=255) tGTTTATGGGCACCGACGGCATCGATCTCAATGCg > 1:3749236/1‑35 (MQ=255) tGTTTATGGGCACCGACGGCATCGATCTCAATGCg > 1:3910147/1‑35 (MQ=255) tGTTTATGGACCCCGACGGCATCGATCTCAATGCg < 1:13718249/35‑1 (MQ=37) tGTTTAGGGGCACCGACGGCATCGATCTCAATGCg < 1:5276307/35‑1 (MQ=255) ttttAGGGGCACCGACGGCATCGATCTCAATGCgg < 1:8548709/34‑1 (MQ=255) gTTTATGGGCACCGACGGCATCGATCTCAATGCgg > 1:9522087/1‑35 (MQ=255) tttATGTGCACCGACGGCATCGATCTCAATGCgtg > 1:1437176/1‑33 (MQ=37) tttATGGGCACCGACGGCATCGATCTCAATGCggg > 1:7872341/1‑35 (MQ=255) tttATGGGCACCGACGGCATCGATCTCAATGCggg > 1:7577906/1‑35 (MQ=255) tttATGGGCACCGACGGCATCGATCTCAATGCggg > 1:20177643/1‑35 (MQ=255) tttATGGGCACCGACGGCATCGATCTCAATGCggg < 1:2882302/35‑1 (MQ=255) ttATGGGCACCGACGGCATCGATCTCAATGCGGGc > 1:17517451/1‑35 (MQ=255) ttAGGGGCCCCGACGGCATCGATCTCAATGCGGGc < 1:18835196/35‑1 (MQ=37) tATGGGCACCGACGGCATCGATCTCAATGCGGGCg > 1:16229046/1‑35 (MQ=255) tATGGGCACCGACGGCATCGATCTCAATGCGGGCg > 1:3275102/1‑35 (MQ=255) tATGGGCACCGACGGCATCGATCTCAATGCGGGCg > 1:18064025/1‑35 (MQ=255) tATGGGCACCGACGGCATCGATCTCAATGCGGGCg > 1:17892015/1‑35 (MQ=255) gtGGGCCCCGACGGCATCGATCTCAATGCGGGCGt < 1:4383321/34‑1 (MQ=255) attggCACCGACGGCATCGATCTCAATGCGGGCGt > 1:2123233/4‑35 (MQ=255) aTGGGCCCCGACGGCATCGATCTCAATGCGGGCGt < 1:17249833/35‑1 (MQ=255) aTGGGCACCGACTGCATCGATCTCAATGCGGGCGt > 1:19339393/1‑35 (MQ=255) aTGGGCACCGACGGCATCGATCTCAATGCGGGCGt > 1:14523774/1‑35 (MQ=255) tGGGCACCGACGGCATCGATCTCAATGCGGGCgca > 1:2289124/1‑33 (MQ=255) tGGGCACCGACGGCATCGATCTCAATGCGGGCGTa > 1:10243953/1‑35 (MQ=255) tGGGCACCGACGGCATCGATCTCAATGCGGGCGTa > 1:11627559/1‑35 (MQ=255) tGGGCACCGACGGCATCGATCTCAATGCGGGCGTa > 1:11163847/1‑35 (MQ=255) tGGGCACCGACGGCATCGATCTCAATGCGGGCGTa > 1:15285990/1‑35 (MQ=255) ggggCCCCGACGGCATCGATCTCAATGCGGGCGTa < 1:6650286/34‑1 (MQ=255) gggCACCGACGGCATCGATCTCAATGCGGGCGTaa > 1:13007088/1‑35 (MQ=255) ggCACCGACGGCATCGATCTCAATGCGGGCGTAac > 1:801130/1‑35 (MQ=255) ggCACCGACGGCATCGATCTCAATGCGGGCGTAac < 1:3133704/35‑1 (MQ=255) cACCGACGGCATCGATCTCAATGCGGGCGTAacca > 1:8231965/1‑35 (MQ=255) cACCGACGGCATCGATCTCAATGCGGGCGTAacca > 1:4560641/1‑35 (MQ=255) cACCGACGGCATCGATCTCAATGCGGGCGTAacca > 1:3620290/1‑35 (MQ=255) aCCGACGGCATCGATCTCAATGCGGGCGTAaccac > 1:19657384/1‑35 (MQ=255) aCCGACGGCATCGATCTCAATGCGGGCGTAaccac > 1:17365023/1‑35 (MQ=255) | CGAGCATTTCACCTTCGATAAATTGTTTATGGGCACCGACGGCATCGATCTCAATGCGGGCGTAACCAC > REL606/2724389‑2724457 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ 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. |