Read alignment evidence... | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|
seq id | position | change | freq | score | reads | annotation | genes | product | ||
* | REL606 | 864,069 | 0 | T→G | 100.0% | 12.7 | 17 | V184G (GTG→GGG) | yliC | predicted peptide transporter subunit: membrane component of ABC superfamily |
Rejected: E-value exceeds prediction threshold. | ||||||||||
Reads supporting (aligned to +/- strand): new base (11/0): ref base (1/5): total (12/5) | ||||||||||
Fisher's exact test strand distribution p-value = 9.70e-04 | ||||||||||
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00 |
TACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCG > REL606/864035‑864103 | tACATTTTACCCTCCCTGACGCTCGGCGCGGCAGt < 1:16144299/35‑1 (MQ=255) aCATTTTACCCTCCCTGACGCTCGGCGCGGCAGTg < 1:758100/35‑1 (MQ=255) aCATTTTACCCTCCCTGACGCTCGGCGCGGCAGTg < 1:9461080/35‑1 (MQ=255) aCATTTTACCCTCCCTGACGCTCGGCGCGGCAGTg < 1:18579117/35‑1 (MQ=255) cATTTTACCCTCCCCGACGCTCGGCGCGGCAGGgg > 1:157771/1‑35 (MQ=37) aTTTTACCCTCCCTGACGCTCGGCGCGGCAGGggg > 1:5404701/1‑34 (MQ=255) aTTTTACCCTCCCTGACGCTCGGCGCGGCAGGggg > 1:9884799/1‑34 (MQ=255) ttttACCCTCCCTGACGCTCGGCGCGGCAGGgggc > 1:6394385/1‑33 (MQ=37) ttttACCCTCCCTGACGCTCGGCGCGGCAGGgggc > 1:16898235/1‑33 (MQ=37) ttttACCCTCCCTGACGCTCGGCGCGGCAGGGgcg > 1:10952454/1‑34 (MQ=255) ttttACCATCCCTGACGCTCGGCGCGGCAGTGgcc < 1:1276403/35‑1 (MQ=255) tttCCCCCCCCTGACGCTCGGCGCGGCAGTGgccg < 1:709797/35‑1 (MQ=37) tttACCCTCCCTGACGCTCGGCGCGGCAGGGggcg > 1:7326549/1‑35 (MQ=37) tttACCCTCCCTGACGCTCGGCGCGGCAGGGggcg > 1:5122166/1‑35 (MQ=37) tttACCCTCCCTGACGCTCGGCGCGGCAGGGggcg > 1:4219233/1‑35 (MQ=37) ttACCCTCCCTGACGCTCGGCGCGGCAGGGggcgg > 1:10896044/1‑34 (MQ=25) tACCCTCCCTGACGCTCGGCGCGGCAGGGggcgcg > 1:18781370/1‑34 (MQ=25) aCCCTCCCTGACGCTCGGCGCGGCAGGGggcgggc > 1:15732436/1‑32 (MQ=25) cctccctGACGCTCGGCGCGGCAGGGggcgcccgg > 1:6130742/1‑32 (MQ=25) cctccctGACGCTCGGCGCGGCAGGGGgcgccggg > 1:3025288/1‑33 (MQ=25) ctccctGACGCTCGGCGCTGCAGTGGccggcgggg > 1:1067847/1‑32 (MQ=21) cggggCGCGGCGGGGCCCGCCGTGATGGCGCGCtt < 1:16048931/32‑1 (MQ=19) cgcgGCGGGGCCCGCCGTGATGGCGCGCTTTAccc < 1:11999410/35‑1 (MQ=21) gggCAGTGCCCGCCGTGATGGCGCGCTTTACCcgc < 1:997639/34‑1 (MQ=255) ggCAGCGCCCGCCGTGATGGCGCGCTTTACCcgcg < 1:2220750/35‑1 (MQ=25) aGTGCCCGCCGTGATGGCGCGCTTTACCCGCGcgt < 1:7234584/35‑1 (MQ=255) tGGCCGCCGTGATGGCGCGCTTTACCCGCGcgtcg > 1:10129760/1‑35 (MQ=255) | TACATTTTACCCTCCCTGACGCTCGGCGCGGCAGTGGCCGCCGTGATGGCGCGCTTTACCCGCGCGTCG > REL606/864035‑864103 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ 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. |