OPENSEQ.org
ss

Genes: A B A+B
Length: 317 153 459
Sequences: 4732 2938 120
Seq/Len: 14.93 19.2 0.26
MirrorTree (Pazo et al. 2001) 0.35
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.03 0.35 0.09
2 0.04 0.39 0.10
5 0.05 0.42 0.14
10 0.05 0.43 0.19
20 0.05 0.44 0.26
100 0.06 0.46 0.40
0.07 0.48 0.61
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
269_A 90_K 1.60 0.44 0.00
258_Y 136_L 1.54 0.41 0.00
285_D 17_S 1.40 0.32 0.00
260_G 136_L 1.34 0.29 0.00
72_N 42_I 1.20 0.22 0.00
15_Y 141_D 1.18 0.21 0.00
191_Q 86_L 1.17 0.21 0.00
96_D 93_L 1.15 0.20 0.00
273_K 72_Y 1.14 0.19 0.00
230_V 42_I 1.12 0.19 0.00
79_T 119_F 1.11 0.18 0.00
16_Q 139_L 1.10 0.18 0.00
35_I 84_D 1.08 0.17 0.00
181_V 82_A 1.08 0.17 0.00
26_P 93_L 1.06 0.16 0.00
85_L 72_Y 1.05 0.16 0.00
18_L 88_Q 1.04 0.16 0.00
207_D 117_D 1.02 0.15 0.00
258_Y 70_G 1.01 0.15 0.00
224_K 6_I 1.00 0.14 0.00
285_D 51_Q 0.99 0.14 0.00
230_V 93_L 0.98 0.14 0.00
15_Y 22_A 0.98 0.14 0.00
270_G 119_F 0.97 0.14 0.00
230_V 129_L 0.97 0.14 0.00
181_V 147_S 0.97 0.13 0.00
226_Y 13_T 0.97 0.13 0.00
22_V 42_I 0.96 0.13 0.00
37_V 5_L 0.96 0.13 0.00
63_Y 139_L 0.96 0.13 0.00
96_D 145_V 0.96 0.13 0.00
205_S 77_E 0.95 0.13 0.00
37_V 127_T 0.95 0.13 0.00
148_V 62_Q 0.94 0.13 0.00
237_G 125_I 0.94 0.13 0.00
18_L 38_G 0.94 0.12 0.00
164_Q 23_L 0.92 0.12 0.00
236_L 117_D 0.92 0.12 0.00
64_A 114_T 0.92 0.12 0.00
69_G 146_I 0.92 0.12 0.00
41_T 6_I 0.92 0.12 0.00
162_V 146_I 0.91 0.12 0.00
85_L 73_E 0.91 0.12 0.00
273_K 119_F 0.90 0.12 0.00
6_V 145_V 0.90 0.12 0.00
75_D 39_Q 0.90 0.11 0.00
186_R 48_S 0.90 0.11 0.00
109_Q 46_Q 0.90 0.11 0.00
67_I 44_I 0.90 0.11 0.00
40_G 75_C 0.90 0.11 0.00
258_Y 69_K 0.90 0.11 0.00
184_Q 140_C 0.89 0.11 0.00
231_V 50_K 0.89 0.11 0.00
99_I 141_D 0.89 0.11 0.00
234_S 69_K 0.89 0.11 0.00
234_S 58_S 0.88 0.11 0.00
240_I 139_L 0.88 0.11 0.00
298_I 116_E 0.88 0.11 0.00
55_L 34_L 0.88 0.11 0.00
298_I 147_S 0.88 0.11 0.00
48_V 86_L 0.87 0.11 0.00
196_S 59_L 0.87 0.11 0.00
236_L 87_G 0.87 0.11 0.00
152_V 111_G 0.87 0.11 0.00
185_G 23_L 0.86 0.10 0.00
161_V 122_P 0.86 0.10 0.00
202_I 127_T 0.86 0.10 0.00
295_A 56_I 0.86 0.10 0.00
55_L 42_I 0.86 0.10 0.00
77_S 84_D 0.86 0.10 0.00
95_Q 48_S 0.86 0.10 0.00
233_E 108_A 0.86 0.10 0.00
267_P 146_I 0.85 0.10 0.00
228_R 74_T 0.84 0.10 0.00
259_Y 81_D 0.84 0.10 0.00
122_I 16_P 0.84 0.10 0.00
181_V 70_G 0.84 0.10 0.00
35_I 93_L 0.83 0.10 0.00
193_T 108_A 0.83 0.10 0.00
214_E 119_F 0.83 0.10 0.00
24_A 94_T 0.83 0.10 0.00
118_D 117_D 0.83 0.10 0.00
67_I 92_F 0.83 0.10 0.00
289_A 114_T 0.83 0.10 0.00
181_V 69_K 0.83 0.10 0.00
18_L 125_I 0.82 0.10 0.00
228_R 118_S 0.82 0.10 0.00
152_V 35_K 0.82 0.10 0.00
163_F 47_A 0.82 0.10 0.00
304_L 48_S 0.82 0.10 0.00
72_N 58_S 0.82 0.10 0.00
235_Y 17_S 0.82 0.10 0.00
201_N 145_V 0.81 0.09 0.00
96_D 25_S 0.81 0.09 0.00
20_E 111_G 0.81 0.09 0.00
251_D 59_L 0.81 0.09 0.00
33_Y 54_K 0.81 0.09 0.00
196_S 82_A 0.81 0.09 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)

WARNING: The input alignment may be corrupted!
  • For sequence B, there is a high ratio (0.44 > 0.4) of paralogs.

ID Seq/Len Name Options I_Prob Status
8229 0.41 pengfei Δgene:(1, ∞) A:(1E-20, 8) B:(1E-20, 8) msa: HHblits (2015_06) 0.00 Done - Shared
8081 0.24 pengfei Δgene:(1, 20) A:(1E-06, 8) B:(1E-06, 8) msa: Jackhmmer (r132) 0.00 Done
8066 0.21 tt Δgene:(1, 20) A:(1E-20, 8) B:(1E-20, 8) msa: HHblits (2015_06) 0.00 Done
8065 0.27 tt Δgene:(1, 20) A:(1E-02, 8) B:(1E-02, 8) msa: Jackhmmer (r132) 0.00 Done - Shared
8064 0.22 test Δgene:(1, 20) A:(1E-02, 8) B:(1E-02, 8) msa: HHblits (2015_06) 0.00 Done
8063 0.26 ss Δgene:(1, 20) A:(1E-04, 8) B:(1E-04, 8) msa: Jackhmmer (r132) 0.00 Done - Shared

Page generated in 0.0388 seconds.