Read alignment evidence... | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|
seq id | position | change | freq | score | reads | annotation | genes | product | ||
* | REL606 | 2,766,958 | 1 | .→C | 1.5% | 5.8 0.0 ‑0.2 | 194 | coding (394/711 nt) | ispD | 2‑C‑methyl‑D‑erythritol 4‑phosphate cytidylyltransferase |
Rejected: Polymorphism prediction not supported by minimum number of reads on both strands. | ||||||||||
Rejected: Polymorphic indel expands or contracts a homopolymer stretch. | ||||||||||
Reads supporting (aligned to +/- strand): new base (2/1): ref base (106/85): total (108/86) | ||||||||||
Fisher's exact test strand distribution p-value = 1.00e+00 | ||||||||||
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.63e-01 |
Only 100 of 239 total aligned reads displayed.
CGTTTCATAGTATCGCGCACTGGTGCGGCGAGGAT‑CCCCCCCGTGCGGCTGGTTTCGCTCAACGCCAAC > REL606/2766924‑2766992 | cGTTTCATAGTATCGCGCACTGGTGCGGCGAGGAt < 1:16286370/35‑1 (MQ=255) gTTTCATAGTATCGCGCACTGGTGCGGCGAGGAT‑c < 1:18062561/35‑1 (MQ=255) gTTTCATAGTATCGCGCACTGGTGCGGCGAGGAT‑c < 1:1579840/35‑1 (MQ=255) gTTTCATAGTATAGCGCACTGGTGCGGCGAGGAT‑c < 1:17687320/35‑1 (MQ=255) tttCATAGTATCGCGCACTGGTGCGGCGAGGAT‑cc < 1:6370505/35‑1 (MQ=255) cATATTATCGCGCACTGGTGCGGCGAGGAT‑ccccc < 1:7781124/35‑1 (MQ=255) cATAGTATCGCGCACTGGTGCGGCGAGGAA‑ccccc > 1:3004864/1‑35 (MQ=255) aTAGTATCGCGCACTGGTGCGGCGAGGAT‑cccccc < 1:19218204/35‑1 (MQ=255) aTAGTATCGCGCACTGGTGCGGCGAGGAT‑cccccc < 1:18132207/35‑1 (MQ=255) aTAGTATCGCGCACTGGTGCGGCGAGGAT‑cccccc < 1:15415644/35‑1 (MQ=255) aTAGTATCGCGCACTGGTGCGGCGAGGAT‑cccccc < 1:11987926/35‑1 (MQ=255) tAGTATCGCGCACTGGTGCGGCGAGGAT‑ccccccc < 1:19174086/35‑1 (MQ=255) tAGTATCGCGCACTGGTGCGGCGAGGAA‑ccccccc > 1:2385676/1‑35 (MQ=255) tAGGATCGCGCACTGGTGCGGCGAGGAT‑ccccccc < 1:316813/35‑1 (MQ=255) aGTATCGCGCACTGGTGCGGCGAGGAT‑CCCCCCCg > 1:20356626/1‑35 (MQ=255) aGTATCGCGCACTGGTGCGGCGAGGAT‑CCCCCCCg < 1:15180614/35‑1 (MQ=255) aGTATCGCGCACTGGTGCGGCGAGGAT‑CCCCCCCg > 1:7625682/1‑35 (MQ=255) aGTATCGCGCACTGGTGCGGCGAGGAT‑CCCCCCCg > 1:113724/1‑35 (MQ=255) aGTATCGCGCACTGGTGCGGCGAGGAT‑CCCCCCCg > 1:10585397/1‑35 (MQ=255) gTATCGCGCACTGGTGCGGCGAGGAT‑CCCCCCCGt > 1:7211183/1‑35 (MQ=255) tATCGCGCACTGGTGCGGCGAGGAT‑CCCCCCCGTg < 1:16472673/35‑1 (MQ=255) aTCGCGCACTGGTGCGGCGAGGAT‑CCCCCCCGTGc > 1:3773285/1‑35 (MQ=255) aTCGCGCACTGGTGCGGCGAGGAT‑CCCCCCCGTGc > 1:2745208/1‑35 (MQ=255) aTCGCGCACTGGTGCGGCGAGGAT‑CCCCCCCGTGc < 1:9628448/35‑1 (MQ=255) tCGCGCACTGGTGCGGCGAGGAT‑CCCCCCCGTGCg < 1:17099643/35‑1 (MQ=255) tCGCGCACTGGTGCGGCGAGGAT‑CCCCCCCGTGCg < 1:17269693/35‑1 (MQ=255) tCGCGCACTGGTGCGGCGAGGAT‑CCCCCCCGTGCg < 1:6804407/35‑1 (MQ=255) cgcgCACTGGTGCGGCGAGGAT‑CCCCCCCGTGCgg > 1:8919144/1‑35 (MQ=255) gcgcACTGGTGCGGCGAGGAT‑CCCCCCCGTGCGGc < 1:12425086/35‑1 (MQ=255) gcgcACTGGTGCGGCGAGGAT‑CCCCCCCGTGCGGc < 1:10571102/35‑1 (MQ=255) cgcACTGGTGCGGCGAGGAT‑CCCCCCCGTGCGGCt < 1:5207385/35‑1 (MQ=255) gcACTGGTGCGGCGAGGAT‑CCCCCCCGTGCGGCTg < 1:7671825/35‑1 (MQ=255) gcACTGGTGCGGCGAGGAT‑CCCCCCCGTGCGGCTg < 1:15940727/35‑1 (MQ=255) gcACTGGTGCGGCGAGGAT‑CCCCCCCGTGCGGCTg > 1:10651071/1‑35 (MQ=255) gcACTGGTGCGGCGAGGAT‑CCCCCCCGTGCGGCTg > 1:1861345/1‑35 (MQ=255) cACTGGTGCGGCGAGGAT‑CCCCCCCGTGCTGCTgg > 1:6721158/1‑35 (MQ=255) cACTGGTGCGGCGAGGAT‑CCCCCCCGTGCGGCtcg > 1:12309459/1‑33 (MQ=255) cACTGGTGCGGCGAGGAT‑CCCCCCCGTGCGGCTgg > 1:6229571/1‑35 (MQ=255) cACTGGTGCGGCGAGGAT‑CCCCCCCGTGCGGCTgg > 1:8696497/1‑35 (MQ=255) cACTGGTGCGGCGAGGAT‑CCCCCCCGTGCGGCTgg < 1:9979467/35‑1 (MQ=255) cACTGGTGCGGCGAGGAT‑CCCCCCCGCGCGGCTgg > 1:20281290/1‑35 (MQ=255) cACTGGTGCGGCGACGAT‑CCCCCCCGTGCGGCTgg < 1:4786713/35‑1 (MQ=255) aCTGGTGCGGCGAGGAT‑CCCCCCCGTGCGGCTggg > 1:16739332/1‑34 (MQ=255) aCTGGTGCGGCGAGGAT‑CCCCCCCGTGCGGCTggg > 1:8095031/1‑34 (MQ=255) aCTGGTGCGGCGAGGAT‑CCCCCCCGTGCGGCTggg > 1:8228553/1‑34 (MQ=255) aCTGGTGCGGCGAGGAT‑CCCCCCCGTGCGGCTggg > 1:9527488/1‑34 (MQ=255) aCTGGTGCGGCGAGGAT‑CCCCCCCGTGCGGCTGGt > 1:1985857/1‑35 (MQ=255) cTGGTGCGGCGAGGAT‑CCCCCCCGTGCGGCTgggt > 1:4324675/1‑33 (MQ=255) cTGGTGCGGCGAGGAT‑CCCCCCCGTGCGGCTgggt > 1:1397442/1‑33 (MQ=255) cTGGTGCGGCGAGGAT‑CCCCCCCGTGCGGCTGGtt < 1:10103735/35‑1 (MQ=255) cTGGTGCGGCGAGGAT‑CCCCCCCGTGCGGCTGGtt < 1:15520766/35‑1 (MQ=255) tGGTGCGGCGAGGAT‑CCCCCCCGTGCGGCTGGttt < 1:15908648/35‑1 (MQ=255) tGGTGCGGCGAGGAT‑CCCCCCCGTGCGGCTGGgtt > 1:4879704/1‑35 (MQ=255) ggTGCGGCGAGGAT‑CCCCCCCGTGCGGCTGGTTTc < 1:13013395/35‑1 (MQ=255) ggTGCGGCGAGGAT‑CCCCCCCGTGCGGCTGGTTTc < 1:4083037/35‑1 (MQ=255) ggTGCGGCGAGGAT‑CCCCCCCGTGCGGCTGGGTTc > 1:1324979/1‑35 (MQ=255) gTGCGGCGAGGAT‑CCCCCCCGTGCGGCTGGgtttc > 1:2098377/1‑33 (MQ=38) tGCGGCGAGGAT‑CCCCCCCGTGCGGCTGGgtttgg > 1:1183903/1‑32 (MQ=38) tGCGGCGAGGAT‑CCCCCCCGTGCGGCTGGgtttgg > 1:5370680/1‑32 (MQ=38) tGCGGCGAGGAT‑CCCCCCCGTGCGGCTGGTTTCGc < 1:5122185/35‑1 (MQ=255) tGCGGCGAGGAT‑CCCCCCCGTGCGGCTGGGTTCGc > 1:8549284/1‑35 (MQ=255) gcggcgAGGAT‑CCCCCCCGTGCGGCTGGTTTCGCt < 1:188814/35‑1 (MQ=255) gcggcgAGGAT‑CCCCCCCGTGCGGCTGGTTTCGCt < 1:9903321/35‑1 (MQ=255) cggcgAGGATCCCCCCCCGTGCGGCTGGGTTCggt > 1:336843/1‑33 (MQ=25) cggcgAGGAT‑CCCCCCCGTGCGGCTGGTTTCGCTc < 1:6096241/35‑1 (MQ=255) cggcgAGGAT‑CCCCCCCGTGCGGCTGGTTTCGCTc < 1:11593708/35‑1 (MQ=255) cggcgAGGAT‑CCCCCCCGTGCGGCTGGTTTCGCTc < 1:4606845/35‑1 (MQ=255) cggcgAGGAT‑CCCCCCCGTGCGGCTGGTTTCGCTc < 1:14795552/35‑1 (MQ=255) cggcgAGGAT‑CCCCCCCGTGCGGCTGGGTTTGCTc > 1:19769286/1‑35 (MQ=37) cggcgAGGAT‑CCCCCCCGTGCGGCTGGGTTCGCTc > 1:8433079/1‑35 (MQ=255) cggcgAGGAT‑CCCCCCCGTGCGGCTGGGTTCGCTc > 1:9093764/1‑35 (MQ=255) cggcgAGGAT‑CCCCCCCGTGCGGCTGGGTTCGCTc > 1:7269014/1‑35 (MQ=255) ggcgAGGAT‑CCCCCCCGTGCGGCTGGTTTCGCTCa > 1:15488984/1‑35 (MQ=255) ggcgAGGAT‑CCCCCCCGTGCGGCTGGGTTCGGTCa > 1:7037499/1‑35 (MQ=38) ggcgAGGAT‑CCCCCCCGTGCGGCTGGGTTCGCTCa > 1:388602/1‑35 (MQ=255) gcgAGGAT‑CCCCCCCGTGCGGCTGGTTTCGCTCaa < 1:4452564/35‑1 (MQ=255) gcgAGGAT‑CCCCCCCGTGCGGCTGGTTTCGCTCaa > 1:10742866/1‑35 (MQ=255) gcgAGGAT‑CCCCCCCGTGCGGCTGGGTTCGCTCaa > 1:5624010/1‑35 (MQ=255) cgAGGAT‑CCCCCCCGTGCGGCTGGTTTCGCTCCAc > 1:18791314/1‑35 (MQ=255) cgAGGAT‑CCCCCCCGTGCGGCTGGTTTCGCTCAAc > 1:18300255/1‑35 (MQ=255) cgAGGAT‑CCCCCCCGTGCGGCTGGGTTTGCTCAAc > 1:13551128/1‑35 (MQ=37) cgAGGAT‑CCCCCCAGTGCGGCTGGTTTCGCTCAAc < 1:1695423/35‑1 (MQ=255) gAGGAT‑CCCCCCCGTGCGGCTGGTTTCGCTCAACg < 1:11073628/35‑1 (MQ=255) gAGGAT‑CCCCCCCGTGCGGCTGGTTTCGCTCAACg < 1:17520929/35‑1 (MQ=255) gAGGAT‑CCCCCCCGTGCGGCTGGTTTCGCTCAACg < 1:1744264/35‑1 (MQ=255) gAGGAT‑CCCCCCCGTGCGGCTGGTTTCGCTCAACg < 1:17200707/35‑1 (MQ=255) gAGGAT‑CCCCCCCGTGCGGCTGGGTTTGCTCAACg > 1:15585026/1‑35 (MQ=37) aGGAT‑CCCCCCCGTGCGGCTGGGTTCGCTCAACGc > 1:14727754/1‑35 (MQ=255) aGGAT‑CCCCCCCGTGCGGCTGGGTTAGCGCAACgg > 1:12202032/1‑34 (MQ=25) gAT‑CCCCCCCGTGCGGCTGGTTTCGCTCAACGCca < 1:10358331/35‑1 (MQ=255) gAT‑CCCCCCCGTGCGGCTGGTTTCGCTCAACGCca < 1:19618225/35‑1 (MQ=255) gAT‑CCCCCCCGTGCGGCTGGTTTCGCTAAACGCca < 1:18435675/35‑1 (MQ=255) gAT‑CCCCCCCGTGCGGCTGGGTTTGCTCCACGCca > 1:14899546/1‑35 (MQ=25) aT‑CCCCCCCGTGCTGCTGGTTTCGCTCAACGCcaa > 1:16206181/1‑35 (MQ=255) aT‑CCCCCCCGTGCGGCTGGTTTCGCTCAACGCcaa < 1:13429087/35‑1 (MQ=255) aT‑CCCCCCCGTGCGGCTGGTTTCGCTCAACGCcaa > 1:9711125/1‑35 (MQ=255) t‑CCCCCCCGTGCGGCTGGTTTCGCTCAACGCcaac < 1:13634383/35‑1 (MQ=255) t‑CCCCCCCGTGCGGCTGGTTTCGCTCAACGCcaac < 1:14144363/35‑1 (MQ=255) t‑CCCCCCCGTGCGGCTGGTTTCGCTCAACGCcaac < 1:10898452/35‑1 (MQ=255) t‑CCCCCCCGTGCGGCTGGTTTCGCTCAACGCcaac < 1:10406808/35‑1 (MQ=255) | CGTTTCATAGTATCGCGCACTGGTGCGGCGAGGAT‑CCCCCCCGTGCGGCTGGTTTCGCTCAACGCCAAC > REL606/2766924‑2766992 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 30 ≤ 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. |