Read alignment evidence... | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|
seq id | position | change | freq | score | reads | annotation | genes | product | ||
* | REL606 | 1,111,095 | 0 | C→G | 100.0% | 8.4 | 17 | G94A (GGC→GCC) | insF‑1 | putative transposase for IS3 |
Rejected: E-value exceeds prediction threshold. | ||||||||||
Reads supporting (aligned to +/- strand): new base (9/0): ref base (7/0): total (17/0) | ||||||||||
Fisher's exact test strand distribution p-value = 1.00e+00 | ||||||||||
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.80e-02 |
CTGCTCCAACAGATTTTCTGACACAGGCAGGCCGTGTGCGCGGTAG > REL606/1111064‑1111109 | cTGCTCCAACATATTTTCTGACACAGGCAGGGCgg > 1:3758104/1‑34 (MQ=1) cTGCTCCAACAGATTTTCTGACACAGGCAGGGCgg > 1:3272260/1‑34 (MQ=0) cTGCTCCAACAGATTTTCTGACACAGGCAGGGCgg > 1:2759660/1‑34 (MQ=0) cTGCTCCAACAGATTTTCTGACACAGGCAGGGCgg > 1:15628429/1‑34 (MQ=0) cTGCTCCAACAGATTTTCTGACACAGGCAGGGCgg > 1:6226455/1‑34 (MQ=0) cTGCTCCAACAGATTTTCTGACACAGGCAGGGCgg > 1:15216934/1‑34 (MQ=0) cTGCTCCAACAGATTTTCTGACACAGGCAGGACgg > 1:7660244/1‑34 (MQ=0) cTGCTCCAACAGATGTTCTGACACAGGCAGGGCgg > 1:1746972/1‑34 (MQ=1) tGCTCCAACAGATTTTCTGACACAGGCAGGGCggg > 1:16962353/1‑33 (MQ=1) tGCTCCAACAGATTTTCTGACACAGGCAGGGCggg > 1:3559215/1‑33 (MQ=1) tGCTCCAACAGATTTTCGGACACAGGCAGGGCggg > 1:2862191/1‑33 (MQ=1) gCTCCAACAGATTTTCTGACACAGGCAGGGCgggg > 1:14867958/1‑32 (MQ=1) gCTCCAACAGATTTTCTGACACAGGCAGGGCgggg > 1:13680243/1‑32 (MQ=1) gCTCCAACAGATTTTCTGACACAGGCAGGGCgggg > 1:15733221/1‑32 (MQ=1) gCTCCAACAGATTTTCTGACACAGGCAGGGCgggg > 1:231342/1‑32 (MQ=1) gCTCCAACAGATTTTCTGACACAGGCAGGGCgggg > 1:12393622/1‑32 (MQ=1) cTCCAACAGATTTTCTGACACAGGCAGGGCgggta > 1:13067323/1‑34 (MQ=1) ccAACAGATTTTCTGACACAGGCAGGCCgggtttg > 1:7478450/1‑32 (MQ=1) ccAACAGATTTTCTGACACAGGCAGGCCGGGGgcg > 1:5149351/1‑35 (MQ=1) ccAACAGATTTTCTGACACAGGCAGGCCGGGGgcg > 1:14935944/1‑35 (MQ=1) ccAACAGATTTTCTGACACAGGCAGGCCGGGGgcg > 1:4661932/1‑35 (MQ=1) aCAGATTTTCTGACACAGGCAGGGCGGGGGGGCgg > 1:2483216/1‑35 (MQ=0) aCAGATTTTCTGACACAGGCAGGGCGGGGGGGCgg > 1:11222377/1‑35 (MQ=0) aCAGATTTTCTGACACAGGCAGGCCGGGGGGGCgg > 1:5833624/1‑35 (MQ=1) aCAGATTTTCTGACACAGGCAGGCCGGGGGGGCgg > 1:9143553/1‑35 (MQ=1) cAGATTTTCTGACACAGGCAGGGCGGGGGGGCggg > 1:602470/1‑34 (MQ=0) cAGATTTTCTGACACAGGCAGGCCGCGGGGGCggg > 1:13327073/1‑34 (MQ=1) aGATTTTCTGACACAGGCAGGCCGGGTGCGGgggg > 1:13730080/1‑33 (MQ=1) gATTTTCTGACACAGGCAGGGCGGGGGCGGggggg > 1:18010159/1‑32 (MQ=0) gATTTTCTGACACAGGCAGGCCGTGGGCGGgggaa > 1:1273468/1‑32 (MQ=1) | CTGCTCCAACAGATTTTCTGACACAGGCAGGCCGTGTGCGCGGTAG > REL606/1111064‑1111109 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 28 ≤ 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. |