Psyllid ID: psy11196

---

Local Sequence Feature Prediction

Residue Number Marker
--------10--------20--------30--------40--------50--------60--------70--------80--------90-------100---
Protein Sequence
MSAYKSGAIVRRKELKMLPREHIYSQYQGVWNLSSDQGNLGVFMATNIRVVCTYYFGEVMDLGFDIRTYYFGEVMDLGFDIRYESFENEVYHELFDSHHVFRR
Secondary Structure (PSIPRED)
ccccccccEEEccEEEcccccEEEEEEccEEEcccccccEEEEEEEEEEEEEEEcccccccEEEEEEEcccccEEEEEEEccHHHHHHHHHHHHHHHHHHHcc
Secondary Structure Prediction (SSPRO)
ccccccccEEEcccEEEccccEEEEEcccEEEEcccccccEEEEEEEEEEEEEHHHHHHHcccccccccccccEEEEEEEEccHHHHHHHHHHHHHHHHHHcc
Coil and Loop (DISEMBL)
MSAYKSGAIVRRKElkmlprehiYSQYQGvwnlssdqgnlgvfMATNIRVVCTYYFGEVMDLGFDIRTYYFGEVMDLGFDIRYESFENEVYHelfdshhvfrr
Flexible Loop (DISEMBL)
msayksgaivrrkelkmlprEHIYSQYQGVWNLSSDQGNLGVFMATNIRVVCTYYFGEVMDLGFDIRTYYFGEVMDLGFDIRYESFENEVYhelfdshhvfrr
Low Complexity Region (SEG)
MSAYKSGAIVRRKELKMLPREHIYSQYQGVWNLSSDQGNLGVFMATNIRVVCTYYFGEVMDLGFDIRTYYFGEVMDLGFDIRYESFENEVYHELFDSHHVFRR
Disordered region (IsUnstruct)
*************ELKMLPREHIYSQYQGVWNLSSDQGNLGVFMATNIRVVCTYYFGEVMDLGFDIRTYYFGEVMDLGFDIRYESFENEVYHELFDSHHV***
Disordered Region (DISOPRED)
******G***RRKELKMLPREHIYSQYQGVWNLSSDQGNLGVFMATNIRVVCTYYFGEVMDLGFDIRTYYFGEVMDLGFDIRYESFENEVYHELFDSHHVFR*
Disordered Region (DISEMBL)
MSAYKSGAIVRRKELKMLPREHIYSQYQGVWNLSSDQGNLGVFMATNIRVVCTYYFGEVMDLGFDIRTYYFGEVMDLGFDIRYESFENEVYHELFDSHHVFRR
Disordered Region (DISPRO)
****KSGAIVRRKELKMLPREHIYSQYQGVWNLSSDQGNLGVFMATNIRVVCTYYFGEVMDLGFDIRTYYFGEVMDLGFDIRYESFENEVYHELFDSHHVFRR
Transmembrane Helix (TMHMM)
ooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo
Transmembrane Helix (HMMTOP)
iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiHHHHHHHHHHHHHHHHHHHHoooooooooooooooooooooooooooooooooooooooooooooo
Transmembrane Helix (MEMSAT)
iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiihhhhhhhhhhhhhhhhhhhhoooooooooooooooooooooooooooooooooooooooooooooo
TM Helix, Signal Peptide (MEMSAT_SVM)
iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiihhhhhhhhhhhhhhhhooooooooooooooooooooooooooooooooooooooooooooooo
TM Helix, Signal Peptide (Phobius)
ooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo
Signal Peptide (SignalP HMM Mode)
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
Signal Peptide (SignalP NN Mode)
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
Coiled Coils (COILS)
MSAYKSGAIVRRKELKMLPREHIYSQYQGVWNLSSDQGNLGVFMATNIRVVCTYYFGEVMDLGFDIRTYYFGEVMDLGFDIRYESFENEVYHELFDSHHVFRR
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	103	2.2.26 [Sep-21-2011]
Q66IS6	365	Bardet-Biedl syndrome 5 p	N/A	N/A	0.436	0.123	0.622	2e-12
Q9CZQ9	341	Bardet-Biedl syndrome 5 p	yes	N/A	0.427	0.129	0.590	1e-11
Q8N3I7	341	Bardet-Biedl syndrome 5 p	yes	N/A	0.427	0.129	0.590	1e-11
Q7ZWB7	342	Bardet-Biedl syndrome 5 p	yes	N/A	0.436	0.131	0.577	1e-11
Q4R649	341	Bardet-Biedl syndrome 5 p	N/A	N/A	0.427	0.129	0.590	1e-11

>sp|Q66IS6|BBS5_XENLA Bardet-Biedl syndrome 5 protein homolog OS=Xenopus laevis GN=bbs5 PE=2 SV=1

 Score = 70.9 bits (172), Expect = 2e-12,   Method: Compositional matrix adjust.
 Identities = 28/45 (62%), Positives = 37/45 (82%)

Query: 7   GAIVRRKELKMLPREHIYSQYQGVWNLSSDQGNLGVFMATNIRVV 51
           GA+++ K+LK+LPRE +Y +  GVWNLSSDQGNLG F  TN+R+V
Sbjct: 145 GALIQNKQLKLLPREQVYDKINGVWNLSSDQGNLGTFFITNVRIV 189

Required for ciliogenesis.
Xenopus laevis (taxid: 8355)

>sp|Q9CZQ9|BBS5_MOUSE Bardet-Biedl syndrome 5 protein homolog OS=Mus musculus GN=Bbs5 PE=2 SV=1

 Score = 68.6 bits (166), Expect = 1e-11,   Method: Composition-based stats.
 Identities = 26/44 (59%), Positives = 37/44 (84%)

Query: 8   AIVRRKELKMLPREHIYSQYQGVWNLSSDQGNLGVFMATNIRVV 51
           A+++ K+L++LP+EH+Y +  GVWNLSSDQGNLG F  TN+R+V
Sbjct: 145 AVIQNKQLRLLPQEHVYDKINGVWNLSSDQGNLGTFFITNVRIV 188

The BBSome complex is required for ciliogenesis but is dispensable for centriolar satellite function. This ciliogenic function is mediated in part by the Rab8 GDP/GTP exchange factor, which localizes to the basal body and contacts the BBSome. Rab8(GTP) enters the primary cilium and promotes extension of the ciliary membrane. Firstly the BBSome associates with the ciliary membrane and binds to RAB3IP/Rabin8, the guanosyl exchange factor (GEF) for Rab8 and then the Rab8-GTP localizes to the cilium and promotes docking and fusion of carrier vesicles to the base of the ciliary membrane.

Mus musculus (taxid: 10090)

>sp|Q8N3I7|BBS5_HUMAN Bardet-Biedl syndrome 5 protein OS=Homo sapiens GN=BBS5 PE=1 SV=1

 Score = 68.2 bits (165), Expect = 1e-11,   Method: Composition-based stats.
 Identities = 26/44 (59%), Positives = 37/44 (84%)

Query: 8   AIVRRKELKMLPREHIYSQYQGVWNLSSDQGNLGVFMATNIRVV 51
           A+++ K+L++LP+EH+Y +  GVWNLSSDQGNLG F  TN+R+V
Sbjct: 145 ALIQNKQLRLLPQEHVYDKINGVWNLSSDQGNLGTFFITNVRIV 188

The BBSome complex is required for ciliogenesis but is dispensable for centriolar satellite function. This ciliogenic function is mediated in part by the Rab8 GDP/GTP exchange factor, which localizes to the basal body and contacts the BBSome. Rab8(GTP) enters the primary cilium and promotes extension of the ciliary membrane. Firstly the BBSome associates with the ciliary membrane and binds to RAB3IP/Rabin8, the guanosyl exchange factor (GEF) for Rab8 and then the Rab8-GTP localizes to the cilium and promotes docking and fusion of carrier vesicles to the base of the ciliary membrane.

Homo sapiens (taxid: 9606)

>sp|Q7ZWB7|BBS5_DANRE Bardet-Biedl syndrome 5 protein homolog OS=Danio rerio GN=bbs5 PE=2 SV=1

 Score = 68.2 bits (165), Expect = 1e-11,   Method: Composition-based stats.
 Identities = 26/45 (57%), Positives = 37/45 (82%)

Query: 7   GAIVRRKELKMLPREHIYSQYQGVWNLSSDQGNLGVFMATNIRVV 51
           GA+++ K+L++LP+E +Y +  GVWNLSSDQGNLG F  TN+R+V
Sbjct: 145 GALIQNKQLRLLPQEQVYDKINGVWNLSSDQGNLGTFFITNVRIV 189

Required for ciliogenesis.

Danio rerio (taxid: 7955)

>sp|Q4R649|BBS5_MACFA Bardet-Biedl syndrome 5 protein homolog OS=Macaca fascicularis GN=BBS5 PE=2 SV=2

 Score = 68.2 bits (165), Expect = 1e-11,   Method: Composition-based stats.
 Identities = 26/44 (59%), Positives = 37/44 (84%)

Query: 8   AIVRRKELKMLPREHIYSQYQGVWNLSSDQGNLGVFMATNIRVV 51
           A+++ K+L++LP+EH+Y +  GVWNLSSDQGNLG F  TN+R+V
Sbjct: 145 ALIQNKQLRLLPQEHVYDKINGVWNLSSDQGNLGTFFITNVRIV 188


  Database: swissprot
    Posted date:  Mar 23, 2013  2:32 AM
  Number of letters in database: 191,569,459
  Number of sequences in database:  539,616
  
Lambda     K      H
   0.327    0.142    0.437 

Lambda     K      H
   0.267   0.0410    0.140 


Matrix: BLOSUM62
Gap Penalties: Existence: 11, Extension: 1
Number of Hits to DB: 39,410,030
Number of Sequences: 539616
Number of extensions: 1488575
Number of successful extensions: 3389
Number of sequences better than 100.0: 7
Number of HSP's better than 100.0 without gapping: 6
Number of HSP's successfully gapped in prelim test: 1
Number of HSP's that attempted gapping in prelim test: 3379
Number of HSP's gapped (non-prelim): 11
length of query: 103
length of database: 191,569,459
effective HSP length: 72
effective length of query: 31
effective length of database: 152,717,107
effective search space: 4734230317
effective search space used: 4734230317
T: 11
A: 40
X1: 15 ( 7.1 bits)
X2: 38 (14.6 bits)
X3: 64 (24.7 bits)
S1: 40 (21.7 bits)
S2: 55 (25.8 bits)

The BBSome complex is required for ciliogenesis but is dispensable for centriolar satellite function. This ciliogenic function is mediated in part by the Rab8 GDP/GTP exchange factor, which localizes to the basal body and contacts the BBSome. Rab8(GTP) enters the primary cilium and promotes extension of the ciliary membrane. Firstly the BBSome associates with the ciliary membrane and binds to RAB3IP/Rabin8, the guanosyl exchange factor (GEF) for Rab8 and then the Rab8-GTP localizes to the cilium and promotes docking and fusion of carrier vesicles to the base of the ciliary membrane.

Macaca fascicularis (taxid: 9541)

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	103
348672493	348	hypothetical protein PHYSODRAFT_515023 [	0.922	0.272	0.368	2e-11
255087046	345	predicted protein [Micromonas sp. RCC299	0.922	0.275	0.368	2e-11
301112689	348	bardet-Biedl syndrome 5 protein [Phytoph	0.922	0.272	0.368	3e-11
195453986	385	GK12823 [Drosophila willistoni] gi|19417	0.922	0.246	0.377	6e-11
148231484	365	Bardet-Biedl syndrome 5 protein homolog	0.436	0.123	0.622	9e-11
432098391	341	Bardet-Biedl syndrome 5 protein like pro	0.427	0.129	0.613	2e-10
348585675	341	PREDICTED: Bardet-Biedl syndrome 5 prote	0.427	0.129	0.613	2e-10
148695091	395	Bardet-Biedl syndrome 5 (human), isoform	0.427	0.111	0.590	2e-10
213510754	342	Bardet-Biedl syndrome 5 protein homolog	0.436	0.131	0.6	2e-10
325182780	345	bardetBiedl syndrome 5 protein putative	0.941	0.281	0.362	2e-10

>gi|348672493|gb|EGZ12313.1| hypothetical protein PHYSODRAFT_515023 [Phytophthora sojae]

 Score = 72.8 bits (177), Expect = 2e-11,   Method: Composition-based stats.
 Identities = 42/114 (36%), Positives = 58/114 (50%), Gaps = 19/114 (16%)

Query: 7   GAIVRRKELKMLPREHIYSQYQGVWNLSSDQGNLGVFMATNIRVVCTYYFGEVMDLGFD- 65
           GAI++ KEL +LP E +Y++  GVWNLSSDQGNLG F  TN+R+V      E  ++    
Sbjct: 150 GAIIKEKELTLLPHEQVYTRLNGVWNLSSDQGNLGTFFVTNVRIVWHANLAENFNVSIPY 209

Query: 66  -------IRTYYFGEVM-----------DLGFDIRYESFENEVYHELFDSHHVF 101
                  IR   FG  +            LGF +  E   +E++ E+   H VF
Sbjct: 210 MQINSARIRGSKFGPALVIETSAASGGYILGFRVDPEERLHELFKEMQSLHSVF 263

Source: Phytophthora sojae

Species: Phytophthora sojae

Genus: Phytophthora

Order: Peronosporales

Superkingdom: Eukaryota

>gi|255087046|ref|XP_002505446.1| predicted protein [Micromonas sp. RCC299] gi|226520716|gb|ACO66704.1| predicted protein [Micromonas sp. RCC299]

 Score = 72.8 bits (177), Expect = 2e-11,   Method: Compositional matrix adjust.
 Identities = 42/114 (36%), Positives = 58/114 (50%), Gaps = 19/114 (16%)

Query: 7   GAIVRRKELKMLPREHIYSQYQGVWNLSSDQGNLGVFMATNIRVVCTYYFGEVMDLGF-- 64
           GAI+  KELK+LP+E +Y++  GVWNLSS+QGNLG F  TN+RVV      E  ++    
Sbjct: 150 GAIIADKELKLLPKEEVYNKVNGVWNLSSEQGNLGTFFITNVRVVWHANMAENFNVSIPY 209

Query: 65  ------DIRTYYFGEVM-----------DLGFDIRYESFENEVYHELFDSHHVF 101
                  +R   FG+ +            LGF +       EV+ E+   H VF
Sbjct: 210 IQMKSVSVRDSKFGQALVVHTTHRSGGYILGFRVDPAEKLKEVFTEIEGLHKVF 263

Source: Micromonas sp. RCC299

Species: Micromonas sp. RCC299

Genus: Micromonas

Order: Mamiellales

Class: Mamiellophyceae

Phylum: Chlorophyta

Superkingdom: Eukaryota

>gi|301112689|ref|XP_002998115.1| bardet-Biedl syndrome 5 protein [Phytophthora infestans T30-4] gi|262112409|gb|EEY70461.1| bardet-Biedl syndrome 5 protein [Phytophthora infestans T30-4]

 Score = 72.4 bits (176), Expect = 3e-11,   Method: Composition-based stats.
 Identities = 42/114 (36%), Positives = 58/114 (50%), Gaps = 19/114 (16%)

Query: 7   GAIVRRKELKMLPREHIYSQYQGVWNLSSDQGNLGVFMATNIRVVCTYYFGEVMDLGFD- 65
           GAI++ KEL +LP E +Y++  GVWNLSSDQGNLG F  TN+R+V      E  ++    
Sbjct: 150 GAIIKEKELTLLPHEQVYTRLNGVWNLSSDQGNLGTFFVTNVRIVWHANLAENFNVSIPY 209

Query: 66  -------IRTYYFGEVM-----------DLGFDIRYESFENEVYHELFDSHHVF 101
                  IR   FG  +            LGF +  E   +E++ E+   H VF
Sbjct: 210 MQINSARIRGSKFGPALVIETSAVSGGYILGFRVDPEDRLHELFKEMQSLHSVF 263

Source: Phytophthora infestans T30-4

Species: Phytophthora infestans

Genus: Phytophthora

Order: Peronosporales

Superkingdom: Eukaryota

>gi|195453986|ref|XP_002074033.1| GK12823 [Drosophila willistoni] gi|194170118|gb|EDW85019.1| GK12823 [Drosophila willistoni]

 Score = 71.6 bits (174), Expect = 6e-11,   Method: Compositional matrix adjust.
 Identities = 43/114 (37%), Positives = 58/114 (50%), Gaps = 19/114 (16%)

Query: 7   GAIVRRKELKMLPREHIYSQYQGVWNLSSDQGNLGVFMATNIRVVCTYYFGEVMDLGFD- 65
           GAI++  +L +LP E +YSQ QGVWNLSSDQGNLG F+ +NIR+V      E  ++    
Sbjct: 157 GAIIQTGQLIILPNEQVYSQVQGVWNLSSDQGNLGSFVVSNIRLVWFADANETFNISLPY 216

Query: 66  -------IRTYYFGEVM-----------DLGFDIRYESFENEVYHELFDSHHVF 101
                  IR   +G  +            LGF I      NE+Y EL   H ++
Sbjct: 217 LQIESLRIRDSKYGPALVIQTAETGGGYVLGFRIDPAERLNELYKELSSLHAIY 270

Source: Drosophila willistoni

Species: Drosophila willistoni

Genus: Drosophila

Family: Drosophilidae

Order: Diptera

Class: Insecta

Phylum: Arthropoda

Superkingdom: Eukaryota

>gi|148231484|ref|NP_001087782.1| Bardet-Biedl syndrome 5 protein homolog [Xenopus laevis] gi|82181251|sp|Q66IS6.1|BBS5_XENLA RecName: Full=Bardet-Biedl syndrome 5 protein homolog gi|51704053|gb|AAH81218.1| Bbs5 protein [Xenopus laevis]

 Score = 70.9 bits (172), Expect = 9e-11,   Method: Compositional matrix adjust.
 Identities = 28/45 (62%), Positives = 37/45 (82%)

Query: 7   GAIVRRKELKMLPREHIYSQYQGVWNLSSDQGNLGVFMATNIRVV 51
           GA+++ K+LK+LPRE +Y +  GVWNLSSDQGNLG F  TN+R+V
Sbjct: 145 GALIQNKQLKLLPREQVYDKINGVWNLSSDQGNLGTFFITNVRIV 189

Source: Xenopus laevis

Species: Xenopus laevis

Genus: Xenopus

Family: Pipidae

Order: Anura

Class: Amphibia

Phylum: Chordata

Superkingdom: Eukaryota

>gi|432098391|gb|ELK28191.1| Bardet-Biedl syndrome 5 protein like protein [Myotis davidii]

 Score = 70.1 bits (170), Expect = 2e-10,   Method: Composition-based stats.
 Identities = 27/44 (61%), Positives = 37/44 (84%)

Query: 8   AIVRRKELKMLPREHIYSQYQGVWNLSSDQGNLGVFMATNIRVV 51
           A+++ KEL++LP+EH+Y +  GVWNLSSDQGNLG F  TN+R+V
Sbjct: 145 ALIQNKELRLLPQEHVYDKINGVWNLSSDQGNLGTFFITNVRIV 188

Source: Myotis davidii

Species: Myotis davidii

Genus: Myotis

Family: Vespertilionidae

Order: Chiroptera

Class: Mammalia

Phylum: Chordata

Superkingdom: Eukaryota

>gi|348585675|ref|XP_003478596.1| PREDICTED: Bardet-Biedl syndrome 5 protein homolog [Cavia porcellus]

 Score = 69.7 bits (169), Expect = 2e-10,   Method: Composition-based stats.
 Identities = 27/44 (61%), Positives = 37/44 (84%)

Query: 8   AIVRRKELKMLPREHIYSQYQGVWNLSSDQGNLGVFMATNIRVV 51
           A++R K+L++LP+EH+Y +  GVWNLSSDQGNLG F  TN+R+V
Sbjct: 145 ALIRNKQLRLLPQEHVYDKISGVWNLSSDQGNLGTFFITNVRIV 188

Source: Cavia porcellus

Species: Cavia porcellus

Genus: Cavia

Family: Caviidae

Order: Rodentia

Class: Mammalia

Phylum: Chordata

Superkingdom: Eukaryota

>gi|148695091|gb|EDL27038.1| Bardet-Biedl syndrome 5 (human), isoform CRA_b [Mus musculus]

 Score = 69.7 bits (169), Expect = 2e-10,   Method: Composition-based stats.
 Identities = 26/44 (59%), Positives = 37/44 (84%)

Query: 8   AIVRRKELKMLPREHIYSQYQGVWNLSSDQGNLGVFMATNIRVV 51
           A+++ K+L++LP+EH+Y +  GVWNLSSDQGNLG F  TN+R+V
Sbjct: 199 AVIQNKQLRLLPQEHVYDKINGVWNLSSDQGNLGTFFITNVRIV 242

Source: Mus musculus

Species: Mus musculus

Genus: Mus

Family: Muridae

Order: Rodentia

Class: Mammalia

Phylum: Chordata

Superkingdom: Eukaryota

>gi|213510754|ref|NP_001134089.1| Bardet-Biedl syndrome 5 protein homolog [Salmo salar] gi|209730652|gb|ACI66195.1| Bardet-Biedl syndrome 5 protein homolog [Salmo salar]

 Score = 69.7 bits (169), Expect = 2e-10,   Method: Composition-based stats.
 Identities = 27/45 (60%), Positives = 37/45 (82%)

Query: 7   GAIVRRKELKMLPREHIYSQYQGVWNLSSDQGNLGVFMATNIRVV 51
           GA+++ K+L++LPRE +Y +  GVWNLSSDQGNLG F  TN+R+V
Sbjct: 145 GALIQNKQLRLLPREEVYDKINGVWNLSSDQGNLGTFFITNVRIV 189

Source: Salmo salar

Species: Salmo salar

Genus: Salmo

Family: Salmonidae

Order: Salmoniformes

Class: Actinopterygii

Phylum: Chordata

Superkingdom: Eukaryota

>gi|325182780|emb|CCA17235.1| bardetBiedl syndrome 5 protein putative [Albugo laibachii Nc14]

 Score = 69.7 bits (169), Expect = 2e-10,   Method: Composition-based stats.
 Identities = 42/116 (36%), Positives = 56/116 (48%), Gaps = 19/116 (16%)

Query: 7   GAIVRRKELKMLPREHIYSQYQGVWNLSSDQGNLGVFMATNIRVVCTYYFGEVMDLGFD- 65
           G I+R KEL +LP+E IY++  GVWNLSSDQGNLG  + TN+R+V      E  ++    
Sbjct: 148 GGIIREKELVLLPQEEIYTRVSGVWNLSSDQGNLGSLLITNVRIVWHATLAENFNVSIPF 207

Query: 66  -------IRTYYFGEVM-----------DLGFDIRYESFENEVYHELFDSHHVFRR 103
                  IR   FG  +            LGF I        V+ EL   H +F +
Sbjct: 208 IQLHSACIRKSKFGLALVLETSAGSGGYVLGFQIEPSEKLESVFQELERLHSIFYK 263

Source: Albugo laibachii Nc14

Species: Albugo laibachii

Genus: Albugo

Family: Albuginaceae

Order: Albuginales

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	103
FB|FBgn0037280	381	CG1126 [Drosophila melanogaste	0.436	0.118	0.644	4.2e-12
MGI|MGI:1919819	341	Bbs5 "Bardet-Biedl syndrome 5	0.427	0.129	0.590	1.1e-11
UNIPROTKB|Q8N3I7	341	BBS5 "Bardet-Biedl syndrome 5	0.427	0.129	0.590	1.4e-11
ZFIN|ZDB-GENE-040426-1083	342	bbs5 "Bardet-Biedl syndrome 5"	0.436	0.131	0.577	2.4e-11
UNIPROTKB|F1S1V8	544	LOC100156983 "Uncharacterized	0.427	0.080	0.590	3.8e-11
UNIPROTKB|A6QLF9	341	BBS5 "Uncharacterized protein"	0.427	0.129	0.568	6.5e-11
UNIPROTKB|E1C569	342	BBS5 "Uncharacterized protein"	0.427	0.128	0.613	6.6e-11
UNIPROTKB|A8HWI3	352	BBS5 "Bardet-biedl syndrome 5	0.436	0.127	0.555	1.9e-09
WB|WBGene00010974	361	bbs-5 [Caenorhabditis elegans	0.582	0.166	0.446	1.9e-08
UNIPROTKB|Q21626	361	bbs-5 "Protein BBS-5" [Caenorh	0.582	0.166	0.446	1.9e-08

FB|FBgn0037280 CG1126 [Drosophila melanogaster (taxid:7227)]

 Score = 151 (58.2 bits), Expect = 4.2e-12, Sum P(2) = 4.2e-12
 Identities = 29/45 (64%), Positives = 36/45 (80%)

Query:     7 GAIVRRKELKMLPREHIYSQYQGVWNLSSDQGNLGVFMATNIRVV 51
             GAIV+  +L +LP E ++S  QGVWNLSSDQGNLG F+ TNIR+V
Sbjct:   162 GAIVQAGQLVILPDEEVFSHVQGVWNLSSDQGNLGSFVVTNIRLV 206

GO:0042384 "cilium assembly" evidence=ISM;IEP

MGI|MGI:1919819 Bbs5 "Bardet-Biedl syndrome 5 (human)" [Mus musculus (taxid:10090)]

 Score = 164 (62.8 bits), Expect = 1.1e-11, P = 1.1e-11
 Identities = 26/44 (59%), Positives = 37/44 (84%)

Query:     8 AIVRRKELKMLPREHIYSQYQGVWNLSSDQGNLGVFMATNIRVV 51
             A+++ K+L++LP+EH+Y +  GVWNLSSDQGNLG F  TN+R+V
Sbjct:   145 AVIQNKQLRLLPQEHVYDKINGVWNLSSDQGNLGTFFITNVRIV 188

GO:0005737 "cytoplasm" evidence=IEA

GO:0005856 "cytoskeleton" evidence=IEA

GO:0005886 "plasma membrane" evidence=IEA

GO:0005929 "cilium" evidence=IEA

GO:0005932 "microtubule basal body" evidence=IDA

GO:0016020 "membrane" evidence=IEA

GO:0032266 "phosphatidylinositol-3-phosphate binding" evidence=ISO

GO:0034464 "BBSome" evidence=ISO;IDA

GO:0042384 "cilium assembly" evidence=ISO

GO:0042995 "cell projection" evidence=IEA

UNIPROTKB|Q8N3I7 BBS5 "Bardet-Biedl syndrome 5 protein" [Homo sapiens (taxid:9606)]

 Score = 163 (62.4 bits), Expect = 1.4e-11, P = 1.4e-11
 Identities = 26/44 (59%), Positives = 37/44 (84%)

Query:     8 AIVRRKELKMLPREHIYSQYQGVWNLSSDQGNLGVFMATNIRVV 51
             A+++ K+L++LP+EH+Y +  GVWNLSSDQGNLG F  TN+R+V
Sbjct:   145 ALIQNKQLRLLPQEHVYDKINGVWNLSSDQGNLGTFFITNVRIV 188

GO:0007601 "visual perception" evidence=IEA

GO:0050896 "response to stimulus" evidence=IEA

GO:0005737 "cytoplasm" evidence=IEA

GO:0060170 "cilium membrane" evidence=IEA

GO:0005515 "protein binding" evidence=IPI

GO:0034464 "BBSome" evidence=IDA

GO:0005622 "intracellular" evidence=IDA

GO:0032402 "melanosome transport" evidence=ISS

GO:0001947 "heart looping" evidence=ISS

GO:0005932 "microtubule basal body" evidence=ISS

GO:0009296 "flagellum assembly" evidence=ISS

GO:0032266 "phosphatidylinositol-3-phosphate binding" evidence=IDA

GO:0042384 "cilium assembly" evidence=IMP

ZFIN|ZDB-GENE-040426-1083 bbs5 "Bardet-Biedl syndrome 5" [Danio rerio (taxid:7955)]

 Score = 161 (61.7 bits), Expect = 2.4e-11, P = 2.4e-11
 Identities = 26/45 (57%), Positives = 37/45 (82%)

Query:     7 GAIVRRKELKMLPREHIYSQYQGVWNLSSDQGNLGVFMATNIRVV 51
             GA+++ K+L++LP+E +Y +  GVWNLSSDQGNLG F  TN+R+V
Sbjct:   145 GALIQNKQLRLLPQEQVYDKINGVWNLSSDQGNLGTFFITNVRIV 189

GO:0001947 "heart looping" evidence=IMP

GO:0032402 "melanosome transport" evidence=IMP

GO:0046907 "intracellular transport" evidence=IMP

GO:0051877 "pigment granule aggregation in cell center" evidence=IMP

GO:0070121 "Kupffer's vesicle development" evidence=IMP

GO:0005622 "intracellular" evidence=IMP

GO:0005886 "plasma membrane" evidence=IEA

GO:0005929 "cilium" evidence=IEA

GO:0042995 "cell projection" evidence=IEA

GO:0016020 "membrane" evidence=IEA

GO:0060170 "cilium membrane" evidence=IEA

UNIPROTKB|F1S1V8 LOC100156983 "Uncharacterized protein" [Sus scrofa (taxid:9823)]

 Score = 163 (62.4 bits), Expect = 3.8e-11, P = 3.8e-11
 Identities = 26/44 (59%), Positives = 37/44 (84%)

Query:     8 AIVRRKELKMLPREHIYSQYQGVWNLSSDQGNLGVFMATNIRVV 51
             A+++ K+L++LP+EH+Y +  GVWNLSSDQGNLG F  TN+R+V
Sbjct:   145 ALIQNKQLRLLPQEHVYDKINGVWNLSSDQGNLGTFFITNVRIV 188

GO:0034464 "BBSome" evidence=IEA

GO:0005932 "microtubule basal body" evidence=IEA

UNIPROTKB|A6QLF9 BBS5 "Uncharacterized protein" [Bos taurus (taxid:9913)]

 Score = 157 (60.3 bits), Expect = 6.5e-11, P = 6.5e-11
 Identities = 25/44 (56%), Positives = 38/44 (86%)

Query:     8 AIVRRKELKMLPREHIYSQYQGVWNLSSDQGNLGVFMATNIRVV 51
             A+++ K+L++LP+E++Y++  GVWNLSSDQGNLG F  TN+R+V
Sbjct:   145 ALIQNKQLRLLPQENVYNKINGVWNLSSDQGNLGTFFITNVRIV 188

GO:0042384 "cilium assembly" evidence=IEA

GO:0034464 "BBSome" evidence=IEA

GO:0032266 "phosphatidylinositol-3-phosphate binding" evidence=IEA

GO:0005932 "microtubule basal body" evidence=IEA

UNIPROTKB|E1C569 BBS5 "Uncharacterized protein" [Gallus gallus (taxid:9031)]

 Score = 157 (60.3 bits), Expect = 6.6e-11, P = 6.6e-11
 Identities = 27/44 (61%), Positives = 36/44 (81%)

Query:     8 AIVRRKELKMLPREHIYSQYQGVWNLSSDQGNLGVFMATNIRVV 51
             A+++ K+L++LP+E IY +  GVWNLSSDQGNLG F  TN+RVV
Sbjct:   146 ALIQNKQLRLLPQEQIYDKINGVWNLSSDQGNLGTFFITNVRVV 189

GO:0005932 "microtubule basal body" evidence=IEA

GO:0034464 "BBSome" evidence=IEA

UNIPROTKB|A8HWI3 BBS5 "Bardet-biedl syndrome 5 protein" [Chlamydomonas reinhardtii (taxid:3055)]

 Score = 144 (55.7 bits), Expect = 1.9e-09, P = 1.9e-09
 Identities = 25/45 (55%), Positives = 35/45 (77%)

Query:     7 GAIVRRKELKMLPREHIYSQYQGVWNLSSDQGNLGVFMATNIRVV 51
             GA++  +EL++LP E  YS+  GVWNLS++QGNLG F  +N+RVV
Sbjct:   153 GALIMDRELRLLPLEQTYSRVNGVWNLSAEQGNLGTFFISNVRVV 197

GO:0009296 "flagellum assembly" evidence=IMP

WB|WBGene00010974 bbs-5 [Caenorhabditis elegans (taxid:6239)]

 Score = 135 (52.6 bits), Expect = 1.9e-08, P = 1.9e-08
 Identities = 29/65 (44%), Positives = 44/65 (67%)

Query:     7 GAIVRRK-ELKMLPREHIYSQYQGVWNLSSDQGNLGVFMATNIRVVCTYYFGEVMDLGFD 65
             G  +R    LK+LP+E I  +  GVWNLS++ G+LGVF+ TNIRVV   ++ E M++G++
Sbjct:   153 GVFIREDGTLKILPQETIIEKVNGVWNLSTETGSLGVFVITNIRVV---WYAE-MNIGYN 208

Query:    66 IRTYY 70
             +   Y
Sbjct:   209 VSVPY 213

GO:0000287 "magnesium ion binding" evidence=IEA

GO:0004427 "inorganic diphosphatase activity" evidence=IEA

GO:0005737 "cytoplasm" evidence=IEA

GO:0006796 "phosphate-containing compound metabolic process" evidence=IEA

GO:0005932 "microtubule basal body" evidence=IDA

GO:0031513 "nonmotile primary cilium" evidence=IDA

UNIPROTKB|Q21626 bbs-5 "Protein BBS-5" [Caenorhabditis elegans (taxid:6239)]

 Score = 135 (52.6 bits), Expect = 1.9e-08, P = 1.9e-08
 Identities = 29/65 (44%), Positives = 44/65 (67%)

Query:     7 GAIVRRK-ELKMLPREHIYSQYQGVWNLSSDQGNLGVFMATNIRVVCTYYFGEVMDLGFD 65
             G  +R    LK+LP+E I  +  GVWNLS++ G+LGVF+ TNIRVV   ++ E M++G++
Sbjct:   153 GVFIREDGTLKILPQETIIEKVNGVWNLSTETGSLGVFVITNIRVV---WYAE-MNIGYN 208

Query:    66 IRTYY 70
             +   Y
Sbjct:   209 VSVPY 213


Parameters:
  V=100
  filter=SEG
  E=0.001

  ctxfactor=1.00

  Query                        -----  As Used  -----    -----  Computed  ----
  Frame  MatID Matrix name     Lambda    K       H      Lambda    K       H
   +0      0   BLOSUM62        0.327   0.142   0.437    same    same    same
               Q=9,R=2         0.244   0.0300  0.180     n/a     n/a     n/a

  Query
  Frame  MatID  Length  Eff.Length     E     S W   T  X   E2     S2
   +0      0      103       103   0.00091  102 3  11 22  0.49    29
                                                     29  0.49    31


Statistics:

  Database:  /share/blast/go-seqdb.fasta
   Title:  go_20130330-seqdb.fasta
   Posted:  5:47:42 AM PDT Apr 1, 2013
   Created:  5:47:42 AM PDT Apr 1, 2013
   Format:  XDF-1
   # of letters in database:  169,044,731
   # of sequences in database:  368,745
   # of database sequences satisfying E:  10
  No. of states in DFA:  563 (60 KB)
  Total size of DFA:  122 KB (2078 KB)
  Time to generate neighborhood:  0.00u 0.00s 0.00t   Elapsed:  00:00:01
  No. of threads or processors used:  24
  Search cpu time:  11.21u 0.07s 11.28t   Elapsed:  00:00:09
  Total cpu time:  11.21u 0.07s 11.28t   Elapsed:  00:00:10
  Start:  Thu Aug 15 12:32:14 2013   End:  Thu Aug 15 12:32:24 2013

GO:0005932 "microtubule basal body" evidence=IDA

---

Prediction of Enzyme Commission (EC) Number

EC Number Prediction by Annotation Transfer from SWISS-PROT Entries

Original result of BLAST against SWISS-PROT

ID	Name	Annotated EC number	Identity	Query coverage	Hit coverage	RBH(Q2H)	RBH(H2Q)
Q8N3I7	BBS5_HUMAN	No assigned EC number	0.5909	0.4271	0.1290	yes	N/A
Q9CZQ9	BBS5_MOUSE	No assigned EC number	0.5909	0.4271	0.1290	yes	N/A
Q7ZWB7	BBS5_DANRE	No assigned EC number	0.5777	0.4368	0.1315	yes	N/A

EC Number Prediction by EFICAz Software

No EC number assignment, probably not an enzyme!

---

Conserved Domains and Related Protein Families

Conserved Domains Detected by RPS-BLAST

Original result of RPS-BLAST against CDD database

ID	Length	Definition	E-value
Query	103
pfam07289	339	pfam07289, DUF1448, Protein of unknown function (D	1e-18
smart00683	55	smart00683, DM16, Repeats in sea squirt COS41	4e-09

>gnl|CDD|219365 pfam07289, DUF1448, Protein of unknown function (DUF1448)

 Score = 78.2 bits (193), Expect = 1e-18
 Identities = 27/45 (60%), Positives = 34/45 (75%)

Query: 7   GAIVRRKELKMLPREHIYSQYQGVWNLSSDQGNLGVFMATNIRVV 51
           GAI++  +L +LP E + S+  GVWNLSSDQGNLG F  TNIR+V
Sbjct: 141 GAIIQDGQLNLLPDEQVISRVNGVWNLSSDQGNLGTFFVTNIRLV 185

This family consists of several eukaryotic proteins of around 375 residues in length. The function of this family is unknown. It appears that this family includes a divergent GRAM domain. Length = 339

>gnl|CDD|128926 smart00683, DM16, Repeats in sea squirt COS41

 Score = 48.0 bits (115), Expect = 4e-09
 Identities = 12/31 (38%), Positives = 20/31 (64%)

Query: 21 EHIYSQYQGVWNLSSDQGNLGVFMATNIRVV 51
          E + ++  GV +   + G+LGVF  TN+R+V
Sbjct: 1  ERVLTRINGVEDTKGNNGDLGVFFVTNLRLV 31

4, worm R01H10.6, fly CG1126 etc. Length = 55

Conserved Domains Detected by HHsearch

Original result of HHsearch against CDD database

ID	Length	Definition	Probability
Query	103
PF07289	339	DUF1448: Protein of unknown function (DUF1448); In	100.0
smart00683	55	DM16 Repeats in sea squirt COS41.4, worm R01H10.6,	99.82
PF07289	339	DUF1448: Protein of unknown function (DUF1448); In	99.57
PF14470	96	bPH_3: Bacterial PH domain	94.21
PF08000	124	bPH_1: Bacterial PH domain; InterPro: IPR012544 Th	92.72
PF11605	89	Vps36_ESCRT-II: Vacuolar protein sorting protein 3	92.6
PF03517	135	Voldacs: Regulator of volume decrease after cellul	91.66
KOG3294|consensus	261		90.1
smart00568	61	GRAM domain in glucosyltransferases, myotubularins	88.16

>PF07289 DUF1448: Protein of unknown function (DUF1448); InterPro: IPR006606 This entry represents the Bardet-Biedl syndrome 5 protein (BBL5)

Probab=100.00  E-value=2.4e-45  Score=300.13  Aligned_cols=96  Identities=46%  Similarity=0.841  Sum_probs=93.6

Q ss_pred             cchhceeecCcccccCCceeeeeecceeeccCCCCceeeEEEEeeeEEeecccccccCccee------------eeeecc
Q psy11196          4 YKSGAIVRRKELKMLPREHIYSQYQGVWNLSSDQGNLGVFMATNIRVVCTYYFGEVMDLGFD------------IRTYYF   71 (103)
Q Consensus         4 ~lR~aii~~~~L~llp~E~i~~~~~gV~nls~dqG~lG~f~vTNiRiV~~~W~a~~~~~~fN------------ir~SKF   71 (103)
                      .||||||++|+|+|||+|+||++++||||||+|||++|+|++||||+|   |||+ +|+.||            ||+|||
T Consensus       138 klR~aii~~~~L~lLp~E~v~~~~~gVwnls~dqGnLGtfivTNvRiV---W~A~-~ne~fNVSiPylqi~~i~ir~SKf  213 (339)
T PF07289_consen  138 KLRGAIIRDGQLKLLPQEQVYSRVNGVWNLSSDQGNLGTFIVTNVRIV---WFAD-MNESFNVSIPYLQIKSIRIRDSKF  213 (339)
T ss_pred             hhheeeeeCCeEeeCCccEEeeccCCEEEcccCCCceeEEEEeeeEEE---EEcc-CCccccccchHhhheeeeeecccc
Confidence            479999999999999999999999999999999999999999999999   9999 999999            799999


Q ss_pred             cceE-----------EEEEecChhhHHHHHHHHHHHhhhhhhC
Q psy11196         72 GEVM-----------DLGFDIRYESFENEVYHELFDSHHVFRR  103 (103)
Q Consensus        72 G~al-----------vLGFridP~e~l~~~~ke~~~L~~~~~~  103 (103)
                      |+||           +|||||||+|||+++||||++||++|+.
T Consensus       214 G~aLVieT~~~sGgYVLGFRvDP~ErL~~l~KEi~sLh~vy~~  256 (339)
T PF07289_consen  214 GPALVIETSESSGGYVLGFRVDPEERLQELFKEIQSLHKVYSA  256 (339)
T ss_pred             ceEEEEEEeccCCcEEEEEEcCHHHHHHHHHHHHHHHHHHHHh
Confidence            9999           9999999999999999999999999973

It consists of eukaryotic proteins of around 375 residues in length.

>smart00683 DM16 Repeats in sea squirt COS41

Probab=99.82  E-value=1.3e-20  Score=120.89  Aligned_cols=43  Identities=33%  Similarity=0.692  Sum_probs=41.8

Q ss_pred             ceeeeeecceeeccCCCCceeeEEEEeeeEEeecccccccCcceeee
Q psy11196         21 EHIYSQYQGVWNLSSDQGNLGVFMATNIRVVCTYYFGEVMDLGFDIR   67 (103)
Q Consensus        21 E~i~~~~~gV~nls~dqG~lG~f~vTNiRiV~~~W~a~~~~~~fNir   67 (103)
                      |+++++++||||+|||||++|+|+|||||+|   |||+ ++++||++
T Consensus         1 E~v~~~~~~Ved~kgn~G~~G~l~VTNlRii---W~s~-~~~~~NlS   43 (55)
T smart00683        1 ERVLTRINGVEDTKGNNGDLGVFFVTNLRLV---WHSD-TNPRFNIS   43 (55)
T ss_pred             CEEEeeecCeEecCCCCCCeeEEEEEeeEEE---EEeC-CCCceEEE
Confidence            8999999999999999999999999999999   9999 99999964

4, worm R01H10.6, fly CG1126 etc.

>PF07289 DUF1448: Protein of unknown function (DUF1448); InterPro: IPR006606 This entry represents the Bardet-Biedl syndrome 5 protein (BBL5)

Probab=99.57  E-value=7e-15  Score=121.18  Aligned_cols=87  Identities=17%  Similarity=0.349  Sum_probs=70.4

Q ss_pred             cCcccccCCceeeeeecceeeccCCCCceeeEEEEeeeEEeecccccccCcceee------------e--ee-ccc--ce
Q psy11196         12 RKELKMLPREHIYSQYQGVWNLSSDQGNLGVFMATNIRVVCTYYFGEVMDLGFDI------------R--TY-YFG--EV   74 (103)
Q Consensus        12 ~~~L~llp~E~i~~~~~gV~nls~dqG~lG~f~vTNiRiV~~~W~a~~~~~~fNi------------r--~S-KFG--~a   74 (103)
                      .+.|+|+|||.+++++++|+|+|||+|+.|+++|||+|++   |||. .++++|+            |  .| ..|  +|
T Consensus        14 ~~~l~lr~GE~~i~~~~~VEDtKGN~G~~G~l~vTNLR~i---W~s~-~~~r~NlSIG~~~i~~i~~~~~~sklrg~tea   89 (339)
T PF07289_consen   14 ASQLKLRPGEFIIDRLDPVEDTKGNNGDRGRLVVTNLRLI---WHSL-KRPRINLSIGYNCITNISTKTVNSKLRGNTEA   89 (339)
T ss_pred             hhhcccccceEEEEeeeceeeccCCCCCeeEEEEEeeeeE---Eecc-CCCceeEEeeceeEEEEEEEEeeccccCceeE
Confidence            4578999999999999999999999999999999999999   9999 9999992            1  22 234  23


Q ss_pred             E-----------EEEEe-cChhhHHHHHHHHHHHhhhhhh
Q psy11196         75 M-----------DLGFD-IRYESFENEVYHELFDSHHVFR  102 (103)
Q Consensus        75 l-----------vLGFr-idP~e~l~~~~ke~~~L~~~~~  102 (103)
                      |           .+=|. +.|...-+.+|..+++.|++|.
T Consensus        90 L~i~~k~~~~rfEFiFt~~~~~~~~~~lf~~v~~v~raY~  129 (339)
T PF07289_consen   90 LYILAKFNNTRFEFIFTNLSPNSPRQRLFTSVQAVYRAYE  129 (339)
T ss_pred             EEEeeecCCceEEEEeccCCCCCccchHHHHHHHHHHHHH
Confidence            3           55565 4444433689999999999985

It consists of eukaryotic proteins of around 375 residues in length.

>PF14470 bPH_3: Bacterial PH domain

Probab=94.21  E-value=0.091  Score=33.67  Aligned_cols=38  Identities=24%  Similarity=0.329  Sum_probs=32.9

Q ss_pred             ccCCceeeeeecceeeccCCCCceeeEEEEeeeEEeeccccc
Q psy11196         17 MLPREHIYSQYQGVWNLSSDQGNLGVFMATNIRVVCTYYFGE   58 (103)
Q Consensus        17 llp~E~i~~~~~gV~nls~dqG~lG~f~vTNiRiV~~~W~a~   58 (103)
                      |.|+|+|.....|..+.. .....|.+++||=|++   -++.
T Consensus         1 L~~~E~I~~~~~~~~~~~-~~~~~g~l~~TnkRli---f~~~   38 (96)
T PF14470_consen    1 LKEDEEIEYVAVGSYNYF-FTSFPGVLVLTNKRLI---FYSK   38 (96)
T ss_pred             CcCCCEEEEEEEEEEeec-ccCceeEEEEeCCEEE---EEEc
Confidence            579999999999998876 6667999999999999   6655

>PF08000 bPH_1: Bacterial PH domain; InterPro: IPR012544 This family contains many bacterial hypothetical proteins

Probab=92.72  E-value=0.1  Score=37.82  Aligned_cols=28  Identities=25%  Similarity=0.299  Sum_probs=22.9

Q ss_pred             ccccCCceeeeeecceeeccCCCCceeeEEEEeeeEE
Q psy11196         15 LKMLPREHIYSQYQGVWNLSSDQGNLGVFMATNIRVV   51 (103)
Q Consensus        15 L~llp~E~i~~~~~gV~nls~dqG~lG~f~vTNiRiV   51 (103)
                      --|+|||+|+..++.+         ...+++||-|+|
T Consensus        25 ~~L~~gE~I~~ayk~i---------RD~~vFTnkRlI   52 (124)
T PF08000_consen   25 PLLLDGEEIEAAYKLI---------RDEIVFTNKRLI   52 (124)
T ss_dssp             GGSSTT--EEEEEEES---------SEEEEEESSEEE
T ss_pred             HhcCCCCeeeeeehhh---------ceeEEEecChhe
Confidence            3578999999999998         567999999999

; PDB: 3DCX_E 3HSA_C 3B77_D.

>PF11605 Vps36_ESCRT-II: Vacuolar protein sorting protein 36 Vps36; InterPro: IPR021648 Vps36 is a subunit of ESCRT-II, a protein involved in driving protein sorting from endosomes to lysosomes

Probab=92.60  E-value=0.15  Score=34.86  Aligned_cols=47  Identities=23%  Similarity=0.250  Sum_probs=30.9

Q ss_pred             cCcccccCCceeeeeecceeeccCC----CCceeeEEEEeeeEEeecccccccCc
Q psy11196         12 RKELKMLPREHIYSQYQGVWNLSSD----QGNLGVFMATNIRVVCTYYFGEVMDL   62 (103)
Q Consensus        12 ~~~L~llp~E~i~~~~~gV~nls~d----qG~lG~f~vTNiRiV~~~W~a~~~~~   62 (103)
                      -++-.|.|+|.++-+-+||-=-.|+    +=.-|++++|+-|++   |... .++
T Consensus         5 s~~p~L~~~E~~~~~q~~V~LYdG~~K~~~~q~G~l~LTsHRli---w~d~-~~~   55 (89)
T PF11605_consen    5 SGRPVLEPNETIVYQQDGVGLYDGDQKTPNFQNGRLYLTSHRLI---WVDD-SDP   55 (89)
T ss_dssp             CS-C---TT--EEEEEEEEEEEETTECSTT-SCEEEEEESSEEE---EEES-SGH
T ss_pred             ccCcccCCCceEEEEecCeeeEcCCccCccccCCEEEEEeeEEE---EEcC-CCC
Confidence            4567789999999999998543344    444899999999999   9866 444

The GLUE domain of Vps36 allows for a tight interaction to occur between the protein and Vps28, a subunit of ESCRT-I. This interaction is critical for ubiquitinated cargo progression from early to late endosomes []. ; PDB: 2HTH_B 2DX5_A 2CAY_B.

>PF03517 Voldacs: Regulator of volume decrease after cellular swelling; InterPro: IPR003521 The nucleotide-sensitive chloride conductance regulatory protein (ICln) is found ubiquitously in mammalian (and other) cell types and is postulated to play a critical role in cell volume regulation

Probab=91.66  E-value=0.37  Score=34.50  Aligned_cols=22  Identities=32%  Similarity=0.652  Sum_probs=17.9

Q ss_pred             eeEEEEeeeEEeecccc-cccCccee
Q psy11196         41 GVFMATNIRVVCTYYFG-EVMDLGFD   65 (103)
Q Consensus        41 G~f~vTNiRiV~~~W~a-~~~~~~fN   65 (103)
                      |+|+|||=||+   |++ .-.+.+|.
T Consensus         1 g~L~Vt~~~l~---w~~~~~~~~G~~   23 (135)
T PF03517_consen    1 GTLYVTESRLI---WFSNEDSSKGFS   23 (135)
T ss_dssp             EEEEEETTEEE---EEET--TTEEEE
T ss_pred             CEEEEecCEEE---EECCCcCCccee
Confidence            89999999999   999 33677777

Initial studies proposed that ICln was itself a swelling-activated anion channel; however, further studies demonstrated that it is localised primarily to the cell cytoplasm. It has therefore been postulated that activation of cell volume regulation may involve reversible translocation of ICln from the cytoplasm, and its insertion into the plasma membrane. It is not resolved whether the anionic channel involved in cell volume regulation after cell-swelling comprises one or more subunits, and if it does, whether ICln is in fact one of them [].; GO: 0006821 chloride transport, 0006884 cell volume homeostasis; PDB: 1ZYI_A.

>KOG3294|consensus

Probab=90.10  E-value=0.52  Score=38.60  Aligned_cols=45  Identities=18%  Similarity=0.258  Sum_probs=32.9

Q ss_pred             ecCc-ccccCCceeeeeecceeec---------cCCCCceeeEEEEeeeEEeeccccc
Q psy11196         11 RRKE-LKMLPREHIYSQYQGVWNL---------SSDQGNLGVFMATNIRVVCTYYFGE   58 (103)
Q Consensus        11 ~~~~-L~llp~E~i~~~~~gV~nl---------s~dqG~lG~f~vTNiRiV~~~W~a~   58 (103)
                      .++. ..+..+|-|+...+||+=+         .+..-..|++.+||.|||   .-+.
T Consensus         9 ~~G~gvl~~n~E~i~~~~dgVeltf~~~~~~~~~~~g~kkGtlyLTs~RiI---Fis~   63 (261)
T KOG3294|consen    9 PNGRGVLLYNGESILLVYDGVELTFDKMPNVSEAFKGTKKGTLYLTSHRII---FISS   63 (261)
T ss_pred             CCCceeEecCCcEEEEeeCCeEEEeeccCCCccccccceeeeEEeecceEE---EecC
Confidence            3454 6778888888888888622         222234899999999999   7776

>smart00568 GRAM domain in glucosyltransferases, myotubularins and other putative membrane-associated proteins

Probab=88.16  E-value=1.1  Score=27.11  Aligned_cols=37  Identities=16%  Similarity=0.304  Sum_probs=30.9

Q ss_pred             ccCCceeeeeecceeeccCCCCceeeEEEEeeeEEeeccccc
Q psy11196         17 MLPREHIYSQYQGVWNLSSDQGNLGVFMATNIRVVCTYYFGE   58 (103)
Q Consensus        17 llp~E~i~~~~~gV~nls~dqG~lG~f~vTNiRiV~~~W~a~   58 (103)
                      |.++|.+++.+.--+.  ++.-.-|++++||=|+.   .+|.
T Consensus         2 l~~~E~l~~~~~C~l~--~~~~~~G~lyiT~~~l~---F~S~   38 (61)
T smart00568        2 LPEEEKLIADYSCYLS--RDGPVQGRLYISNYRLC---FRSD   38 (61)
T ss_pred             cCCCcEEEEEEEeEEC--CCccccEEEEEECCEEE---EEcc
Confidence            6789999999875555  55667999999999999   8886

---

Homologous Structure Templates

Structure Templates Detected by BLAST

Original result of BLAST against Protein Data Bank

No homologous structure with e-value below 0.005

Structure Templates Detected by RPS-BLAST

Original result of RPS-BLAST against PDB70 database

No hit with e-value below 0.005

Structure Templates Detected by HHsearch

Original result of HHsearch against PDB70 database

ID	Length	Definition	Probability
Query	103
1zsq_A	528	Myotubularin-related protein 2; protein-phospholip	95.73
1lw3_A	657	Myotubularin-related protein 2; protein-phosphate	95.41
2cay_A	145	VPS36 PH domain, vacuolar protein sorting protein	93.54
3hsa_A	126	Pleckstrin homology domain; YP_926556.1, structura	92.29
2hth_B	140	Vacuolar protein sorting protein 36; GLUE domain,	88.79
1zyi_A	171	Methylosome subunit picln; PH domain, cell volume	85.28
3b77_A	193	Uncharacterized protein; pleckstrin-homology domai	82.86

>1zsq_A Myotubularin-related protein 2; protein-phospholipid complex, hydrolase; HET: PIB; 1.82A {Homo sapiens} SCOP: b.55.1.8 c.45.1.3 PDB: 1zvr_A*

Probab=95.73  E-value=0.0068  Score=51.49  Aligned_cols=45  Identities=22%  Similarity=0.229  Sum_probs=36.0

Q ss_pred             ecCcccccCCceeeeeecceeeccCC-CCceeeEEEEeeeEEeeccccc
Q psy11196         11 RRKELKMLPREHIYSQYQGVWNLSSD-QGNLGVFMATNIRVVCTYYFGE   58 (103)
Q Consensus        11 ~~~~L~llp~E~i~~~~~gV~nls~d-qG~lG~f~vTNiRiV~~~W~a~   58 (103)
                      +...+.+||||.|....++|.-+... ....|++.+||-|++   +.+.
T Consensus         3 ~~~~~~lLpGE~i~~~~~~V~~~~~~~~~~~GtL~lTn~rLi---F~~~   48 (528)
T 1zsq_A            3 SMEEPPLLPGENIKDMAKDVTYICPFTGAVRGTLTVTNYRLY---FKSM   48 (528)
T ss_dssp             -CCCCCCCTTCCEEEEEEEEEEEETTTEEEEEEEEEESSEEE---EEES
T ss_pred             CCCCCccCCCCeEEEecCCEEEeCCCCCceeEEEEEECCEEE---EEcc
Confidence            34578899999999999999866433 245899999999999   7765

>1lw3_A Myotubularin-related protein 2; protein-phosphate complex, hydrolase; 2.30A {Homo sapiens} SCOP: b.55.1.8 c.45.1.3 PDB: 1m7r_A

Probab=95.41  E-value=0.066  Score=46.82  Aligned_cols=43  Identities=26%  Similarity=0.309  Sum_probs=34.8

Q ss_pred             CcccccCCceeeeeecceeeccCCCC-ceeeEEEEeeeEEeeccccc
Q psy11196         13 KELKMLPREHIYSQYQGVWNLSSDQG-NLGVFMATNIRVVCTYYFGE   58 (103)
Q Consensus        13 ~~L~llp~E~i~~~~~gV~nls~dqG-~lG~f~vTNiRiV~~~W~a~   58 (103)
                      ..+.+||||.++...++|--+....| ..|++++||-|++   +-+.
T Consensus        77 ~~~~LLPGE~v~~~~~~V~~~~p~~~~~~GtL~lTnyrLi---F~~~  120 (657)
T 1lw3_A           77 EEPPLLPGENIKDMAKDVTYICPFTGAVRGTLTVTNYRLY---FKSM  120 (657)
T ss_dssp             CCCCCCTTCCEEEEEEEEEEEETTTEEEEEEEEEESSEEE---EEES
T ss_pred             CCCccCCCeeEEEEeCCEEEecCCCCceeEEEEEECCEEE---EEec
Confidence            36889999999999999975543333 5799999999999   7764

>2cay_A VPS36 PH domain, vacuolar protein sorting protein 36; transport, ESCRT-II, lipid-binding, multivesicular bodies, membrane; 1.9A {Saccharomyces cerevisiae} SCOP: b.55.1.12

Probab=93.54  E-value=0.051  Score=38.91  Aligned_cols=50  Identities=16%  Similarity=0.065  Sum_probs=38.6

Q ss_pred             eecCcccccCCceeeeeecceeeccC----CCCceeeEEEEeeeEEeecccccccCcc
Q psy11196         10 VRRKELKMLPREHIYSQYQGVWNLSS----DQGNLGVFMATNIRVVCTYYFGEVMDLG   63 (103)
Q Consensus        10 i~~~~L~llp~E~i~~~~~gV~nls~----dqG~lG~f~vTNiRiV~~~W~a~~~~~~   63 (103)
                      -..++..|+|+|.++-..++|-=--|    .+-.-|++++||-|||   |-.. .++.
T Consensus        18 t~s~rp~L~~~E~il~~~d~V~Ly~g~~k~~~~k~G~lyLTshRII---yv~~-~~~~   71 (145)
T 2cay_A           18 TSSGQPLLREGEKDIFIDQSVGLYHGKSKILQRQRGRIFLTSQRII---YIDD-AKPT   71 (145)
T ss_dssp             CTTSCBCCCTTCCEEEEEEEEEEEETTEECTTCSSEEEEEESSEEE---EEES-SCHH
T ss_pred             ccCCCcccCCCCEEEEEcCCcceecCCCCCcCCCccEEEEEcCeEE---EecC-CCCc
Confidence            34677889999999999998853333    3445899999999999   9766 5554

>3hsa_A Pleckstrin homology domain; YP_926556.1, structural genomics center for structural genomics, JCSG, protein structure INI PSI-2; HET: MSE MLY; 1.99A {Shewanella amazonensis SB2B} SCOP: b.55.1.13 PDB: 3dcx_A*

Probab=92.29  E-value=0.085  Score=37.63  Aligned_cols=29  Identities=21%  Similarity=0.256  Sum_probs=25.6

Q ss_pred             cccccCCceeeeeecceeeccCCCCceeeEEEEeeeEE
Q psy11196         14 ELKMLPREHIYSQYQGVWNLSSDQGNLGVFMATNIRVV   51 (103)
Q Consensus        14 ~L~llp~E~i~~~~~gV~nls~dqG~lG~f~vTNiRiV   51 (103)
                      +--|+|||+|+..+..+         ...+++||-|+|
T Consensus        25 ~~~L~~gE~I~~ayk~i---------RD~~iFTnkRLI   53 (126)
T 3hsa_A           25 SPILGDNEELQLAYXMV---------RDLFVFTSXRLI   53 (126)
T ss_dssp             GGGSCTTCCEEEEEECS---------SEEEEEESSEEE
T ss_pred             HHhcCCCCEEEEEehhc---------CceEEEccchhe
Confidence            34689999999999988         567999999999

>2hth_B Vacuolar protein sorting protein 36; GLUE domain, PH domain, viral budding, ubiquitin complex, protein transport; 2.70A {Homo sapiens} SCOP: b.55.1.12 PDB: 2dx5_A

Probab=88.79  E-value=0.48  Score=33.46  Aligned_cols=43  Identities=14%  Similarity=0.203  Sum_probs=33.7

Q ss_pred             CcccccCCceeeeeecceeeccCC---CCceeeEEEEeeeEEeeccccc
Q psy11196         13 KELKMLPREHIYSQYQGVWNLSSD---QGNLGVFMATNIRVVCTYYFGE   58 (103)
Q Consensus        13 ~~L~llp~E~i~~~~~gV~nls~d---qG~lG~f~vTNiRiV~~~W~a~   58 (103)
                      ++-.|.|+|.++-..++|-=--||   +=.-|++++|+-|||   |...
T Consensus         9 ~~p~L~~~E~il~~~~~V~Ly~g~~k~~~k~G~l~LTshRli---~~d~   54 (140)
T 2hth_B            9 TSGLLEINETLVIQQRGVRIYDGEEKIKFDAGTLLLSTHRLI---WRDQ   54 (140)
T ss_dssp             CCSCCCTTCCEEEEEEEEEEEETTCSSCCCCEEEEEESSEEE---EEET
T ss_pred             CCCccCCCceEEEEcCCEEEecCCCCCCCcCCEEEEecCEEE---EecC
Confidence            556689999999999988533333   224899999999999   9765

>1zyi_A Methylosome subunit picln; PH domain, cell volume regulation, RNA splicing, translation; NMR {Canis lupus familiaris}

Probab=85.28  E-value=0.64  Score=34.33  Aligned_cols=42  Identities=21%  Similarity=0.357  Sum_probs=26.0

Q ss_pred             CC--ceeeeeecceeeccC--CCCceeeEEEEeeeEEeecccccccCccee
Q psy11196         19 PR--EHIYSQYQGVWNLSS--DQGNLGVFMATNIRVVCTYYFGEVMDLGFD   65 (103)
Q Consensus        19 p~--E~i~~~~~gV~nls~--dqG~lG~f~vTNiRiV~~~W~a~~~~~~fN   65 (103)
                      |.  |.|.-+-++|.=.=+  +.| .|+|+||+=|++   |++. ...+|.
T Consensus        19 p~~~E~v~~~~~~v~~vl~~~~~g-~GtLyVt~~~l~---W~~~-~~~G~s   64 (171)
T 1zyi_A           19 PGSAEGLRQQQPETEAVLNGKGLG-TGTLYIAESRLS---WLDG-SGLGFS   64 (171)
T ss_dssp             --------CBCSSEEEEBTTBCCC-EEEEEEETTEEE---EEET-TTEEEE
T ss_pred             CCCccceEEEcCCEEEEECCCccC-ceEEEEEcCEEE---EEcC-CCcEEE
Confidence            56  666655554432221  223 799999999999   9999 888888

>3b77_A Uncharacterized protein; pleckstrin-homology domain, structural genomics, joint cente structural genomics, JCSG; 2.42A {Exiguobacterium sibiricum} SCOP: b.55.1.13

Probab=82.86  E-value=0.73  Score=35.12  Aligned_cols=28  Identities=25%  Similarity=0.250  Sum_probs=24.5

Q ss_pred             ccccC--CceeeeeecceeeccCCCCceeeEEEEeeeEE
Q psy11196         15 LKMLP--REHIYSQYQGVWNLSSDQGNLGVFMATNIRVV   51 (103)
Q Consensus        15 L~llp--~E~i~~~~~gV~nls~dqG~lG~f~vTNiRiV   51 (103)
                      =-|+|  ||+|+..+..+-|         .+++||-|+|
T Consensus        20 ~~Ll~EdgE~I~~ayk~irD---------~~iFTNkRLI   49 (193)
T 3b77_A           20 DYVLHEDGEKIYFLIKSKTD---------EYCFTNLALV   49 (193)
T ss_dssp             GGSCGGGTCCEEEEEECSSE---------EEEEESSEEE
T ss_pred             hhccCCCCCEEeeeeeecce---------EEEEecchhh
Confidence            35789  9999999998855         7999999999

---

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	103
d1zsqa1	125	Myotubularin-related protein 2, N-terminal domain	97.46
d2caya1	130	Vacuolar protein sorting protein 36, VPS36 {Baker'	94.77
d2hthb1	129	Vacuolar protein sorting protein 36, VPS36 {Human	93.38
d3dcxa1	116	Uncharacterized protein Shew0819 {Shewanella loihi	91.26

>d1zsqa1 b.55.1.8 (A:74-198) Myotubularin-related protein 2, N-terminal domain {Human (Homo sapiens) [TaxId: 9606]}

class: All beta proteins
fold: PH domain-like barrel
superfamily: PH domain-like
family: GRAM domain
domain: Myotubularin-related protein 2, N-terminal domain
species: Human (Homo sapiens) [TaxId: 9606]

Probab=97.46  E-value=0.00022  Score=47.97  Aligned_cols=81  Identities=23%  Similarity=0.288  Sum_probs=53.1

Q ss_pred             cccccCCceeeeeecceeeccCCCC-ceeeEEEEeeeEEeecccccccCccee-------ee-e---------e--cccc
Q psy11196         14 ELKMLPREHIYSQYQGVWNLSSDQG-NLGVFMATNIRVVCTYYFGEVMDLGFD-------IR-T---------Y--YFGE   73 (103)
Q Consensus        14 ~L~llp~E~i~~~~~gV~nls~dqG-~lG~f~vTNiRiV~~~W~a~~~~~~fN-------ir-~---------S--KFG~   73 (103)
                      +..|||||.++...++|-=...-.| --|++.+||-|++   --+...++.+.       |. -         +  .+|.
T Consensus         2 ~~~llpGE~i~~~~~~V~~~~p~~~~~~G~L~lTnyrli---F~~~~~~~~~~~~ipl~~I~~v~k~~~~~~~~~~~~~L   78 (125)
T d1zsqa1           2 EPPLLPGENIKDMAKDVTYICPFTGAVRGTLTVTNYRLY---FKSMERDPPFVLDASLGVINRVEKIGGASSRGENSYGL   78 (125)
T ss_dssp             CCCCCTTCCEEEEEEEEEEEETTTEEEEEEEEEESSEEE---EEESSSSSCEEEEEEGGGEEEEEEECCTTCCSTTCSEE
T ss_pred             CCCcCCCcEEEEeeccEEEEccCCCccceEEEEEeeEEE---EEcCCCCccEEEEeccceeeeeeecccccccCCccccE
Confidence            4579999999999999955543333 4699999999999   65542343332       11 1         1  1233


Q ss_pred             eE------EEEEecChh-hHHHHHHHHHHHh
Q psy11196         74 VM------DLGFDIRYE-SFENEVYHELFDS   97 (103)
Q Consensus        74 al------vLGFridP~-e~l~~~~ke~~~L   97 (103)
                      -+      ++-|..+|+ +.-+++++-++.+
T Consensus        79 ~I~CKDfr~~~f~f~~e~~~~~~v~~~L~~~  109 (125)
T d1zsqa1          79 ETVCKDIRNLRFAHKPEGRTRRSIFENLMKY  109 (125)
T ss_dssp             EEEETTTEEEEEECCGGGCCHHHHHHHHHHH
T ss_pred             EEEeccCeEEEEEEcCCccHHHHHHHHHHHh
Confidence            22      788888776 4456777777665

>d2caya1 b.55.1.12 (A:1-99,A:252-282) Vacuolar protein sorting protein 36, VPS36 {Baker's yeast (Saccharomyces cerevisiae) [TaxId: 4932]}

class: All beta proteins
fold: PH domain-like barrel
superfamily: PH domain-like
family: VPS36 N-terminal domain-like
domain: Vacuolar protein sorting protein 36, VPS36
species: Baker's yeast (Saccharomyces cerevisiae) [TaxId: 4932]

Probab=94.77  E-value=0.009  Score=40.62  Aligned_cols=53  Identities=13%  Similarity=0.068  Sum_probs=40.1

Q ss_pred             eeecCcccccCCceeeeeecceeeccCCC----CceeeEEEEeeeEEeecccccccCccee
Q psy11196          9 IVRRKELKMLPREHIYSQYQGVWNLSSDQ----GNLGVFMATNIRVVCTYYFGEVMDLGFD   65 (103)
Q Consensus         9 ii~~~~L~llp~E~i~~~~~gV~nls~dq----G~lG~f~vTNiRiV~~~W~a~~~~~~fN   65 (103)
                      +...|+-.|.|||.++-..+||-=--+++    -..|++++||-|||   |... .++.++
T Consensus         9 ~t~sg~p~L~~~E~~l~~qd~V~Ly~g~~k~~~~~~G~l~LTshRiI---~v~~-~~~~~~   65 (130)
T d2caya1           9 TTSSGQPLLREGEKDIFIDQSVGLYHGKSKILQRQRGRIFLTSQRII---YIDD-AKPTQN   65 (130)
T ss_dssp             BCTTSCBCCCTTCCEEEEEEEEEEEETTEECTTCSSEEEEEESSEEE---EEES-SCHHHH
T ss_pred             cccCCCeEccCCeEEEEEeCCEEEEECCccCccccCceEEEEeeEEE---EecC-CCCcee
Confidence            34578889999999999999885333322    23799999999999   9887 555443

>d2hthb1 b.55.1.12 (B:3-131) Vacuolar protein sorting protein 36, VPS36 {Human (Homo sapiens) [TaxId: 9606]}

class: All beta proteins
fold: PH domain-like barrel
superfamily: PH domain-like
family: VPS36 N-terminal domain-like
domain: Vacuolar protein sorting protein 36, VPS36
species: Human (Homo sapiens) [TaxId: 9606]

Probab=93.38  E-value=0.035  Score=37.06  Aligned_cols=39  Identities=15%  Similarity=0.213  Sum_probs=30.2

Q ss_pred             cccCCceeeeeeccee----eccCCCCceeeEEEEeeeEEeeccccc
Q psy11196         16 KMLPREHIYSQYQGVW----NLSSDQGNLGVFMATNIRVVCTYYFGE   58 (103)
Q Consensus        16 ~llp~E~i~~~~~gV~----nls~dqG~lG~f~vTNiRiV~~~W~a~   58 (103)
                      .|.|+|.++-+-+||-    +-+.+ =.-|++++||-|+|   |-..
T Consensus         8 ll~~~E~~l~~~~~V~lydG~~k~~-~~~G~l~LTthRli---~~~~   50 (129)
T d2hthb1           8 LLEINETLVIQQRGVRIYDGEEKIK-FDAGTLLLSTHRLI---WRDQ   50 (129)
T ss_dssp             CCCTTCCEEEEEEEEEEEETTCSSC-CCCEEEEEESSEEE---EEET
T ss_pred             CCCCCcEEEEEeCCEEeecCCCccc-ccceEEEEEeceEE---EecC
Confidence            3789999999999996    33332 13599999999999   9644

>d3dcxa1 b.55.1.13 (A:9-124) Uncharacterized protein Shew0819 {Shewanella loihica [TaxId: 359303]}

class: All beta proteins
fold: PH domain-like barrel
superfamily: PH domain-like
family: BPHL domain
domain: Uncharacterized protein Shew0819
species: Shewanella loihica [TaxId: 359303]

Probab=91.26  E-value=0.05  Score=36.67  Aligned_cols=29  Identities=24%  Similarity=0.301  Sum_probs=24.8

Q ss_pred             cccccCCceeeeeecceeeccCCCCceeeEEEEeeeEE
Q psy11196         14 ELKMLPREHIYSQYQGVWNLSSDQGNLGVFMATNIRVV   51 (103)
Q Consensus        14 ~L~llp~E~i~~~~~gV~nls~dqG~lG~f~vTNiRiV   51 (103)
                      +--|+|||+|+..+..+         ...+++||-|+|
T Consensus        16 ~~lL~~gE~i~~~yk~~---------rD~~~FTnkRlI   44 (116)
T d3dcxa1          16 GPIMGDNEQLALAYRVI---------RDMFVFTNKRLI   44 (116)
T ss_dssp             GGGSCTTCCEEEEEEET---------TEEEEEESSEEE
T ss_pred             HHhcCCCCEEEEEEEee---------cceEEEecceEE
Confidence            34589999999999888         556889999999