Citrus Sinensis ID: 040703

---

Local Sequence Feature Prediction

Residue Number Marker
--------10--------20--------30--------40--------50--------60--------70--------80--------
Protein Sequence
VLYTLPFRDSISLLSCIDTICHDPTQIVISVFYQVDPPDMRNQNGNFGDSFLKLEERLKENTEKLQSWRNALKEFMPSLSGFHSLNIT
Secondary Structure (PSIPRED)
cEEcccccHHHHHHHHHHHcccccccEEEEEEccccccccccccccHHHHHHHHHHHHcccHHHHHHHHHHHHHHccccccccccccc
Secondary Structure Prediction (SSPRO)
ccccccccccHHHHHHHHHHHHccccEEEEEEEEccHHHHHHccccHHHHHHHHHHHHHHcHHHHHHHHHHHHHHHHHcccccccccc
Coil and Loop (DISEMBL)
vlytlpfrdsisllscidtichdptQIVISVFyqvdppdmrnqngnfgdSFLKLEERLKENTEKLQSWRNALKEfmpslsgfhslnit
Flexible Loop (DISEMBL)
vlytlpfrdsiSLLSCIDTICHDPTQIVISVFYQVDPPDMRNQNGNFGDSFLKLEERLKENTEKLQSWRNALKefmpslsgfhslnit
Low Complexity Region (SEG)
VLYTLPFRDSISLLSCIDTICHDPTQIVISVFYQVDPPDMRNQNGNFGDSFLKLEERLKENTEKLQSWRNALKEFMPSLSGFHSLNIT
Disordered region (IsUnstruct)
**YTLPFRDSISLLSCIDTICHDPTQIVISVFYQVD****************************************************
Disordered Region (DISOPRED)
VLYTLPFRDSISLLSCIDTICHDPTQIVISVFYQVDPPDMRNQNGNFGDSFLKLEER*****EKLQSWRNALKEFMPSLSGFHSLNI*
Disordered Region (DISEMBL)
VLYTLPFRDSISLLSCIDTICHDPTQIVISVFYQVDPPDMRNQNGNFGDSFLKLEERLKENTEKLQSWRNALKEFMPSLSGFHSLNIT
Disordered Region (DISPRO)
VLYTLPFRDSISLLSCIDTICHDPTQIVISVFYQVDPPDMRNQNGNFGDSFLKLEERLKENTEKLQSWRNALKEFMPSLS********
Transmembrane Helix (TMHMM)
oooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo
Transmembrane Helix (HMMTOP)
oooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo
Transmembrane Helix (MEMSAT)
ooooooooooooooooooooooooooooooooooooooooooooooooooooooohhhhhhhhhhhhhhhhhhhhiiiiiiiiiiiii
TM Helix, Signal Peptide (MEMSAT_SVM)
oooooooooooooooooooooooohhhhhhhhhhhhhhhhiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiii
TM Helix, Signal Peptide (Phobius)
oooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo
Signal Peptide (SignalP HMM Mode)
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
Signal Peptide (SignalP NN Mode)
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
Coiled Coils (COILS)
VLYTLPFRDSISLLSCIDTICHDPTQIVISVFYQVDPPDMRNQNGNFGDSxxxxxxxxxxxxxxxxxxxxxLKEFMPSLSGFHSLNIT
Positional Conservation
no confident homologs detected

---

Close Homologs for Annotation Transfer

Close Homologs in SWISS-PROT Database Detected by BLAST

Original result of BLAST against SWISS-PROT Database

ID	Length	Definition	RBH(Q2H)	RBH(H2Q)	Q cover	H cover	Identity	E-value
Query	88	2.2.26 [Sep-21-2011]
Q40392	1144	TMV resistance protein N	N/A	no	0.556	0.042	0.448	5e-06
Q9FHE9	354	Protein PHLOEM PROTEIN 2-	no	no	0.602	0.149	0.388	6e-05
O82500	1095	Putative disease resistan	no	no	0.488	0.039	0.488	0.0001

>sp|Q40392|TMVRN_NICGU TMV resistance protein N OS=Nicotiana glutinosa GN=N PE=1 SV=1

 Score = 49.7 bits (117), Expect = 5e-06,   Method: Compositional matrix adjust.
 Identities = 22/49 (44%), Positives = 31/49 (63%)

Query: 26  QIVISVFYQVDPPDMRNQNGNFGDSFLKLEERLKENTEKLQSWRNALKE 74
           Q VI +FY VDP  +RNQ  +F  +F + E + K++ E +Q WR AL E
Sbjct: 99  QTVIPIFYDVDPSHVRNQKESFAKAFEEHETKYKDDVEGIQRWRIALNE 147

Disease resistance protein. Resistance proteins guard the plant against pathogens that contain an appropriate avirulence protein via a direct or indirect interaction with this avirulence protein. That triggers a defense system including the hypersensitive response, which restricts the pathogen growth.
Nicotiana glutinosa (taxid: 35889)

>sp|Q9FHE9|P2A08_ARATH Protein PHLOEM PROTEIN 2-LIKE A8 OS=Arabidopsis thaliana GN=PP2A8 PE=2 SV=1

 Score = 46.2 bits (108), Expect = 6e-05,   Method: Composition-based stats.
 Identities = 21/54 (38%), Positives = 32/54 (59%), Gaps = 1/54 (1%)

Query: 29  ISVFYQVDPPDMRNQNGNFGDSFLKLEERLKENTEKLQSWRNALKEFMPSLSGF 82
           I +FY++ P  +    G FGD+F  L+E+ K + E+ Q W+ AL E +P L G 
Sbjct: 102 IPIFYKLAPSSVLELKGGFGDTFRVLKEKYKNDPERTQKWQEAL-ESIPKLKGL 154

Arabidopsis thaliana (taxid: 3702)

>sp|O82500|Y4117_ARATH Putative disease resistance protein At4g11170 OS=Arabidopsis thaliana GN=At4g11170 PE=2 SV=1

 Score = 45.1 bits (105), Expect = 1e-04,   Method: Compositional matrix adjust.
 Identities = 22/45 (48%), Positives = 28/45 (62%), Gaps = 2/45 (4%)

Query: 28  VISVFYQVDPPDMRNQNGNFGDSFLKLEERLKENTEKLQSWRNAL 72
           V+ VFY+VDP D+R Q G FG SF  LE    +  E+  +WR AL
Sbjct: 99  VMPVFYKVDPSDIRKQTGKFGMSF--LETCCGKTEERQHNWRRAL 141

Arabidopsis thaliana (taxid: 3702)

Close Homologs in the Non-Redundant Database Detected by BLAST

Original result of BLAST against Nonredundant Database

GI	Length	Definition	Q cover	H cover	Identity	E-value
Query	88
351727669	711	TIR-NBS-LRR type disease resistance prot	0.647	0.080	0.568	3e-11
148807171	176	putative TIR-NBS-LRR class disease resis	0.795	0.397	0.447	5e-11
255555357	1094	leucine-rich repeat-containing protein,	0.784	0.063	0.486	3e-10
357462141	1545	Tir-nbs-lrr resistance protein [Medicago	0.715	0.040	0.484	3e-10
255561496	876	TMV resistance protein N, putative [Rici	0.670	0.067	0.583	4e-10
27764539	895	R 4 protein [Glycine max]	0.761	0.074	0.438	4e-10
22037313	158	disease resistance-like protein GS0-1 [G	0.761	0.424	0.465	4e-10
356560037	1083	PREDICTED: TMV resistance protein N-like	0.761	0.061	0.438	5e-10
224111328	188	predicted protein [Populus trichocarpa]	0.738	0.345	0.472	6e-10
105923279	1021	TIR-NBS-LRR type disease resistance prot	0.75	0.064	0.465	7e-10

>gi|351727669|ref|NP_001235633.1| TIR-NBS-LRR type disease resistance protein [Glycine max] gi|223452586|gb|ACM89620.1| TIR-NBS-LRR type disease resistance protein [Glycine max]

 Score = 72.4 bits (176), Expect = 3e-11,   Method: Compositional matrix adjust.
 Identities = 33/58 (56%), Positives = 45/58 (77%), Gaps = 1/58 (1%)

Query: 26 QIVISVFYQVDPPDMRNQNGNFGDSFLKLEERLKENTEKLQSWRNALKEFMPSLSGFH 83
          QI++ VF+ VDP D+R Q+G +GD+  K EE+LKEN  K+QSWR+ALK+   +LSGFH
Sbjct: 10 QILVPVFFNVDPSDVRQQHGEYGDALAKHEEKLKENMFKVQSWRSALKK-AANLSGFH 66

Source: Glycine max

Species: Glycine max

Genus: Glycine

Family: Fabaceae

Order: Fabales

Phylum: Streptophyta

Superkingdom: Eukaryota

>gi|148807171|gb|ABR13295.1| putative TIR-NBS-LRR class disease resistance protein [Prunus dulcis]

 Score = 72.0 bits (175), Expect = 5e-11,   Method: Compositional matrix adjust.
 Identities = 34/76 (44%), Positives = 48/76 (63%), Gaps = 6/76 (7%)

Query: 16  CIDTI-----CHDPTQIVISVFYQVDPPDMRNQNGNFGDSFLKLEERLKENTEKLQSWRN 70
           C+D +     C      ++ VFY V+P D+R Q+G+F D F KLEER +E+ E+++SWR 
Sbjct: 36  CLDELTEILQCMKSKGTILPVFYHVEPSDVRKQSGSFADVFAKLEERFREDIERVKSWRA 95

Query: 71  ALKEFMPSLSGFHSLN 86
           AL E + SLSG  S N
Sbjct: 96  ALTE-VASLSGLDSKN 110

Source: Prunus dulcis

Species: Prunus dulcis

Genus: Prunus

Family: Rosaceae

Order: Rosales

Phylum: Streptophyta

Superkingdom: Eukaryota

>gi|255555357|ref|XP_002518715.1| leucine-rich repeat-containing protein, putative [Ricinus communis] gi|223542096|gb|EEF43640.1| leucine-rich repeat-containing protein, putative [Ricinus communis]

 Score = 69.3 bits (168), Expect = 3e-10,   Method: Compositional matrix adjust.
 Identities = 36/74 (48%), Positives = 50/74 (67%), Gaps = 5/74 (6%)

Query: 11  ISLLSCIDTICHDPTQIVISVFYQVDPPDMRNQNGNFGDSFLKLEERLKENTEKLQSWRN 70
           + +L C+     +  Q+VI VFY+VDP  +RNQ G+F D F + EE L  + EK+QSWR 
Sbjct: 98  VKILECM----KNRGQMVIPVFYRVDPSHVRNQTGSFEDVFARHEESLSVSKEKVQSWRA 153

Query: 71  ALKEFMPSLSGFHS 84
           ALKE + +LSG+HS
Sbjct: 154 ALKE-VANLSGWHS 166

Source: Ricinus communis

Species: Ricinus communis

Genus: Ricinus

Family: Euphorbiaceae

Order: Malpighiales

Phylum: Streptophyta

Superkingdom: Eukaryota

>gi|357462141|ref|XP_003601352.1| Tir-nbs-lrr resistance protein [Medicago truncatula] gi|355490400|gb|AES71603.1| Tir-nbs-lrr resistance protein [Medicago truncatula]

 Score = 69.3 bits (168), Expect = 3e-10,   Method: Compositional matrix adjust.
 Identities = 31/64 (48%), Positives = 45/64 (70%), Gaps = 1/64 (1%)

Query: 23  DPTQIVISVFYQVDPPDMRNQNGNFGDSFLKLEERLKENTEKLQSWRNALKEFMPSLSGF 82
           D  QIVI +FY++DP  +RNQNG++G +F K    LK+N E L+ W++AL E   +L+G+
Sbjct: 540 DQGQIVIPIFYEIDPSHVRNQNGSYGQAFAKHARDLKQNKEMLKKWKDALTE-AANLAGW 598

Query: 83  HSLN 86
           HS N
Sbjct: 599 HSQN 602

Source: Medicago truncatula

Species: Medicago truncatula

Genus: Medicago

Family: Fabaceae

Order: Fabales

Phylum: Streptophyta

Superkingdom: Eukaryota

>gi|255561496|ref|XP_002521758.1| TMV resistance protein N, putative [Ricinus communis] gi|223538971|gb|EEF40568.1| TMV resistance protein N, putative [Ricinus communis]

 Score = 68.9 bits (167), Expect = 4e-10,   Method: Compositional matrix adjust.
 Identities = 35/60 (58%), Positives = 45/60 (75%), Gaps = 1/60 (1%)

Query: 26  QIVISVFYQVDPPDMRNQNGNFGDSFLKLEERLKENTEKLQSWRNALKEFMPSLSGFHSL 85
           QIVI VFY+VDP  +RNQ G+F  +F K EE LKE  EK++SWR+AL E   S+SG++SL
Sbjct: 97  QIVIPVFYRVDPVHVRNQRGSFACAFAKHEETLKERMEKVESWRSALNE-AGSISGWNSL 155

Source: Ricinus communis

Species: Ricinus communis

Genus: Ricinus

Family: Euphorbiaceae

Order: Malpighiales

Phylum: Streptophyta

Superkingdom: Eukaryota

>gi|27764539|gb|AAO23069.1| R 4 protein [Glycine max]

 Score = 68.9 bits (167), Expect = 4e-10,   Method: Compositional matrix adjust.
 Identities = 32/73 (43%), Positives = 47/73 (64%), Gaps = 6/73 (8%)

Query: 16  CIDTI-----CHDPTQIVISVFYQVDPPDMRNQNGNFGDSFLKLEERLKENTEKLQSWRN 70
           C+D +     C     +VI VFY+VDP D+R+Q G++G++  K ++R K   EKLQ WR 
Sbjct: 83  CLDELVTILHCKSEGLLVIPVFYKVDPSDVRHQKGSYGEAMAKHQKRFKAKKEKLQKWRM 142

Query: 71  ALKEFMPSLSGFH 83
           ALK+ +  LSG+H
Sbjct: 143 ALKQ-VADLSGYH 154

Source: Glycine max

Species: Glycine max

Genus: Glycine

Family: Fabaceae

Order: Fabales

Phylum: Streptophyta

Superkingdom: Eukaryota

>gi|22037313|gb|AAM89998.1|AF403250_1 disease resistance-like protein GS0-1 [Glycine max] gi|22037315|gb|AAM89999.1|AF403251_1 disease resistance-like protein GS0-2 [Glycine max]

 Score = 68.9 bits (167), Expect = 4e-10,   Method: Compositional matrix adjust.
 Identities = 34/73 (46%), Positives = 49/73 (67%), Gaps = 6/73 (8%)

Query: 16  CID---TICHDPTQ--IVISVFYQVDPPDMRNQNGNFGDSFLKLEERLKENTEKLQSWRN 70
           C+D   TI H  +Q  +VI VFY+VDP  +R+Q G++G++  K ++R K N EKLQ WR 
Sbjct: 83  CLDELVTILHCKSQGLLVIPVFYKVDPSHVRHQKGSYGEAMAKHQKRFKANKEKLQKWRM 142

Query: 71  ALKEFMPSLSGFH 83
           AL + +  LSG+H
Sbjct: 143 ALHQ-VADLSGYH 154

Source: Glycine max

Species: Glycine max

Genus: Glycine

Family: Fabaceae

Order: Fabales

Phylum: Streptophyta

Superkingdom: Eukaryota

>gi|356560037|ref|XP_003548302.1| PREDICTED: TMV resistance protein N-like [Glycine max]

 Score = 68.6 bits (166), Expect = 5e-10,   Method: Compositional matrix adjust.
 Identities = 32/73 (43%), Positives = 47/73 (64%), Gaps = 6/73 (8%)

Query: 16  CIDTI-----CHDPTQIVISVFYQVDPPDMRNQNGNFGDSFLKLEERLKENTEKLQSWRN 70
           C+D +     C     +VI VFY+VDP D+R+Q G++G++  K ++R K   EKLQ WR 
Sbjct: 83  CLDELVTILHCKSEGLLVIPVFYKVDPSDVRHQKGSYGEAMAKHQKRFKAKKEKLQKWRM 142

Query: 71  ALKEFMPSLSGFH 83
           ALK+ +  LSG+H
Sbjct: 143 ALKQ-VADLSGYH 154

Source: Glycine max

Species: Glycine max

Genus: Glycine

Family: Fabaceae

Order: Fabales

Phylum: Streptophyta

Superkingdom: Eukaryota

>gi|224111328|ref|XP_002332938.1| predicted protein [Populus trichocarpa] gi|222834192|gb|EEE72669.1| predicted protein [Populus trichocarpa]

 Score = 68.2 bits (165), Expect = 6e-10,   Method: Compositional matrix adjust.
 Identities = 34/72 (47%), Positives = 46/72 (63%), Gaps = 7/72 (9%)

Query: 16  CIDTI-----CHDPT-QIVISVFYQVDPPDMRNQNGNFGDSFLKLEERLKENTEKLQSWR 69
           C+D +     C  PT QIV+ VFY +DP D+R Q G+F ++F + EER KE  EK+Q WR
Sbjct: 86  CLDELEKILDCRQPTGQIVLPVFYDIDPSDIRKQTGSFAEAFDRHEERFKEEMEKVQKWR 145

Query: 70  NALKEFMPSLSG 81
            AL E   ++SG
Sbjct: 146 KALVE-AANISG 156

Source: Populus trichocarpa

Species: Populus trichocarpa

Genus: Populus

Family: Salicaceae

Order: Malpighiales

Phylum: Streptophyta

Superkingdom: Eukaryota

>gi|105923279|gb|ABF81468.1| TIR-NBS-LRR type disease resistance protein [Populus trichocarpa]

 Score = 67.8 bits (164), Expect = 7e-10,   Method: Compositional matrix adjust.
 Identities = 34/73 (46%), Positives = 46/73 (63%), Gaps = 7/73 (9%)

Query: 16  CIDTI-----CHDPT-QIVISVFYQVDPPDMRNQNGNFGDSFLKLEERLKENTEKLQSWR 69
           C+D +     C  PT QIV+ VFY +DP D+R Q G+F ++F + EER KE  EK+Q WR
Sbjct: 110 CLDELEKILDCRQPTGQIVLPVFYDIDPSDIRKQTGSFAEAFDRHEERFKEEMEKVQKWR 169

Query: 70  NALKEFMPSLSGF 82
            AL E   ++SG 
Sbjct: 170 KALVE-AANISGL 181

Source: Populus trichocarpa

Species: Populus trichocarpa

Genus: Populus

Family: Salicaceae

Order: Malpighiales

Phylum: Streptophyta

Superkingdom: Eukaryota

---

Prediction of Gene Ontology (GO) Terms

Close Homologs with Gene Ontology terms Detected by BLAST

Original result of BLAST against Gene Ontology (AMIGO)

ID	Length	Definition	Q cover	H cover	Identity	E-value
Query	88
TAIR|locus:2206455	172	AT1G57630 "AT1G57630" [Arabido	0.772	0.395	0.376	9.7e-09
TAIR|locus:2160472	1038	AT5G41540 [Arabidopsis thalian	0.806	0.068	0.397	2e-07
TAIR|locus:2175075	1068	AT5G41750 [Arabidopsis thalian	0.784	0.064	0.394	3.3e-07
TAIR|locus:2160487	1085	AT5G41550 [Arabidopsis thalian	0.784	0.063	0.381	7e-07
TAIR|locus:1006230150	313	AT2G20142 [Arabidopsis thalian	0.875	0.246	0.364	1.1e-06
TAIR|locus:2076043	1194	RPP1 "recognition of peronospo	0.625	0.046	0.416	1.3e-06
TAIR|locus:2100182	1007	AT3G44400 [Arabidopsis thalian	0.738	0.064	0.355	1.3e-06
TAIR|locus:2098110	1219	AT3G44670 [Arabidopsis thalian	0.647	0.046	0.419	1.7e-06
TAIR|locus:2094498	1981	AT3G25510 [Arabidopsis thalian	0.681	0.030	0.430	1.8e-06
UNIPROTKB|Q40392	1144	N "TMV resistance protein N" [	0.556	0.042	0.448	2e-06

TAIR|locus:2206455 AT1G57630 "AT1G57630" [Arabidopsis thaliana (taxid:3702)]

 Score = 131 (51.2 bits), Expect = 9.7e-09, P = 9.7e-09
 Identities = 29/77 (37%), Positives = 49/77 (63%)

Query:    16 CIDTICH----DPTQIVISVFYQVDPPDMRNQNGNFGDSFLKLEERLKENTE-KLQSWRN 70
             C+D +      +  Q VI++FY+VDP +++ Q G+FG +F   +E  +  TE K+Q+WR 
Sbjct:    93 CLDELVEIMKKESGQTVITIFYEVDPNEVKKQTGDFGKAF---KETCQGKTEEKVQTWRK 149

Query:    71 ALKEF-MPSLSGFHSLN 86
             AL+   + +++GFHS N
Sbjct:   150 ALEGVGLATIAGFHSSN 166

GO:0006952 "defense response" evidence=ISS

GO:0007165 "signal transduction" evidence=IEA

GO:0009507 "chloroplast" evidence=ISM

GO:0071456 "cellular response to hypoxia" evidence=IEP

GO:0000165 "MAPK cascade" evidence=RCA

GO:0002237 "response to molecule of bacterial origin" evidence=RCA

GO:0006612 "protein targeting to membrane" evidence=RCA

GO:0009697 "salicylic acid biosynthetic process" evidence=RCA

GO:0009814 "defense response, incompatible interaction" evidence=RCA

GO:0009863 "salicylic acid mediated signaling pathway" evidence=RCA

GO:0009867 "jasmonic acid mediated signaling pathway" evidence=RCA

GO:0010363 "regulation of plant-type hypersensitive response" evidence=RCA

GO:0031348 "negative regulation of defense response" evidence=RCA

GO:0043090 "amino acid import" evidence=RCA

GO:0045087 "innate immune response" evidence=RCA

TAIR|locus:2160472 AT5G41540 [Arabidopsis thaliana (taxid:3702)]

 Score = 132 (51.5 bits), Expect = 2.0e-07, P = 2.0e-07
 Identities = 31/78 (39%), Positives = 46/78 (58%)

Query:    11 ISLLSCIDTICHDPTQIVISVFYQVDPPDMRNQNGNFGDSFLKLEERLKENTEKLQSWRN 70
             + +L C +    D  QIV+ +FY++DP D+R Q+G+FG +F K    + +  E  Q W N
Sbjct:    86 VEILKCKE----DQEQIVMPIFYEIDPSDVRKQSGDFGKAFGKT--CVGKTKEVKQRWTN 139

Query:    71 ALKEFMPSLSGFHSLNIT 88
             AL E   ++ G HSLN T
Sbjct:   140 ALTE-AANIGGEHSLNWT 156

GO:0000166 "nucleotide binding" evidence=IEA

GO:0005622 "intracellular" evidence=IEA

GO:0005737 "cytoplasm" evidence=ISM

GO:0006952 "defense response" evidence=IEA;ISS

GO:0007165 "signal transduction" evidence=IEA

GO:0017111 "nucleoside-triphosphatase activity" evidence=IEA

GO:0043531 "ADP binding" evidence=IEA

TAIR|locus:2175075 AT5G41750 [Arabidopsis thaliana (taxid:3702)]

 Score = 130 (50.8 bits), Expect = 3.3e-07, P = 3.3e-07
 Identities = 30/76 (39%), Positives = 49/76 (64%)

Query:    11 ISLLSCIDTICHDPTQIVISVFYQVDPPDMRNQNGNFGDSFLKLEERLKENTEKLQSWRN 70
             + +L C + +     QIV++VFY+VDP D++ Q+G FG++F K  +   E   K++ WRN
Sbjct:    88 VEILKCKEAL----GQIVMTVFYEVDPSDVKKQSGVFGEAFEKTCQGKNEEV-KIR-WRN 141

Query:    71 ALKEFMPSLSGFHSLN 86
             AL   + +++G HSLN
Sbjct:   142 ALAH-VATIAGEHSLN 156

GO:0000166 "nucleotide binding" evidence=IEA

GO:0005622 "intracellular" evidence=IEA

GO:0005634 "nucleus" evidence=ISM

GO:0006952 "defense response" evidence=IEA;ISS

GO:0007165 "signal transduction" evidence=IEA

GO:0017111 "nucleoside-triphosphatase activity" evidence=IEA

GO:0043531 "ADP binding" evidence=IEA

TAIR|locus:2160487 AT5G41550 [Arabidopsis thaliana (taxid:3702)]

 Score = 127 (49.8 bits), Expect = 7.0e-07, P = 7.0e-07
 Identities = 29/76 (38%), Positives = 43/76 (56%)

Query:    11 ISLLSCIDTICHDPTQIVISVFYQVDPPDMRNQNGNFGDSFLKLEERLKENTEKLQSWRN 70
             + +L C +    D  QIV+++FY VDP  +R Q G+FG +F+K  E   E  +  Q W  
Sbjct:    88 VEILKCKE----DQGQIVMTIFYDVDPSSVRKQKGDFGSTFMKTCEGKSEEVK--QRWTK 141

Query:    71 ALKEFMPSLSGFHSLN 86
             AL   + ++ G HSLN
Sbjct:   142 ALTH-VANIKGEHSLN 156

GO:0000166 "nucleotide binding" evidence=IEA

GO:0005622 "intracellular" evidence=IEA

GO:0006952 "defense response" evidence=IEA;ISS

GO:0007165 "signal transduction" evidence=IEA

GO:0017111 "nucleoside-triphosphatase activity" evidence=IEA

GO:0043531 "ADP binding" evidence=IEA

TAIR|locus:1006230150 AT2G20142 [Arabidopsis thaliana (taxid:3702)]

 Score = 118 (46.6 bits), Expect = 1.1e-06, P = 1.1e-06
 Identities = 31/85 (36%), Positives = 49/85 (57%)

Query:     1 VLYTLPFRDSISLLSCIDTI--CHDPT-QIVISVFYQVDPPDMRNQNGNFGDSFLKLEER 57
             ++++  F DS   L  +  I  C +   QIVI V Y VDP D+ NQ G+FGD+F K  ++
Sbjct:    66 IVFSENFGDSRWCLDEVVAILKCKEKFGQIVIPVLYHVDPLDIENQTGSFGDAFAKRRDK 125

Query:    58 LKENTEKLQSWRNALKEFMPSLSGF 82
                  E+LQ W+++  E + +L G+
Sbjct:   126 ----AEQLQEWKDSFTEAI-NLPGW 145

GO:0005622 "intracellular" evidence=IEA

GO:0005737 "cytoplasm" evidence=ISM

GO:0007165 "signal transduction" evidence=IEA

TAIR|locus:2076043 RPP1 "recognition of peronospora parasitica 1" [Arabidopsis thaliana (taxid:3702)]

 Score = 125 (49.1 bits), Expect = 1.3e-06, P = 1.3e-06
 Identities = 25/60 (41%), Positives = 41/60 (68%)

Query:    26 QIVISVFYQVDPPDMRNQNGNFGDSFLKL-EERLKENTEKLQSWRNALKEFMPSLSGFHS 84
             QIV+++FY+VDP D++ Q G FG +F K    + KE  E+   WR AL++ + +++G+HS
Sbjct:   184 QIVMTIFYEVDPTDIKKQTGEFGKAFTKTCRGKPKEQVER---WRKALED-VATIAGYHS 239

GO:0005737 "cytoplasm" evidence=ISM

GO:0006952 "defense response" evidence=IEA;ISS;TAS

GO:0007165 "signal transduction" evidence=IEA

GO:0043531 "ADP binding" evidence=IEA

GO:0002239 "response to oomycetes" evidence=IDA

GO:0000139 "Golgi membrane" evidence=IDA

GO:0005789 "endoplasmic reticulum membrane" evidence=IDA

GO:0005886 "plasma membrane" evidence=IDA

GO:0009817 "defense response to fungus, incompatible interaction" evidence=IMP;TAS

GO:0030275 "LRR domain binding" evidence=ISS

TAIR|locus:2100182 AT3G44400 [Arabidopsis thaliana (taxid:3702)]

 Score = 124 (48.7 bits), Expect = 1.3e-06, P = 1.3e-06
 Identities = 27/76 (35%), Positives = 48/76 (63%)

Query:    16 CIDTI-----CHDPT-QIVISVFYQVDPPDMRNQNGNFGDSFLKL-EERLKENTEKLQSW 68
             C+D +     C +   QIV+++FY+V+P D++ Q G FG +F K    + KE+ E+   W
Sbjct:   118 CLDELAEIMKCREMVGQIVMTIFYEVEPTDIKKQTGEFGKAFTKTCRGKTKEHIER---W 174

Query:    69 RNALKEFMPSLSGFHS 84
             R AL++ + +++G+HS
Sbjct:   175 RKALED-VATIAGYHS 189

GO:0000166 "nucleotide binding" evidence=IEA

GO:0005622 "intracellular" evidence=IEA

GO:0005737 "cytoplasm" evidence=ISM

GO:0006952 "defense response" evidence=IEA;ISS

GO:0007165 "signal transduction" evidence=IEA

GO:0017111 "nucleoside-triphosphatase activity" evidence=IEA

GO:0043531 "ADP binding" evidence=IEA

TAIR|locus:2098110 AT3G44670 [Arabidopsis thaliana (taxid:3702)]

 Score = 124 (48.7 bits), Expect = 1.7e-06, P = 1.7e-06
 Identities = 26/62 (41%), Positives = 42/62 (67%)

Query:    26 QIVISVFYQVDPPDMRNQNGNFGDSFLKL-EERLKENTEKLQSWRNALKEFMPSLSGFHS 84
             QIV+++FY+VDP D++ Q G FG +F K  + + KE  E+   WR AL++ + +++G HS
Sbjct:   180 QIVMTIFYEVDPTDIKKQTGEFGKAFTKTCKGKTKEYVER---WRKALED-VATIAGEHS 235

Query:    85 LN 86
              N
Sbjct:   236 RN 237

GO:0000166 "nucleotide binding" evidence=IEA

GO:0005622 "intracellular" evidence=IEA

GO:0005737 "cytoplasm" evidence=ISM

GO:0006952 "defense response" evidence=IEA;ISS;TAS

GO:0007165 "signal transduction" evidence=IEA

GO:0017111 "nucleoside-triphosphatase activity" evidence=IEA

GO:0043531 "ADP binding" evidence=IEA

GO:0042742 "defense response to bacterium" evidence=ISS

TAIR|locus:2094498 AT3G25510 [Arabidopsis thaliana (taxid:3702)]

 Score = 126 (49.4 bits), Expect = 1.8e-06, P = 1.8e-06
 Identities = 28/65 (43%), Positives = 43/65 (66%)

Query:    23 DPTQIVISVFYQVDPPDMRNQNGNFGDSFLKLEERLKENT-EKLQSWRNALKEFMPSLSG 81
             D  QIVI++FY+VDP  ++ Q G+FG  F   +E  K  T E+++ WR AL E + +++G
Sbjct:   136 DLGQIVITIFYEVDPTHIKKQTGDFGKVF---KETCKGKTKEEIKRWRKAL-EGVATIAG 191

Query:    82 FHSLN 86
             +HS N
Sbjct:   192 YHSSN 196

GO:0000166 "nucleotide binding" evidence=IEA

GO:0005622 "intracellular" evidence=IEA

GO:0006952 "defense response" evidence=IEA;ISS

GO:0007165 "signal transduction" evidence=IEA

GO:0017111 "nucleoside-triphosphatase activity" evidence=IEA

GO:0043531 "ADP binding" evidence=IEA

UNIPROTKB|Q40392 N "TMV resistance protein N" [Nicotiana glutinosa (taxid:35889)]

 Score = 113 (44.8 bits), Expect = 2.0e-06, Sum P(2) = 2.0e-06
 Identities = 22/49 (44%), Positives = 31/49 (63%)

Query:    26 QIVISVFYQVDPPDMRNQNGNFGDSFLKLEERLKENTEKLQSWRNALKE 74
             Q VI +FY VDP  +RNQ  +F  +F + E + K++ E +Q WR AL E
Sbjct:    99 QTVIPIFYDVDPSHVRNQKESFAKAFEEHETKYKDDVEGIQRWRIALNE 147

GO:0005515 "protein binding" evidence=IPI

---

Prediction of Enzyme Commission (EC) Number

EC Number Prediction by Annotation Transfer from SWISS-PROT Entries

Original result of BLAST against SWISS-PROT

No confident hit for EC number transfering in SWISSPROT detected by BLAST

EC Number Prediction by Ezypred Server

Original result from Ezypred Server

Fail to connect to Ezypred Server

EC Number Prediction by EFICAz Software

No EC number assignment, probably not an enzyme!

---

Prediction of Functionally Associated Proteins

Functionally Associated Proteins Detected by STRING

Original result from the STRING server

Your Input: GSVIVG00021873001 SubName- Full=Chromosome chr18 scaffold_24, whole genome shotgun sequence; (890 aa) (Vitis vinifera)

Predicted Functional Partners:

Sorry, there are no predicted associations at the current settings.

---

Conserved Domains and Related Protein Families

Conserved Domains Detected by RPS-BLAST

Original result of RPS-BLAST against CDD database part I

ID	Length	Definition	E-value
Query	88
pfam01582	135	pfam01582, TIR, TIR domain	5e-09
PLN03210	1153	PLN03210, PLN03210, Resistant to P	2e-07
smart00255	140	smart00255, TIR, Toll - interleukin 1 - resistance	7e-05

>gnl|CDD|216585 pfam01582, TIR, TIR domain

 Score = 49.2 bits (118), Expect = 5e-09
 Identities = 18/48 (37%), Positives = 33/48 (68%)

Query: 26  QIVISVFYQVDPPDMRNQNGNFGDSFLKLEERLKENTEKLQSWRNALK 73
           ++++ +FY+VDP D+R Q+G FG +FLK  +   +  +K++ W+ AL 
Sbjct: 88  KVILPIFYKVDPSDVRPQSGKFGKAFLKTLKWSGDKEDKIRFWKKALY 135

The Toll/interleukin-1 receptor (TIR) homology domain is an intracellular signalling domain found in MyD88, interleukin 1 receptor and the Toll receptor. It contains three highly-conserved regions, and mediates protein-protein interactions between the Toll-like receptors (TLRs) and signal-transduction components. TIR-like motifs are also found in plant proteins thought to be involved in resistance to disease. When activated, TIR domains recruit cytoplasmic adaptor proteins MyD88 and TOLLIP (Toll interacting protein). In turn, these associate with various kinases to set off signalling cascades. Length = 135

>gnl|CDD|215633 PLN03210, PLN03210, Resistant to P

 Score = 46.4 bits (110), Expect = 2e-07
 Identities = 23/61 (37%), Positives = 38/61 (62%), Gaps = 3/61 (4%)

Query: 26  QIVISVFYQVDPPDMRNQNGNFGDSFLKLEERLKENTEKLQSWRNALKEFMPSLSGFHSL 85
           Q+VI VFY +DP  +R Q G+FG++F K  +   E+ +    W+ AL + + ++ G+HS 
Sbjct: 99  QLVIPVFYGLDPSHVRKQTGDFGEAFEKTCQNKTEDEKIQ--WKQALTD-VANILGYHSQ 155

Query: 86  N 86
           N
Sbjct: 156 N 156

syringae 6; Provisional. Length = 1153

>gnl|CDD|214587 smart00255, TIR, Toll - interleukin 1 - resistance

 Score = 38.5 bits (90), Expect = 7e-05
 Identities = 22/59 (37%), Positives = 29/59 (49%), Gaps = 2/59 (3%)

Query: 21  CHDPTQIVISVFYQVDPPDMRNQNGNFGDSFLKLEERLKENTEKLQSWRNALKEFMPSL 79
             +    VI +FY+V P D+R Q G F   F K   +  E  EK Q W+ AL   +PS 
Sbjct: 84  LEEGGLRVIPIFYEVIPSDVRKQPGKFRKVFKKNYLKWPE-DEKEQFWKKALY-AVPSK 140

Length = 140

Conserved Domains Detected by HHsearch

Original result of HHsearch against CDD database

ID	Length	Definition	Probability
Query	88
PLN03210	1153	Resistant to P. syringae 6; Provisional	99.97
PLN03194	187	putative disease resistance protein; Provisional	99.95
PF01582	141	TIR: TIR domain; InterPro: IPR000157 In Drosophila	99.56
smart00255	140	TIR Toll - interleukin 1 - resistance.	99.39
PF13676	102	TIR_2: TIR domain; PDB: 3H16_B 3UB4_A 2Y92_A 3UB3_	97.77

>PLN03210 Resistant to P

Probab=99.97  E-value=1.5e-32  Score=230.52  Aligned_cols=84  Identities=30%  Similarity=0.615  Sum_probs=78.4

Q ss_pred             EEEEecCCCCCchhhHHHH-----hhc-CCCeEEEEeeecCCCcccccccchHHHHHHHHHhhccChHHHHHHHHHHHhh
Q 040703            2 LYTLPFRDSISLLSCIDTI-----CHD-PTQIVISVFYQVDPPDMRNQNGNFGDSFLKLEERLKENTEKLQSWRNALKEF   75 (88)
Q Consensus         2 ~~vifS~~Ya~S~wCLdEL-----~~~-~~~~ViPVFY~V~psdVr~q~g~fg~af~~~~~~~~~~~e~v~~Wr~AL~~v   75 (88)
                      -+||||+|||+|+||||||     |++ .+++|+||||+|||||||+|+|.||++|.+++++.  +.+++++||+||++|
T Consensus        69 ~ivv~s~~ya~s~wcl~el~~i~~~~~~~~~~v~pvfy~v~p~~v~~~~g~f~~~f~~~~~~~--~~~~~~~w~~al~~~  146 (1153)
T PLN03210         69 AVVVFSKNYASSSWCLNELLEIVRCKEELGQLVIPVFYGLDPSHVRKQTGDFGEAFEKTCQNK--TEDEKIQWKQALTDV  146 (1153)
T ss_pred             EEEEecCCcccchHHHHHHHHHHHhhhhcCceEEEEEecccHHHHhhccchHHHHHHHHhccc--chhHHHHHHHHHHHH
Confidence            5899999999999999999     788 89999999999999999999999999999998764  578999999999999


Q ss_pred             cCCceeeeecCCC
Q 040703           76 MPSLSGFHSLNIT   88 (88)
Q Consensus        76 ~~~l~G~~~~~~~   88 (88)
                       ++++|+++++++
T Consensus       147 -~~~~g~~~~~~~  158 (1153)
T PLN03210        147 -ANILGYHSQNWP  158 (1153)
T ss_pred             -hCcCceecCCCC
Confidence             999999998763

syringae 6; Provisional

>PLN03194 putative disease resistance protein; Provisional

Probab=99.95  E-value=6.8e-29  Score=175.55  Aligned_cols=72  Identities=21%  Similarity=0.278  Sum_probs=63.0

Q ss_pred             CEEEEecCCCCCchhhHHHH--hhcCCCeEEEEeeecCCCccccc-ccchHHHHHHHHHhhccChHHHHHHHHHHHhhcC
Q 040703            1 VLYTLPFRDSISLLSCIDTI--CHDPTQIVISVFYQVDPPDMRNQ-NGNFGDSFLKLEERLKENTEKLQSWRNALKEFMP   77 (88)
Q Consensus         1 ~~~vifS~~Ya~S~wCLdEL--~~~~~~~ViPVFY~V~psdVr~q-~g~fg~af~~~~~~~~~~~e~v~~Wr~AL~~v~~   77 (88)
                      +.+||||+|||+|.||||||  .++.+..||||||+|+|+|||+| .|.             .+.+++++||+||++| +
T Consensus        83 i~IvVfS~~Ya~S~WCLdEL~~I~e~~~~ViPIFY~VdPsdVr~q~~~~-------------~~~e~v~~Wr~AL~~v-a  148 (187)
T PLN03194         83 VGVAVFSPRYCESYFCLHELALIMESKKRVIPIFCDVKPSQLRVVDNGT-------------CPDEEIRRFNWALEEA-K  148 (187)
T ss_pred             EEEEEECCCcccchhHHHHHHHHHHcCCEEEEEEecCCHHHhhccccCC-------------CCHHHHHHHHHHHHHH-h
Confidence            36899999999999999999  66656789999999999999997 443             1468999999999999 9


Q ss_pred             CceeeeecC
Q 040703           78 SLSGFHSLN   86 (88)
Q Consensus        78 ~l~G~~~~~   86 (88)
                      +++|++++.
T Consensus       149 ~l~G~~~~~  157 (187)
T PLN03194        149 YTVGLTFDS  157 (187)
T ss_pred             ccccccCCC
Confidence            999998754

>PF01582 TIR: TIR domain; InterPro: IPR000157 In Drosophila melanogaster the Toll protein is involved in establishment of dorso-ventral polarity in the embryo

Probab=99.56  E-value=2e-16  Score=105.24  Aligned_cols=72  Identities=21%  Similarity=0.340  Sum_probs=61.3

Q ss_pred             EEEEecCCCCCchhhHHHH-----hhc-C--CCeEEEEeeecCCCccc-ccccchHHHHHHHHHhhccC--hHHHHHHHH
Q 040703            2 LYTLPFRDSISLLSCIDTI-----CHD-P--TQIVISVFYQVDPPDMR-NQNGNFGDSFLKLEERLKEN--TEKLQSWRN   70 (88)
Q Consensus         2 ~~vifS~~Ya~S~wCLdEL-----~~~-~--~~~ViPVFY~V~psdVr-~q~g~fg~af~~~~~~~~~~--~e~v~~Wr~   70 (88)
                      +++|||++|++|.||++||     ++. .  ..+|+||||+|.|++++ .|.+.|+..|..+....+++  .++...|++
T Consensus        58 ~I~VlS~~y~~s~wc~~el~~a~~~~~~~~~~~~Il~v~~~v~~~~~~~~~~~~~~~~~~~~~~w~~~~~~~~~~~fW~~  137 (141)
T PF01582_consen   58 TIVVLSRNYLSSEWCLFELQEALERLLEEGRDKLILPVFYDVSPSDVRPDQSLRFLLRFLTYLRWPDDDSREDRSWFWKK  137 (141)
T ss_dssp             EEEEESHHHHHHTHHHHHHHHHHHHHHCSTCTTEEEEESSSS-CHHCHTHHHHHHHHHCTHCEETSSSGGGGGHHHHHHH
T ss_pred             eEEEeecccccccchhhhhhhhhhhccccccccceeeEeccCChhhcChhhhHHHHHHhhhheeCCCCCCccHHHHHHHH
Confidence            6899999999999999999     333 3  59999999999999999 79999999998887765543  578999999


Q ss_pred             HHH
Q 040703           71 ALK   73 (88)
Q Consensus        71 AL~   73 (88)
                      ++.
T Consensus       138 l~~  140 (141)
T PF01582_consen  138 LRY  140 (141)
T ss_dssp             HHH
T ss_pred             Hhc
Confidence            975

In addition, members of the Toll family play a key role in innate antibacterial and antifungal immunity in insects as well as in mammals. These proteins are type-I transmembrane receptors that share an intracellular 200 residue domain with the interleukin-1 receptor (IL-1R), the Toll/IL-1R homologous region (TIR). The similarity between Toll-like receptors (LTRs) and IL-1R is not restricted to sequence homology since these proteins also share a similar signalling pathway. They both induce the activation of a Rel type transcription factor via an adaptor protein and a protein kinase []. Interestingly, MyD88, a cytoplasmic adaptor protein found in mammals, contains a TIR domain associated to a DEATH domain (see IPR000488 from INTERPRO) [, , ]. Besides the mammalian and Drosophila melanogaster proteins, a TIR domain is also found in a number of plant proteins implicated in host defence []. As MyD88, these proteins are cytoplasmic. Site directed mutagenesis and deletion analysis have shown that the TIR domain is essential for Toll and IL-1R activities. Sequence analysis have revealed the presence of three highly conserved regions among the different members of the family: box 1 (FDAFISY), box 2 (GYKLC-RD-PG), and box 3 (a conserved W surrounded by basic residues). It has been proposed that boxes 1 and 2 are involved in the binding of proteins involved in signalling, whereas box 3 is primarily involved in directing localization of receptor, perhaps through interactions with cytoskeletal elements [].; GO: 0005515 protein binding, 0007165 signal transduction, 0005622 intracellular; PDB: 3J0A_A 2J67_B 3JRN_A 1FYV_A 1O77_D 1FYX_A 1FYW_A 3OZI_B 1T3G_B 2JS7_A ....

>smart00255 TIR Toll - interleukin 1 - resistance

Probab=99.39  E-value=1.8e-12  Score=84.70  Aligned_cols=76  Identities=26%  Similarity=0.439  Sum_probs=62.9

Q ss_pred             CEEEEecCCCCCchhhHHHH-----hh-c-CCCeEEEEeeecCCCcccccccchHHHHHHHHHhhccChHHHHHHHHHHH
Q 040703            1 VLYTLPFRDSISLLSCIDTI-----CH-D-PTQIVISVFYQVDPPDMRNQNGNFGDSFLKLEERLKENTEKLQSWRNALK   73 (88)
Q Consensus         1 ~~~vifS~~Ya~S~wCLdEL-----~~-~-~~~~ViPVFY~V~psdVr~q~g~fg~af~~~~~~~~~~~e~v~~Wr~AL~   73 (88)
                      .+++|+|++|..|.||..|+     ++ + .+..||||+|+..|+++..+.+.++..+..+..+..++..+ +.|++++.
T Consensus        57 ~~i~vlS~~~~~S~w~~~E~~~a~~~~~~~~~~~iIPI~~~~~~~~~~~~~~~l~~~~~~~~~~w~~~~~~-~fW~~~~~  135 (140)
T smart00255       57 IAIVVLSPNYAESEWCLDELVAALENALEEGGLRVIPIFYEVIPSDVRKQPGKFRKVLKKNYLKWPEDEKE-RFWKKALY  135 (140)
T ss_pred             EEEEEECcccccChhHHHHHHHHHHHHHHcCCCeEEEEEEecChHHHHhcccHHHHHHHHHHhhcCCchhH-HHHHHHHH
Confidence            36899999999999999999     22 3 56899999999999999999999999998875555444443 78999998


Q ss_pred             hhcCC
Q 040703           74 EFMPS   78 (88)
Q Consensus        74 ~v~~~   78 (88)
                      .+ ++
T Consensus       136 ~l-~~  139 (140)
T smart00255      136 AV-PS  139 (140)
T ss_pred             Hh-cc
Confidence            87 53

>PF13676 TIR_2: TIR domain; PDB: 3H16_B 3UB4_A 2Y92_A 3UB3_A 3UB2_A

Probab=97.77  E-value=8.8e-06  Score=50.72  Aligned_cols=32  Identities=13%  Similarity=0.013  Sum_probs=27.8

Q ss_pred             EEEEecCCCCCchhhHHHH-hhc-CCCeEEEEee
Q 040703            2 LYTLPFRDSISLLSCIDTI-CHD-PTQIVISVFY   33 (88)
Q Consensus         2 ~~vifS~~Ya~S~wCLdEL-~~~-~~~~ViPVFY   33 (88)
                      +++++|++|..|.||..|+ ..+ .+..|+||..
T Consensus        53 ~i~~~S~~~~~s~~~~~E~~~a~~~~~~iipv~~   86 (102)
T PF13676_consen   53 VIVLLSPNYLKSPWCRFELGAAWKRGKPIIPVRL   86 (102)
T ss_dssp             EEEEEEHHHHCTHHHHHHHHHHHCTSESEEEEEC
T ss_pred             EEEEECcccccChHHHHHHHHHHHCCCEEEEEEE
Confidence            6889999999999999999 444 7779999983

---

Homologous Structure Templates

Structure Templates Detected by BLAST

Original result of BLAST against Protein Data Bank

ID	Length	Definition	E-value
Query	88
3ozi_A	204	Crystal Structure Of The Tir Domain From The Flax D		4e-07

>pdb|3OZI|A Chain A, Crystal Structure Of The Tir Domain From The Flax Disease Resistance Protein L6 Length = 204

Iteration: 1
 Score = 49.7 bits (117), Expect = 4e-07,   Method: Compositional matrix adjust.
 Identities = 23/61 (37%), Positives = 39/61 (63%), Gaps = 3/61 (4%)
Query: 23  DPTQIVISVFYQVDPPDMRNQNGNFGDSFLKLEERLKENTEKLQSWRNALKEFMPSLSGF 82
           DP +I++ +FY VDP D+R+Q G +  +F K   +    T  +Q+W++ALK+ +  L G+
Sbjct: 121 DPRRIILPIFYMVDPSDVRHQTGCYKKAFRKHANKFDGQT--IQNWKDALKK-VGDLKGW 177
Query: 83  H 83
           H
Sbjct: 178 H 178

Structure Templates Detected by RPS-BLAST

Original result of RPS-BLAST against PDB70 database

ID	Length	Definition	E-value
Query	88
3jrn_A	176	AT1G72930 protein; TIR domain arabidopsis thaliana	4e-20
3ozi_A	204	L6TR; plant TIR domain, plant protein; 2.30A {Linu	9e-20
3h16_A	154	TIR protein; bacteria TIR domain, signaling protei	5e-11

>3jrn_A AT1G72930 protein; TIR domain arabidopsis thaliana, plant protein; 2.00A {Arabidopsis thaliana} Length = 176

 Score = 78.0 bits (193), Expect = 4e-20
 Identities = 22/59 (37%), Positives = 29/59 (49%), Gaps = 3/59 (5%)

Query: 26  QIVISVFYQVDPPDMRNQNGNFGDSFLKLEERLKENTEKLQSWRNALKEFMPSLSGFHS 84
             V+ +FY V+P  +R Q G   + F K   R  E+ EK+  WR AL      LSG  S
Sbjct: 96  ITVMPIFYGVEPNHVRWQTGVLAEQFKKHASR--EDPEKVLKWRQALTN-FAQLSGDCS 151

>3ozi_A L6TR; plant TIR domain, plant protein; 2.30A {Linum usitatissimum} Length = 204

 Score = 78.1 bits (193), Expect = 9e-20
 Identities = 20/61 (32%), Positives = 37/61 (60%), Gaps = 3/61 (4%)

Query: 26  QIVISVFYQVDPPDMRNQNGNFGDSFLKLEERLKENTEKLQSWRNALKEFMPSLSGFHSL 85
           +I++ +FY VDP D+R+Q G +  +F K   +     + +Q+W++ALK+ +  L G+H  
Sbjct: 124 RIILPIFYMVDPSDVRHQTGCYKKAFRKHANKFDG--QTIQNWKDALKK-VGDLKGWHIG 180

Query: 86  N 86
            
Sbjct: 181 K 181

>3h16_A TIR protein; bacteria TIR domain, signaling protein; 2.50A {Paracoccus denitrificans PD1222} Length = 154

 Score = 54.2 bits (131), Expect = 5e-11
 Identities = 3/47 (6%), Positives = 20/47 (42%), Gaps = 2/47 (4%)

Query: 26  QIVISVFYQVDPPDMRNQNGNFGDSFLKLEERLKENTEKLQSWRNAL 72
             ++ ++++V   ++ + +    D          ++ +++ +   A+
Sbjct: 107 SRILPIWHKVSKDEVASFSPTMADKLAFNTST--KSVDEIVADLMAI 151

Structure Templates Detected by HHsearch

Original result of HHsearch against PDB70 database

ID	Length	Definition	Probability
Query	88
3jrn_A	176	AT1G72930 protein; TIR domain arabidopsis thaliana	100.0
3ozi_A	204	L6TR; plant TIR domain, plant protein; 2.30A {Linu	100.0
3h16_A	154	TIR protein; bacteria TIR domain, signaling protei	99.73
3ub2_A	146	TOLL/interleukin-1 receptor domain-containing ADA	99.29
1fyx_A	149	TOLL-like receptor 2; beta-alpha-beta fold, signal	98.69
2js7_A	160	Myeloid differentiation primary response protein M	98.67
2j67_A	178	TOLL like receptor 10; TIR, IL-1, TLR10, membrane,	98.66
1t3g_A	159	X-linked interleukin-1 receptor accessory protein-	98.44
3j0a_A	844	TOLL-like receptor 5; membrane protein, leucine-ri	97.65
1eiw_A	111	Hypothetical protein MTH538; CHEY-like fold, flavo	91.17

>3jrn_A AT1G72930 protein; TIR domain arabidopsis thaliana, plant protein; 2.00A {Arabidopsis thaliana}

Probab=100.00  E-value=4.1e-37  Score=214.60  Aligned_cols=81  Identities=31%  Similarity=0.507  Sum_probs=77.1

Q ss_pred             EEEEecCCCCCchhhHHHH-----hhc-CCCeEEEEeeecCCCcccccccchHHHHHHHHHhhccChHHHHHHHHHHHhh
Q 040703            2 LYTLPFRDSISLLSCIDTI-----CHD-PTQIVISVFYQVDPPDMRNQNGNFGDSFLKLEERLKENTEKLQSWRNALKEF   75 (88)
Q Consensus         2 ~~vifS~~Ya~S~wCLdEL-----~~~-~~~~ViPVFY~V~psdVr~q~g~fg~af~~~~~~~~~~~e~v~~Wr~AL~~v   75 (88)
                      .+||||+||++|.|||+||     |++ .+++||||||+|+|||||+|+|.||++|.+|+++  .+.+++++||.||++|
T Consensus        66 ~IvV~S~ny~~S~WCl~EL~~i~~~~~~~~~~ViPIfy~V~ps~Vr~q~g~fg~af~~~~~~--~~~~~~~~Wr~AL~~v  143 (176)
T 3jrn_A           66 AVVVVSENYAASSWCLDELVTIMDFEKKGSITVMPIFYGVEPNHVRWQTGVLAEQFKKHASR--EDPEKVLKWRQALTNF  143 (176)
T ss_dssp             EEEEECTTTTTCHHHHHHHHHHHHHHHTTSCEEEEEECSSCHHHHHHTCTHHHHHHHHHHTT--SCHHHHHHHHHHHHHH
T ss_pred             EEEEecCCcCCChhHHHHHHHHHhhhccCCCEEEEEEecCCHHHhhhccCcHHHHHHHHHhc--cCHHHHHHHHHHHHHH
Confidence            6899999999999999999     777 8999999999999999999999999999999987  4789999999999999


Q ss_pred             cCCceeeeec
Q 040703           76 MPSLSGFHSL   85 (88)
Q Consensus        76 ~~~l~G~~~~   85 (88)
                       |+++||++.
T Consensus       144 -a~~~G~~~~  152 (176)
T 3jrn_A          144 -AQLSGDCSG  152 (176)
T ss_dssp             -TTSCCEECC
T ss_pred             -hcccceecC
Confidence             999999984

>3ozi_A L6TR; plant TIR domain, plant protein; 2.30A {Linum usitatissimum}

Probab=100.00  E-value=9.2e-37  Score=217.16  Aligned_cols=84  Identities=26%  Similarity=0.578  Sum_probs=78.2

Q ss_pred             CEEEEecCCCCCchhhHHHH-----hhc--CCCeEEEEeeecCCCcccccccchHHHHHHHHHhhccChHHHHHHHHHHH
Q 040703            1 VLYTLPFRDSISLLSCIDTI-----CHD--PTQIVISVFYQVDPPDMRNQNGNFGDSFLKLEERLKENTEKLQSWRNALK   73 (88)
Q Consensus         1 ~~~vifS~~Ya~S~wCLdEL-----~~~--~~~~ViPVFY~V~psdVr~q~g~fg~af~~~~~~~~~~~e~v~~Wr~AL~   73 (88)
                      +.+||||+||++|.|||+||     |++  .+++||||||+|+||+||+|+|.||++|.+|++++  +.+++++||.||+
T Consensus        92 i~IvV~S~nYa~S~WCl~EL~~I~e~~~~~~~~~ViPIFY~VdPs~Vr~q~g~fg~af~~~~~~~--~~~~v~~Wr~AL~  169 (204)
T 3ozi_A           92 IYVPIISSGYADSKWCLMELAEIVRRQEEDPRRIILPIFYMVDPSDVRHQTGCYKKAFRKHANKF--DGQTIQNWKDALK  169 (204)
T ss_dssp             EEEEEECTTGGGCHHHHHHHHHHHHHHHHCTTSEECCEEESSCHHHHHHTCTTHHHHHHHHTTTS--CHHHHHHHHHHHH
T ss_pred             EeeEEEEcccccCcHHHHHHHHHHHHHHhcCCeeeEEEEeecCHHHHHhccccHHHHHHHHHHhh--CHHHHHHHHHHHH
Confidence            36899999999999999999     664  57999999999999999999999999999999987  4789999999999


Q ss_pred             hhcCCceeeeecCC
Q 040703           74 EFMPSLSGFHSLNI   87 (88)
Q Consensus        74 ~v~~~l~G~~~~~~   87 (88)
                      +| |+++||++.++
T Consensus       170 ~v-a~lsG~~~~~~  182 (204)
T 3ozi_A          170 KV-GDLKGWHIGKN  182 (204)
T ss_dssp             HH-HTSCBEEECTT
T ss_pred             HH-hccCceecCCC
Confidence            99 99999999875

>3h16_A TIR protein; bacteria TIR domain, signaling protein; 2.50A {Paracoccus denitrificans PD1222}

Probab=99.73  E-value=1.4e-19  Score=121.66  Aligned_cols=70  Identities=6%  Similarity=0.001  Sum_probs=55.3

Q ss_pred             CEEEEecCCCCCchhhHHHH-----hhc-CCCeEEEEeeecCCCcccccccchHHHHHHHHHhhccChHHHHHHHHHHHh
Q 040703            1 VLYTLPFRDSISLLSCIDTI-----CHD-PTQIVISVFYQVDPPDMRNQNGNFGDSFLKLEERLKENTEKLQSWRNALKE   74 (88)
Q Consensus         1 ~~~vifS~~Ya~S~wCLdEL-----~~~-~~~~ViPVFY~V~psdVr~q~g~fg~af~~~~~~~~~~~e~v~~Wr~AL~~   74 (88)
                      ++++|+|++|++|.||++||     |++ .++.||||||+|+|++||+|.|.||++|.....     ...+.++-+.|.+
T Consensus        76 ~~i~v~S~~y~~S~wc~~El~~~~~~~~~~~~~iiPV~~~v~p~~v~~~~~~~~~~~~~~~~-----~~~~~~ia~~l~~  150 (154)
T 3h16_A           76 FGIVVLSTHFFKKEWPQKELDGLFQLESSGRSRILPIWHKVSKDEVASFSPTMADKLAFNTS-----TKSVDEIVADLMA  150 (154)
T ss_dssp             EEEEEEEHHHHTTCCCHHHHHHHTCCCTTSCCCEEEEEESCCTGGGTTTCCCCCSSCCEETT-----TSCHHHHHHHHHH
T ss_pred             EEEEEeCcchhcChHHHHHHHHHHHHHhcCCCEEEEEEecCCHHHHhhCCccHHHHHhhhcC-----cccHHHHHHHHHH
Confidence            36899999999999999999     445 678999999999999999999999998764332     2334555555554


Q ss_pred             h
Q 040703           75 F   75 (88)
Q Consensus        75 v   75 (88)
                      +
T Consensus       151 l  151 (154)
T 3h16_A          151 I  151 (154)
T ss_dssp             H
T ss_pred             H
Confidence            4

>3ub2_A TOLL/interleukin-1 receptor domain-containing ADA protein; TIR domain, TLRS adaptor, immune system; 2.40A {Homo sapiens} PDB: 3ub3_A 3ub4_A 2y92_A

Probab=99.29  E-value=2.3e-13  Score=91.30  Aligned_cols=68  Identities=13%  Similarity=0.109  Sum_probs=47.7

Q ss_pred             EEEEecCCCCCchhhHHHH-----hh-cCCCeEEEEeeecC----CCcccccccchHHHHHHHHHhhccChHHHHHHHHH
Q 040703            2 LYTLPFRDSISLLSCIDTI-----CH-DPTQIVISVFYQVD----PPDMRNQNGNFGDSFLKLEERLKENTEKLQSWRNA   71 (88)
Q Consensus         2 ~~vifS~~Ya~S~wCLdEL-----~~-~~~~~ViPVFY~V~----psdVr~q~g~fg~af~~~~~~~~~~~e~v~~Wr~A   71 (88)
                      +++|+|+||++|.||+.||     +. .....||||||+|+    |..+|.......   ..++..+....+.|.+|++|
T Consensus        68 ~I~VlS~~y~~S~wc~~El~~al~~~~~~~~~vIpv~~~v~~~~lp~~Lr~~~~id~---~~~d~~f~~l~~~v~~~~~~  144 (146)
T 3ub2_A           68 RVLLITPGFLQDPWCKYQMLQALTEAPGAEGCTIPLLSGLSRAAYPPELRFMYYVDG---RGPDGGFRQVKEAVMRYLQT  144 (146)
T ss_dssp             EEEEECHHHHHCHHHHHHHHHHHHTSSSSSSEEEEEECSCCGGGSCGGGGGSCCEET---TSGGGGHHHHHHHHHHHHTT
T ss_pred             EEEEECcccccCHHHHHHHHHHHHHHhhcCCcEEEEEcCCChhhCCHHHhCeeeeec---cChHhhHHHHHHHHHHHHHh
Confidence            6899999999999999999     32 22236779999999    445566543222   12334444456889999988


Q ss_pred             H
Q 040703           72 L   72 (88)
Q Consensus        72 L   72 (88)
                      |
T Consensus       145 ~  145 (146)
T 3ub2_A          145 L  145 (146)
T ss_dssp             C
T ss_pred             c
Confidence            6

>1fyx_A TOLL-like receptor 2; beta-alpha-beta fold, signaling protein; 2.80A {Homo sapiens} SCOP: c.23.2.1 PDB: 1fyw_A 1o77_A

Probab=98.69  E-value=1.3e-09  Score=72.83  Aligned_cols=39  Identities=5%  Similarity=0.027  Sum_probs=31.9

Q ss_pred             EEEEecCCCCCchhhHHHH------hhc-CCCeEEEEeee-cCCCcc
Q 040703            2 LYTLPFRDSISLLSCIDTI------CHD-PTQIVISVFYQ-VDPPDM   40 (88)
Q Consensus         2 ~~vifS~~Ya~S~wCLdEL------~~~-~~~~ViPVFY~-V~psdV   40 (88)
                      +++|+|++|++|.||+.||      +.+ .++.||||||+ +++..+
T Consensus        64 ~I~VlS~~y~~S~wc~~El~~a~~~~~~~~~~~vIpv~~~~i~~~~~  110 (149)
T 1fyx_A           64 TVFVLSENFVKSEWXKYELDFSHFRLFDENNDAAILILLEPIEKKAI  110 (149)
T ss_dssp             EEEEECHHHHHHHTHHHHSCCSCCTTCGGGTTCCEEEESSCCCTTTS
T ss_pred             EEEEeCcchhccchHHHHHHHHHHHHHhcCCCEEEEEEecCCChhhc
Confidence            6899999999999999999      335 67889999996 555443

>2js7_A Myeloid differentiation primary response protein MYD88; MYD88_human, TIR domain, TOLL like receptor adaptor domain, innate immune signaling; NMR {Homo sapiens} PDB: 2z5v_A

Probab=98.67  E-value=1.5e-09  Score=73.43  Aligned_cols=36  Identities=11%  Similarity=0.106  Sum_probs=30.2

Q ss_pred             EEEEecCCCCCchhhHHHH-----hh-c-CCCeEEEEeeecCC
Q 040703            2 LYTLPFRDSISLLSCIDTI-----CH-D-PTQIVISVFYQVDP   37 (88)
Q Consensus         2 ~~vifS~~Ya~S~wCLdEL-----~~-~-~~~~ViPVFY~V~p   37 (88)
                      +++|+|++|++|.||+.||     +. + .++.||||||+.-+
T Consensus        75 ~IvVlS~~y~~S~wc~~El~~a~~~~~~~~~~~vIpV~~~~~~  117 (160)
T 2js7_A           75 MVVVVSDDYLQSKECDFQTKFALSLSPGAHQKRLIPIKYKAMK  117 (160)
T ss_dssp             EEEECCHHHHHSHHHHHHHHHHHHHCTTHHHHTEEEEESSCCC
T ss_pred             EEEEECcchhcCHHHHHHHHHHHHHHHccCCCEEEEEEEcccc
Confidence            6899999999999999999     22 3 46789999998653

>2j67_A TOLL like receptor 10; TIR, IL-1, TLR10, membrane, innate immunity, immune response, leucine-rich repeat, glycoprotein, transmembrane; 2.20A {Homo sapiens} PDB: 1fyv_A

Probab=98.66  E-value=2e-09  Score=74.24  Aligned_cols=34  Identities=6%  Similarity=0.001  Sum_probs=26.2

Q ss_pred             EEEEecCCCCCchhhHHHH------hhc-CCCeEEEEeeec
Q 040703            2 LYTLPFRDSISLLSCIDTI------CHD-PTQIVISVFYQV   35 (88)
Q Consensus         2 ~~vifS~~Ya~S~wCLdEL------~~~-~~~~ViPVFY~V   35 (88)
                      +++|+|+||++|.||+.||      +.+ .++.||||||+-
T Consensus        94 ~IvVlS~~yl~S~wc~~El~~a~~~~~~~~~~~vIpV~~~~  134 (178)
T 2j67_A           94 SIFVLSPNFVQNEWCHYEFYFAHHNLFHENSDHIILILLEP  134 (178)
T ss_dssp             EEEEECHHHHHHTGGGTHHHHTTCC-------CEEEEESSC
T ss_pred             EEEEecccccccchHHHHHHHHHHHHHhcCCCEEEEEEecC
Confidence            6899999999999999999      234 578899999973

>1t3g_A X-linked interleukin-1 receptor accessory protein-like 1; TIR, IL-1RAPL, IL-1R, TLR, membrane protein; 2.30A {Homo sapiens}

Probab=98.44  E-value=3.7e-08  Score=66.52  Aligned_cols=37  Identities=11%  Similarity=-0.189  Sum_probs=31.2

Q ss_pred             EEEEecCCCC-CchhhHHHH-----hh-c-CCCeEEEEeeecCCC
Q 040703            2 LYTLPFRDSI-SLLSCIDTI-----CH-D-PTQIVISVFYQVDPP   38 (88)
Q Consensus         2 ~~vifS~~Ya-~S~wCLdEL-----~~-~-~~~~ViPVFY~V~ps   38 (88)
                      +++|+|+||. +|.||+.||     ++ + .++.||||||.-.|.
T Consensus        71 ~IvVlS~~y~~~S~wc~~El~~a~~~~~~~~~~~vI~I~~~~~~~  115 (159)
T 1t3g_A           71 LIIVMTPNYVVRRGWSIFELETRLRNMLVTGEIKVILIECSELRG  115 (159)
T ss_dssp             EEEEECHHHHHTTTTHHHHHSHHHHHHHHTTSSEEEEEECSCCCS
T ss_pred             EEEEEccchhhcChHHHHHHHHHHHHHHhcCCCEEEEEEeccccc
Confidence            5899999996 999999999     33 4 678999999986654

>3j0a_A TOLL-like receptor 5; membrane protein, leucine-rich repeat, asymmetric homodimer, glycoprotein, immune system; HET: NAG FUC; 26.00A {Homo sapiens}

Probab=97.65  E-value=1.1e-05  Score=64.74  Aligned_cols=37  Identities=16%  Similarity=0.192  Sum_probs=31.0

Q ss_pred             EEEEecCCCCCchhhHHHH------hhc-CCCeEEEEeeecCCC
Q 040703            2 LYTLPFRDSISLLSCIDTI------CHD-PTQIVISVFYQVDPP   38 (88)
Q Consensus         2 ~~vifS~~Ya~S~wCLdEL------~~~-~~~~ViPVFY~V~ps   38 (88)
                      +++|+|++|+.|.||..|+      +.+ .+..||||||+--|.
T Consensus       732 ~i~vls~~~~~s~wc~~e~~~a~~~~~~~~~~~~i~i~~~~~~~  775 (844)
T 3j0a_A          732 IVCLVSRHFLRDGWCLEAFSYAQGRCLSDLNSALIMVVVGSLSQ  775 (844)
T ss_dssp             EEEEECTTHHHHTSTTHHHHHHHSCCCCSSCTTEEEEESSCCCS
T ss_pred             EEEEeccccccChHHHHHHHHHHHHHHHhcCCcEEEEEeccCCh
Confidence            5889999999999999999      224 678999999985553

>1eiw_A Hypothetical protein MTH538; CHEY-like fold, flavodoxin-like fold, (A/B)5 doubly wound fold, parallel beta sheet; NMR {Methanothermobacterthermautotrophicus} SCOP: c.23.3.1

Probab=91.17  E-value=0.2  Score=31.75  Aligned_cols=40  Identities=15%  Similarity=0.262  Sum_probs=28.5

Q ss_pred             EEEEecCCCCCchhhHHHH-hhc-CCCeEEEEe-eecC--CCccc
Q 040703            2 LYTLPFRDSISLLSCIDTI-CHD-PTQIVISVF-YQVD--PPDMR   41 (88)
Q Consensus         2 ~~vifS~~Ya~S~wCLdEL-~~~-~~~~ViPVF-Y~V~--psdVr   41 (88)
                      ++|+.+++-++|.||-.|+ +.+ .|..++=|. |+.+  |+.+.
T Consensus        42 vIvL~G~~t~~s~wv~~EI~~A~~~gkpIigV~~~g~~~~P~~l~   86 (111)
T 1eiw_A           42 VIVLAGLWGTRRDEILGAVDLARKSSKPIITVRPYGLENVPPELE   86 (111)
T ss_dssp             EEEEGGGTTTSHHHHHHHHHHHTTTTCCEEEECCSSSSCCCTTHH
T ss_pred             EEEEeCCCcCCChHHHHHHHHHHHcCCCEEEEEcCCCCcCCHHHH
Confidence            4555666667999999999 666 777777776 5544  55554

---

Homologous Structure Domains

Structure Domains Detected by RPS-BLAST

Original result of RPS-BLAST against SCOP70(version1.75) database

No hit with e-value below 0.005

Homologous Domains Detected by HHsearch

Original result of HHsearch against SCOP70(version1.75) database

ID	Length	Definition	Probability
Query	88
d1fyxa_	149	Toll-like receptor 2, TLR2 {Human (Homo sapiens) [	97.99
d1fyva_	161	Toll-like receptor 1, TLR1 {Human (Homo sapiens) [	97.61

>d1fyxa_ c.23.2.1 (A:) Toll-like receptor 2, TLR2 {Human (Homo sapiens) [TaxId: 9606]}

class: Alpha and beta proteins (a/b)
fold: Flavodoxin-like
superfamily: Toll/Interleukin receptor TIR domain
family: Toll/Interleukin receptor TIR domain
domain: Toll-like receptor 2, TLR2
species: Human (Homo sapiens) [TaxId: 9606]

Probab=97.99  E-value=1.3e-06  Score=56.24  Aligned_cols=38  Identities=8%  Similarity=0.016  Sum_probs=30.8

Q ss_pred             EEEEecCCCCCchhhHHHH------hhc-CCCeEEEEeeecCCCc
Q 040703            2 LYTLPFRDSISLLSCIDTI------CHD-PTQIVISVFYQVDPPD   39 (88)
Q Consensus         2 ~~vifS~~Ya~S~wCLdEL------~~~-~~~~ViPVFY~V~psd   39 (88)
                      +++|+|++|..|.||..|+      ..+ ++..+|||+++--|.+
T Consensus        64 ~i~vlS~~~l~s~wc~~E~~~a~~~~~~~~~~~iIpV~l~~~~~~  108 (149)
T d1fyxa_          64 TVFVLSENFVKSEWCKYELDFSHFRLFDENNDAAILILLEPIEKK  108 (149)
T ss_dssp             EEEEECHHHHHHHTHHHHSCCSCCTTCGGGTTCCEEEESSCCCTT
T ss_pred             EEEEEcCccccccchHHHHHHHHHHHHHcCCceEEEEEeccCchh
Confidence            6899999999999999999      224 6788999999744433

>d1fyva_ c.23.2.1 (A:) Toll-like receptor 1, TLR1 {Human (Homo sapiens) [TaxId: 9606]}

class: Alpha and beta proteins (a/b)
fold: Flavodoxin-like
superfamily: Toll/Interleukin receptor TIR domain
family: Toll/Interleukin receptor TIR domain
domain: Toll-like receptor 1, TLR1
species: Human (Homo sapiens) [TaxId: 9606]

Probab=97.61  E-value=4.2e-06  Score=54.39  Aligned_cols=35  Identities=6%  Similarity=-0.007  Sum_probs=29.3

Q ss_pred             EEEEecCCCCCchhhHHHH--h----hc-CCCeEEEEeeecC
Q 040703            2 LYTLPFRDSISLLSCIDTI--C----HD-PTQIVISVFYQVD   36 (88)
Q Consensus         2 ~~vifS~~Ya~S~wCLdEL--~----~~-~~~~ViPVFY~V~   36 (88)
                      +++|+|++|..|.||..|+  .    .+ ++..+|||+++--
T Consensus        69 ~i~vlS~~~~~s~w~~~E~~~a~~~~~~~~~~~lIpV~l~~~  110 (161)
T d1fyva_          69 SIFVLSPNFVQSEWCHYELYFAHHNLFHEGSNSLILILLEPI  110 (161)
T ss_dssp             EEEEESHHHHHHTSHHHHHHTTSCCCSCSSCSSEEEEESSCC
T ss_pred             EEEEecCccccChHHHHHHHHHHHHHHHcCCCceeEEEEecC
Confidence            6899999999999999999  2    23 5678999999743