OPENSEQ.org
tt

Genes: A B A+B
Length: 317 153 459
Sequences: 4641 2959 122
Seq/Len: 14.64 19.34 0.27
MirrorTree (Pazo et al. 2001) 0.33
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.20
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.50 0.39 0.00
258_Y 136_L 1.40 0.33 0.00
273_K 72_Y 1.31 0.28 0.00
285_D 17_S 1.27 0.26 0.00
37_V 127_T 1.15 0.20 0.00
26_P 93_L 1.15 0.20 0.00
85_L 72_Y 1.14 0.20 0.00
230_V 42_I 1.12 0.19 0.00
16_Q 139_L 1.11 0.19 0.00
96_D 93_L 1.06 0.17 0.00
82_S 146_I 1.06 0.17 0.00
77_S 84_D 1.06 0.17 0.00
72_N 42_I 1.06 0.17 0.00
15_Y 141_D 1.06 0.17 0.00
168_L 95_S 1.04 0.16 0.00
181_V 82_A 1.04 0.16 0.00
224_K 6_I 1.03 0.16 0.00
164_Q 23_L 1.03 0.16 0.00
181_V 147_S 1.03 0.16 0.00
37_V 5_L 0.98 0.14 0.00
99_I 110_D 0.98 0.14 0.00
23_A 123_P 0.97 0.14 0.00
233_E 108_A 0.97 0.14 0.00
15_Y 22_A 0.97 0.14 0.00
15_Y 123_P 0.97 0.14 0.00
35_I 84_D 0.97 0.14 0.00
15_Y 44_I 0.96 0.14 0.00
48_V 45_A 0.96 0.14 0.00
118_D 89_A 0.96 0.14 0.00
165_K 5_L 0.96 0.13 0.00
205_S 77_E 0.96 0.13 0.00
40_G 75_C 0.96 0.13 0.00
152_V 111_G 0.95 0.13 0.00
79_T 119_F 0.95 0.13 0.00
82_S 73_E 0.94 0.13 0.00
280_Y 102_N 0.94 0.13 0.00
99_I 141_D 0.93 0.13 0.00
18_L 88_Q 0.92 0.12 0.00
33_Y 117_D 0.92 0.12 0.00
258_Y 119_F 0.92 0.12 0.00
285_D 21_Q 0.92 0.12 0.00
103_A 13_T 0.92 0.12 0.00
41_T 6_I 0.91 0.12 0.00
63_Y 139_L 0.91 0.12 0.00
236_L 117_D 0.91 0.12 0.00
55_L 34_L 0.91 0.12 0.00
269_A 25_S 0.91 0.12 0.00
260_G 136_L 0.91 0.12 0.00
157_G 10_C 0.91 0.12 0.00
157_G 19_C 0.91 0.12 0.00
22_V 59_L 0.90 0.12 0.00
64_A 114_T 0.90 0.12 0.00
157_G 75_C 0.90 0.12 0.00
15_Y 92_F 0.89 0.12 0.00
184_Q 140_C 0.89 0.12 0.00
170_A 6_I 0.89 0.12 0.00
69_G 146_I 0.89 0.12 0.00
228_R 118_S 0.89 0.11 0.00
285_D 51_Q 0.89 0.11 0.00
186_R 48_S 0.89 0.11 0.00
196_S 59_L 0.88 0.11 0.00
24_A 94_T 0.88 0.11 0.00
191_Q 86_L 0.88 0.11 0.00
170_A 105_A 0.88 0.11 0.00
237_G 125_I 0.88 0.11 0.00
193_T 108_A 0.88 0.11 0.00
157_G 115_C 0.87 0.11 0.00
297_L 22_A 0.87 0.11 0.00
292_F 88_Q 0.86 0.11 0.00
35_I 59_L 0.86 0.11 0.00
128_D 101_L 0.86 0.11 0.00
258_Y 69_K 0.86 0.11 0.00
196_S 139_L 0.85 0.10 0.00
95_Q 48_S 0.85 0.10 0.00
95_Q 89_A 0.85 0.10 0.00
289_A 114_T 0.84 0.10 0.00
200_C 8_E 0.84 0.10 0.00
27_D 38_G 0.84 0.10 0.00
12_T 142_I 0.84 0.10 0.00
95_Q 71_R 0.84 0.10 0.00
251_D 59_L 0.84 0.10 0.00
240_I 34_L 0.84 0.10 0.00
22_V 42_I 0.83 0.10 0.00
259_Y 129_L 0.83 0.10 0.00
186_R 101_L 0.83 0.10 0.00
259_Y 81_D 0.83 0.10 0.00
78_T 150_L 0.83 0.10 0.00
65_T 13_T 0.82 0.10 0.00
152_V 94_T 0.82 0.10 0.00
96_D 145_V 0.82 0.10 0.00
18_L 38_G 0.82 0.10 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.0326 seconds.