Read alignment evidence... | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|
seq id | position | change | freq | score | reads | annotation | genes | product | ||
* | REL606 | 2,721,831 | 0 | C→T | 8.1% | 25.8 ‑42.0 ‑0.0 | 233 | S194F (TCT→TTT) | srlE | glucitol/sorbitol‑specific enzyme IIB component of PTS |
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 (19/0): ref base (15/198): total (34/199) | ||||||||||
Fisher's exact test strand distribution p-value = 5.47e-19 | ||||||||||
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.97e-01 |
Only 100 of 248 total aligned reads displayed.
TTCGTCTCGGCGCTCATTGGCATCATTATGGCTTCTGGCCTTGGTGACTGGATTGCCCACGGCCTTGCT > REL606/2721797‑2721865 | tCGTCTCGGCGCTCATTGGCATCATTATGGCTTCt < 1:16522782/35‑1 (MQ=255) tCGTCTCGGCGCTCATTGGCATCATTATGGCTTCt < 1:1585676/35‑1 (MQ=255) tCGTCTCGGCGCTCATTGGCATCATTATGGCTTCt < 1:4155697/35‑1 (MQ=255) tCGTCTCGGCGCTCATTGGCATCATTATGGCTTCt < 1:5776457/35‑1 (MQ=255) cGTCTCGGCGCTCATTGGCATCATTATGGCTTCTg < 1:17697498/35‑1 (MQ=255) cGTCTCGGCGCTCATTGGCATCATTATGGCTTCTg < 1:3617979/35‑1 (MQ=255) cGTCTCGGCGCTCATTGGCATCATTATGGCTTCTg < 1:4740266/35‑1 (MQ=255) gTCTCGGCGCTCATTGGCATCATTATGGCTTCTgg < 1:16940916/35‑1 (MQ=255) gTCTCGGCGCTCATTGGCATCATTATGGCTTCTgg < 1:7999051/35‑1 (MQ=255) gTCTCGGCGCTCATTGGCATCATTATGGCTTCTgg < 1:9730548/35‑1 (MQ=255) tctcGGCGCTCATTGGCATCATTATGGCTTCTGGc < 1:850770/35‑1 (MQ=255) tctcGGCGCTCATTGGCATCATTATGGCTTCTGGc < 1:10253179/35‑1 (MQ=255) ctcGGCGCTCATTGGCATCATTATGGCTTCTGGcc < 1:15684439/35‑1 (MQ=255) ctcGGCGCTCATTGGCATCATTATGGCTTCTGGcc < 1:20373034/35‑1 (MQ=255) ctcGGCGCTCATTAGCATCATTATGGCTTTTggtc > 1:7290535/1‑33 (MQ=25) tcGGCGCTCATTGGCATCATTATGGCTTTTTGCCt > 1:11097531/1‑35 (MQ=25) tcGGCGCTCATTGGCATCATTATGGCTTCTGGCCt < 1:16115750/35‑1 (MQ=255) tcGGCGCTCATTGGCATCATTATGGCTTCTGGCCt < 1:8426705/35‑1 (MQ=255) tcGGCGCTCATTGGCATCATTATGGCTTCTGGCCt < 1:952214/35‑1 (MQ=255) tcGGCGCTCATTGGCATCATTATGGCTTCTGGCCt < 1:19787271/35‑1 (MQ=255) tcGGAGCTCATTGGCATCATTATGGCTTCTGGCCt < 1:3880046/35‑1 (MQ=255) cGGCGCTCATTGGCATCATTATGGCTTCTGGCCtt < 1:10379926/35‑1 (MQ=255) cGGCGCTCATTGGCATCATTATGGCTTCTGGCCtt < 1:6880924/35‑1 (MQ=255) cGGCGCTCATTGGCATCATTATGGCTTCTGGCCtt < 1:7617623/35‑1 (MQ=255) cGGCGCTCATTGGCATCATTATGGCTTCTGGCCtt < 1:19249570/35‑1 (MQ=255) ggCGCTCATTGGCATCATTATGGCTTCTGGCCTTg < 1:6058828/35‑1 (MQ=255) ggCGCTCATTGGCATCATTATGGCTTCTGGCCTTg < 1:11973403/35‑1 (MQ=255) ggCGCCCATTGGCATCATTATGGCTTCTGGCCTTg < 1:2127065/35‑1 (MQ=255) ggCACTCATTGGCATCATTATGGCTTCTGGCCTTg < 1:1650613/35‑1 (MQ=255) gcgcTCATTGGCATCATTATGGCTTTTTTCCTTgt > 1:8264990/1‑34 (MQ=21) gcgcTCATTGGCATCATTATGGCTTCTGGCCTTgg < 1:7511504/35‑1 (MQ=255) cgcTCATTTGCATCATTATGACTTCTGGCCTTGGt > 1:18697506/1‑35 (MQ=37) cgcTCATTGGCATCATTATGGCTTTTTGTCTTGGt > 1:9156774/1‑35 (MQ=25) cgcTCATTGGCATCATTATGGCTTTTGGCTTTTGt > 1:17498398/1‑35 (MQ=25) cgcTCATTGGCATCATTATGGCTTCTGGCCTTGGt < 1:9238110/35‑1 (MQ=255) cgcTCATTGGCATCATTATGGCTTCTGGCCTTGGt < 1:11265568/35‑1 (MQ=255) cgcTCATTGGCATCATTATGGCTTCTGGCCTTGGt < 1:18895271/35‑1 (MQ=255) gcTCATTGGCATCATTATGGCTTTTTTGCTTGGTg > 1:5365318/1‑35 (MQ=16) gcTCATTGGCATCATTATGGCTTTTTGTCTTGGTg > 1:16794915/1‑35 (MQ=25) gcTCATTGGCATCATTATGGCTTCTGGCCTTGGTg < 1:7018638/35‑1 (MQ=255) gcTCATTGGCATCATTATGGCTTCTGGCATTGGTg < 1:11015207/35‑1 (MQ=255) cTCATTGGCATCATTATGGCTTTTTTGCTTGGTgg > 1:8984643/1‑34 (MQ=12) cTCATTGGCATCATTATGGCTTTTTTCCTTTGTgg > 1:18418136/1‑34 (MQ=12) cTCATTGGCATCATTATGGCTTTTTGTCTTGGTgg > 1:6218229/1‑34 (MQ=21) cTCATTGGCATCATTATGGCTTTTTGTCTTGGTgg > 1:19639015/1‑34 (MQ=21) cTCATTGGCATCATTATGGCTTCTGGCCTTGGTGa < 1:10145221/35‑1 (MQ=255) tCATTGGCATCATTATGGCTTTTTTCCTTTGTggc > 1:10023283/1‑33 (MQ=12) tCATTGGCATCATTATGGCTTCTGGCCTTGGTGAc < 1:7722391/35‑1 (MQ=255) tCATTGGCATCATTATGGCTTCTGGCCTTGGTGAc < 1:5864297/35‑1 (MQ=255) cATTGGCATCATTATGGCTTCTGGCCTTGGTGACt < 1:13852819/35‑1 (MQ=255) cATTGGCATCATTATGGCTTCTGGCCTTGGTGACt < 1:529766/35‑1 (MQ=255) aTTGGCATCATTATGGCTTCTGGCCTTGGTGACTg < 1:10654723/35‑1 (MQ=255) aTTGGCATCATTATGGCTTCTGGCCTTGGTGACTg < 1:6660529/35‑1 (MQ=255) aTTGGCATCATTATGGCTTCTGGCCTTGGTGACTg < 1:12908146/35‑1 (MQ=255) aTTGGCATCATTATGGCTTCTGGCCTTGGTGACTg < 1:7675790/35‑1 (MQ=255) ttGGCATCATTATGGCTTGTGGCCTTGGTGACTgg < 1:1337674/35‑1 (MQ=255) ttGGCATCATTATGGCTTCTGGCCTTGGTGACTgg < 1:6181735/35‑1 (MQ=255) ttGGCATCATTATGGCTTCTGGCCTTGGTGACTgg < 1:18870857/35‑1 (MQ=255) tGGCATCATTATGGCTTCTGGCCTTGGTGACTGGa < 1:8856489/35‑1 (MQ=255) tGGCATCATTATGGCTTCTGGCCTTGGTGACTGGa < 1:19602222/35‑1 (MQ=255) tGGCATCATTATGGCTTCTGGCCTTGGTGACTGGa < 1:12421375/35‑1 (MQ=255) gCATCATTATGGCTTCTGGCCTTGGTGACTGGAtt < 1:860316/35‑1 (MQ=255) gCATCATTATGGCTTCTGGCCTTGGTGACTGGAtt < 1:17280627/35‑1 (MQ=255) catcatTATGGCTTCTGGCCTTGGTGACTGGATTg < 1:20195101/35‑1 (MQ=255) catcatTATGGCTTCTGGCCTTGGTGACTGGATTg < 1:14435560/35‑1 (MQ=255) atcatTATGGCTTCTGGCCTTGGTGACTGGATTGc < 1:11211089/35‑1 (MQ=255) tcatTATGGCTTCTGGCCTTGGTGACTGGATTGcc < 1:3244670/35‑1 (MQ=255) tcatTATGGCTTCTGGCCTTGGTGACTGGATTGcc < 1:13476664/35‑1 (MQ=255) catTATGGCTTCTGGCCTTGGTGACTGGATTGccc < 1:12361639/35‑1 (MQ=255) catTATGGCTTCTGGCCTTGGTGACTGGATTGccc < 1:2279130/35‑1 (MQ=255) catTATGGCTTCTGGCCTTGGTGACTGGATTGccc < 1:14815886/35‑1 (MQ=255) ttATGGCTTCTGGCCTTGGTGACTGGATTGCCCAc < 1:18482040/35‑1 (MQ=255) ttATGGCTTCTGGCCTTGGTGACTGGATTGCCCAc < 1:20538499/35‑1 (MQ=255) ttATGGCTTCTGGCCTTGGTGACTGGATTGACCAc < 1:3726412/35‑1 (MQ=255) tATGGCTTCTGGCCTTGGTGACTGGATTGCCCACg < 1:12493144/35‑1 (MQ=255) tATGGCTTCTGGCCTTGGTGACTGGATTGCCCACg < 1:17538316/35‑1 (MQ=255) tATGGCTTCTGGCCTTGGTGACTGGATTGCCCACg < 1:7833842/35‑1 (MQ=255) aTGGCTTCTGGCCTTGGTGACTGGATTGCCCACgg < 1:17859702/35‑1 (MQ=255) aTGGCTTCTGGCCTTGGTGACTGGATTGCCCACgg < 1:2478219/35‑1 (MQ=255) tGGCTTCTGGCCTTGGTGACTGGATTGCCCACGGc < 1:18074797/35‑1 (MQ=255) tGGCTTCTGGCCTTGGTGACTGGATTGCCCACGGc < 1:11427101/35‑1 (MQ=255) tGGCTTCTGGCCTTGGTGACTGGATTGCCCACGGc < 1:3048668/35‑1 (MQ=255) gCTTCTGGCCTTGGTGACTGGATTGCCCACGGCCt > 1:17800035/1‑35 (MQ=255) gCTTCTGGCCTTGGTGACTGGATTGCCCACGGCCt < 1:11715720/35‑1 (MQ=255) gCTTCTGGCCTTGGTGACTGGATTGCCCACGGCCt < 1:4459313/35‑1 (MQ=255) gCTTCTGGCCTTGGTGACTGGATTGCCCACGGCCt < 1:5876735/35‑1 (MQ=255) gCTTCTGGCCTTGGTGACTGGATTGCCCACGGCCt > 1:13983359/1‑35 (MQ=255) gCTTCTGGCCTTGGTGACTGGATTGCCAACGGCCt < 1:11932991/35‑1 (MQ=255) cTTCTGGCCTTGGTGACTGGATTGCCCACGGCCtt < 1:3487921/35‑1 (MQ=255) ttCTGGCCTTGGTGACTGGATTGCCCACGGCCTTg > 1:16200483/1‑35 (MQ=255) tCTGGCCTTGGTGACTGGATTGCCCACGGCCTTGc < 1:15114823/35‑1 (MQ=255) tCTGGCCTTGGTGACTGGATTGCCCACGGCCTTGc > 1:1544898/1‑35 (MQ=255) tCTGGCCTTGGTGACTGGATTGCCCACGGCCTTGc < 1:5646340/35‑1 (MQ=255) tCTGGCCTTGGTGACTGGATTGCCCACGGCATTGc < 1:11581785/35‑1 (MQ=255) cTGGCCTTTGTGACTGGATTGCCCACGGCCTTGCt > 1:1421663/1‑35 (MQ=255) cTGGCCTTGGTGACTGGATTGCCCACGGCCTTGCt < 1:7591368/35‑1 (MQ=255) cTGGCCTTGGTGACTGGATTGCCCACGGCCTTGCt > 1:12070199/1‑35 (MQ=255) cTGGCCTTGGTGACTGGATTGCCCACGGCCTTGCt < 1:734532/35‑1 (MQ=255) cTGGCCTTGGTGACTGGATTGCCCACGGCCTTGCt < 1:17341292/35‑1 (MQ=255) cTGGCCTTGGTGACGGGATTGCCCACGGCCTTGCt < 1:2817052/35‑1 (MQ=255) | TTCGTCTCGGCGCTCATTGGCATCATTATGGCTTCTGGCCTTGGTGACTGGATTGCCCACGGCCTTGCT > REL606/2721797‑2721865 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ 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. |