Read alignment evidence... | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|
seq id | position | change | freq | score | reads | annotation | genes | product | ||
* | REL606 | 4,557,629 | 0 | T→G | 100.0% | 21.9 | 27 | T53P (ACC→CCC) | yjiW | hypothetical protein |
Rejected: E-value exceeds prediction threshold. | ||||||||||
Reads supporting (aligned to +/- strand): new base (17/0): ref base (1/9): total (18/9) | ||||||||||
Fisher's exact test strand distribution p-value = 2.13e-06 | ||||||||||
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00 |
ATCACTTTGACATCTACCGCCGTGCCGGTGGTAAAACCGGCGGCTTCCAGCCACT > REL606/4557598‑4557652 | aTCACTTTGACATCTACCGCCGTGCCGGTGGGaaa > 1:18659954/1‑35 (MQ=255) aTCACTTTGACATCTACCGCCGTGCCGGTGGGaaa > 1:6093991/1‑35 (MQ=255) aTCACTTTGACATCTACCGCCGTGCCGGTGGGaaa > 1:19163292/1‑35 (MQ=255) tCACTTTGACATCTACCGCCGTGCCGGTGGGaaaa > 1:5136784/1‑35 (MQ=255) tCACTTTGACATCTACCGCCGTGCCGGTGGGaaaa > 1:12839008/1‑35 (MQ=255) tCACTTTGACATCTACCGCCGTGCCGGTGGGaaaa > 1:13827665/1‑35 (MQ=255) cACTTTGACATCTACCGCCGTGCCGGTGGTAAAAc > 1:12214194/1‑35 (MQ=255) cACTTTGACATCTACCGCCGTGCCGGTGGGaaaaa > 1:2015973/1‑34 (MQ=255) cACTTTGACATCTACCGCCGTGCCGGTGGGAAAAc > 1:11711235/1‑35 (MQ=255) cACTTTGACATCTACCGCCGTGCCGGTGGGAAAAc > 1:20237183/1‑35 (MQ=255) aCTTTGACATCTACCGCCGTGCCGGTGGGaaaaac > 1:17265330/1‑33 (MQ=37) aCTTTGACATCTACCGCCGTGCCGGTGGGaaaaac > 1:15264169/1‑33 (MQ=37) aCTTTGACATCTACCGCCGTGCCGGTGGGaaaaac > 1:19591200/1‑33 (MQ=37) aCTTTGACATCTACCGCCGTGCCGGTGGGAAAAcc > 1:13347954/1‑35 (MQ=255) tttGACATCTACCGCCGTGCCGGTGGGAAAACcgg > 1:16077790/1‑35 (MQ=255) ttGACATCTACCGCCGTGCCGGTGGGAAAAcccgg > 1:10879095/1‑32 (MQ=25) ttGACATCTACCGCCGTGCCGGTGGGAAAACcggc > 1:1402512/1‑35 (MQ=255) ttGACATCTACCGCCGTGCCGGTGGGAAAACcggc > 1:18694735/1‑35 (MQ=255) tGACATCTACCGCCGTGCCGGTGGGAAAACcgggg > 1:18670817/1‑33 (MQ=37) gACATCTACCGCCGTGCCGGTGGGGAAAAcggggg > 1:4657888/1‑35 (MQ=16) gACATCTACCGCCGTGCCGGTGGGAAAACccgcgg > 1:17393501/1‑35 (MQ=25) cATCTACCGCCGTGCCGGTGGGAAAACcggcgggg > 1:8332214/1‑33 (MQ=37) aTCTACCGCCGTGCCGGTGGGGAAAAccgcggggt > 1:19234587/1‑32 (MQ=12) cccccgccgTGCCGGTGGTAAAACCGGCGGCTTcc < 1:7134782/32‑1 (MQ=255) ccccgccgTGCCGGTGGTAAAACCGGCGGCTTCCa < 1:15484695/33‑1 (MQ=255) ccccgccgTGCCGGTGGTAAAACCGGCGGCTTCCa < 1:15108456/33‑1 (MQ=255) ccccgccgTGCCGGTGGTAAAACCGGCGGCTTCCa < 1:769560/33‑1 (MQ=255) ccccgccgTGCCGGTGATAAAACCGGCGGCTTCCa < 1:14430849/33‑1 (MQ=37) gccgccgTGCCGGGGGTAAAACCGGCGGCTTCCAg < 1:9889178/34‑1 (MQ=255) cccgccgTGCCGGTGGTAAAACCGGCGGCTTCCAg < 1:5457324/34‑1 (MQ=255) ccgccgTGCCGGTGGTAAAACCGGCGGCTTCCAGc < 1:11292922/35‑1 (MQ=255) cgccgTGCCGGTGGTAAACCCGGCGGCTTCCAGcc < 1:12170104/35‑1 (MQ=255) cgTGCCGGTGTTAAAACCGGCGGGCTCCCGCCCCt > 1:15309437/1‑35 (MQ=12) | ATCACTTTGACATCTACCGCCGTGCCGGTGGTAAAACCGGCGGCTTCCAGCCACT > REL606/4557598‑4557652 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 28 ≤ 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. |