May 4, 2021 - We are working on upgrading the webserver, some pages may not work.
OPENSEQ.org
ACRA_TOLC
Genes: A B A+B
Length: 397 493 806
Sequences: 17232 14400 7503
Seq/Len: 43.41 29.21 9.31
MirrorTree (Pazo et al. 2001) 0.57
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.07 0.05 4.25
2 0.08 0.06 7.74
5 0.10 0.08 8.43
10 0.11 0.10 8.58
20 0.13 0.12 8.73
100 0.17 0.16 9.43
∞ 0.22 0.21 10.35
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
135_T 382_A 2.46 1.00 1.00
131_K 386_V 1.92 1.00 0.98
136_Q 173_I 1.81 1.00 0.98
135_T 164_Q 1.75 1.00 0.97
136_Q 166_F 1.47 1.00 0.92
137_Y 382_A 1.29 0.99 0.85
131_K 168_V 1.27 0.98 0.84
129_Y 388_T 1.23 0.98 0.81
139_S 172_A 1.17 0.97 0.77
135_T 378_D 1.14 0.96 0.74
129_Y 170_L 1.09 0.95 0.69
137_Y 164_Q 1.03 0.93 0.63
141_Q 170_L 0.94 0.88 0.53
136_Q 385_S 0.84 0.78 0.40
137_Y 379_A 0.81 0.75 0.37
137_Y 161_Q 0.77 0.70 0.32
136_Q 157_R 0.73 0.63 0.27
135_T 160_D 0.71 0.61 0.25
142_E 383_G 0.71 0.60 0.25
132_L 170_L 0.61 0.44 0.15
265_V 276_T 0.60 0.42 0.15
135_T 167_N 0.60 0.42 0.14
139_S 171_V 0.54 0.32 0.10
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
5917 6.61 acra-tolc Δgene:(1, 20) A:(1E-20, 8) B:(1E-20, 8) msa: HHblits (2015_06) 0.99 Done - Shared
4292 2.61 ACRA_TOLC Δgene:(1, 5) A:(1E-80, 4) B:(1E-40, 4) msa: Jackhmmer (2015_06) 1.00 Done
4290 0.35 ACRA_TOLC Δgene:(1, 1) A:(1E-80, 4) B:(1E-40, 4) msa: Jackhmmer (2015_06) 0.01 Done
4288 0.03 ACRA_TOLC Δgene:(1, 1) A:(1E-80, 4) B:(1E-60, 4) msa: Jackhmmer (2015_06) Killed
4286 0.02 ACRA_TOLC Δgene:(1, 1) A:(1E-80, 4) B:(1E-80, 4) msa: Jackhmmer (2015_06) Killed
4283 4.75 ACRA_TOLC Δgene:(1, 1) A:(1E-20, 8) B:(1E-20, 8) msa: Jackhmmer (2015_06) 0.99 Done
4276 9.31 ACRA_TOLC Δgene:(1, 5) A:(1E-20, 8) B:(1E-20, 8) msa: Jackhmmer (2015_06) 1.00 Done
0669 4.67 A-C Δgene:(1, 20) A:(1E-20, 8) B:(1E-20, 8) msa: HHblits (2013_03) 0.99 Done
Page generated in 2.1759 seconds.