GREMLIN.BAKERLAB.org
Powered by OPENSEQ.org
MFD_MYCTU_1-450-UVRA_MYCTU_89-515

Genes: A B A+B
Length: 480 427 874
Sequences: 9529 7181 3650
Seq/Len: 19.85 16.82 4.18
MirrorTree (Pazo et al. 2001) 0.09
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.00 1.13
2 0.00 0.00 1.31
5 0.00 0.00 1.55
10 0.00 0.00 1.74
20 0.00 0.01 1.87
100 0.01 0.01 2.29
0.09 0.02 4.02
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
161_E 123_K 1.65 0.99 0.92
194_I 88_S 1.15 0.87 0.61
212_G 121_T 1.09 0.82 0.55
211_W 119_R 0.97 0.71 0.42
216_T 119_R 0.96 0.70 0.41
193_G 119_R 0.93 0.66 0.37
161_E 58_G 0.88 0.60 0.32
211_W 89_R 0.86 0.57 0.29
193_G 89_R 0.79 0.48 0.23
212_G 119_R 0.78 0.47 0.22
447_E 361_S 0.76 0.44 0.20
165_D 128_R 0.74 0.41 0.18
216_T 89_R 0.73 0.40 0.17
213_D 123_K 0.69 0.35 0.14
105_V 147_L 0.69 0.34 0.14
462_K 195_P 0.67 0.32 0.13
438_A 405_S 0.66 0.32 0.13
231_E 100_T 0.65 0.30 0.12
190_V 132_D 0.65 0.30 0.12
214_E 119_R 0.64 0.29 0.11
216_T 108_Q 0.63 0.28 0.11
224_A 135_E 0.63 0.28 0.11
216_T 58_G 0.62 0.27 0.10
374_A 67_V 0.62 0.27 0.10
462_K 181_E 0.61 0.26 0.10
24_P 65_A 0.60 0.25 0.09
212_G 89_R 0.60 0.25 0.09
241_A 254_P 0.60 0.25 0.09
64_V 209_E 0.60 0.25 0.09
217_E 209_E 0.60 0.25 0.09
269_V 14_I 0.59 0.24 0.09
108_V 217_R 0.59 0.24 0.09
330_I 116_V 0.59 0.24 0.09
303_L 302_K 0.59 0.24 0.09
452_A 266_I 0.59 0.24 0.09
458_G 425_T 0.59 0.23 0.09
239_A 384_Q 0.58 0.23 0.08
212_G 120_L 0.58 0.23 0.08
197_I 412_T 0.58 0.23 0.08
373_A 63_V 0.58 0.23 0.08
222_S 293_E 0.58 0.23 0.08
314_V 161_Q 0.58 0.23 0.08
227_R 83_N 0.58 0.23 0.08
114_A 87_Y 0.58 0.23 0.08
34_G 382_A 0.57 0.22 0.08
315_L 224_V 0.57 0.22 0.08
188_F 404_L 0.57 0.22 0.08
180_D 305_Q 0.57 0.21 0.08
315_L 120_L 0.56 0.21 0.07
459_Q 96_V 0.56 0.21 0.07
287_D 5_R 0.56 0.21 0.07
314_V 127_K 0.56 0.21 0.07
443_E 47_Q 0.56 0.21 0.07
240_F 131_T 0.56 0.21 0.07
214_E 89_R 0.56 0.21 0.07
307_Q 299_L 0.56 0.21 0.07
30_M 48_I 0.55 0.20 0.07
33_A 189_S 0.55 0.20 0.07
153_P 131_T 0.55 0.20 0.07
196_D 385_V 0.55 0.20 0.07
170_R 385_V 0.55 0.20 0.07
208_V 94_G 0.55 0.20 0.07
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
7440 1.81 MFD_MYCTU_1-450-UVRA_MYCTU_89-515 Δgene:(1, 10) A:(1E-20, 8) B:(1E-20, 8) msa: Jackhmmer (r132) 0.10 Done - Shared
7439 1.17 MFD_MYCTU_1-450-UVRA_MYCTU_89-515 Δgene:(1, 1) A:(1E-20, 8) B:(1E-20, 8) msa: Jackhmmer (r132) 0.00 Done - Shared
7375 4.18 MFD_MYCTU_1-450-UVRA_MYCTU_89-515 Δgene:(1, ∞) A:(1E-20, 8) B:(1E-20, 8) msa: Jackhmmer (r132) 0.92 Done - Shared

Page generated in 0.0371 seconds.