May 4, 2021 - We are working on upgrading the webserver, some pages may not work.
OPENSEQ.org
cI_B_cIII_cytb
Genes: A B A+B
Length: 613 440 982
Sequences: 12009 3168 3548
Seq/Len: 19.59 7.2 3.61
MirrorTree (Pazo et al. 2001) 0.71
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.39 0.00 0.26
2 0.43 0.00 1.42
5 0.60 0.00 2.45
10 0.60 0.00 2.63
20 0.60 0.00 2.69
100 0.60 0.00 2.79
∞ 0.60 0.01 3.37
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
465_T 60_A 1.13 0.82 0.08
556_I 60_A 0.98 0.68 0.05
181_L 70_T 0.91 0.59 0.04
464_S 364_K 0.90 0.59 0.04
13_I 169_I 0.90 0.58 0.04
13_I 334_I 0.86 0.53 0.03
464_S 176_I 0.82 0.48 0.03
415_M 364_K 0.81 0.47 0.03
413_N 261_V 0.81 0.46 0.03
462_I 248_F 0.80 0.46 0.03
421_G 334_I 0.80 0.45 0.03
62_Q 197_L 0.79 0.44 0.03
606_L 110_I 0.79 0.44 0.03
216_M 211_V 0.79 0.43 0.03
122_L 253_L 0.78 0.42 0.02
180_F 383_M 0.78 0.42 0.02
294_A 357_G 0.77 0.42 0.02
89_L 69_Y 0.77 0.42 0.02
544_L 207_A 0.77 0.41 0.02
600_I 260_L 0.77 0.41 0.02
355_A 96_L 0.76 0.41 0.02
379_F 121_Y 0.76 0.40 0.02
108_R 121_Y 0.75 0.39 0.02
430_F 37_T 0.75 0.39 0.02
571_I 207_A 0.75 0.39 0.02
272_L 145_M 0.75 0.39 0.02
507_L 206_A 0.74 0.38 0.02
181_L 356_S 0.74 0.38 0.02
602_A 338_G 0.74 0.38 0.02
606_L 395_A 0.74 0.37 0.02
273_F 131_V 0.74 0.37 0.02
295_G 41_L 0.73 0.37 0.02
217_W 137_L 0.73 0.37 0.02
286_V 112_I 0.73 0.36 0.02
369_R 71_P 0.72 0.36 0.02
474_P 286_L 0.72 0.35 0.02
273_F 136_Y 0.72 0.35 0.02
141_L 397_S 0.72 0.35 0.02
406_A 172_V 0.72 0.35 0.02
125_A 22_R 0.71 0.34 0.02
435_I 284_N 0.71 0.34 0.02
562_R 116_L 0.71 0.34 0.02
201_V 206_A 0.70 0.33 0.02
606_L 398_A 0.70 0.33 0.02
608_L 406_I 0.70 0.33 0.02
188_L 243_P 0.70 0.33 0.02
384_A 214_W 0.70 0.33 0.02
490_M 207_A 0.70 0.33 0.02
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 A, there is a high ratio (0.60 > 0.4) of paralogs.
ID Seq/Len Name Options I_Prob Status
2276 3.48 cI_B_cIII_cytb Δgene:(1, ∞) A:(1E-80, 8) B:(1E-20, 8) msa: Jackhmmer (2015_03) 0.02 Done - Shared
2275 3.58 cI_B_cIII_cytb Δgene:(1, ∞) A:(1E-60, 8) B:(1E-20, 8) msa: Jackhmmer (2015_03) 0.01 Done - Shared
2270 3.61 cI_B_cIII_cytb Δgene:(1, ∞) A:(1E-40, 8) B:(1E-20, 8) msa: Jackhmmer (2015_03) 0.08 Done - Shared
2268 3.63 cI_B_cIII_cytb Δgene:(1, ∞) A:(1E-20, 8) B:(1E-20, 8) msa: Jackhmmer (2015_03) 0.01 Done - Shared
Page generated in 0.2035 seconds.