analysis of sequence from tem16.fa
~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~

>gi|15987493|gb|AAL11993.1|AF378756_1 tumor endothelial marker 6 [Homo sapiens]
MVKCYHKKYRSATRDVIFRLQFHTGAVQGYGLVFGKEDLDNASKDDRFPDYGKVELVFSATPEKIQGSEHLYNDHGVIVD
YNTTDPLIRWDSYENLSADGEVLHTQGPVDGSLYAKVRKKSSSDPGIPGGPQAIPATNSPDHSDHTLSVSSDSGHSTASA
RTDKTEERLAPGTRRGLSAQEKAELDQLLSGFGLEDPGSSLKEMTDARSKYSGTRHVVPAQVHVNGDAALKDRETDILDD
EMPHHDLHSVDSLGTLSSSEGPQSAHLGPFTCHKSSQNSLLSDGFGSNVGEDPQGTLVPDLGLGMDGPYERERTFGSREP
KQPQPLLRKPSVSAQMQAYGQSSYSTQTWVRQQQMVVAHQYSFAPDGEARLVSRCPADNPGLVQAQPRVPLTPTRGTSSR
VAVQRGVGSGPHPPDTQQPSPSKAFKPRFPGDQVVNGAGPELSTGPSPGSPTLDIDQSIEQLNRLILELDPTFEPIPTHM
NALGSQANGSVSPDSVGGGLRASSRLPDTGEGPSRATGRQGSSAEQPLGGRLRKLSLGQYDNDAGGQLPFSKCAWGKAGV
DYAPNLPPFPSPADVKETMTPGYPQDLDIIDGRILSSKESMCSTPAFPVSPETPYVKTALRHPPFSPPEPPLSSPASQHK
GGREPRSCPETLTHAVGMSESPIGPKSTMLRADASSTPSFQQAFASSCTISSNGPGQRRESSSSAERQWVESSPKPMVSL
LGSGRPTGSPLSAEFSGTRKDSPVLSCFPPSELQAPFHSHELSLAEPPDSLAPPSSQAFLGFGTAPVGSGLPPEEDLGAL
LANSHGASPTPSIPLTATGAADNGFLSHNFLTVAPGHSSHHSPGLQGQGVTLPGQPPLPEKKRASEGDRSLGSVSPSSSG
FSSPHSGSTISIPFPNVLPDFSKASEAASPLPDSPGDKLVIVKFVQDTSKFWYKADISREQAIAMLKDKEPGSFIVRDSH
SFRGAYGLAMKVATPPPSVLQLNKKAGDLANELVRHFLIECTPKGVRLKGCSNEPYFGSLTALVCQHSITPLALPCKLLI
PERDPLEEIAESSPQTAANSAAELLKQGAACNVWYLNSVEMESLTGHQAIQKALSITLVQEPPPVSTVVHFKVSAQGITL
TDNQRKLFFRRHYPVNSVIFCALDPQDRKWIKDGPSSKVFGFVARKQGSATDNVCHLFAEHDPEQPASAIVNFVSKVMIG
SPKKV
~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~

sec.str. with predator

> gi|15987493|gb|AAL11993.1|AF378756_1
              .         .         .         .         .
1    MVKCYHKKYRSATRDVIFRLQFHTGAVQGYGLVFGKEDLDNASKDDRFPD   50
     ___________HHHHHH__HHHHHEEEE__EEEEE_______________

              .         .         .         .         .
51   YGKVELVFSATPEKIQGSEHLYNDHGVIVDYNTTDPLIRWDSYENLSADG  100
     __EEEEEE__________________EEEE______EEE___________

              .         .         .         .         .
101  EVLHTQGPVDGSLYAKVRKKSSSDPGIPGGPQAIPATNSPDHSDHTLSVS  150
     EEEEE__________HHHHH__________________________EEEE

              .         .         .         .         .
151  SDSGHSTASARTDKTEERLAPGTRRGLSAQEKAELDQLLSGFGLEDPGSS  200
     _______________EEE_______HHHHHHHHHHHHHHH__________

              .         .         .         .         .
201  LKEMTDARSKYSGTRHVVPAQVHVNGDAALKDRETDILDDEMPHHDLHSV  250
     HHHHHHHHHH__________EEE____HHHHHHHH_______________

              .         .         .         .         .
251  DSLGTLSSSEGPQSAHLGPFTCHKSSQNSLLSDGFGSNVGEDPQGTLVPD  300
     __________________________________________________

              .         .         .         .         .
301  LGLGMDGPYERERTFGSREPKQPQPLLRKPSVSAQMQAYGQSSYSTQTWV  350
     ________________________________HHHHHHH_______HHHH

              .         .         .         .         .
351  RQQQMVVAHQYSFAPDGEARLVSRCPADNPGLVQAQPRVPLTPTRGTSSR  400
     HHHHHHHHHHH________EEEE________EEE________________

              .         .         .         .         .
401  VAVQRGVGSGPHPPDTQQPSPSKAFKPRFPGDQVVNGAGPELSTGPSPGS  450
     EEEEEE____________________________________________

              .         .         .         .         .
451  PTLDIDQSIEQLNRLILELDPTFEPIPTHMNALGSQANGSVSPDSVGGGL  500
     ___HHHHHHHHHHHHHHH________________________________

              .         .         .         .         .
501  RASSRLPDTGEGPSRATGRQGSSAEQPLGGRLRKLSLGQYDNDAGGQLPF  550
     __________________________________________________

              .         .         .         .         .
551  SKCAWGKAGVDYAPNLPPFPSPADVKETMTPGYPQDLDIIDGRILSSKES  600
     ___________________________________________EEE____

              .         .         .         .         .
601  MCSTPAFPVSPETPYVKTALRHPPFSPPEPPLSSPASQHKGGREPRSCPE  650
     ______________EEEE________________________________

              .         .         .         .         .
651  TLTHAVGMSESPIGPKSTMLRADASSTPSFQQAFASSCTISSNGPGQRRE  700
     ___EEE___________HHHHHH______HHHHHHH___________EEE

              .         .         .         .         .
701  SSSSAERQWVESSPKPMVSLLGSGRPTGSPLSAEFSGTRKDSPVLSCFPP  750
     ____HHHHHHH_____EEEE______________________________

              .         .         .         .         .
751  SELQAPFHSHELSLAEPPDSLAPPSSQAFLGFGTAPVGSGLPPEEDLGAL  800
     ___________EEEE___________EEEE_______________HHHHH

              .         .         .         .         .
801  LANSHGASPTPSIPLTATGAADNGFLSHNFLTVAPGHSSHHSPGLQGQGV  850
     HHH__________EEEE________________________________E

              .         .         .         .         .
851  TLPGQPPLPEKKRASEGDRSLGSVSPSSSGFSSPHSGSTISIPFPNVLPD  900
     EE____________________________________EEEEE_______

              .         .         .         .         .
901  FSKASEAASPLPDSPGDKLVIVKFVQDTSKFWYKADISREQAIAMLKDKE  950
     HHHHHHH__________EEEEEEEE_____EEE___HHHHHHHHHHHH__

              .         .         .         .         .
951  PGSFIVRDSHSFRGAYGLAMKVATPPPSVLQLNKKAGDLANELVRHFLIE 1000
     __EEEEE____________EEEE_____HHHHH___HHHHHHHHHHHHEE

              .         .         .         .         .
1001 CTPKGVRLKGCSNEPYFGSLTALVCQHSITPLALPCKLLIPERDPLEEIA 1050
     E____EEEE_________HHHHHHH__________EEEEEE____HHHHH

              .         .         .         .         .
1051 ESSPQTAANSAAELLKQGAACNVWYLNSVEMESLTGHQAIQKALSITLVQ 1100
     H____HHHHHHHHHHHH___EEEEEE_________HHHHHHHHHHHH___

              .         .         .         .         .
1101 EPPPVSTVVHFKVSAQGITLTDNQRKLFFRRHYPVNSVIFCALDPQDRKW 1150
     _____EEEEEEEE__________HHHHHHH______EEEEEE______EE

              .         .         .         .         .
1151 IKDGPSSKVFGFVARKQGSATDNVCHLFAEHDPEQPASAIVNFVSKVMIG 1200
     EE_______EEEEE__________EEEEE________EEEHHHHHHH___

          
1201 SPKKV                                              1205
     _____


~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~


method         :         1
alpha-contents :      11.3 %
beta-contents  :       3.4 %
coil-contents  :      85.4 %
class          : irregular


method         :         2
alpha-contents :       0.0 %
beta-contents  :       0.5 %
coil-contents  :      99.5 %
class          : irregular


~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~

GPI: learning from metazoa
-18.02  -5.20  -2.80  -0.25  -4.00   0.00   0.00   0.00  -0.15  -2.17  -0.58   0.00 -12.00   0.00   0.00   0.00  -45.17
  0.56   0.00   0.00   0.00   0.00   0.00 -16.00  -0.53  -0.36  -2.68  -0.58 -12.00   0.00   0.00   0.00   0.00  -31.59
ID: gi|15987493|gb|AAL11993.1|AF378756_1	AC: xxx Len: 1200 1:I  1187 Sc:  -31.59 Pv: 8.232018e-02 NO_GPI_SITE
GPI: learning from protozoa
-13.82  -0.26  -0.42  -0.20  -4.00   0.00 -12.00  -2.74  -0.03  -2.46  -3.20 -12.00 -12.00   0.00   0.00   0.00  -63.13
 -6.84   0.00   0.00  -0.01   0.00   0.00 -16.00  -1.18  -0.04  -2.46  -3.20 -12.00 -12.00   0.00   0.00   0.00  -53.73
ID: gi|15987493|gb|AAL11993.1|AF378756_1	AC: xxx Len: 1200 1:I  1187 Sc:  -53.73 Pv: 1.678742e-01 NO_GPI_SITE

~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~

# SignalP euk predictions
# name       Cmax  pos ?  Ymax  pos ?  Smax  pos ?  Smean ?
gi|15987493  0.270  68 N  0.137 1036 N  0.346 1027 N  0.056 N
# SignalP gram- predictions
# name       Cmax  pos ?  Ymax  pos ?  Smax  pos ?  Smean ?
gi|15987493  0.460 808 N  0.193 1034 N  0.906 1134 Y  0.117 N
# SignalP gram+ predictions
# name       Cmax  pos ?  Ymax  pos ?  Smax  pos ?  Smean ?
gi|15987493  0.552 161 Y  0.346 905 Y  0.961   5 Y  0.152 N

~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~

low complexity regions: SEG 12 2.2 2.5
>gi|15987493|gb|AAL11993.1|AF378756_1 tumor endothelial marker 6 [Homo sapiens]

                                  1-124  MVKCYHKKYRSATRDVIFRLQFHTGAVQGY
                                         GLVFGKEDLDNASKDDRFPDYGKVELVFSA
                                         TPEKIQGSEHLYNDHGVIVDYNTTDPLIRW
                                         DSYENLSADGEVLHTQGPVDGSLYAKVRKK
                                         SSSD
                  pgipggpqaipa  125-136  
                                137-622  TNSPDHSDHTLSVSSDSGHSTASARTDKTE
                                         ERLAPGTRRGLSAQEKAELDQLLSGFGLED
                                         PGSSLKEMTDARSKYSGTRHVVPAQVHVNG
                                         DAALKDRETDILDDEMPHHDLHSVDSLGTL
                                         SSSEGPQSAHLGPFTCHKSSQNSLLSDGFG
                                         SNVGEDPQGTLVPDLGLGMDGPYERERTFG
                                         SREPKQPQPLLRKPSVSAQMQAYGQSSYST
                                         QTWVRQQQMVVAHQYSFAPDGEARLVSRCP
                                         ADNPGLVQAQPRVPLTPTRGTSSRVAVQRG
                                         VGSGPHPPDTQQPSPSKAFKPRFPGDQVVN
                                         GAGPELSTGPSPGSPTLDIDQSIEQLNRLI
                                         LELDPTFEPIPTHMNALGSQANGSVSPDSV
                                         GGGLRASSRLPDTGEGPSRATGRQGSSAEQ
                                         PLGGRLRKLSLGQYDNDAGGQLPFSKCAWG
                                         KAGVDYAPNLPPFPSPADVKETMTPGYPQD
                                         LDIIDGRILSSKESMCSTPAFPVSPETPYV
                                         KTALRH
               ppfsppepplsspas  623-637  
                                638-700  QHKGGREPRSCPETLTHAVGMSESPIGPKS
                                         TMLRADASSTPSFQQAFASSCTISSNGPGQ
                                         RRE
                          ssss  701-704  
                                705-869  AERQWVESSPKPMVSLLGSGRPTGSPLSAE
                                         FSGTRKDSPVLSCFPPSELQAPFHSHELSL
                                         AEPPDSLAPPSSQAFLGFGTAPVGSGLPPE
                                         EDLGALLANSHGASPTPSIPLTATGAADNG
                                         FLSHNFLTVAPGHSSHHSPGLQGQGVTLPG
                                         QPPLPEKKRASEGDR
      slgsvspsssgfssphsgstisip  870-893  
                                894-1205 FPNVLPDFSKASEAASPLPDSPGDKLVIVK
                                         FVQDTSKFWYKADISREQAIAMLKDKEPGS
                                         FIVRDSHSFRGAYGLAMKVATPPPSVLQLN
                                         KKAGDLANELVRHFLIECTPKGVRLKGCSN
                                         EPYFGSLTALVCQHSITPLALPCKLLIPER
                                         DPLEEIAESSPQTAANSAAELLKQGAACNV
                                         WYLNSVEMESLTGHQAIQKALSITLVQEPP
                                         PVSTVVHFKVSAQGITLTDNQRKLFFRRHY
                                         PVNSVIFCALDPQDRKWIKDGPSSKVFGFV
                                         ARKQGSATDNVCHLFAEHDPEQPASAIVNF
                                         VSKVMIGSPKKV

low complexity regions: SEG 25 3.0 3.3
>gi|15987493|gb|AAL11993.1|AF378756_1 tumor endothelial marker 6 [Homo sapiens]

                                  1-135  MVKCYHKKYRSATRDVIFRLQFHTGAVQGY
                                         GLVFGKEDLDNASKDDRFPDYGKVELVFSA
                                         TPEKIQGSEHLYNDHGVIVDYNTTDPLIRW
                                         DSYENLSADGEVLHTQGPVDGSLYAKVRKK
                                         SSSDPGIPGGPQAIP
     atnspdhsdhtlsvssdsghstasa  136-160  
                                161-330  RTDKTEERLAPGTRRGLSAQEKAELDQLLS
                                         GFGLEDPGSSLKEMTDARSKYSGTRHVVPA
                                         QVHVNGDAALKDRETDILDDEMPHHDLHSV
                                         DSLGTLSSSEGPQSAHLGPFTCHKSSQNSL
                                         LSDGFGSNVGEDPQGTLVPDLGLGMDGPYE
                                         RERTFGSREPKQPQPLLRKP
svsaqmqaygqssystqtwvrqqqmvvahq  331-362  
                            ys
                                363-379  FAPDGEARLVSRCPADN
pglvqaqprvpltptrgtssrvavqrgvgs  380-422  
                 gphppdtqqpsps
                                423-602  KAFKPRFPGDQVVNGAGPELSTGPSPGSPT
                                         LDIDQSIEQLNRLILELDPTFEPIPTHMNA
                                         LGSQANGSVSPDSVGGGLRASSRLPDTGEG
                                         PSRATGRQGSSAEQPLGGRLRKLSLGQYDN
                                         DAGGQLPFSKCAWGKAGVDYAPNLPPFPSP
                                         ADVKETMTPGYPQDLDIIDGRILSSKESMC
stpafpvspetpyvktalrhppfsppeppl  603-649  
             sspasqhkggreprscp
                                650-741  ETLTHAVGMSESPIGPKSTMLRADASSTPS
                                         FQQAFASSCTISSNGPGQRRESSSSAERQW
                                         VESSPKPMVSLLGSGRPTGSPLSAEFSGTR
                                         KD
spvlscfppselqapfhshelslaeppdsl  742-827  
appssqaflgfgtapvgsglppeedlgall
    anshgasptpsipltatgaadngfls
                                828-834  HNFLTVA
pghsshhspglqgqgvtlpgqpplpekkra  835-895  
segdrslgsvspsssgfssphsgstisipf
                             p
                                896-1205 NVLPDFSKASEAASPLPDSPGDKLVIVKFV
                                         QDTSKFWYKADISREQAIAMLKDKEPGSFI
                                         VRDSHSFRGAYGLAMKVATPPPSVLQLNKK
                                         AGDLANELVRHFLIECTPKGVRLKGCSNEP
                                         YFGSLTALVCQHSITPLALPCKLLIPERDP
                                         LEEIAESSPQTAANSAAELLKQGAACNVWY
                                         LNSVEMESLTGHQAIQKALSITLVQEPPPV
                                         STVVHFKVSAQGITLTDNQRKLFFRRHYPV
                                         NSVIFCALDPQDRKWIKDGPSSKVFGFVAR
                                         KQGSATDNVCHLFAEHDPEQPASAIVNFVS
                                         KVMIGSPKKV

low complexity regions: SEG 45 3.4 3.75
>gi|15987493|gb|AAL11993.1|AF378756_1 tumor endothelial marker 6 [Homo sapiens]

                                  1-104  MVKCYHKKYRSATRDVIFRLQFHTGAVQGY
                                         GLVFGKEDLDNASKDDRFPDYGKVELVFSA
                                         TPEKIQGSEHLYNDHGVIVDYNTTDPLIRW
                                         DSYENLSADGEVLH
tqgpvdgslyakvrkksssdpgipggpqai  105-217  
patnspdhsdhtlsvssdsghstasartdk
teerlapgtrrglsaqekaeldqllsgfgl
       edpgsslkemtdarskysgtrhv
                                218-363  VPAQVHVNGDAALKDRETDILDDEMPHHDL
                                         HSVDSLGTLSSSEGPQSAHLGPFTCHKSSQ
                                         NSLLSDGFGSNVGEDPQGTLVPDLGLGMDG
                                         PYERERTFGSREPKQPQPLLRKPSVSAQMQ
                                         AYGQSSYSTQTWVRQQQMVVAHQYSF
apdgearlvsrcpadnpglvqaqprvpltp  364-454  
trgtssrvavqrgvgsgphppdtqqpspsk
afkprfpgdqvvngagpelstgpspgsptl
                             d
                                455-481  IDQSIEQLNRLILELDPTFEPIPTHMN
algsqangsvspdsvggglrassrlpdtge  482-574  
gpsratgrqgssaeqplggrlrklslgqyd
ndaggqlpfskcawgkagvdyapnlppfps
                           pad
                                575-658  VKETMTPGYPQDLDIIDGRILSSKESMCST
                                         PAFPVSPETPYVKTALRHPPFSPPEPPLSS
                                         PASQHKGGREPRSCPETLTHAVGM
sespigpkstmlradasstpsfqqafassc  659-916  
tissngpgqrressssaerqwvesspkpmv
sllgsgrptgsplsaefsgtrkdspvlscf
ppselqapfhshelslaeppdslappssqa
flgfgtapvgsglppeedlgallanshgas
ptpsipltatgaadngflshnfltvapghs
shhspglqgqgvtlpgqpplpekkrasegd
rslgsvspsssgfssphsgstisipfpnvl
            pdfskaseaasplpdspg
                                917-1205 DKLVIVKFVQDTSKFWYKADISREQAIAML
                                         KDKEPGSFIVRDSHSFRGAYGLAMKVATPP
                                         PSVLQLNKKAGDLANELVRHFLIECTPKGV
                                         RLKGCSNEPYFGSLTALVCQHSITPLALPC
                                         KLLIPERDPLEEIAESSPQTAANSAAELLK
                                         QGAACNVWYLNSVEMESLTGHQAIQKALSI
                                         TLVQEPPPVSTVVHFKVSAQGITLTDNQRK
                                         LFFRRHYPVNSVIFCALDPQDRKWIKDGPS
                                         SKVFGFVARKQGSATDNVCHLFAEHDPEQP
                                         ASAIVNFVSKVMIGSPKKV


low complexity regions: XNU
# Score cutoff = 21, Search from offsets 1 to 4
# both members of each repeat flagged
# lambda = 0.347, K = 0.200, H = 0.664
>gi|15987493|gb|AAL11993.1|AF378756_1 tumor endothelial marker 6 [Homo sapiens]
MVKCYHKKYRSATRDVIFRLQFHTGAVQGYGLVFGKEDLDNASKDDRFPDYGKVELVFSA
TPEKIQGSEHLYNDHGVIVDYNTTDPLIRWDSYENLSADGEVLHTQGPVDGSLYAKVRKK
SSSDpgipggpqaipatnspdhsdhtlsvssdsghstasartdkteERLAPGTRRGLSAQ
EKAELDQLLSGFGLEDPGSSLKEMTDARSKYSGTRHVVPAQVHVNGDAALKDRETDILDD
EMPHhdlhsvdslgtlsssegpqsahLGPFTCHKSSQNSLLSDGFGSNVGEDPQGTLVPD
LGLGMDGPYERERTFGSREPKQPQPLLRKPSVSAQMQAYGQSSYSTQTWVRQQQMVVAHQ
YSFAPDGEARLVSRCPADNPGLVQAQPRVPLTPTRGTSSRVAVQRGVGSGPHPPDTQQPS
PSKAFKPRFPGDQVVNGAGPELSTGPSPGSPTLDIDQSIEQLNRLILELDPTFEPIPTHM
NALGSQANGSVSPDSVGGGLRASSRLPDTGEGPSRATGRQGSSAEQPLGGRLRKLSLGQY
DNDAGGQLPFSKCAWGKAGVDYAPNLPPFPSPADVKETMTPGYPQDLDIIDGRILSSKES
MCSTPAFPVSPETPYVKTALRHPPFSPPEPPLSSPASQHKGGREPRSCPETLTHAVGMSE
SPIGPKSTMLRADASSTPSFQQAFASSCTISSNGPGQRRESSSSAERQWVESSPKPMVSL
LGSGRPTGSPLSAEFSGTRKDSPVLSCFPPSELQAPFHSHELSLAEPPDSLAPPSSQAFL
GFGTAPVGSGLPPEEDLGALLANSHGASPTPSIPLTATGAADNGFLSHNFLTVAPGHSSH
HSPGLQGQGVTLPGQPPLPEKKRASEGDRSLGSVSPSSSGFSSPHSGSTISIPFPNVLPD
FSKASEAASPLPDSPGDKLVIVKFVQDTSKFWYKADISREQAIAMLKDKEPGSFIVRDSH
SFRGAYGLAMKVATPPPSVLQLNKKAGDLANELVRHFLIECTPKGVRLKGCSNEPYFGSL
TALVCQHSITPLALPCKLLIPERDPLEEIAESSPQTAANSAAELLKQGAACNVWYLNSVE
MESLTGHQAIQKALSITLVQEPPPVSTVVHFKVSAQGITLTDNQRKLFFRRHYPVNSVIF
CALDPQDRKWIKDGPSSKVFGFVARKQGSATDNVCHLFAEHDPEQPASAIVNFVSKVMIG
SPKKV
    1 -  124 MVKCYHKKYR SATRDVIFRL QFHTGAVQGY GLVFGKEDLD NASKDDRFPD YGKVELVFSA 
             TPEKIQGSEH LYNDHGVIVD YNTTDPLIRW DSYENLSADG EVLHTQGPVD GSLYAKVRKK 
             SSSD
  125 -  166   pgipgg pqaipatnsp dhsdhtlsvs sdsghstasa rtdkte
  167 -  244 ERLA PGTRRGLSAQ EKAELDQLLS GFGLEDPGSS LKEMTDARSK YSGTRHVVPA QVHVNG
             DAAL KDRETDILDD EMPH
  245 -  266   hdlhsv dslgtlssse gpqsah
  267 - 1205 LGPF TCHKSSQNSL LSDGFGSNVG EDPQGTLVPD LGLGMDGPYE RERTFGSREP KQPQPL
             LRKP SVSAQMQAYG QSSYSTQTWV RQQQMVVAHQ YSFAPDGEAR LVSRCPADNP GLVQAQ
             PRVP LTPTRGTSSR VAVQRGVGSG PHPPDTQQPS PSKAFKPRFP GDQVVNGAGP ELSTGP
             SPGS PTLDIDQSIE QLNRLILELD PTFEPIPTHM NALGSQANGS VSPDSVGGGL RASSRL
             PDTG EGPSRATGRQ GSSAEQPLGG RLRKLSLGQY DNDAGGQLPF SKCAWGKAGV DYAPNL
             PPFP SPADVKETMT PGYPQDLDII DGRILSSKES MCSTPAFPVS PETPYVKTAL RHPPFS
             PPEP PLSSPASQHK GGREPRSCPE TLTHAVGMSE SPIGPKSTML RADASSTPSF QQAFAS
             SCTI SSNGPGQRRE SSSSAERQWV ESSPKPMVSL LGSGRPTGSP LSAEFSGTRK DSPVLS
             CFPP SELQAPFHSH ELSLAEPPDS LAPPSSQAFL GFGTAPVGSG LPPEEDLGAL LANSHG
             ASPT PSIPLTATGA ADNGFLSHNF LTVAPGHSSH HSPGLQGQGV TLPGQPPLPE KKRASE
             GDRS LGSVSPSSSG FSSPHSGSTI SIPFPNVLPD FSKASEAASP LPDSPGDKLV IVKFVQ
             DTSK FWYKADISRE QAIAMLKDKE PGSFIVRDSH SFRGAYGLAM KVATPPPSVL QLNKKA
             GDLA NELVRHFLIE CTPKGVRLKG CSNEPYFGSL TALVCQHSIT PLALPCKLLI PERDPL
             EEIA ESSPQTAANS AAELLKQGAA CNVWYLNSVE MESLTGHQAI QKALSITLVQ EPPPVS
             TVVH FKVSAQGITL TDNQRKLFFR RHYPVNSVIF CALDPQDRKW IKDGPSSKVF GFVARK
             QGSA TDNVCHLFAE HDPEQPASAI VNFVSKVMIG SPKKV

low complexity regions: DUST
>gi|15987493|gb|AAL11993.1|AF378756_1 tumor endothelial marker 6 [Homo sapiens]
MVKCYHKKYRSATRDVIFRLQFHTGAVQGYGLVFGKEDLDNASKDDRFPDYGKVELVFSA
TPEKIQGSEHLYNDHGVIVDYNTTDPLIRWDSYENLSADGEVLHTQGPVDGSLYAKVRKK
SSSDPGIPGGPQAIPATNSPDHSDHTLSVSSDSGHSTASARTDKTEERLAPGTRRGLSAQ
EKAELDQLLSGFGLEDPGSSLKEMTDARSKYSGTRHVVPAQVHVNGDAALKDRETDILDD
EMPHHDLHSVDSLGTLSSSEGPQSAHLGPFTCHKSSQNSLLSDGFGSNVGEDPQGTLVPD
LGLGMDGPYERERTFGSREPKQPQPLLRKPSVSAQMQAYGQSSYSTQTWVRQQQMVVAHQ
YSFAPDGEARLVSRCPADNPGLVQAQPRVPLTPTRGTSSRVAVQRGVGSGPHPPDTQQPS
PSKAFKPRFPGDQVVNGAGPELSTGPSPGSPTLDIDQSIEQLNRLILELDPTFEPIPTHM
NALGSQANGSVSPDSVGGGLRASSRLPDTGEGPSRATGRQGSSAEQPLGGRLRKLSLGQY
DNDAGGQLPFSKCAWGKAGVDYAPNLPPFPSPADVKETMTPGYPQDLDIIDGRILSSKES
MCSTPAFPVSPETPYVKTALRHPPFSPPEPPLSSPASQHKGGREPRSCPETLTHAVGMSE
SPIGPKSTMLRADASSTPSFQQAFASSCTISSNGPGQRRESSSSAERQWVESSPKPMVSL
LGSGRPTGSPLSAEFSGTRKDSPVLSCFPPSELQAPFHSHELSLAEPPDSLAPPSSQAFL
GFGTAPVGSGLPPEEDLGALLANSHGASPTPSIPLTATGAADNGFLSHNFLTVAPGHSSH
HSPGLQGQGVTLPGQPPLPEKKRASEGDRSLGSVSPSSSGFSSPHSGSTISIPFPNVLPD
FSKASEAASPLPDSPGDKLVIVKFVQDTSKFWYKADISREQAIAMLKDKEPGSFIVRDSH
SFRGAYGLAMKVATPPPSVLQLNKKAGDLANELVRHFLIECTPKGVRLKGCSNEPYFGSL
TALVCQHSITPLALPCKLLIPERDPLEEIAESSPQTAANSAAELLKQGAACNVWYLNSVE
MESLTGHQAIQKALSITLVQEPPPVSTVVHFKVSAQGITLTDNQRKLFFRRHYPVNSVIF
CALDPQDRKWIKDGPSSKVFGFVARKQGSATDNVCHLFAEHDPEQPASAIVNFVSKVMIG
SPKKV

~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~

coiled coil prediction for gi|15987493|gb|AAL11993.1|AF378756_1
sequence: 1200 amino acids, 0 residue(s) in coiled coil state

    .    |     .    |     .    |     .    |     .    |     .   60
MVKCYHKKYR SATRDVIFRL QFHTGAVQGY GLVFGKEDLD NASKDDRFPD YGKVELVFSA
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 M'95 -w border
---------- ---------- ---------- ---------- ---------- ---------- * 21 M'95 -w register
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 M'95 +w polar
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 MTK  -w class
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 28 M'95 -w signif.
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 14 M'95 -w local

    .    |     .    |     .    |     .    |     .    |     .  120
TPEKIQGSEH LYNDHGVIVD YNTTDPLIRW DSYENLSADG EVLHTQGPVD GSLYAKVRKK
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 M'95 -w border
---------- ---------- ---------- ---------- ---------- ---------- * 21 M'95 -w register
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 M'95 +w polar
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 MTK  -w class
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 28 M'95 -w signif.
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 14 M'95 -w local

    .    |     .    |     .    |     .    |     .    |     .  180
SSSDPGIPGG PQAIPATNSP DHSDHTLSVS SDSGHSTASA RTDKTEERLA PGTRRGLSAQ
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 M'95 -w border
---------- ---------- ---------- ---------- ---------- ---------- * 21 M'95 -w register
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 M'95 +w polar
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 MTK  -w class
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 28 M'95 -w signif.
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 14 M'95 -w local

    .    |     .    |     .    |     .    |     .    |     .  240
EKAELDQLLS GFGLEDPGSS LKEMTDARSK YSGTRHVVPA QVHVNGDAAL KDRETDILDD
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 M'95 -w border
---------- ---------- ---------- ---------- ---------- ---------- * 21 M'95 -w register
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 M'95 +w polar
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 MTK  -w class
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 28 M'95 -w signif.
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 14 M'95 -w local

    .    |     .    |     .    |     .    |     .    |     .  300
EMPHHDLHSV DSLGTLSSSE GPQSAHLGPF TCHKSSQNSL LSDGFGSNVG EDPQGTLVPD
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 M'95 -w border
---------- ---------- ---------- ---------- ---------- ---------- * 21 M'95 -w register
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 M'95 +w polar
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 MTK  -w class
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 28 M'95 -w signif.
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 14 M'95 -w local

    .    |     .    |     .    |     .    |     .    |     .  360
LGLGMDGPYE RERTFGSREP KQPQPLLRKP SVSAQMQAYG QSSYSTQTWV RQQQMVVAHQ
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 M'95 -w border
---------- ---------- ---------- ---------- ---------- ---------- * 21 M'95 -w register
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 M'95 +w polar
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 MTK  -w class
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 28 M'95 -w signif.
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 14 M'95 -w local

    .    |     .    |     .    |     .    |     .    |     .  420
YSFAPDGEAR LVSRCPADNP GLVQAQPRVP LTPTRGTSSR VAVQRGVGSG PHPPDTQQPS
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 M'95 -w border
---------- ---------- ---------- ---------- ---------- ---------- * 21 M'95 -w register
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 M'95 +w polar
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 MTK  -w class
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 28 M'95 -w signif.
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 14 M'95 -w local

    .    |     .    |     .    |     .    |     .    |     .  480
PSKAFKPRFP GDQVVNGAGP ELSTGPSPGS PTLDIDQSIE QLNRLILELD PTFEPIPTHM
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 M'95 -w border
---------- ---------- ---------- ---------- ---------- ---------- * 21 M'95 -w register
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 M'95 +w polar
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 MTK  -w class
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 28 M'95 -w signif.
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~222244444 4444444444 ~~~~~~~~~~ * 14 M'95 -w local

    .    |     .    |     .    |     .    |     .    |     .  540
NALGSQANGS VSPDSVGGGL RASSRLPDTG EGPSRATGRQ GSSAEQPLGG RLRKLSLGQY
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 M'95 -w border
---------- ---------- ---------- ---------- ---------- ---------- * 21 M'95 -w register
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 M'95 +w polar
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 MTK  -w class
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 28 M'95 -w signif.
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 14 M'95 -w local

    .    |     .    |     .    |     .    |     .    |     .  600
DNDAGGQLPF SKCAWGKAGV DYAPNLPPFP SPADVKETMT PGYPQDLDII DGRILSSKES
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 M'95 -w border
---------- ---------- ---------- ---------- ---------- ---------- * 21 M'95 -w register
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 M'95 +w polar
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 MTK  -w class
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 28 M'95 -w signif.
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 14 M'95 -w local

    .    |     .    |     .    |     .    |     .    |     .  660
MCSTPAFPVS PETPYVKTAL RHPPFSPPEP PLSSPASQHK GGREPRSCPE TLTHAVGMSE
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 M'95 -w border
---------- ---------- ---------- ---------- ---------- ---------- * 21 M'95 -w register
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 M'95 +w polar
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 MTK  -w class
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 28 M'95 -w signif.
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 14 M'95 -w local

    .    |     .    |     .    |     .    |     .    |     .  720
SPIGPKSTML RADASSTPSF QQAFASSCTI SSNGPGQRRE SSSSAERQWV ESSPKPMVSL
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 M'95 -w border
---------- ---------- ---------- ---------- ---------- ---------- * 21 M'95 -w register
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 M'95 +w polar
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 MTK  -w class
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 28 M'95 -w signif.
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 14 M'95 -w local

    .    |     .    |     .    |     .    |     .    |     .  780
LGSGRPTGSP LSAEFSGTRK DSPVLSCFPP SELQAPFHSH ELSLAEPPDS LAPPSSQAFL
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 M'95 -w border
---------- ---------- ---------- ---------- ---------- ---------- * 21 M'95 -w register
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 M'95 +w polar
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 MTK  -w class
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 28 M'95 -w signif.
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 14 M'95 -w local

    .    |     .    |     .    |     .    |     .    |     .  840
GFGTAPVGSG LPPEEDLGAL LANSHGASPT PSIPLTATGA ADNGFLSHNF LTVAPGHSSH
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 M'95 -w border
---------- ---------- ---------- ---------- ---------- ---------- * 21 M'95 -w register
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 M'95 +w polar
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 MTK  -w class
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 28 M'95 -w signif.
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 14 M'95 -w local

    .    |     .    |     .    |     .    |     .    |     .  900
HSPGLQGQGV TLPGQPPLPE KKRASEGDRS LGSVSPSSSG FSSPHSGSTI SIPFPNVLPD
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 M'95 -w border
---------- ---------- ---------- ---------- ---------- ---------- * 21 M'95 -w register
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 M'95 +w polar
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 MTK  -w class
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 28 M'95 -w signif.
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 14 M'95 -w local

    .    |     .    |     .    |     .    |     .    |     .  960
FSKASEAASP LPDSPGDKLV IVKFVQDTSK FWYKADISRE QAIAMLKDKE PGSFIVRDSH
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 M'95 -w border
---------- ---------- ---------- ---------- ---------- ---------- * 21 M'95 -w register
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 M'95 +w polar
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 MTK  -w class
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 28 M'95 -w signif.
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 14 M'95 -w local

    .    |     .    |     .    |     .    |     .    |     . 1020
SFRGAYGLAM KVATPPPSVL QLNKKAGDLA NELVRHFLIE CTPKGVRLKG CSNEPYFGSL
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 M'95 -w border
---------- ---------- ---------- ---------- ---------- ---------- * 21 M'95 -w register
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 M'95 +w polar
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 MTK  -w class
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 28 M'95 -w signif.
~~~~~~~~~~ ~~~~~~~~~2 3333333333 33332~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 14 M'95 -w local

    .    |     .    |     .    |     .    |     .    |     . 1080
TALVCQHSIT PLALPCKLLI PERDPLEEIA ESSPQTAANS AAELLKQGAA CNVWYLNSVE
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 M'95 -w border
---------- ---------- ---------- ---------- ---------- ---------- * 21 M'95 -w register
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 M'95 +w polar
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 MTK  -w class
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 28 M'95 -w signif.
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 14 M'95 -w local

    .    |     .    |     .    |     .    |     .    |     . 1140
MESLTGHQAI QKALSITLVQ EPPPVSTVVH FKVSAQGITL TDNQRKLFFR RHYPVNSVIF
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 M'95 -w border
---------- ---------- ---------- ---------- ---------- ---------- * 21 M'95 -w register
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 M'95 +w polar
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 MTK  -w class
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 28 M'95 -w signif.
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 14 M'95 -w local

    .    |     .    |     .    |     .    |     .    |     . 1200
CALDPQDRKW IKDGPSSKVF GFVARKQGSA TDNVCHLFAE HDPEQPASAI VNFVSKVMIG
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 M'95 -w border
---------- ---------- ---------- ---------- ---------- ---------- * 21 M'95 -w register
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 M'95 +w polar
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 MTK  -w class
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 28 M'95 -w signif.
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 14 M'95 -w local



~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~

prediction of transmembrane regions with toppred2

     ***********************************
     *TOPPREDM with eukaryotic function*
     ***********************************

tem16.fa.___inter___ is a single sequence
Using hydrophobicity file: /bio_software/2D/toppredm/lib/Engelman-Steitz.scale
Using cyt/ext file: /bio_software/2D/toppredm/lib/Cyt-Ext.prok
Using sequence file: tem16.fa.___inter___

 (1 sequences)
MVKCYHKKYRSATRDVIFRLQFHTGAVQGYGLVFGKEDLDNASKDDRFPD
YGKVELVFSATPEKIQGSEHLYNDHGVIVDYNTTDPLIRWDSYENLSADG
EVLHTQGPVDGSLYAKVRKKSSSDPGIPGGPQAIPATNSPDHSDHTLSVS
SDSGHSTASARTDKTEERLAPGTRRGLSAQEKAELDQLLSGFGLEDPGSS
LKEMTDARSKYSGTRHVVPAQVHVNGDAALKDRETDILDDEMPHHDLHSV
DSLGTLSSSEGPQSAHLGPFTCHKSSQNSLLSDGFGSNVGEDPQGTLVPD
LGLGMDGPYERERTFGSREPKQPQPLLRKPSVSAQMQAYGQSSYSTQTWV
RQQQMVVAHQYSFAPDGEARLVSRCPADNPGLVQAQPRVPLTPTRGTSSR
VAVQRGVGSGPHPPDTQQPSPSKAFKPRFPGDQVVNGAGPELSTGPSPGS
PTLDIDQSIEQLNRLILELDPTFEPIPTHMNALGSQANGSVSPDSVGGGL
RASSRLPDTGEGPSRATGRQGSSAEQPLGGRLRKLSLGQYDNDAGGQLPF
SKCAWGKAGVDYAPNLPPFPSPADVKETMTPGYPQDLDIIDGRILSSKES
MCSTPAFPVSPETPYVKTALRHPPFSPPEPPLSSPASQHKGGREPRSCPE
TLTHAVGMSESPIGPKSTMLRADASSTPSFQQAFASSCTISSNGPGQRRE
SSSSAERQWVESSPKPMVSLLGSGRPTGSPLSAEFSGTRKDSPVLSCFPP
SELQAPFHSHELSLAEPPDSLAPPSSQAFLGFGTAPVGSGLPPEEDLGAL
LANSHGASPTPSIPLTATGAADNGFLSHNFLTVAPGHSSHHSPGLQGQGV
TLPGQPPLPEKKRASEGDRSLGSVSPSSSGFSSPHSGSTISIPFPNVLPD
FSKASEAASPLPDSPGDKLVIVKFVQDTSKFWYKADISREQAIAMLKDKE
PGSFIVRDSHSFRGAYGLAMKVATPPPSVLQLNKKAGDLANELVRHFLIE
CTPKGVRLKGCSNEPYFGSLTALVCQHSITPLALPCKLLIPERDPLEEIA
ESSPQTAANSAAELLKQGAACNVWYLNSVEMESLTGHQAIQKALSITLVQ
EPPPVSTVVHFKVSAQGITLTDNQRKLFFRRHYPVNSVIFCALDPQDRKW
IKDGPSSKVFGFVARKQGSATDNVCHLFAEHDPEQPASAIVNFVSKVMIG
SPKKV


(p)rokaryotic or (e)ukaryotic: e


Charge-pair energy: 0

Length of full window (odd number!): 21

Length of core window (odd number!): 11

Number of residues to add to each end of helix: 1

Critical length: 60

Upper cutoff for candidates: 1

Lower cutoff for candidates: 0.6
Total of 8 structures are to be tested


Candidate membrane-spanning segments:

 Helix Begin   End   Score Certainity
     1   674   694   0.666 Putative
     2   773   793   1.228 Certain
     3   801   821   0.715 Putative
     4   875   895   0.881 Putative
     5  1016  1036   1.083 Certain

----------------------------------------------------------------------
Structure 1

Transmembrane segments included in this structure:
     Segment       1     2     3     4     5
 Loop length   673    78     7    53   120   169
 K+R profile     +        0.00           +      
                       +        4.00           +      
CYT-EXT prof  0.95           -        0.97      
                    0.67           -        1.16      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: -4.00
Tm probability: 0.03
-> Orientation: N-out

Charge-difference over N-terminal Tm (+-15 residues): -1.00
 (NEG-POS)/(NEG+POS): 0.1096
                 NEG: 81.0000
                 POS: 65.0000
-> Orientation: N-out

CYT-EXT difference:   0.09
-> Orientation: N-out

----------------------------------------------------------------------
Structure 2

Transmembrane segments included in this structure:
     Segment       2     3     4     5
 Loop length   772     7    53   120   169
 K+R profile     +        4.00           +      
                    0.00           +      
CYT-EXT prof  0.94           -        1.16      
                       -        0.97      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: 4.00
Tm probability: 0.20
-> Orientation: N-in

Charge-difference over N-terminal Tm (+-15 residues): 0.00
 (NEG-POS)/(NEG+POS): 0.1111
                 NEG: 90.0000
                 POS: 72.0000
-> Orientation: undecided

CYT-EXT difference:   1.13
-> Orientation: N-out

----------------------------------------------------------------------
Structure 3

Transmembrane segments included in this structure:
     Segment       1     2     4     5
 Loop length   673    78    81   120   169
 K+R profile     +           +           +      
                       +           +      
CYT-EXT prof  0.95        0.79        1.16      
                    0.67        0.97      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: 0.00
Tm probability: 0.12
-> Orientation: undecided

Charge-difference over N-terminal Tm (+-15 residues): -1.00
 (NEG-POS)/(NEG+POS): 0.1096
                 NEG: 81.0000
                 POS: 65.0000
-> Orientation: N-out

CYT-EXT difference:   1.26
-> Orientation: N-out

----------------------------------------------------------------------
Structure 4

Transmembrane segments included in this structure:
     Segment       2     4     5
 Loop length   772    81   120   169
 K+R profile     +           +      
                       +           +      
CYT-EXT prof  0.94        0.97      
                    0.79        1.16      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: 0.00
Tm probability: 0.70
-> Orientation: undecided

Charge-difference over N-terminal Tm (+-15 residues): 0.00
 (NEG-POS)/(NEG+POS): 0.1111
                 NEG: 90.0000
                 POS: 72.0000
-> Orientation: undecided

CYT-EXT difference:  -0.05
-> Orientation: N-in

----------------------------------------------------------------------
Structure 5

Transmembrane segments included in this structure:
     Segment       1     2     3     5
 Loop length   673    78     7   194   169
 K+R profile     +        0.00           +      
                       +           +      
CYT-EXT prof  0.95           -        1.16      
                    0.67        1.26      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: 0.00
Tm probability: 0.05
-> Orientation: undecided

Charge-difference over N-terminal Tm (+-15 residues): -1.00
 (NEG-POS)/(NEG+POS): 0.1096
                 NEG: 81.0000
                 POS: 65.0000
-> Orientation: N-out

CYT-EXT difference:   0.19
-> Orientation: N-out

----------------------------------------------------------------------
Structure 6

Transmembrane segments included in this structure:
     Segment       2     3     5
 Loop length   772     7   194   169
 K+R profile     +           +      
                    0.00           +      
CYT-EXT prof  0.94        1.26      
                       -        1.16      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: 0.00
Tm probability: 0.29
-> Orientation: undecided

Charge-difference over N-terminal Tm (+-15 residues): 0.00
 (NEG-POS)/(NEG+POS): 0.1111
                 NEG: 90.0000
                 POS: 72.0000
-> Orientation: undecided

CYT-EXT difference:   1.03
-> Orientation: N-out

----------------------------------------------------------------------
Structure 7

Transmembrane segments included in this structure:
     Segment       1     2     5
 Loop length   673    78   222   169
 K+R profile     +           +      
                       +           +      
CYT-EXT prof  0.95        1.13      
                    0.67        1.16      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: 0.00
Tm probability: 0.16
-> Orientation: undecided

Charge-difference over N-terminal Tm (+-15 residues): -1.00
 (NEG-POS)/(NEG+POS): 0.1096
                 NEG: 81.0000
                 POS: 65.0000
-> Orientation: N-out

CYT-EXT difference:   0.25
-> Orientation: N-out

----------------------------------------------------------------------
Structure 8

Transmembrane segments included in this structure:
     Segment       2     5
 Loop length   772   222   169
 K+R profile     +           +      
                       +      
CYT-EXT prof  0.94        1.16      
                    1.13      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: 0.00
Tm probability: 1.00
-> Orientation: undecided

Charge-difference over N-terminal Tm (+-15 residues): 0.00
 (NEG-POS)/(NEG+POS): 0.1111
                 NEG: 90.0000
                 POS: 72.0000
-> Orientation: undecided

CYT-EXT difference:   0.97
-> Orientation: N-out

----------------------------------------------------------------------

"tem16" 1205 
 674 694 #f 0.665625
 773 793 #t 1.22812
 801 821 #f 0.714583
 875 895 #f 0.88125
 1016 1036 #t 1.08333



     ************************************
     *TOPPREDM with prokaryotic function*
     ************************************

tem16.fa.___inter___ is a single sequence
Using hydrophobicity file: /bio_software/2D/toppredm/lib/Engelman-Steitz.scale
Using cyt/ext file: /bio_software/2D/toppredm/lib/Cyt-Ext.prok
Using sequence file: tem16.fa.___inter___

 (1 sequences)
MVKCYHKKYRSATRDVIFRLQFHTGAVQGYGLVFGKEDLDNASKDDRFPD
YGKVELVFSATPEKIQGSEHLYNDHGVIVDYNTTDPLIRWDSYENLSADG
EVLHTQGPVDGSLYAKVRKKSSSDPGIPGGPQAIPATNSPDHSDHTLSVS
SDSGHSTASARTDKTEERLAPGTRRGLSAQEKAELDQLLSGFGLEDPGSS
LKEMTDARSKYSGTRHVVPAQVHVNGDAALKDRETDILDDEMPHHDLHSV
DSLGTLSSSEGPQSAHLGPFTCHKSSQNSLLSDGFGSNVGEDPQGTLVPD
LGLGMDGPYERERTFGSREPKQPQPLLRKPSVSAQMQAYGQSSYSTQTWV
RQQQMVVAHQYSFAPDGEARLVSRCPADNPGLVQAQPRVPLTPTRGTSSR
VAVQRGVGSGPHPPDTQQPSPSKAFKPRFPGDQVVNGAGPELSTGPSPGS
PTLDIDQSIEQLNRLILELDPTFEPIPTHMNALGSQANGSVSPDSVGGGL
RASSRLPDTGEGPSRATGRQGSSAEQPLGGRLRKLSLGQYDNDAGGQLPF
SKCAWGKAGVDYAPNLPPFPSPADVKETMTPGYPQDLDIIDGRILSSKES
MCSTPAFPVSPETPYVKTALRHPPFSPPEPPLSSPASQHKGGREPRSCPE
TLTHAVGMSESPIGPKSTMLRADASSTPSFQQAFASSCTISSNGPGQRRE
SSSSAERQWVESSPKPMVSLLGSGRPTGSPLSAEFSGTRKDSPVLSCFPP
SELQAPFHSHELSLAEPPDSLAPPSSQAFLGFGTAPVGSGLPPEEDLGAL
LANSHGASPTPSIPLTATGAADNGFLSHNFLTVAPGHSSHHSPGLQGQGV
TLPGQPPLPEKKRASEGDRSLGSVSPSSSGFSSPHSGSTISIPFPNVLPD
FSKASEAASPLPDSPGDKLVIVKFVQDTSKFWYKADISREQAIAMLKDKE
PGSFIVRDSHSFRGAYGLAMKVATPPPSVLQLNKKAGDLANELVRHFLIE
CTPKGVRLKGCSNEPYFGSLTALVCQHSITPLALPCKLLIPERDPLEEIA
ESSPQTAANSAAELLKQGAACNVWYLNSVEMESLTGHQAIQKALSITLVQ
EPPPVSTVVHFKVSAQGITLTDNQRKLFFRRHYPVNSVIFCALDPQDRKW
IKDGPSSKVFGFVARKQGSATDNVCHLFAEHDPEQPASAIVNFVSKVMIG
SPKKV


(p)rokaryotic or (e)ukaryotic: p


Charge-pair energy: 0

Length of full window (odd number!): 21

Length of core window (odd number!): 11

Number of residues to add to each end of helix: 1

Critical length: 60

Upper cutoff for candidates: 1

Lower cutoff for candidates: 0.6
Total of 8 structures are to be tested


Candidate membrane-spanning segments:

 Helix Begin   End   Score Certainity
     1   674   694   0.666 Putative
     2   773   793   1.228 Certain
     3   801   821   0.715 Putative
     4   875   895   0.881 Putative
     5  1016  1036   1.083 Certain

----------------------------------------------------------------------
Structure 1

Transmembrane segments included in this structure:
     Segment       1     2     3     4     5
 Loop length   673    78     7    53   120   169
 K+R profile     +        0.00           +      
                       +        4.00           +      
CYT-EXT prof  0.95           -        0.97      
                    0.67           -        1.16      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: -4.00
Tm probability: 0.03
-> Orientation: N-out

Charge-difference over N-terminal Tm (+-15 residues): -1.00
 (NEG-POS)/(NEG+POS): 0.1096
                 NEG: 81.0000
                 POS: 65.0000
-> Orientation: N-out

CYT-EXT difference:   0.09
-> Orientation: N-out

----------------------------------------------------------------------
Structure 2

Transmembrane segments included in this structure:
     Segment       2     3     4     5
 Loop length   772     7    53   120   169
 K+R profile     +        4.00           +      
                    0.00           +      
CYT-EXT prof  0.94           -        1.16      
                       -        0.97      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: 4.00
Tm probability: 0.20
-> Orientation: N-in

Charge-difference over N-terminal Tm (+-15 residues): 0.00
 (NEG-POS)/(NEG+POS): 0.1111
                 NEG: 90.0000
                 POS: 72.0000
-> Orientation: undecided

CYT-EXT difference:   1.13
-> Orientation: N-out

----------------------------------------------------------------------
Structure 3

Transmembrane segments included in this structure:
     Segment       1     2     4     5
 Loop length   673    78    81   120   169
 K+R profile     +           +           +      
                       +           +      
CYT-EXT prof  0.95        0.79        1.16      
                    0.67        0.97      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: 0.00
Tm probability: 0.12
-> Orientation: undecided

Charge-difference over N-terminal Tm (+-15 residues): -1.00
 (NEG-POS)/(NEG+POS): 0.1096
                 NEG: 81.0000
                 POS: 65.0000
-> Orientation: N-out

CYT-EXT difference:   1.26
-> Orientation: N-out

----------------------------------------------------------------------
Structure 4

Transmembrane segments included in this structure:
     Segment       2     4     5
 Loop length   772    81   120   169
 K+R profile     +           +      
                       +           +      
CYT-EXT prof  0.94        0.97      
                    0.79        1.16      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: 0.00
Tm probability: 0.70
-> Orientation: undecided

Charge-difference over N-terminal Tm (+-15 residues): 0.00
 (NEG-POS)/(NEG+POS): 0.1111
                 NEG: 90.0000
                 POS: 72.0000
-> Orientation: undecided

CYT-EXT difference:  -0.05
-> Orientation: N-in

----------------------------------------------------------------------
Structure 5

Transmembrane segments included in this structure:
     Segment       1     2     3     5
 Loop length   673    78     7   194   169
 K+R profile     +        0.00           +      
                       +           +      
CYT-EXT prof  0.95           -        1.16      
                    0.67        1.26      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: 0.00
Tm probability: 0.05
-> Orientation: undecided

Charge-difference over N-terminal Tm (+-15 residues): -1.00
 (NEG-POS)/(NEG+POS): 0.1096
                 NEG: 81.0000
                 POS: 65.0000
-> Orientation: N-out

CYT-EXT difference:   0.19
-> Orientation: N-out

----------------------------------------------------------------------
Structure 6

Transmembrane segments included in this structure:
     Segment       2     3     5
 Loop length   772     7   194   169
 K+R profile     +           +      
                    0.00           +      
CYT-EXT prof  0.94        1.26      
                       -        1.16      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: 0.00
Tm probability: 0.29
-> Orientation: undecided

Charge-difference over N-terminal Tm (+-15 residues): 0.00
 (NEG-POS)/(NEG+POS): 0.1111
                 NEG: 90.0000
                 POS: 72.0000
-> Orientation: undecided

CYT-EXT difference:   1.03
-> Orientation: N-out

----------------------------------------------------------------------
Structure 7

Transmembrane segments included in this structure:
     Segment       1     2     5
 Loop length   673    78   222   169
 K+R profile     +           +      
                       +           +      
CYT-EXT prof  0.95        1.13      
                    0.67        1.16      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: 0.00
Tm probability: 0.16
-> Orientation: undecided

Charge-difference over N-terminal Tm (+-15 residues): -1.00
 (NEG-POS)/(NEG+POS): 0.1096
                 NEG: 81.0000
                 POS: 65.0000
-> Orientation: N-out

CYT-EXT difference:   0.25
-> Orientation: N-out

----------------------------------------------------------------------
Structure 8

Transmembrane segments included in this structure:
     Segment       2     5
 Loop length   772   222   169
 K+R profile     +           +      
                       +      
CYT-EXT prof  0.94        1.16      
                    1.13      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: 0.00
Tm probability: 1.00
-> Orientation: undecided

Charge-difference over N-terminal Tm (+-15 residues): 0.00
 (NEG-POS)/(NEG+POS): 0.1111
                 NEG: 90.0000
                 POS: 72.0000
-> Orientation: undecided

CYT-EXT difference:   0.97
-> Orientation: N-out

----------------------------------------------------------------------

"tem16" 1205 
 674 694 #f 0.665625
 773 793 #t 1.22812
 801 821 #f 0.714583
 875 895 #f 0.88125
 1016 1036 #t 1.08333



~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~

SAPS.  Version of April 11, 1996.
Date run: Wed Nov 14 12:20:46 2001

File: /people/maria/tem16.fa.___saps___
ID   gi|15987493|gb|AAL11993.1|AF378756_1
DE   tumor endothelial marker 6 [Homo sapiens]

number of residues: 1205;   molecular weight: 128.4 kdal
 
         1  MVKCYHKKYR SATRDVIFRL QFHTGAVQGY GLVFGKEDLD NASKDDRFPD YGKVELVFSA 
        61  TPEKIQGSEH LYNDHGVIVD YNTTDPLIRW DSYENLSADG EVLHTQGPVD GSLYAKVRKK 
       121  SSSDPGIPGG PQAIPATNSP DHSDHTLSVS SDSGHSTASA RTDKTEERLA PGTRRGLSAQ 
       181  EKAELDQLLS GFGLEDPGSS LKEMTDARSK YSGTRHVVPA QVHVNGDAAL KDRETDILDD 
       241  EMPHHDLHSV DSLGTLSSSE GPQSAHLGPF TCHKSSQNSL LSDGFGSNVG EDPQGTLVPD 
       301  LGLGMDGPYE RERTFGSREP KQPQPLLRKP SVSAQMQAYG QSSYSTQTWV RQQQMVVAHQ 
       361  YSFAPDGEAR LVSRCPADNP GLVQAQPRVP LTPTRGTSSR VAVQRGVGSG PHPPDTQQPS 
       421  PSKAFKPRFP GDQVVNGAGP ELSTGPSPGS PTLDIDQSIE QLNRLILELD PTFEPIPTHM 
       481  NALGSQANGS VSPDSVGGGL RASSRLPDTG EGPSRATGRQ GSSAEQPLGG RLRKLSLGQY 
       541  DNDAGGQLPF SKCAWGKAGV DYAPNLPPFP SPADVKETMT PGYPQDLDII DGRILSSKES 
       601  MCSTPAFPVS PETPYVKTAL RHPPFSPPEP PLSSPASQHK GGREPRSCPE TLTHAVGMSE 
       661  SPIGPKSTML RADASSTPSF QQAFASSCTI SSNGPGQRRE SSSSAERQWV ESSPKPMVSL 
       721  LGSGRPTGSP LSAEFSGTRK DSPVLSCFPP SELQAPFHSH ELSLAEPPDS LAPPSSQAFL 
       781  GFGTAPVGSG LPPEEDLGAL LANSHGASPT PSIPLTATGA ADNGFLSHNF LTVAPGHSSH 
       841  HSPGLQGQGV TLPGQPPLPE KKRASEGDRS LGSVSPSSSG FSSPHSGSTI SIPFPNVLPD 
       901  FSKASEAASP LPDSPGDKLV IVKFVQDTSK FWYKADISRE QAIAMLKDKE PGSFIVRDSH 
       961  SFRGAYGLAM KVATPPPSVL QLNKKAGDLA NELVRHFLIE CTPKGVRLKG CSNEPYFGSL 
      1021  TALVCQHSIT PLALPCKLLI PERDPLEEIA ESSPQTAANS AAELLKQGAA CNVWYLNSVE 
      1081  MESLTGHQAI QKALSITLVQ EPPPVSTVVH FKVSAQGITL TDNQRKLFFR RHYPVNSVIF 
      1141  CALDPQDRKW IKDGPSSKVF GFVARKQGSA TDNVCHLFAE HDPEQPASAI VNFVSKVMIG 
      1201  SPKKV

--------------------------------------------------------------------------------
COMPOSITIONAL ANALYSIS (extremes relative to: swp23s)

A  : 90( 7.5%); C  : 15( 1.2%); D  : 68( 5.6%); E  : 60( 5.0%); F  : 43( 3.6%)
G  :104( 8.6%); H  : 37( 3.1%); I  : 34( 2.8%); K  : 56( 4.6%); L  :102( 8.5%)
M  : 16( 1.3%); N  : 30( 2.5%); P+ :123(10.2%); Q  : 59( 4.9%); R  : 55( 4.6%)
S+ :151(12.5%); T  : 62( 5.1%); V  : 71( 5.9%); W  :  7( 0.6%); Y  : 22( 1.8%)

KR      :  111 (  9.2%);   ED      :  128 ( 10.6%);   AGP     :  317 ( 26.3%);
KRED    :  239 ( 19.8%);   KR-ED   :  -17 ( -1.4%);   FIKMNY  :  201 ( 16.7%);
LVIFM   :  266 ( 22.1%);   ST      :  213 ( 17.7%).

--------------------------------------------------------------------------------
CHARGE DISTRIBUTIONAL ANALYSIS
 
         1  00+000++0+ 000+-000+0 0000000000 00000+--0- 000+--+00- 00+0-00000 
        61  00-+0000-0 000-00000- 0000-000+0 -00-0000-0 -00000000- 00000+0+++ 
       121  000-000000 0000000000 -00-000000 0-00000000 +0-+0--+00 000++00000 
       181  -+0-0-0000 0000--0000 0+-00-0+0+ 0000+00000 000000-000 +-+-0-00-- 
       241  -0000-0000 -00000000- 0000000000 000+000000 00-0000000 --0000000- 
       301  00000-000- +-+0000+-0 +000000++0 0000000000 0000000000 +000000000 
       361  00000-0-0+ 000+000-00 0000000+00 0000+0000+ 0000+00000 0000-00000 
       421  00+00+0+00 0-00000000 -000000000 000-0-000- 000+000-0- 000-000000 
       481  0000000000 000-000000 +000+00-00 -000+000+0 0000-00000 +0++000000 
       541  -0-0000000 0+0000+000 -000000000 000-0+-000 00000-0-00 -0+0000+-0 
       601  0000000000 0-0000+000 +0000000-0 000000000+ 00+-0+000- 000000000- 
       661  00000+0000 +0-0000000 0000000000 0000000++- 00000-+000 -000+00000 
       721  0000+00000 000-0000++ -000000000 0-00000000 -0000-00-0 0000000000 
       781  0000000000 000---0000 0000000000 0000000000 0-00000000 0000000000 
       841  0000000000 000000000- +++00-0-+0 0000000000 0000000000 000000000- 
       901  00+00-0000 00-000-+00 00+000-00+ 000+0-00+- 000000+-+- 000000+-00 
       961  00+0000000 +000000000 000++00-00 0-00+0000- 000+00+0+0 000-000000 
      1021  0000000000 000000+000 0-+-00--00 -000000000 00-00+0000 000000000- 
      1081  0-00000000 0+00000000 -000000000 0+00000000 0-00++000+ +000000000 
      1141  000-00-++0 0+-0000+00 0000++0000 0-0000000- 0-0-000000 00000+0000 
      1201  00++0

A. CHARGE CLUSTERS.


Positive charge clusters (cmin =  9/30 or 12/45 or 15/60):  none


Negative charge clusters (cmin = 10/30 or 14/45 or 17/60):  none


Mixed charge clusters (cmin = 15/30 or 20/45 or 25/60):  none


B. HIGH SCORING (UN)CHARGED SEGMENTS.

There are no high scoring positive charge segments.
There are no high scoring negative charge segments.
There are no high scoring mixed charge segments.

________________________________
High scoring uncharged segments:

score=   1.00 frequency=   0.802  ( LAGSVTIPNFQYHMCW )
score=   0.00 frequency=   0.000  ( BZX )
score=  -8.00 frequency=   0.198  ( KEDR )

 Expected score/letter:  -0.785
 M_0.01=  54.79; M_0.05=  45.08

 1) From  797 to  859:  length= 63, score=54.00  * 
     797  LGALLANSHG ASPTPSIPLT ATGAADNGFL SHNFLTVAPG HSSHHSPGLQ 
     847  GQGVTLPGQP PLP
    L:  9(14.3%);  A:  7(11.1%);  G:  9(14.3%);  S:  7(11.1%);
    P:  9(14.3%);


C. CHARGE RUNS AND PATTERNS.

pattern  (+)|  (-)|  (*)|  (0)| (+0)| (-0)| (*0)|(+00)|(-00)|(*00)| (H.)|(H..)|
lmin0     5 |   5 |   7 |  46 |  10 |  10 |  13 |  12 |  13 |  16 |   8 |  10 | 
lmin1     6 |   7 |   9 |  55 |  12 |  12 |  16 |  15 |  15 |  19 |   9 |  12 | 
lmin2     7 |   8 |  10 |  61 |  13 |  14 |  18 |  16 |  17 |  22 |  11 |  14 | 
 (Significance level: 0.010000; Minimal displayed length:  6)

  (0)  62(1,0,0); at  797- 859:   see sequence above
      (3. quartile)               
      L:  9 (14.3%);  A:  7 (11.1%);  G:  9 (14.3%);  S:  7 (11.1%);
      P:  9 (14.3%);  ST: 12 (19.0%);


Run count statistics:

  +  runs >=   3:   2, at  118;  861;
  -  runs >=   3:   2, at  239;  794;
  *  runs >=   5:   0
  0  runs >=  30:   2, at  823;  870;

--------------------------------------------------------------------------------
DISTRIBUTION OF OTHER AMINO ACID TYPES

1. HIGH SCORING SEGMENTS.
There are no high scoring hydrophobic segments.
There are no high scoring transmembrane segments.


2. SPACINGS OF C.


H2N-3-C-267-C-102-C-177-C-48-C-45-C-39-C-58-C-253-C-9-C-13-C-10-C-34-C-69-C-33-C-30-COOH


2*. SPACINGS OF C and H. (additional deluxe function for ALEX)


H2N-3-C-1-H-16-H-46-H-4-H-28-H-37-H-2-H-9-H-60-H-6-H-20-H-H-2-H-17-H-5-C-H-85-H-15-C-36-H-66-H-73-C-48-C-19-H-16-H-8-C-5-H-33-C-58-C-10-H-1-H-44-H-22-H-8-H-2-H-H-43-H-74-H-35-H-4-C-9-C-13-C-1-H-8-C-34-C-15-H-22-H-21-H-8-C-33-C-H-4-H-24-COOH

--------------------------------------------------------------------------------
REPETITIVE STRUCTURES.

A. SEPARATED, TANDEM, AND PERIODIC REPEATS: amino acid alphabet.
Repeat core block length:  5

Aligned matching blocks:


[ 489- 493]   GSVSP
[ 872- 876]   GSVSP


B. SEPARATED AND TANDEM REPEATS: 11-letter reduced alphabet.
   (i= LVIF; += KR; -= ED; s= AG; o= ST; n= NQ; a= YW; p= P; h= H; m= M; c= C)
Repeat core block length:  9

--------------------------------------------------------------------------------

MULTIPLETS.

A. AMINO ACID ALPHABET.

1. Total number of amino acid multiplets:  77  (Expected range:  50--102)

2. Histogram of spacings between consecutive amino acid multiplets:
   (1-5) 28   (6-10) 15   (11-20) 18   (>=21) 17

3. Clusters of amino acid multiplets (cmin = 11/30 or 14/45 or 17/60):  none


B. CHARGE ALPHABET.

1. Total number of charge multiplets:  22  (Expected range:   7-- 36)
   14 +plets (f+: 9.2%), 8 -plets (f-: 10.6%)
   Total number of charge altplets: 27 (Critical number: 40)

2. Histogram of spacings between consecutive charge multiplets:
   (1-5) 2   (6-10) 3   (11-20) 3   (>=21) 15

--------------------------------------------------------------------------------
PERIODICITY ANALYSIS.

A. AMINO ACID ALPHABET (core:  4; !-core: 5)

Location	Period	Element		Copies	Core	Errors
 150- 161	 3	S..       	 4	 4  	 0
 335- 358	 6	Q.....    	 4	 4  	 0
 605- 616	 3	P..       	 4	 4  	 0
 623- 638	 4	P...      	 4	 4  	 0
 701- 704	 1	S         	 4	 4  	 0
 753- 788	 9	L........ 	 4	 4  	 0
 873- 880	 2	S.        	 4	 4  	 0
 873- 892	 5	S....     	 4	 4  	 0


B. CHARGE ALPHABET ({+= KR; -= ED; 0}; core:  5; !-core: 6)
   and HYDROPHOBICITY ALPHABET ({*= KRED; i= LVIF; 0}; core:  6; !-core: 9)

Location	Period	Element		Copies	Core	Errors
 735- 788	 9	i000.000. 	 6	 6  	/0/1/2/0/./1/1/2/./


--------------------------------------------------------------------------------
SPACING ANALYSIS.

Location (Quartile) Spacing     Rank       P-value   Interpretation

 740- 861  (3.)     +( 121)+     1 of 112   0.0005   large maximal spacing



~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~

Start with Pfam (from /data/patterns/pfam)
hmmpfam - search a single seq against HMM database
HMMER 2.1.1 (Dec 1998)
Copyright (C) 1992-1998 Washington University School of Medicine
HMMER is freely distributed under the GNU General Public License (GPL).
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
HMM file:                 /data/patterns/pfam/Pfam
Sequence file:            tem16.fa
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query:  gi|15987493|gb|AAL11993.1|AF378756_1  tumor endothelial marker 6 [Homo sapiens]

Scores for sequence family classification (score includes all domains):
Model        Description                                Score    E-value  N 
--------     -----------                                -----    ------- ---
SH2          SH2 domain                                  40.5      9e-13   1
MR_MLE       Mandelate racemase / muconate lactonizin    -0.6         66   1
Peptidase_M7 Streptomyces extracellular neutral prote   -33.7         51   1
Keratin_B2   Keratin, high sulfur B2 protein            -97.6         75   1
Peptidase_A4 Peptidase A4 family                       -127.7        7.3   1

Parsed for domains:
Model        Domain  seq-f seq-t    hmm-f hmm-t      score  E-value
--------     ------- ----- -----    ----- -----      -----  -------
MR_MLE         1/1     213   237 ..     1    29 [.    -0.6       66
Peptidase_A4   1/1      77   262 ..     1   230 []  -127.7      7.3
Keratin_B2     1/1     644   751 ..     1   177 []   -97.6       75
Peptidase_M7   1/1     667   804 ..     1   133 []   -33.7       51
SH2            1/1     932  1027 ..     1    79 []    40.5    9e-13

Alignments of top-scoring domains:
MR_MLE: domain 1 of 1, from 213 to 237: score -0.6, E = 66
                   *->rvRdgipiAytlasGtdvdelLDDtardi<-*
                        R ++  A+ ++ G   d +L D++ di   
  gi|1598749   213    GTRHVV-PAQVHVNG---DAALKDRETDI    237  

Peptidase_A4: domain 1 of 1, from 77 to 262: score -127.7, E = 7.3
                   *->veysSNWAGAVlIGsGydkVTGtfvVPspsapSgGSSGyGGGyGyWk
                      v ++ N        + y+  ++   V    +p  GS        y k
  gi|1598749    77    VIVDYNTTDPLIRWDSYENLSADGEVLHTQGPVDGSL-------YAK 116  

                   nkrqteeYaASAWVGIDGDTCqtAILQTGvDFYaeggsvsYDAWYEWYPd
                    ++++     S   GI G     AI  T                    Pd
  gi|1598749   117 VRKKS-----SSDPGIPG--GPQAIPAT------NS------------PD 141  

                   yAYDFSGisISeGDsIkvTVtATSktSGsAtlENlTtGqkVThsFSgees
                    ++    +s+S  Ds   T++A +++    t E l  G+   +  S++e 
  gi|1598749   142 HSDH--TLSVSS-DSGHSTASARTDK----TEERLAPGTR--RGLSAQEK 182  

                   geLCetNAEWIVEDFEeCnssGsseLVpFAd....FgPtVtFtdasaTss
                   +eL  +   +  ED       Gss L    d +++++ t     a + + 
  gi|1598749   183 AELDQLLSGFGLED------PGSS-LKEMTDarskYSGTRHVVPAQVHVN 225  

                   GetVsldD..AtivDIeqngevLTecSVSGtTtVTveYV<-*
                   G    l D++  i D e    +L ++   Gt   + e     
  gi|1598749   226 G-DAALKDreTDILDDEMPHHDLHSVDSLGTL-SSSEGP    262  

Keratin_B2: domain 1 of 1, from 644 to 751: score -97.6, E = 75
                   *->qtScCGfptCStlgtrPsCGsscCQPsCCeSCCQpsCcqpSCCqPtc
                       +  C  +     g       s   P         +  +       +
  gi|1598749   644    EPRSCPETLTHAVGM----SESPIGP-------KSTMLRAD----AS 675  

                   sqtscCqPtcfqssCCrPsCcqTSCCqPtCcqsssCqtgCgigGsiGyGQ
                   s+ s+ q  +f ssC+ +s  +      +   sss  +            
  gi|1598749   676 STPSFQQ--AFASSCTISSNGP-----GQRRESSSSAERQ---------- 708  

                   eGsSGAvScrirWCRPdCrvegtClPpCCvvsCtaPTCCqpvsaQasCCR
                                      ve++ +P +      +PT   p sa+ s  R
  gi|1598749   709 ------------------WVESSPKPMVSLLGSGRPT-GSPLSAEFSGTR 739  

                   PsCqPyCgqsCCRPaCccsvtCtrTccePc<-*
                            s   P+ +c++        P+   
  gi|1598749   740 -------KDS---PVLSCFP--------PS    751  

Peptidase_M7: domain 1 of 1, from 667 to 804: score -33.7, E = 51
                   *->avtvtYdasnAPsFrsqiarsaqiWnssvsnvrLqegssgaDFtYrE
                      +     das+ PsF++  a+s  i  ss    + +e+ss a+  + E
  gi|1598749   667    STMLRADASSTPSFQQAFASSCTI--SSNGPGQRRESSSSAERQWVE 711  

                   Gn...........DsRGsYAstdgHGrGYiFLDYrqnqqYdstRvtaHEt
                     +++  +  +++ + Gs  s +  G      D      +      a   
  gi|1598749   712 SSpkpmvsllgsgRPTGSPLSAEFSGTRK---DSPVLSCFPPSELQAPFH 758  

                   GHvLGL...PDhYsGPCseLmsGGGPGPsCtnayPnsaERsRvnqLWAnG
                    H L L ++PD    P s+   G G  P      P+        +L An 
  gi|1598749   759 SHELSLaepPDSLAPPSSQAFLGFGTAPVGSGLPPEED----LGALLANS 804  

                   <-*
                      
  gi|1598749     -     -    

SH2: domain 1 of 1, from 932 to 1027: score 40.5, E = 9e-13
                   *->WyhGkIsRkeAErlLlnpggkdGtFLvReSestpGdYtLSVrddgp.
                      Wy  +IsR++A  +L+   +++G+F+vR+S s  G Y L  +   p+
  gi|1598749   932    WYKADISREQAIAMLKD--KEPGSFIVRDSHSFRGAYGLAMKVATPp 976  

                   ................rVkHyrIqsldng..gyyitggvkFssLpeLveh
                   ++  + +++ ++  ++ V+H+ I +++ g +   ++ +++F+sL  Lv  
  gi|1598749   977 psvlqlnkkagdlaneLVRHFLIECTPKGvrLKGCSNEPYFGSLTALVCQ 1026 

                   Y<-*
                   +   
  gi|1598749  1027 H    1027 

//

Start with PfamFrag (from /data/patterns/pfam)
hmmpfam - search a single seq against HMM database
HMMER 2.1.1 (Dec 1998)
Copyright (C) 1992-1998 Washington University School of Medicine
HMMER is freely distributed under the GNU General Public License (GPL).
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
HMM file:                 /data/patterns/pfam/PfamFrag
Sequence file:            tem16.fa
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query:  gi|15987493|gb|AAL11993.1|AF378756_1  tumor endothelial marker 6 [Homo sapiens]

Scores for sequence family classification (score includes all domains):
Model           Description                             Score    E-value  N 
--------        -----------                             -----    ------- ---
SH2             SH2 domain                               40.5      9e-13   1
FA_desaturase   Fatty acid desaturase                     3.3        4.7   1
Ndr             Ndr family                                2.2        6.9   1
Peptidase_M7    Streptomyces extracellular neutral pr     1.2         30   1
DUF63           Membrane protein of unknown function     -0.2         62   1
3A              3A/RNA2 movement protein family          -0.3         88   1
MR_MLE          Mandelate racemase / muconate lactoni    -0.6         66   1
Thymidylate_kin Thymidylate kinase                       -0.6         61   1
Prismane        Prismane                                 -1.4         33   1

Parsed for domains:
Model           Domain  seq-f seq-t    hmm-f hmm-t      score  E-value
--------        ------- ----- -----    ----- -----      -----  -------
FA_desaturase     1/1      63    88 ..   325   350 .]     3.3      4.7
MR_MLE            1/1     213   237 ..     1    29 [.    -0.6       66
Thymidylate_kin   1/1     452   466 ..   183   197 .]    -0.6       61
3A                1/1     468   498 ..   262   295 .]    -0.3       88
Peptidase_M7      1/1     667   690 ..     1    24 [.     1.2       30
Prismane          1/1     793   802 ..   556   565 .]    -1.4       33
DUF63             1/1     822   836 ..    74    90 ..    -0.2       62
Ndr               1/1     845   886 ..   299   339 .]     2.2      6.9
SH2               1/1     932  1027 ..     1    79 []    40.5    9e-13

Alignments of top-scoring domains:
FA_desaturase: domain 1 of 1, from 63 to 88: score 3.3, E = 4.7
                   *->eraqkrahlvsdtgdivsyifdrkvk<-*
                      e++q+  hl+ d g+iv y ++ ++    
  gi|1598749    63    EKIQGSEHLYNDHGVIVDYNTTDPLI    88   

MR_MLE: domain 1 of 1, from 213 to 237: score -0.6, E = 66
                   *->rvRdgipiAytlasGtdvdelLDDtardi<-*
                        R ++  A+ ++ G   d +L D++ di   
  gi|1598749   213    GTRHVV-PAQVHVNG---DAALKDRETDI    237  

Thymidylate_kin: domain 1 of 1, from 452 to 466: score -0.6, E = 61
                   *->vidAsnsieevheeI<-*
                      ++d ++sie++ + I   
  gi|1598749   452    TLDIDQSIEQLNRLI    466  

3A: domain 1 of 1, from 468 to 498: score -0.3, E = 88
                   *->krepsleessesPiavgestssrsksvtasvvaG<-*
                      + +p+ e+    P++ +  +s++  sv++ +v G   
  gi|1598749   468    ELDPTFEPI---PTHMNALGSQANGSVSPDSVGG    498  

Peptidase_M7: domain 1 of 1, from 667 to 690: score 1.2, E = 30
                   *->avtvtYdasnAPsFrsqiarsaqi<-*
                      +     das+ PsF++  a+s  i   
  gi|1598749   667    STMLRADASSTPSFQQAFASSCTI    690  

Prismane: domain 1 of 1, from 793 to 802: score -1.4, E = 33
                   *->peeDlkalLg<-*
                      peeDl alL+   
  gi|1598749   793    PEEDLGALLA    802  

DUF63: domain 1 of 1, from 822 to 836: score -0.2, E = 62
                   *->DaGvlprPisyltvTPG<-*
                      D+G+l    ++ltv PG   
  gi|1598749   822    DNGFLSH--NFLTVAPG    836  

Ndr: domain 1 of 1, from 845 to 886: score 2.2, E = 6.9
                   *->lqGmGylpsvrmprlpRsrtasless.tssssvdGsrlrsrs<-*
                      lqG+G+  + + p  +++r++  + s +s s ++ ++++++s   
  gi|1598749   845    LQGQGVTLPGQPPLPEKKRASEGDRSlGSVSPSSSGFSSPHS    886  

SH2: domain 1 of 1, from 932 to 1027: score 40.5, E = 9e-13
                   *->WyhGkIsRkeAErlLlnpggkdGtFLvReSestpGdYtLSVrddgp.
                      Wy  +IsR++A  +L+   +++G+F+vR+S s  G Y L  +   p+
  gi|1598749   932    WYKADISREQAIAMLKD--KEPGSFIVRDSHSFRGAYGLAMKVATPp 976  

                   ................rVkHyrIqsldng..gyyitggvkFssLpeLveh
                   ++  + +++ ++  ++ V+H+ I +++ g +   ++ +++F+sL  Lv  
  gi|1598749   977 psvlqlnkkagdlaneLVRHFLIECTPKGvrLKGCSNEPYFGSLTALVCQ 1026 

                   Y<-*
                   +   
  gi|1598749  1027 H    1027 

//

Start with Repeat Library (from /data/patterns/repeats-Miguel-Andrade/hmm)
hmmpfam - search a single seq against HMM database
HMMER 2.1.1 (Dec 1998)
Copyright (C) 1992-1998 Washington University School of Medicine
HMMER is freely distributed under the GNU General Public License (GPL).
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
HMM file:                 /data/patterns/repeats-Miguel-Andrade/hmm/repeats.hmm-lib
Sequence file:            tem16.fa
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query:  gi|15987493|gb|AAL11993.1|AF378756_1  tumor endothelial marker 6 [Homo sapiens]

Scores for sequence family classification (score includes all domains):
Model    Description                                    Score    E-value  N 
-------- -----------                                    -----    ------- ---
	[no hits above thresholds]

Parsed for domains:
Model    Domain  seq-f seq-t    hmm-f hmm-t      score  E-value
-------- ------- ----- -----    ----- -----      -----  -------
	[no hits above thresholds]

Alignments of top-scoring domains:
	[no hits above thresholds]
//

~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~

Start with Prosite
---------------------------------------------------------
|          ppsearch (c) 1994 EMBL Data Library          |
|       based on MacPattern (c) 1990-1994 R. Fuchs      |
---------------------------------------------------------

PROSITE pattern search started: Wed Nov 14 12:23:50 2001

Sequence file: tem16.fa

----------------------------------------
Sequence gi|15987493|gb|AAL11993.1|AF378756_1 (1205 residues):

Matching pattern PS00001 ASN_GLYCOSYLATION:
   41: NASK
   82: NTTD
   95: NLSA
  488: NGSV
Total matches: 4

Matching pattern PS00004 CAMP_PHOSPHO_SITE:
  118: RKKS
  119: KKSS
  328: RKPS
  533: RKLS
  698: RRES
  739: RKDS
  862: KRAS
Total matches: 7

Matching pattern PS00005 PKC_PHOSPHO_SITE:
  159: SAR
  162: TDK
  173: TRR
  200: SLK
  398: SSR
  503: SSR
  517: TGR
  596: SSK
  713: SPK
  723: SGR
  738: TRK
  928: TSK
  961: SFR
 1002: TPK
 1156: SSK
 1201: SPK
Total matches: 16

Matching pattern PS00006 CK2_PHOSPHO_SITE:
   43: SKDD
  121: SSSD
  178: SAQE
  200: SLKE
  257: SSSE
  522: SSAE
  571: SPAD
  596: SSKE
  626: SPPE
  647: SCPE
  703: SSAE
  738: TRKD
  763: SLAE
  865: SEGD
  914: SPGD
 1060: SAAE
 1119: TLTD
 1169: SATD
Total matches: 18

Matching pattern PS00008 MYRISTYL:
   25: GAVQGY
   31: GLVFGK
  111: GSLYAK
  129: GGPQAI
  172: GTRRGL
  254: GTLSSS
  286: GSNVGE
  381: GLVQAQ
  484: GSQANG
  498: GGLRAS
  499: GLRASS
  538: GQYDND
  728: GSPLSA
  798: GALLAN
  819: GAADNG
  847: GQGVTL
  887: GSTISI
 1018: GSLTAL
 1068: GAACNV
 1117: GITLTD
Total matches: 20

Matching pattern PS00343 GRAM_POS_ANCHORING:
  506: LPDTGE
Total matches: 1

Total no of hits in this sequence: 66

========================================

1314 pattern(s) searched in 1 sequence(s), 1205 residues.
Total no of hits in all sequences: 66.
Search time: 00:00 min

~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~

Start with Profile Search
L=0  16.1287  1061 pos.    932 -  1042 PS50001|SH2 Src homology 2 (SH2) domain profile.
#
# P       1 WYHGKISREEAERLLRNKGBEGTFLVRESESSPGDYTLSVRXDG----------------     -49
# S     932 WYKADISREQAIAMLKDK-EPGSFIVRDSHSFRGAYGLAMKVATpppsvlqlnkkagdla    -216
#
# P      44 --QVKHYRIQKTDNG-KYYISER-XRFNSLPELIEHYQQNSXGLCXRLRYPV      -1 
# S     991 neLVRHFLIECTPKGvRLKGCSNePYFGSLTALVCQHSITPLALPCKLLIPE    -164 
#

~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~

Start with motif search against own library
     ***** bioMotif : Version V41a DB, 1999 Nov 11 *****
          SeqTyp=2 : PROTEIN  search; 


>APC D-Box is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>ER-GOLGI-traffic signal is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>INTRA-SIGNAL-M minimal SH3 binding  is the MOTIF name

>gi|15987493|gb|AAL11993.1|AF378756_1 tumor endothelial marker 6 [Homo sapiens] ;LENGTH=1205; DIRECT_SEQUENCE
n 23 solutions 
m %_PXXP 125-128
f
m %_PXXP 128-131
f
m %_PXXP 320-323
f
m %_PXXP 387-390
f
m %_PXXP 390-393
f
m %_PXXP 411-414
f
m %_PXXP 427-430
f
m %_PXXP 448-451
f
m %_PXXP 564-567
f
m %_PXXP 567-570
f
m %_PXXP 581-584
f
m %_PXXP 605-608
f
m %_PXXP 608-611
f
m %_PXXP 611-614
f
m %_PXXP 624-627
f
m %_PXXP 627-630
f
m %_PXXP 628-631
f
m %_PXXP 662-665
f
m %_PXXP 811-814
f
m %_PXXP 853-856
f
m %_PXXP 856-859
f
m %_PXXP 912-915
f
m %_PXXP 1183-1186
f

>STATISTICS Total   : 23 solutions in 1 sequences, 1205 units;  out of 1 sequences, 1205 units

>INTRA-SIGNAL-M deubiquitinating enzyme SH3 domain binding motif (Kato, 2000) is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>INTRA-SIGNAL-M minimal class I consensus-SH3 binding motif  is the MOTIF name

>gi|15987493|gb|AAL11993.1|AF378756_1 tumor endothelial marker 6 [Homo sapiens] ;LENGTH=1205; DIRECT_SEQUENCE
n 2 solutions 
m %_R 318-318 %_XPXXP 319-323
f
m %_R 388-388 %_XPXXP 389-393
f

>STATISTICS Total   : 2 solutions in 1 sequences, 1205 units;  out of 1 sequences, 1205 units

>INTRA-SIGNAL-M minimal class II consensus-SH3 binding motif  is the MOTIF name

>gi|15987493|gb|AAL11993.1|AF378756_1 tumor endothelial marker 6 [Homo sapiens] ;LENGTH=1205; DIRECT_SEQUENCE
n 1 solutions 
m %_XPXXPXR 389-395
f

>STATISTICS Total   : 1 solutions in 1 sequences, 1205 units;  out of 1 sequences, 1205 units

>INTRA-SIGNAL-M exact 14-3-3 binding consensus (Muslin 1996 Cell 84 889) is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>INTRA-SIGNAL-M 14-3-3 binding motif in RAF and others (Muslin 1996 Cell 84 889) is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>INTRA-SIGNAL-M WW domain binding motif in formins (Bedford 1997) is the MOTIF name

>gi|15987493|gb|AAL11993.1|AF378756_1 tumor endothelial marker 6 [Homo sapiens] ;LENGTH=1205; DIRECT_SEQUENCE
n 1 solutions 
m %_PPLP 856-859
f

>STATISTICS Total   : 1 solutions in 1 sequences, 1205 units;  out of 1 sequences, 1205 units

>INTRA-SIGNAL-M PY motif for WW domain is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>TM-CYTOPLASMIC-M di-hydrophobic endocytosis motifs for internalized transmembrane proteins is the MOTIF name

>gi|15987493|gb|AAL11993.1|AF378756_1 tumor endothelial marker 6 [Homo sapiens] ;LENGTH=1205; DIRECT_SEQUENCE
n 3 solutions 
m %_E 184-184 %_XXXL 185-188 %_L 189-189
f
m %_D 378-378 %_XXXL 379-382 %_V 383-383
f
m %_D 796-796 %_XXXL 797-800 %_L 801-801
f

>STATISTICS Total   : 3 solutions in 1 sequences, 1205 units;  out of 1 sequences, 1205 units

>TM-CYTOPLASMIC-M tyrosine-based endocytosis motif for internalized transmembrane proteins is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>TM-EXTRACELL-M Endocytosis signal for internalized transmembrane proteins is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>EXTRACELL-M minimal furin protease cleavage site motif  is the MOTIF name

>gi|15987493|gb|AAL11993.1|AF378756_1 tumor endothelial marker 6 [Homo sapiens] ;LENGTH=1205; DIRECT_SEQUENCE
n 1 solutions 
m %_RXXR 643-646
f

>STATISTICS Total   : 1 solutions in 1 sequences, 1205 units;  out of 1 sequences, 1205 units

>EXTRACELL-M extended furin protease cleavage site motif  is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>EXTRACELL-M  zinc binding motif in MMPs is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>EXTRACELL-M g alpha binding go loco is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS PDX-1 experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS QKI-5 experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS HCDA experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS SV40 LrgT experimentally determined  is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS H2B experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS v-Rel experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS Amida experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS Amida experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS RanBP3 experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS Pho4p experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS DNAhelicaseQ1 experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS LEF-1 experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS TCF-1 experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR p53-NLS1 NLS experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS hum-Ku70 experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS GAL4 experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS act/inh betaA experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS BDV-P experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS BDV-P experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS BDV-P experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS TR2 experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS THOV NP experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS polyomaVP1 experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS HIV-1 Tat experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS HIV-1 Rev experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS Rex experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS SRY experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS SRY experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS SOX9 experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS SOX9 experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS NS5A experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS DNAse EBV experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS DNAse EBV experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS adenovE1a experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS ystDNApolalpha experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS hVDR experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS CPV capsid experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS hGlu.cort.experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS cFOS experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS cJUN experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS hDNApolalpha experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS  hDNAtopoII experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS  hDNAtopoII experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS hBLM experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS hARNT experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS influenzaNP experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS influenzaNP experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS p54 experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS hProTalpha experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS Tst1/Oct6 experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS protHsc9 experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS protHsci experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS protHsc3 experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS Ta alpha experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS Pax-QNR experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS Hunt.Dis.pro experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS MyoD experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS MyoD experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS opaque2 experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS CTP experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS HCV experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS HCV experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS p110RB1 experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS VirD2-Nterm experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS VirD2-Cterm experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS Nucloplasmin experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS Nucleolin experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS ICP-8 experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS Nab2 experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS M9 experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS lscMyc experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS humKprotein experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS FluA experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS Mat-alpha experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS polyoma Lrg-T experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS polyoma Lrg-T experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS SV40 VP1 experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS SV40 VP2 experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS polyoma VP2 experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS c-myb experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS N-myc experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS p53 experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS c-erb-A experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS yeast SKI3 experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS L29 experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS L29 experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS Max experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS L3 experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>NUCLEAR NLS dyskerin experimentally determined NLS is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>PDZ domain binding motif science 278_2075_pawson is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units

>WW domain binding motif science 278_2075_pawson is the MOTIF name

>STATISTICS Total   : 0 solutions in 0 sequences, 0 units;  out of 1 sequences, 1205 units


~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~ ~~~

Start with HMM-search search against own library
hmmpfam - search a single seq against HMM database
HMMER 2.1.1 (Dec 1998)
Copyright (C) 1992-1998 Washington University School of Medicine
HMMER is freely distributed under the GNU General Public License (GPL).
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
HMM file:                 /data/patterns/own/own-hmm.lib
Sequence file:            tem16.fa
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query:  gi|15987493|gb|AAL11993.1|AF378756_1  tumor endothelial marker 6 [Homo sapiens]

Scores for sequence family classification (score includes all domains):
Model    Description                                    Score    E-value  N 
-------- -----------                                    -----    ------- ---
	[no hits above thresholds]

Parsed for domains:
Model    Domain  seq-f seq-t    hmm-f hmm-t      score  E-value
-------- ------- ----- -----    ----- -----      -----  -------
	[no hits above thresholds]

Alignments of top-scoring domains:
	[no hits above thresholds]
//
hmmpfam - search a single seq against HMM database
HMMER 2.1.1 (Dec 1998)
Copyright (C) 1992-1998 Washington University School of Medicine
HMMER is freely distributed under the GNU General Public License (GPL).
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
HMM file:                 /data/patterns/own/own-hmm-f.lib
Sequence file:            tem16.fa
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query:  gi|15987493|gb|AAL11993.1|AF378756_1  tumor endothelial marker 6 [Homo sapiens]

Scores for sequence family classification (score includes all domains):
Model       Description                                 Score    E-value  N 
--------    -----------                                 -----    ------- ---
kreil-Nterm                                               1.0         44   1

Parsed for domains:
Model       Domain  seq-f seq-t    hmm-f hmm-t      score  E-value
--------    ------- ----- -----    ----- -----      -----  -------
kreil-Nterm   1/1    1121  1134 ..   172   191 .]     1.0       44

Alignments of top-scoring domains:
kreil-Nterm: domain 1 of 1, from 1121 to 1134: score 1.0, E = 44
                   *->fanqdtdytPLPIfFRhyri<-*
                      + nq+ +     +fFR+++    
  gi|1598749  1121    TDNQR-K-----LFFRRHYP    1134 

//

~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~

L. Aravind's signalling DB+ PSSM from other authors
IMPALA version 1.1 [20-December-1999]


Reference: Alejandro A. Schaffer, Yuri I. Wolf, Chris P. Ponting, 
Eugene V. Koonin, L. Aravind, Stephen F. Altschul (1999), 
"IMPALA: Matching a Protein Sequence Against a Collection of 
"PSI-BLAST-Constructed Position-Specific Score Matrices",
Bioinformatics 15:1000-1011.

Query= gi|15987493|gb|AAL11993.1|AF378756_1 tumor endothelial marker 6
[Homo sapiens]
         (1205 letters)

Searching..................................done
Results from profile search


                                                                   Score     E
Sequences producing significant alignments:                        (bits)  Value

PTB Phospho-Tyrosine Binding domain                                82  4e-18
SH2 Src Homology domain 2                                          64  9e-13
INSL Insulinase like Metallo protease domain                       25  0.42
CNMP Cyclic nucleotide binding domain                              25  0.62
DHHC Novel zinc finger domain with DHHC signature                  24  1.7
PH  Pleckstrin homology domain (lipid and protein interactio...    23  1.8
DSP Dual specificity protein phosphatase                           23  1.8
ACTIN Actin ATPase/ Cytoskeletal ATPase domain                     23  1.9
UBHYD  Ubiquitin C-terminal hydrolase domain                       23  2.8
AP2  A plant specific DNA binding domain (Apetala 2 like)          22  3.5
SNARE Alpha helical domains which are involved in vesicle fu...    21  6.4
HSF Heat Shock Factor Domain (Winged helix type DNA binding ...    21  6.5
CBS cystathionine beta -synthase domain (A  predicted ligand...    21  6.5
BROMO Bromodomain (A chromatin associated domain)                  21  7.7
FHA Forkhead Associated domain (Phosphoserine/threonine bind...    21  9.8

>PTB Phospho-Tyrosine Binding domain 
          Length = 138

 Score = 82.1 bits (201), Expect = 4e-18
 Identities = 24/137 (17%), Positives = 24/137 (17%), Gaps = 24/137 (17%)

Query: 1072 NVWYLNSVEMESLTGHQAIQKALSITLV---------QEPPPVSTVVHFKVSAQGITLTD 1122
                                                                        
Sbjct: 8    VARFLGCVETPKANGSDVAREAIHAIRFQRDLKRSEQTRETAKLQKVEIRISIDNVIIAD 67

Query: 1123 -NQRKLFFRRHYPVNSVIFCALDPQDRKWIKDGPSSKVFGFVARKQGSATDNVCHLFAEH 1181
                                                                        
Sbjct: 68   IKTKAPMY--TFPLGRISFCADDKDDK---------RMFSFIARAEGASGKPSCYAF--- 113

Query: 1182 DPEQPASAIVNFVSKVM 1198
                             
Sbjct: 114  TSEKLAEDITLTIGEAF 130


>SH2 Src Homology domain 2 
          Length = 119

 Score = 64.2 bits (155), Expect = 9e-13
 Identities = 20/118 (16%), Positives = 20/118 (16%), Gaps = 17/118 (14%)

Query: 932  WYKADISREQAIAMLKDK--EPGSFIVRDSHSFRGAYGLAMKVATPPPSVLQLNKKAGDL 989
                                                                        
Sbjct: 7    WFFKDITRKDAERQLLAPGNSAGAFLIRESETLKGSFSLSVRD-------------YDPM 53

Query: 990  ANELVRHFLIECTPKGVRLKGCSNEPYFGSLTALVCQHSITPLALPCKLLIPERDPLE 1047
                                                                      
Sbjct: 54   HGDVIKHYKIRSLDNGGYY--ISPRITFPCISDMIKHYQKQSDGLCRRLEKACISPKP 109


>INSL Insulinase like Metallo protease domain 
          Length = 433

 Score = 25.3 bits (55), Expect = 0.42
 Identities = 18/168 (10%), Positives = 18/168 (10%), Gaps = 39/168 (23%)

Query: 1   MVKCYHKKYRSATRDVIFRLQFHTGAVQGYGLVFGKEDLDNASK--DDRFPDY---GKVE 55
                                                                       
Sbjct: 189 KLLKFYKSFYQPRNMAVV--------------IVGKVNPKEVEEEVMKTFGKEEGRPVPK 234

Query: 56  LVFSATPEKIQGSEHLYNDHG-----VIVDYNTTDPLIRWDSYENLSADGEVLHT--QGP 108
                                                                       
Sbjct: 235 VQIPTEPEQIGIRFKKLKDPRIEKAYWIIGWRVPA----IGKTDYKGLL--VFSEILCGG 288

Query: 109 VDGSLYAKVRKK-------SSSDPGIPGGPQAIPATNSPDHSDHTLSV 149
                                                           
Sbjct: 289 RISVFYRELREKGLVYSYSCGDMGRPRDNIFIITATFPPENYEKVKKR 336


>CNMP Cyclic nucleotide binding domain 
          Length = 118

 Score = 24.8 bits (54), Expect = 0.62
 Identities = 4/42 (9%), Positives = 4/42 (9%), Gaps = 2/42 (4%)

Query: 1016 YFGSLTALVCQ-HSITPLAL-PCKLLIPERDPLEEIAESSPQ 1055
                                                      
Sbjct: 67   MAGEIGVVFNIPQPFTVRTRRLSQVIRIGHHKFKEMVQSDND 108


>DHHC Novel zinc finger domain with DHHC signature 
          Length = 217

 Score = 23.6 bits (50), Expect = 1.7
 Identities = 16/48 (33%), Positives = 16/48 (33%), Gaps = 3/48 (6%)

Query: 598 KESMCSTPAFPVSPETPYVKTALRHPPFSPPEPPLSSPASQHKGGREP 645
                                                           
Sbjct: 3   KMNICNKPSNKTAPEKS-VWTAPAQP--SGPSPELQGQRSRRNGWSWP 47


>PH  Pleckstrin homology domain (lipid and protein interaction domain) 
          Length = 138

 Score = 23.5 bits (50), Expect = 1.8
 Identities = 19/126 (15%), Positives = 19/126 (15%), Gaps = 17/126 (13%)

Query: 1041 PERDPLEEIAESSPQTAANSAAELLKQGAACNVW-----YLNSVEMESLTGHQAIQKALS 1095
                                                                        
Sbjct: 1    PEDDG-NDLTHTFFNPDREG--WLLKLGGRVKTWKRRWFILTDNCLYYFEYTTDKEPRGI 57

Query: 1096 ITL----VQEPPPVSTVVHFKV---SAQGITLTDNQRKLFFRRHYP-VNSVIFCALDPQD 1147
                                                                        
Sbjct: 58   IPLENLSIREVDDPRKPNCFELYIPNNKGQLIKACKTEADGRVVEGNHMVYRISAPTQEE 117

Query: 1148 RK-WIK 1152
                  
Sbjct: 118  KDEWIK 123


>DSP Dual specificity protein phosphatase 
          Length = 185

 Score = 23.4 bits (50), Expect = 1.8
 Identities = 7/23 (30%), Positives = 7/23 (30%)

Query: 594 ILSSKESMCSTPAFPVSPETPYV 616
                                  
Sbjct: 15  LLSDGSGCYSLPSQPCNEVTPRI 37


>ACTIN Actin ATPase/ Cytoskeletal ATPase domain 
          Length = 376

 Score = 23.3 bits (50), Expect = 1.9
 Identities = 22/94 (23%), Positives = 22/94 (23%), Gaps = 11/94 (11%)

Query: 74  DHGVIVDYNTTDPLIRWDSYENL---SADGEVLHTQGPVDGSLYAKVRKKSSS----DPG 126
                                                                       
Sbjct: 74  EHGIVSNWDDMEKIWHHTFYNELRVVSEEPXVLSXEAP----LNPKVNRESAQIMFETFN 129

Query: 127 IPGGPQAIPATNSPDHSDHTLSVSSDSGHSTASA 160
                                             
Sbjct: 130 VPAMYVAIQAVLSLYASGRTTGIVLDSGDGVSHT 163


>UBHYD  Ubiquitin C-terminal hydrolase domain 
          Length = 884

 Score = 22.7 bits (48), Expect = 2.8
 Identities = 11/47 (23%), Positives = 11/47 (23%), Gaps = 9/47 (19%)

Query: 928 TSKFWYKAD------ISREQAIAMLKDKEPGSFIVRDSHSFRGAYGL 968
                                                          
Sbjct: 255 GDGKWCKFDDDVVSRCTKEEAIEHNYGGHDDDLSVRHCTN---AYML 298


>AP2  A plant specific DNA binding domain (Apetala 2 like) 
          Length = 218

 Score = 22.4 bits (47), Expect = 3.5
 Identities = 17/103 (16%), Positives = 17/103 (16%), Gaps = 14/103 (13%)

Query: 483 LGSQANGSVSPDSVGGGLRASSRLPDTGEGPS--------RATGRQGSSAEQPLGGRLRK 534
                                                                       
Sbjct: 2   IASESTKSWEASAVRQENEEEKKKPVKDSGKHPVYRGVRKRNWGKWVSEIREPRKK--SR 59

Query: 535 LSLGQYDNDAGGQLPFSKCA---WG-KAGVDYAPNLPPFPSPA 573
                                                      
Sbjct: 60  IWLGTFPSPEMAARAHDVAALSIKGASAILNFPDLAGSFPRPS 102


>SNARE Alpha helical domains which are involved in vesicle fusion 
          Length = 254

 Score = 21.4 bits (45), Expect = 6.4
 Identities = 3/25 (12%), Positives = 3/25 (12%)

Query: 256 LSSSEGPQSAHLGPFTCHKSSQNSL 280
                                    
Sbjct: 44  LAQFPDRCTLKTGRHNINFISSLGV 68


>HSF Heat Shock Factor Domain (Winged helix type DNA binding domain) 
          Length = 112

 Score = 21.4 bits (45), Expect = 6.5
 Identities = 5/18 (27%), Positives = 5/18 (27%)

Query: 78 IVDYNTTDPLIRWDSYEN 95
                            
Sbjct: 11 LVDDHRTDDVVSWNEEGT 28


>CBS cystathionine beta -synthase domain (A  predicted ligand binding domain) 
          Length = 214

 Score = 21.5 bits (45), Expect = 6.5
 Identities = 11/39 (28%), Positives = 11/39 (28%)

Query: 920 VIVKFVQDTSKFWYKADISREQAIAMLKDKEPGSFIVRD 958
                                                  
Sbjct: 1   MKVRDLMDKNFAKIYVDETVEDAINLLKKKKRFSAPIVD 39


>BROMO Bromodomain (A chromatin associated domain) 
          Length = 94

 Score = 21.3 bits (45), Expect = 7.7
 Identities = 7/37 (18%), Positives = 7/37 (18%), Gaps = 8/37 (21%)

Query: 568 PFPSPADVKETMTPGY------PQDLDIIDGRILSSK 598
                                                
Sbjct: 22  IFSEPVPLSE--VPDYLDHIKKPMDFFTMKQNLEAYR 56


>FHA Forkhead Associated domain (Phosphoserine/threonine binding domain) 
          Length = 89

 Score = 21.0 bits (44), Expect = 9.8
 Identities = 11/30 (36%), Positives = 11/30 (36%), Gaps = 2/30 (6%)

Query: 937 ISREQAIAMLKDKEPGSFIVRDSHSFRGAY 966
                                         
Sbjct: 46  ISRNHAI--LQATEEDTFLLIDLGSRNGTF 73


Underlying Matrix: BLOSUM62
Number of sequences tested against query: 105
Number of sequences better than 10.0: 15 
Number of calls to ALIGN: 15 
Length of query: 1205 
Total length of test sequences: 20182  
Effective length of test sequences: 16136.0
Effective search space size: 18820320.7
Initial X dropoff for ALIGN: 25.0 bits

Y. Wolf's SCOP PSSM
IMPALA version 1.1 [20-December-1999]


Reference: Alejandro A. Schaffer, Yuri I. Wolf, Chris P. Ponting, 
Eugene V. Koonin, L. Aravind, Stephen F. Altschul (1999), 
"IMPALA: Matching a Protein Sequence Against a Collection of 
"PSI-BLAST-Constructed Position-Specific Score Matrices",
Bioinformatics 15:1000-1011.

Query= gi|15987493|gb|AAL11993.1|AF378756_1 tumor endothelial marker 6
[Homo sapiens]
         (1205 letters)

Searching.................................................done
Results from profile search


                                                                   Score     E
Sequences producing significant alignments:                        (bits)  Value

gi|642068 [217..364] Double-stranded beta-helix                    28  0.84
gi|84873 [335..1022] P-loop containing nucleotide triphospha...    27  1.9
gi|2120618 [16..149] Double-stranded beta-helix                    27  2.1
gi|1172934 [123..294] ATP-grasp                                    27  2.1
gi|65848 [1..155] Cupredoxins                                      27  2.2
gi|951110 [5..279] ConA-like lectins/glucanases                    27  2.3
gi|68234 [34..414] Citrate synthase                                26  3.8
gi|2500693 [47..427] Periplasmic binding protein-like II           26  4.7
gi|120356 [34..141] Double-stranded beta-helix                     26  4.9
gi|1078939 [54..180] C-type lectin-like                            25  7.5
gi|153548 [125..455] beta-Lactamase/D-ala carboxypeptidase         25  8.0
gi|117484 [2..138] Double-stranded beta-helix                      25  8.6
gi|1361422 [77..475] Periplasmic binding protein-like II           25  8.6

>gi|642068 [217..364] Double-stranded beta-helix 
          Length = 148

 Score = 28.0 bits (62), Expect = 0.84
 Identities = 8/54 (14%), Positives = 8/54 (14%), Gaps = 3/54 (5%)

Query: 1016 YFGSLTALVCQH-SITPLAL-PCKLLIPERDPLEEIAESSPQTAANSAAELLKQ 1067
                                                                  
Sbjct: 94   YFGELALLTNKPRAASVYAVEDVDCLLD-VNAFERLLGPCMDIMKRNIEHYEQE 146


>gi|84873 [335..1022] P-loop containing nucleotide triphosphate hydrolases 
          Length = 688

 Score = 26.8 bits (59), Expect = 1.9
 Identities = 6/64 (9%), Positives = 6/64 (9%)

Query: 274 KSSQNSLLSDGFGSNVGEDPQGTLVPDLGLGMDGPYERERTFGSREPKQPQPLLRKPSVS 333
                                                                       
Sbjct: 512 RSSLDESIMLMFTNQLTKAGNLTMPFEAVQHKDESERKSYALNTLSAGCISQVNNLRTLA 571

Query: 334 AQMQ 337
               
Sbjct: 572 ANFR 575


>gi|2120618 [16..149] Double-stranded beta-helix 
          Length = 134

 Score = 26.9 bits (59), Expect = 2.1
 Identities = 8/55 (14%), Positives = 8/55 (14%), Gaps = 3/55 (5%)

Query: 1016 YFGSLTAL--VCQHSITPLAL-PCKLLIPERDPLEEIAESSPQTAANSAAELLKQ 1067
                                                                   
Sbjct: 68   TFAEAMMFMDTPNYVATAQAVVPSQLFRFSNKAYLRQLQDNTPLALALLAKLSTR 122


>gi|1172934 [123..294] ATP-grasp 
          Length = 172

 Score = 26.9 bits (59), Expect = 2.1
 Identities = 13/133 (9%), Positives = 13/133 (9%), Gaps = 18/133 (13%)

Query: 885  HSGSTISIPFPNVLPDFSKASEAASPLPDSPGDKLVIVKFVQDTSKFWYKADISREQAIA 944
                                                                        
Sbjct: 29   NGSQGIGVILAEKPQSAVSIMEAFK----QTNISMLQQDFIEEAGNADIRCFVIGDQVVA 84

Query: 945  MLKDKEPGSFIVRDSHSFRGAYGLAMKVATPPPSVLQLNKKAGDLANELVRHFLIECTPK 1004
                                                                        
Sbjct: 85   TMQRIGQNGEFRANCH---------RGGKTEKI---TLSDEEKQIAIQATKAI--GLDVA 130

Query: 1005 GVRLKGCSNEPYF 1017
                         
Sbjct: 131  GVDLIRSKKGLLV 143


>gi|65848 [1..155] Cupredoxins 
          Length = 155

 Score = 26.7 bits (58), Expect = 2.2
 Identities = 23/87 (26%), Positives = 23/87 (26%), Gaps = 10/87 (11%)

Query: 858 LPEKKRASEGDRSLGSVSPSSSGFSSPHSGSTISIPFPNVLPDFSKAS----EAASPLPD 913
                                                                       
Sbjct: 12  LPQVKAMLEKD------TGKVSGDTVTYSGKTVHVVAAAVLPGFPFPSFEVHDKKNPTLE 65

Query: 914 SPGDKLVIVKFVQDTSKFWYKADISRE 940
                                      
Sbjct: 66  IPAGATVDVTFINTNKGFGHSFDITKK 92


>gi|951110 [5..279] ConA-like lectins/glucanases 
          Length = 275

 Score = 26.8 bits (59), Expect = 2.3
 Identities = 14/85 (16%), Positives = 14/85 (16%), Gaps = 9/85 (10%)

Query: 90  WDSYENLSADGEVLHT-------QGPVDGSLYAKVRKKSSSDPGIPGGPQAI--PATNSP 140
                                                                       
Sbjct: 83  WEKSTNRLTNFQAQFSFVIKSPIDNGADGIAFFIAAPDSEIPKNSAGGTLGLFDPSTAQN 142

Query: 141 DHSDHTLSVSSDSGHSTASARTDKT 165
                                    
Sbjct: 143 PSANQVLAVEFDTFYAQDSNGWDPN 167


 Score = 26.4 bits (58), Expect = 2.8
 Identities = 22/121 (18%), Positives = 22/121 (18%), Gaps = 17/121 (14%)

Query: 881 FSSPHSGSTISIPFPNVLPDFSK---ASEAASPLPDSPGDKLVIVKFVQDTSKFWYKADI 937
                                                                       
Sbjct: 17  LNKAHSLDSLSFSYNNFEQDDERNLILQGDATFSASK---GIQLTKVDDNGTPAKSTVGR 73

Query: 938 SREQAIAMLKDKEPG---------SFIVRDSHSFRGAYGLAMKVATPPPSVLQLNKKAGD 988
                                                                       
Sbjct: 74  VLHSTQVRLWEKSTNRLTNFQAQFSFVIK-SPIDNGADGIAFFIA-APDSEIPKNSAGGT 131

Query: 989 L 989
            
Sbjct: 132 L 132


>gi|68234 [34..414] Citrate synthase 
          Length = 381

 Score = 26.0 bits (57), Expect = 3.8
 Identities = 5/22 (22%), Positives = 5/22 (22%)

Query: 142 HSDHTLSVSSDSGHSTASARTD 163
                                 
Sbjct: 194 HADHEQNASTSTVRLAGSSGAN 215


>gi|2500693 [47..427] Periplasmic binding protein-like II 
          Length = 381

 Score = 25.7 bits (56), Expect = 4.7
 Identities = 10/109 (9%), Positives = 10/109 (9%), Gaps = 23/109 (21%)

Query: 896 NVLPDFSKASEAASPLPDSPGDKLVIVKFVQDTSKFWYKADISREQAIAMLK----DKEP 951
                                                                       
Sbjct: 274 AAYGNPQYFKPVTSMFGNTA-------AVTNDENTGWFKPGGDPEKAKQLFKEAGYAGEK 326

Query: 952 GSFIVRDSHSFR-----------GAYGLAMKVATPPPSVLQLNKKAGDL 989
                                                            
Sbjct: 327 IVILQATDWAEQSNASQVVAAKLREIGVNAELAPSDWGG-LVSRRSKRD 374


>gi|120356 [34..141] Double-stranded beta-helix 
          Length = 108

 Score = 25.6 bits (56), Expect = 4.9
 Identities = 13/97 (13%), Positives = 13/97 (13%), Gaps = 19/97 (19%)

Query: 985  KAGDLANELVRHFLIEC--------TPKG----VRLKGCSNEPYFGSLTAL-VCQHSITP 1031
                                                                        
Sbjct: 6    ISPDTAAHL---LVVAHGKLKTYQLATNGKEQLLRVDGVGD--YEGEAGLLNIANPNVYT 60

Query: 1032 LAL-PCKLLIPERDPLEEIAESSPQTAANSAAELLKQ 1067
                                                 
Sbjct: 61   ETLTAATVCTISATDFQQLLLKQPQISLQLLTENARK 97


>gi|1078939 [54..180] C-type lectin-like 
          Length = 127

 Score = 24.8 bits (53), Expect = 7.5
 Identities = 3/41 (7%), Positives = 3/41 (7%), Gaps = 1/41 (2%)

Query: 1127 LFFRRHYPVNSVIFCALDPQDRKWIKDGPSSKVFGFVARKQ 1167
                                                     
Sbjct: 87   PRWSPGQPDGQGECLAVYTHSGFL-DDKQCETQYNYICKTK 126


>gi|153548 [125..455] beta-Lactamase/D-ala carboxypeptidase 
          Length = 331

 Score = 24.9 bits (54), Expect = 8.0
 Identities = 15/67 (22%), Positives = 15/67 (22%), Gaps = 8/67 (11%)

Query: 679 SFQQAFASSCTISS-----NGPGQRRESSSSAERQWVESSPKPMVSLL---GSGRPTGSP 730
                                                                       
Sbjct: 117 TGVQALEYSSNAYMVQTALGLMGQTYQPNMFVGTSNLESAMEKLRSTFGEYGLGAATGID 176

Query: 731 LSAEFSG 737
                  
Sbjct: 177 LPDESTG 183


>gi|117484 [2..138] Double-stranded beta-helix 
          Length = 137

 Score = 24.9 bits (54), Expect = 8.6
 Identities = 7/62 (11%), Positives = 7/62 (11%), Gaps = 3/62 (4%)

Query: 1016 YFGSLTAL--VCQHSITPLAL-PCKLLIPERDPLEEIAESSPQTAANSAAELLKQGAACN 1072
                                                                        
Sbjct: 69   FIGELGLFEEGQERSAWVRAKTACEVAEISYKKFRQLIQVNPDILMRLSAQMARRLQVTS 128

Query: 1073 VW 1074
              
Sbjct: 129  EK 130


>gi|1361422 [77..475] Periplasmic binding protein-like II 
          Length = 399

 Score = 24.5 bits (53), Expect = 8.6
 Identities = 9/75 (12%), Positives = 9/75 (12%), Gaps = 4/75 (5%)

Query: 915 PGDKLVIVKFVQDTSKFWYKADISREQAIAMLKDKEPGSFIVRDSHSFRGAYGLAMKVAT 974
                                                                       
Sbjct: 142 RGELAVLEPN----PHWDPETDPVRVQRASRIEVHLGKDPHEVDRMLLAGEAHVDLAGFG 197

Query: 975 PPPSVLQLNKKAGDL 989
                          
Sbjct: 198 VQPAAQERILAEPEL 212


Underlying Matrix: BLOSUM62
Number of sequences tested against query: 1187
Number of sequences better than 10.0: 13 
Number of calls to ALIGN: 14 
Length of query: 1205 
Total length of test sequences: 256703  
Effective length of test sequences: 202633.0
Effective search space size: 235012464.0
Initial X dropoff for ALIGN: 25.0 bits

~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~

calculation of internal repeats with prospero
***** PROSPERO v1.3  Wed Nov 14 12:24:46 2001 *****

Copyright 2000, Richard Mott, Wellcome Trust Centre for Human Genetics, University of Oxford
For help see http://www.well.ox.ac.uk/ariadne  For usage use -help
using gap penalty 11+1k
using matrix BLOSUM62
printing all alignments with eval < 0.100000
using sequence1 gi|15987493|gb|AAL11993.1|AF378756_1
using self-comparison


~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~

TIGRFAM
hmmpfam - search a single seq against HMM database
HMMER 2.1.1 (Dec 1998)
Copyright (C) 1992-1998 Washington University School of Medicine
HMMER is freely distributed under the GNU General Public License (GPL).
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
HMM file:                 /data/patterns/tigrfam/tigrfam.hmm
Sequence file:            tem16.fa
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query:  gi|15987493|gb|AAL11993.1|AF378756_1  tumor endothelial marker 6 [Homo sapiens]

Scores for sequence family classification (score includes all domains):
Model    Description                                    Score    E-value  N 
-------- -----------                                    -----    ------- ---
	[no hits above thresholds]

Parsed for domains:
Model    Domain  seq-f seq-t    hmm-f hmm-t      score  E-value
-------- ------- ----- -----    ----- -----      -----  -------
	[no hits above thresholds]

Alignments of top-scoring domains:
	[no hits above thresholds]
//
hmmpfam - search a single seq against HMM database
HMMER 2.1.1 (Dec 1998)
Copyright (C) 1992-1998 Washington University School of Medicine
HMMER is freely distributed under the GNU General Public License (GPL).
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
HMM file:                 /data/patterns/tigrfam/tigrfam.hmm-f
Sequence file:            tem16.fa
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query:  gi|15987493|gb|AAL11993.1|AF378756_1  tumor endothelial marker 6 [Homo sapiens]

Scores for sequence family classification (score includes all domains):
Model     Description                                   Score    E-value  N 
--------  -----------                                   -----    ------- ---
TIGR00041 DTMP_kinase: thymidylate kinase                 2.4         12   1
TIGR01006 8a0302s01: capsular polysaccharide transloc     0.2         35   1
TIGR00748 fabH_rel: condensing enzyme, putative, FabH    -0.4         39   1
TIGR00704 NaPi_cotrn_rel: Na/Pi cotransporter II-rela    -0.5         40   1
TIGR00687 pyridox_kin: pyridoxal kinase                  -0.8         85   1

Parsed for domains:
Model     Domain  seq-f seq-t    hmm-f hmm-t      score  E-value
--------  ------- ----- -----    ----- -----      -----  -------
TIGR00748   1/1     109   123 ..   333   347 .]    -0.4       39
TIGR00041   1/1     452   466 ..   196   211 .]     2.4       12
TIGR01006   1/1     703   723 ..   166   186 ..     0.2       35
TIGR00704   1/1    1016  1043 ..   295   322 .]    -0.5       40
TIGR00687   1/1    1041  1087 ..     1    47 [.    -0.8       85

Alignments of top-scoring domains:
TIGR00748: domain 1 of 1, from 109 to 123: score -0.4, E = 39
                   *->iDYAlYAKfrkKlki<-*
                      +D  lYAK+rkK      
  gi|1598749   109    VDGSLYAKVRKKSSS    123  

TIGR00041: domain 1 of 1, from 452 to 466: score 2.4, E = 12
                   *->viDatnsgveeveqdI<-*
                      ++D+++s +e++ + I   
  gi|1598749   452    TLDIDQS-IEQLNRLI    466  

TIGR01006: domain 1 of 1, from 703 to 723: score 0.2, E = 35
                   *->eeAkpattpsSPkiklnlliG<-*
                      + A  + + sSPk++++ll++   
  gi|1598749   703    SSAERQWVESSPKPMVSLLGS    723  

TIGR00704: domain 1 of 1, from 1016 to 1043: score -0.5, E = 40
                   *->lfnlvncllflPftdplarfvewlirpk<-*
                      +f  +++l+    + pla  +++li+++   
  gi|1598749  1016    YFGSLTALVCQHSITPLALPCKLLIPER    1043 

TIGR00687: domain 1 of 1, from 1041 to 1087: score -0.8, E = 85
                   *->lpknvlsiqshvvyGhaGnkaatfPlqrlGfevwalntvqfsnhtgy
                        ++ l   +  +   a n+aa++  q    +vw ln+v++ + tg+
  gi|1598749  1041    PERDPLEEIAESSPQTAANSAAELLKQGAACNVWYLNSVEMESLTGH 1087 

                   <-*
                      
  gi|1598749     -     -    

//
SMART
hmmpfam - search a single seq against HMM database
HMMER 2.1.1 (Dec 1998)
Copyright (C) 1992-1998 Washington University School of Medicine
HMMER is freely distributed under the GNU General Public License (GPL).
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
HMM file:                 /data/patterns/iprscan/data/smart.HMMs
Sequence file:            tem16.fa
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query:  gi|15987493|gb|AAL11993.1|AF378756_1  tumor endothelial marker 6 [Homo sapiens]

Scores for sequence family classification (score includes all domains):
Model    Description                                    Score    E-value  N 
-------- -----------                                    -----    ------- ---
PTB      Phosphotyrosine-binding domain, phosphotyros    92.1    1.1e-23   1
SH2      Src homology 2 domains                          77.9    2.2e-19   1

Parsed for domains:
Model    Domain  seq-f seq-t    hmm-f hmm-t      score  E-value
-------- ------- ----- -----    ----- -----      -----  -------
SH2        1/1     930  1033 ..     1    87 []    77.9  2.2e-19
PTB        1/1    1066  1203 ..     1   156 []    92.1  1.1e-23

Alignments of top-scoring domains:
SH2: domain 1 of 1, from 930 to 1033: score 77.9, E = 2.2e-19
                   *->epWYHGkisReeAEklLkneggpdGtFLvRdSesspgdyvLSvrvk.
                        WY  +isRe+A ++Lk   +++G+F+vRdS s++g+y L ++v +
  gi|1598749   930    KFWYKADISREQAIAMLKD--KEPGSFIVRDSHSFRGAYGLAMKVAt 974  

                   .................gkvkHyrIrrtddgkfylggt.prrkFpsL.eL
                   ++++  + +++ ++  ++ v+H++I++t+ ++++l g++++++F sL+ L
  gi|1598749   975 pppsvlqlnkkagdlanELVRHFLIECTP-KGVRLKGCsNEPYFGSLtAL 1023 

                   vehYqknslg<-*
                   v  ++ ++l    
  gi|1598749  1024 VCQHSITPLA    1033 

PTB: domain 1 of 1, from 1066 to 1203: score 92.1, E = 1.1e-23
                   *->segvsfeVkYLGsvevpearkSlapgkglqvlqeairklr..sekkk
                       +g++++V+YL sve+++ +       g q++q+a++    ++  + 
  gi|1598749  1066    KQGAACNVWYLNSVEMESLT-------GHQAIQKALSITLvqEPPPV 1105 

                   pqkvilsvSsdgvklidektk.qvlaehplrrIsfcaadpq.r.wnt.Gg
                    + v+++vS +g++l+d+++k ++++++p+++++fca dpq+r+w ++G 
  gi|1598749  1106 STVVHFKVSAQGITLTDNQRKlFFRRHYPVNSVIFCALDPQdRkWIKdGP 1155 

                   dgstdddrvFgyiardpgssgSnrgtsrfaCHVFrcekaaaedialaigq
                         ++vFg++ar++gs      ++ + CH+F   +++   ++++++ 
  gi|1598749  1156 -----SSKVFGFVARKQGS------ATDNVCHLFAEHDPE---QPASAIV 1191 

                   aFqvayelklkarse<-*
                   +F   +++ +  ++    
  gi|1598749  1192 NF---VSKVMIGSPK    1203 

//
COG
hmmpfam - search a single seq against HMM database
HMMER 2.1.1 (Dec 1998)
Copyright (C) 1992-1998 Washington University School of Medicine
HMMER is freely distributed under the GNU General Public License (GPL).
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
HMM file:                 /data/patterns/cogs/cogs.hmm
Sequence file:            tem16.fa
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query:  gi|15987493|gb|AAL11993.1|AF378756_1  tumor endothelial marker 6 [Homo sapiens]

Scores for sequence family classification (score includes all domains):
Model    Description                                    Score    E-value  N 
-------- -----------                                    -----    ------- ---
COG3193                                                 -82.9         82   1

Parsed for domains:
Model    Domain  seq-f seq-t    hmm-f hmm-t      score  E-value
-------- ------- ----- -----    ----- -----      -----  -------
COG3193    1/1     324   472 ..     1   181 []   -82.9       82

Alignments of top-scoring domains:
COG3193: domain 1 of 1, from 324 to 472: score -82.9, E = 82
                   *->vlKffFtPlLvlIaLGLSlaLvmktlAletyyvLslelAlriiaaAi
                            +PlL   +        m++     y + s ++   + ++  
  gi|1598749   324    ------QPLLRKPSV----SAQMQA-----YGQSSYSTQTWVRQQQM 355  

                   aeaQrqgyaP......VtiAVVDagGhlLAliRmDdaaphtiqiafeKAr
                    +a + + aP+++ + V+   +D  G+++A  R    +p+t+   +   r
  gi|1598749   356 VVAHQYSFAPdgearlVSRCPADNPGLVQAQPR----VPLTPTRGTS-SR 400  

                   sAvllrtnitldrssgalaamanpGpslisafvp.lpllplTGlsLfpGG
                    Av  +        sg    + ++  +   af p+ p+++        G 
  gi|1598749   401 VAVQRG------VGSGPHPPDTQQPSPS-KAFKPrFPGDQV-----VNGA 438  

                   VpiivdGaqvVGAVGVS.GataGaqDeqcAikGveavklrld<-*
                   +p ++ G         S+G++  ++D+++    ++ + l++    
  gi|1598749   439 GPELSTG--------PSpGSPTLDIDQSIEQLNRLILELDPT    472  

//
hmmpfam - search a single seq against HMM database
HMMER 2.1.1 (Dec 1998)
Copyright (C) 1992-1998 Washington University School of Medicine
HMMER is freely distributed under the GNU General Public License (GPL).
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
HMM file:                 /data/patterns/cogs/cogs.hmm-f
Sequence file:            tem16.fa
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query:  gi|15987493|gb|AAL11993.1|AF378756_1  tumor endothelial marker 6 [Homo sapiens]

Scores for sequence family classification (score includes all domains):
Model    Description                                    Score    E-value  N 
-------- -----------                                    -----    ------- ---
COG2229                                                   1.0         46   1
COG1776                                                   0.9         41   1
COG0131                                                  -0.2         50   1
COG2175                                                  -0.6         58   1
COG0647                                                  -0.8         74   1
COG1967                                                  -0.9         98   1
COG2189                                                  -0.9         30   1
COG1507                                                  -1.0         93   1
COG0358                                                  -1.4         93   1
COG1222                                                  -2.1         79   1
COG1372                                                  -2.5         96   1
COG2225                                                  -2.7         93   1

Parsed for domains:
Model    Domain  seq-f seq-t    hmm-f hmm-t      score  E-value
-------- ------- ----- -----    ----- -----      -----  -------
COG2189    1/1     219   231 ..   901   913 .]    -0.9       30
COG0647    1/1     238   260 ..   286   308 .]    -0.8       74
COG1507    1/1     352   390 ..     1    43 [.    -1.0       93
COG1776    1/1     468   487 ..   229   248 .]     0.9       41
COG2175    1/1     505   526 ..   314   335 .]    -0.6       58
COG1222    1/1     658   681 ..     1    24 [.    -2.1       79
COG1372    1/1     823   829 ..   554   560 .]    -2.5       96
COG1967    1/1     822   836 ..   108   124 ..    -0.9       98
COG2225    1/1     982   998 ..   159   175 ..    -2.7       93
COG0131    1/1     984   998 ..   128   144 ..    -0.2       50
COG0358    1/1    1113  1141 ..   300   328 ..    -1.4       93
COG2229    1/1    1159  1187 ..   166   195 .]     1.0       46

Alignments of top-scoring domains:
COG2189: domain 1 of 1, from 219 to 231: score -0.9, E = 30
                   *->PDEiYinGDkaVs<-*
                      P  +++nGD+a++   
  gi|1598749   219    PAQVHVNGDAALK    231  

COG0647: domain 1 of 1, from 238 to 260: score -0.8, E = 74
                   *->lDalkPdyvldSlkelvkllkel<-*
                      lD ++P++ l+S+  l +l  ++   
  gi|1598749   238    LDDEMPHHDLHSVDSLGTLSSSE    260  

COG1507: domain 1 of 1, from 352 to 390: score -1.0, E = 93
                   *->mGDPRDEKvVtrQLGRkirnvlkiayrcpyGePvVVkeaPRlk<-*
                              vV  Q   ++ +  + + rcp   P+ V+  PR++   
  gi|1598749   352    ----QQQMVVAHQYSFAPDGEARLVSRCPADNPGLVQAQPRVP    390  

COG1776: domain 1 of 1, from 468 to 487: score 0.9, E = 41
                   *->llipdpesfskllnkLgaql<-*
                      +l+p++e +++++n+Lg q+   
  gi|1598749   468    ELDPTFEPIPTHMNALGSQA    487  

COG2175: domain 1 of 1, from 505 to 526: score -0.6, E = 58
                   *->RvtilgdrPfrnpgllsdqiel<-*
                      R+++ g++P r+ g+++  +e+   
  gi|1598749   505    RLPDTGEGPSRATGRQGSSAEQ    526  

COG1222: domain 1 of 1, from 658 to 681: score -2.1, E = 79
                   *->mseeadlemaelgevdsqtkpske<-*
                      mse++ ++  ++   d++++ps++   
  gi|1598749   658    MSESPIGPKSTMLRADASSTPSFQ    681  

COG1372: domain 1 of 1, from 823 to 829: score -2.5, E = 96
                   *->NgivvhN<-*
                      Ng+++hN   
  gi|1598749   823    NGFLSHN    829  

COG1967: domain 1 of 1, from 822 to 836: score -0.9, E = 98
                   *->DaGvltPprtyltvTPG<-*
                      D+G+l    ++ltv PG   
  gi|1598749   822    DNGFL--SHNFLTVAPG    836  

COG2225: domain 1 of 1, from 982 to 998: score -2.7, E = 93
                   *->ltpeAlefviqLhrrFl<-*
                      l+++A++++ +L+r+Fl   
  gi|1598749   982    LNKKAGDLANELVRHFL    998  

COG0131: domain 1 of 1, from 984 to 998: score -0.2, E = 50
                   *->ekvgvGdlstelieHFf<-*
                      +k  +Gdl  el++HF+   
  gi|1598749   984    KK--AGDLANELVRHFL    998  

COG0358: domain 1 of 1, from 1113 to 1141: score -1.4, E = 93
                   *->VAvlGTalteehlkLlkrfk.ygtkkiilc<-*
                      V ++G +lt+ + kL++r  +y ++++i+c   
  gi|1598749  1113    VSAQGITLTDNQRKLFFR-RhYPVNSVIFC    1141 

COG2229: domain 1 of 1, from 1159 to 1187: score 1.0, E = 46
                   *->vigvsareedsvkdaLdaLLelklksetea<-*
                      v+g+ ar+ +s++d + +L     ++e +a   
  gi|1598749  1159    VFGFVARKQGSATDNVCHLFA-EHDPEQPA    1187 

//