May 4, 2021 - We are working on upgrading the webserver, some pages may not work.
OPENSEQ.org

10vs49_-04

Genes: A B A+B
Length: 86 424 462
Sequences: 196 1369 68
Seq/Len: 2.28 3.23 0.15
MirrorTree (Pazo et al. 2001) 0.24
Download Alignment
We filter this alignment to remove positions that have > 75% gaps before running GREMLIN.

[Show Jackhmmer results] - Maybe useful in deciding what regions to trim.

Δgene Ratio of paralogs Joined
A B Seq/Len
1 0.00 0.04 0.01
2 0.01 0.04 0.01
5 0.01 0.05 0.01
10 0.01 0.05 0.01
20 0.01 0.05 0.01
100 0.01 0.07 0.04
0.01 0.10 0.13
Paired alignment generation
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
Sequence A E-value:
Iterations:
Sequence B E-value:
Iterations:
Δgene MSA
Figure 1: The effect of Δgene on ratio of paralogs in the genome for each gene, and the number of sequences in the resulting joined alignment. Figure 2: Distribution of Δgene across all genomes. Figure 3: Current parameters. You can use the form to adjust the parameters and resubmit the job!

GREMLIN results (Scaled_score > 1):
Figure 4: The darker and larger the blue dots, the higher strength in coevolution.
Inter Residue pairs sorted by strength in coevolution signal:
i j Scaled Score Prob I_Prob
17_K 165_Q 1.27 0.18 0.00
41_C 154_A 1.23 0.16 0.00
50_L 33_F 1.15 0.14 0.00
37_R 124_G 1.14 0.14 0.00
31_W 124_G 1.12 0.13 0.00
45_G 124_G 1.11 0.13 0.00
44_M 295_P 1.10 0.13 0.00
84_E 336_R 1.10 0.13 0.00
44_M 288_H 1.09 0.13 0.00
64_V 265_P 1.07 0.12 0.00
82_K 262_P 1.05 0.12 0.00
58_N 320_G 1.04 0.11 0.00
12_E 158_A 1.03 0.11 0.00
35_Q 353_P 1.03 0.11 0.00
58_N 126_I 0.99 0.10 0.00
69_M 61_F 0.98 0.10 0.00
72_M 120_F 0.97 0.10 0.00
40_Y 279_G 0.97 0.10 0.00
58_N 165_Q 0.97 0.10 0.00
55_I 308_Q 0.97 0.10 0.00
27_T 328_G 0.96 0.10 0.00
1_M 153_D 0.96 0.10 0.00
39_S 272_G 0.96 0.10 0.00
44_M 230_S 0.94 0.09 0.00
58_N 313_H 0.94 0.09 0.00
19_I 127_L 0.94 0.09 0.00
37_R 158_A 0.94 0.09 0.00
80_A 239_P 0.93 0.09 0.00
18_Y 345_G 0.93 0.09 0.00
82_K 123_F 0.93 0.09 0.00
82_K 30_W 0.91 0.09 0.00
49_L 321_H 0.91 0.09 0.00
31_W 116_L 0.90 0.09 0.00
68_L 48_R 0.90 0.09 0.00
84_E 266_M 0.89 0.08 0.00
77_G 120_F 0.88 0.08 0.00
55_I 145_A 0.87 0.08 0.00
54_A 319_L 0.87 0.08 0.00
76_C 331_G 0.87 0.08 0.00
51_N 268_P 0.86 0.08 0.00
84_E 169_N 0.86 0.08 0.00
39_S 33_F 0.86 0.08 0.00
49_L 244_P 0.85 0.08 0.00
33_V 272_G 0.85 0.08 0.00
39_S 154_A 0.85 0.08 0.00
11_L 146_F 0.85 0.08 0.00
50_L 37_G 0.85 0.08 0.00
24_A 120_F 0.85 0.08 0.00
41_C 27_P 0.84 0.08 0.00
2_T 29_L 0.84 0.08 0.00
58_N 337_P 0.84 0.08 0.00
63_R 131_K 0.84 0.08 0.00
65_R 57_G 0.84 0.08 0.00
62_A 306_M 0.84 0.07 0.00
81_D 298_P 0.83 0.07 0.00
12_E 120_F 0.83 0.07 0.00
70_E 265_P 0.83 0.07 0.00
80_A 102_I 0.83 0.07 0.00
24_A 123_F 0.83 0.07 0.00
20_G 149_F 0.83 0.07 0.00
38_D 149_F 0.83 0.07 0.00
46_H 149_F 0.83 0.07 0.00
55_I 344_P 0.82 0.07 0.00
52_Y 247_V 0.82 0.07 0.00
3_D 99_G 0.81 0.07 0.00
Figure 5: The i (protein A) and j (protein B) are positions as given in the query sequences.

Scaled Score = raw_score/average(raw_scores)
Prob = P(contact | scaled_score, seq/len)
I_Prob = P(contact | scaled_score, seq/len, top_inter_score)

ID Seq/Len Name Options I_Prob Status
3911 0.16 10vs49-10 Δgene:(1, ∞) A:(1E-10, 8) B:(1E-10, 8) msa: Jackhmmer (2015_06) 0.00 Done
3909 0.15 10vs49_-04 Δgene:(1, ∞) A:(1E-04, 8) B:(1E-04, 8) msa: Jackhmmer (2015_06) 0.00 Done
3908 0.14 10vs49_-02 Δgene:(1, ∞) A:(1E-02, 8) B:(1E-02, 8) msa: Jackhmmer (2015_06) 0.00 Done - Shared
3883 0.16 10vs49 Δgene:(1, ∞) A:(1E-20, 8) B:(1E-20, 8) msa: Jackhmmer (2015_06) 0.00 Done

Page generated in 1.7664 seconds.