IMP GPI Biosynthesis Report IMP-Bioinformatics
GPI Biosynthesis Main Page  
GPI Site Motif   GPI Site Prediction   Home Page B.E.  

GPI Anchor Biosynthesis Report: T38153 (PIG-P family, Schizosaccharomyces pombe)




BLASTP 2.1.1 [Aug-8-2000]

Reference:
Altschul, Stephen F., Thomas L. Madden, Alejandro A. Schäffer, 
Jinghui Zhang, Zheng Zhang, Webb Miller, and David J. Lipman (1997), 
"Gapped BLAST and PSI-BLAST: a new generation of protein database search
programs",  Nucleic Acids Res. 25:3389-3402.

Query= 
         (157 letters)

Database: nr
           887,402 sequences; 277,845,442 total letters

Searching..................................................

Converged !!!


Results of PSI-Blast iteration 2

Distribution of 15 Blast Hits on the Query Sequence


Sequences with E-value BETTER than threshold (0.002)
Sequences producing significant alignments:
pir||T38153 hypothetical protein SPAC22A12.13 - fission yeast (... 262 3e-69
emb|CAB16583.2| (Z99295) human downs syndrome critical region-li... 195 2e-49
ref|NP_062416.1| (NM_019543) Down syndrome critical region prote... 147 1e-34
sp|P57054|DSR5_HUMAN Down syndrome critical region protein 5 (Do... 146 1e-34
dbj|BAA95636.1| (AB035745) DSCR5b [Homo sapiens] >gi|9963906|gb|... 146 2e-34
pir||JC7302 Down syndrome critical region protein 5 beta - human 146 2e-34
pir||JC7301 Down syndrome critical region protein 5 alpha - huma... 145 4e-34
ref|NP_176323.1| (NM_104809) hypothetical protein [Arabidopsis t... 126 2e-28
ref|NP_057514.1| (NM_016430) Down syndrome critical region prote... 115 4e-25
dbj|BAA95512.1| (AP001727) Down syndrome critical region protein... 114 6e-25
ref|NP_010725.1| (NC_001136) Protein required for cell viability... 109 3e-23
ref|NP_496845.1| (NM_064444) BR-1 protein-like [Caenorhabditis e... 99 6e-20
pir||T27023 hypothetical protein Y48E1B.2a - Caenorhabditis elegans 97 2e-19
ref|NP_181476.1| (NM_129502) hypothetical protein [Arabidopsis t... 64 1e-09
Sequences with E-value WORSE than threshold

gb|AAF56563.1| (AE003754) CG14550 gene product [Drosophila melan... 39 0.060
Alignments
>pir||T38153 hypothetical protein SPAC22A12.13 - fission yeast
           (Schizosaccharomyces pombe)
          Length = 157

 Score =  262 bits (663), Expect = 3e-69
 Identities = 157/157 (100%), Positives = 157/157 (100%)

Query: 1   MLVPPSPSLTTMFRTRSTTPDEEGTAIAEIDQQDWDIMVKVPTYEYYGFVMYLVSMLGFG 60
           MLVPPSPSLTTMFRTRSTTPDEEGTAIAEIDQQDWDIMVKVPTYEYYGFVMYLVSMLGFG
Sbjct: 1   MLVPPSPSLTTMFRTRSTTPDEEGTAIAEIDQQDWDIMVKVPTYEYYGFVMYLVSMLGFG 60

Query: 61  VYIVWALTPAPVLKFFEIHYYLSRWWALAIPTWLFVLVIYIHVVLNAYNTEVLTKPFSSL 120
           VYIVWALTPAPVLKFFEIHYYLSRWWALAIPTWLFVLVIYIHVVLNAYNTEVLTKPFSSL
Sbjct: 61  VYIVWALTPAPVLKFFEIHYYLSRWWALAIPTWLFVLVIYIHVVLNAYNTEVLTKPFSSL 120

Query: 121 ECIVDQYALVGEEDGAAHGRVVDLRLCDVNKQQLEET 157
           ECIVDQYALVGEEDGAAHGRVVDLRLCDVNKQQLEET
Sbjct: 121 ECIVDQYALVGEEDGAAHGRVVDLRLCDVNKQQLEET 157
>emb|CAB16583.2| (Z99295) human downs syndrome critical region-like
           [Schizosaccharomyces pombe]
          Length = 120

 Score =  195 bits (493), Expect = 2e-49
 Identities = 120/120 (100%), Positives = 120/120 (100%)

Query: 38  MVKVPTYEYYGFVMYLVSMLGFGVYIVWALTPAPVLKFFEIHYYLSRWWALAIPTWLFVL 97
           MVKVPTYEYYGFVMYLVSMLGFGVYIVWALTPAPVLKFFEIHYYLSRWWALAIPTWLFVL
Sbjct: 1   MVKVPTYEYYGFVMYLVSMLGFGVYIVWALTPAPVLKFFEIHYYLSRWWALAIPTWLFVL 60

Query: 98  VIYIHVVLNAYNTEVLTKPFSSLECIVDQYALVGEEDGAAHGRVVDLRLCDVNKQQLEET 157
           VIYIHVVLNAYNTEVLTKPFSSLECIVDQYALVGEEDGAAHGRVVDLRLCDVNKQQLEET
Sbjct: 61  VIYIHVVLNAYNTEVLTKPFSSLECIVDQYALVGEEDGAAHGRVVDLRLCDVNKQQLEET 120
>ref|NP_062416.1| (NM_019543) Down syndrome critical region protein c [Mus musculus]
 gb|AAF32294.1|AF216306_1 (AF216306) DCRC-1 [Mus musculus]
 gb|AAF32295.1|AF216307_1 (AF216307) DCRC-2 [Mus musculus]
 gb|AAF32296.1|AF216308_1 (AF216308) DCRC-3 [Mus musculus]
 gb|AAF32297.1|AF216309_1 (AF216309) DCRC-4 [Mus musculus]
 dbj|BAB23204.1| (AK004169) Down syndrome critical region gene c~data source:MGD,
           source key:MGI:1860433, evidence:ISS~putative [Mus
           musculus]
          Length = 132

 Score =  147 bits (369), Expect = 1e-34
 Identities = 32/88 (36%), Positives = 50/88 (56%), Gaps = 1/88 (1%)

Query: 41  VPTYEYYGFVMYLVSMLGFGVYIVWALTPAPVLKFFEIHYYLSRWWALAIPTWLFVLVIY 100
           +P    YGFV++L S  GF +Y+VWA  P   L    + Y+  ++WA+A+P +L + V+ 
Sbjct: 9   LPERAIYGFVLFLSSQFGFILYLVWAFVPESWLNSLGLTYWPQKYWAVALPVYLLITVVI 68

Query: 101 IHVVLNAYNTEVLTKPFSSLECIVDQYA 128
            +V+L   N  + T P  S+  I D YA
Sbjct: 69  GYVLLFGINM-MSTSPLDSIHTITDNYA 95
>sp|P57054|DSR5_HUMAN Down syndrome critical region protein 5 (Down syndrome critical
           region protein C)
 dbj|BAA96872.1| (AB037163) DSCR5b [Homo sapiens]
          Length = 158

 Score =  146 bits (367), Expect = 1e-34
 Identities = 31/88 (35%), Positives = 50/88 (56%), Gaps = 1/88 (1%)

Query: 41  VPTYEYYGFVMYLVSMLGFGVYIVWALTPAPVLKFFEIHYYLSRWWALAIPTWLFVLVIY 100
           +P    YGFV++L S  GF +Y+VWA  P   L    + Y+  ++WA+A+P +L + ++ 
Sbjct: 33  LPERAIYGFVLFLSSQFGFILYLVWAFIPESWLNSLGLTYWPQKYWAVALPVYLLIAIVI 92

Query: 101 IHVVLNAYNTEVLTKPFSSLECIVDQYA 128
            +V+L   N  + T P  S+  I D YA
Sbjct: 93  GYVLLFGINM-MSTSPLDSIHTITDNYA 119
>dbj|BAA95636.1| (AB035745) DSCR5b [Homo sapiens]
 gb|AAG09757.1|AF237812_1 (AF237812) NPD010 [Homo sapiens]
          Length = 158

 Score =  146 bits (367), Expect = 2e-34
 Identities = 31/88 (35%), Positives = 50/88 (56%), Gaps = 1/88 (1%)

Query: 41  VPTYEYYGFVMYLVSMLGFGVYIVWALTPAPVLKFFEIHYYLSRWWALAIPTWLFVLVIY 100
           +P    YGFV++L S  GF +Y+VWA  P   L    + Y+  ++WA+A+P +L + ++ 
Sbjct: 33  LPERAIYGFVLFLSSQFGFILYLVWAFIPESWLNSLGLTYWPQKYWAVALPVYLLIAIVI 92

Query: 101 IHVVLNAYNTEVLTKPFSSLECIVDQYA 128
            +V+L   N  + T P  S+  I D YA
Sbjct: 93  GYVLLFGINM-MSTSPLDSIHTITDNYA 119
>pir||JC7302 Down syndrome critical region protein 5 beta - human
          Length = 158

 Score =  146 bits (367), Expect = 2e-34
 Identities = 31/88 (35%), Positives = 50/88 (56%), Gaps = 1/88 (1%)

Query: 41  VPTYEYYGFVMYLVSMLGFGVYIVWALTPAPVLKFFEIHYYLSRWWALAIPTWLFVLVIY 100
           +P    YGFV++L S  GF +Y+VWA  P   L    + Y+  ++WA+A+P +L + ++ 
Sbjct: 33  LPERAIYGFVLFLSSQFGFILYLVWAFIPESWLNSLGLTYWPQKYWAVALPVYLLIAIVI 92

Query: 101 IHVVLNAYNTEVLTKPFSSLECIVDQYA 128
            +V+L   N  + T P  S+  I D YA
Sbjct: 93  GYVLLFGLNM-MSTSPLDSIHTITDNYA 119
>pir||JC7301 Down syndrome critical region protein 5 alpha - human
 dbj|BAA95633.1| (AB035742) DSCR5a1 [Homo sapiens]
 dbj|BAA95634.1| (AB035743) DSCR5a2 [Homo sapiens]
 dbj|BAA95635.1| (AB035744) DSCR5a3 [Homo sapiens]
 dbj|BAA96871.1| (AB037162) DSCR5a [Homo sapiens]
 dbj|BAB12395.1| (AB039659) PIG-P [Homo sapiens]
 gb|AAH05180.1|AAH05180 (BC005180) Unknown (protein for MGC:12168) [Homo sapiens]
 gb|AAH11007.1|AAH11007 (BC011007) Unknown (protein for MGC:13327) [Homo sapiens]
          Length = 134

 Score =  145 bits (364), Expect = 4e-34
 Identities = 31/88 (35%), Positives = 50/88 (56%), Gaps = 1/88 (1%)

Query: 41  VPTYEYYGFVMYLVSMLGFGVYIVWALTPAPVLKFFEIHYYLSRWWALAIPTWLFVLVIY 100
           +P    YGFV++L S  GF +Y+VWA  P   L    + Y+  ++WA+A+P +L + ++ 
Sbjct: 9   LPERAIYGFVLFLSSQFGFILYLVWAFIPESWLNSLGLTYWPQKYWAVALPVYLLIAIVI 68

Query: 101 IHVVLNAYNTEVLTKPFSSLECIVDQYA 128
            +V+L   N  + T P  S+  I D YA
Sbjct: 69  GYVLLFGINM-MSTSPLDSIHTITDNYA 95
>ref|NP_176323.1| (NM_104809) hypothetical protein [Arabidopsis thaliana]
 sp|O64792|DSR5_ARATH Down syndrome critical region protein 5 homolog
 gb|AAC13913.1|AAC13913 (AC004255) T1F9.23 [Arabidopsis thaliana]
          Length = 137

 Score =  126 bits (315), Expect = 2e-28
 Identities = 26/86 (30%), Positives = 46/86 (53%), Gaps = 1/86 (1%)

Query: 43  TYEYYGFVMYLVSMLGFGVYIVWALTPAPVLKFFEIHYYLSRWWALAIPTWLFVLVIYIH 102
           T E YGFV  +  ++   ++++W   P   L+   I+YY S++WA+A+P +  V ++   
Sbjct: 13  TSEVYGFVGSISIVVATVIFLIWGYVPDKFLESIGIYYYPSKYWAMAMPMYSMVTLLVAL 72

Query: 103 VVLNAYNTEVLTKPFSSLECIVDQYA 128
           V     N    +KP +SL  + D Y+
Sbjct: 73  VFYIGLNFMSTSKP-TSLNTLFDDYS 97
>ref|NP_057514.1| (NM_016430) Down syndrome critical region protein 5; Down syndrome
           critical region protein C short form; Down syndrome
           critical region protein C [Homo sapiens]
 gb|AAF32289.1|AF216305_1 (AF216305) Down syndrome critical region protein C short form [Homo
           sapiens]
 dbj|BAA96873.1| (AB037164) DSCR5c [Homo sapiens]
          Length = 108

 Score =  115 bits (287), Expect = 4e-25
 Identities = 23/69 (33%), Positives = 38/69 (54%), Gaps = 1/69 (1%)

Query: 60  GVYIVWALTPAPVLKFFEIHYYLSRWWALAIPTWLFVLVIYIHVVLNAYNTEVLTKPFSS 119
            +Y+VWA  P   L    + Y+  ++WA+A+P +L + V+  +V+L   N  + T P  S
Sbjct: 2   VLYLVWAFIPESWLNSLGLTYWPQKYWAVALPVYLLIAVVIGYVLLFGINM-MSTSPLDS 60

Query: 120 LECIVDQYA 128
           +  I D YA
Sbjct: 61  IHTITDNYA 69
>dbj|BAA95512.1| (AP001727) Down syndrome critical region protein C [Homo sapiens]
          Length = 108

 Score =  114 bits (285), Expect = 6e-25
 Identities = 22/69 (31%), Positives = 38/69 (54%), Gaps = 1/69 (1%)

Query: 60  GVYIVWALTPAPVLKFFEIHYYLSRWWALAIPTWLFVLVIYIHVVLNAYNTEVLTKPFSS 119
            +Y+VWA  P   L    + Y+  ++WA+A+P +L + ++  +V+L   N  + T P  S
Sbjct: 2   VLYLVWAFIPESWLNSLGLTYWPQKYWAVALPVYLLIAIVIGYVLLFGINM-MSTSPLDS 60

Query: 120 LECIVDQYA 128
           +  I D YA
Sbjct: 61  IHTITDNYA 69
>ref|NP_010725.1| (NC_001136) Protein required for cell viability; Ydr437wp
           [Saccharomyces cerevisiae]
 pir||S69717 hypothetical protein YDR437w - yeast (Saccharomyces cerevisiae)
 gb|AAB64854.1| (U33007) Ydr437wp; CAI: 0.10 [Saccharomyces cerevisiae]
          Length = 140

 Score =  109 bits (271), Expect = 3e-23
 Identities = 20/91 (21%), Positives = 41/91 (44%), Gaps = 9/91 (9%)

Query: 43  TYEYYGFVMYLVSMLGFGVYIVWALTPAP--------VLKFFEIHYYLSRWWALAIPTWL 94
           T EYY F  Y++      + I+W++ P+          +    +  +  R W + + + +
Sbjct: 3   TKEYYWFSQYMIITSTLVLTIIWSILPSSLGEAAPKQFINTL-LDIFPQRRWIITLESIM 61

Query: 95  FVLVIYIHVVLNAYNTEVLTKPFSSLECIVD 125
            + ++  ++ L  YN + LT P  SL  + D
Sbjct: 62  LMGMLCTYIGLLMYNEDTLTPPLDSLSTVTD 92
>ref|NP_496845.1| (NM_064444) BR-1 protein-like [Caenorhabditis elegans]
 emb|CAB07689.2| (Z93393) contains similarity to Pfam domain: PF01926 (GTPase of
           unknown function), E-value=0.00012,~cDNA EST
           EMBL:AU111138 comes from this gene~cDNA EST
           EMBL:AU114976 comes from this gene [Caenorhabditis
           elegans]
          Length = 890

 Score = 98.6 bits (243), Expect = 6e-20
 Identities = 25/62 (40%), Positives = 36/62 (57%), Gaps = 1/62 (1%)

Query: 42  PTYEYYGFVMYLVSMLGFGVYIVWALTPAPVLKFFEIHYYLSRWWALAIPTWLFVLVIYI 101
           P    YGF +Y+VS   F +Y++WA+TP P+L    I Y  S+ WALAI  + F     +
Sbjct: 30  PARGIYGFALYIVSWTLFVIYLIWAITPVPILYRLGITYIPSKLWALAIGIF-FPTAACL 88

Query: 102 HV 103
           +V
Sbjct: 89  YV 90
>pir||T27023 hypothetical protein Y48E1B.2a - Caenorhabditis elegans
          Length = 1768

 Score = 96.6 bits (238), Expect = 2e-19
 Identities = 25/62 (40%), Positives = 36/62 (57%), Gaps = 1/62 (1%)

Query: 42  PTYEYYGFVMYLVSMLGFGVYIVWALTPAPVLKFFEIHYYLSRWWALAIPTWLFVLVIYI 101
           P    YGF +Y+VS   F +Y++WA+TP P+L    I Y  S+ WALAI  + F     +
Sbjct: 30  PARGIYGFALYIVSWTLFVIYLIWAITPVPILYRLGITYIPSKLWALAIGIF-FPTAACL 88

Query: 102 HV 103
           +V
Sbjct: 89  YV 90
>ref|NP_181476.1| (NM_129502) hypothetical protein [Arabidopsis thaliana]
 pir||T00554 hypothetical protein F12L6.10 - Arabidopsis thaliana
 gb|AAC27835.1| (AC004218) hypothetical protein [Arabidopsis thaliana]
          Length = 773

 Score = 64.3 bits (155), Expect = 1e-09
 Identities = 11/39 (28%), Positives = 20/39 (51%), Gaps = 1/39 (2%)

Query: 90  IPTWLFVLVIYIHVVLNAYNTEVLTKPFSSLECIVDQYA 128
           +PT+L V ++         N  + T P +SL  + D+Y+
Sbjct: 1   MPTYLMVTLMLALAFYIGLNF-IATPPPTSLNTLFDEYS 38
>gb|AAF56563.1| (AE003754) CG14550 gene product [Drosophila melanogaster]
          Length = 138

 Score = 38.6 bits (89), Expect = 0.060
 Identities = 18/75 (24%), Positives = 35/75 (46%), Gaps = 7/75 (9%)

Query: 53  LVSMLGFGVYIVWALTPAPVLKFFEIH-YYLSRWWALAIPTWLFV-LVIYIHVVLNAYNT 110
           ++  + F VY+ WAL P        +H Y   +++A+ +P  + V    +  ++  A N 
Sbjct: 1   MLFTVLFLVYVTWALLPVEF----GLHSYLPDKYFAVFVPFLVLVFAWFFAFLIYPAINL 56

Query: 111 EVLTKPFSSLECIVD 125
            + T    S+  +VD
Sbjct: 57  SM-TVDVDSIASVVD 70
CPU time:    44.37 user secs.	    1.73 sys. secs	   46.10 total secs.

  Database: nr
    Posted date:  Apr 21, 2002  2:19 PM
  Number of letters in database: 277,845,442
  Number of sequences in database:  887,402
  
Lambda     K      H
   0.319    0.160    0.492 

Gapped
Lambda     K      H
   0.270   0.0560    0.230 


Matrix: BLOSUM62
Gap Penalties: Existence: 11, Extension: 1
Number of Hits to DB: 106723385
Number of Sequences: 887402
Number of extensions: 5100822
Number of successful extensions: 18670
Number of sequences better than 10.0: 49
Number of HSP's better than 10.0 without gapping: 16
Number of HSP's successfully gapped in prelim test: 33
Number of HSP's that attempted gapping in prelim test: 18637
Number of HSP's gapped (non-prelim): 53
length of query: 157
length of database: 277,845,442
effective HSP length: 49
effective length of query: 108
effective length of database: 234,362,744
effective search space: 25311176352
effective search space used: 25311176352
T: 11
A: 40
X1: 15 ( 6.9 bits)
X2: 38 (14.8 bits)
X3: 64 (24.9 bits)
S1: 40 (21.1 bits)
S2: 70 (31.4 bits)