analysis of sequence from tem5
~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~

>gi|2224615|dbj|BAA20795.1| KIAA0337 [Homo sapiens]
MARRLPRTSALKSSSSELLLTGPGAEEDPLPLIVQDQYVQEARQVFEKIQRMGAQQDDGSDAPPGSPDWA
GDVTRGQRSQEELSGPESSLTDEGIGADPEPPVAAFCGLGTTGMWRPLSSSSAQTNHHGPGTEDSLGGWA
LVSPETPPTPGALRRRRKVPPSGSGGSEFSNGEAGEAYRSLSDPIPQRHRAATSEEPTGFSVDSNLLGSL
SPKTGLPATSAMDEGLTSGHSDWSVGSEESKGYQEVIQSIVQGPGTLGRVVDDRIAGKAPKKKSLSDPSR
RGELAGPGFEGPGGEPIREVEPMLPPSSSEPILVEQRAEPEEPGATRSRAQSERALPEALPPPATAHRNF
HLDPKLADILSPRLIRRGSKKRPARSSHQELRRDEGSQDQTGSLSRARPSSRHVRHASVPATFMPIVVPE
PPTSVGPPVAVPEPIGFPTRAHPTLQAPSLEDVTKQYMLNLHSGEVPAPVPVDMPCLPLAAPPSAEAKPP
EAARPADEPTPASKCCSKPQVDMRKHVAMTLLDTEQSYVESLRTLMQGYMQPLKQPENSVLCDPSLVDEI
FDQIPELLEHHEQFLEQVRHCMQTWHAQQKVGALLVQSFSKDVLVNIYSAYIDNFLNAKDAVRVAKEARP
AFLKFLEQSMRENKEKQALSDLMIKPVQRIPRYELLVKDLLKHTPEDHPDHPLLLEAQRNIKQVAERINK
GVRSAEEAERHARVLQEIEAHIEGMEDLQAPLRRFLRQEMVIEVKAIGGKKDRSLFLFTDLIVCTTLKRK
SGSLRRSSMSLYTAASVIDTASKYKMLWKLPLEDADIIKGASQATNRENIQKAISRLDEDLTTLGQMSKL
SESLGFPHQSLDDALRDLSAAMHRDLSEKQALCYALSFPPTKLELCATRPEGTDSYIFEFPHPDARLGFE
QAFDEAKRKLASSKSCLDPEFLKAIPIMKTRSGMQFSCAAPTLNSCPEPSPEVWVCNSDGYVGQVCLLSL
RAEPDVEACIAVCSARILCIGAVPGLQPRCHREPPPSLRSPPETAPEPAGPELDVEAAADEEAATLAEPG
PQPCLHISIAGSGLEMTPGLGEGDPRPELVPFDSDSDDESSPSPSGTLQSQASRSTISSSFGNEETPSSK
EATAETTSSEEEQEPGFLPLSGSFGPGGPCGTSPMDGRALRRSSHGSFTRGSLEDLLSVDPEAYQSSVWL
GTEDGCVHVYQSSDSIRDRRNSMKLQHAASVTCILYLNNQVFVSLANGELVVYQREAGHFWDPQNFKSVT
LGTQGSPITKMVSVGGRLWCGCQNRVLVLSPDTLQLEHMFYVGQDSSRCVACMVDSSLGVWVTLKGSAHV
CLYHPDTFEQLAEVDVTPPVHRMLAGSDAIIRQHKAACLRITALLVCEELLWVGTSAGVVLTMPTSPGTV
SCPRAPLSPTGLGQGHTGHVRFLAAVQLPDGFNLLCPTPPPPPDTGPEKLPSLEHRDSPWHRGPAPARPK
MLVISGGDGYEDFRLSSGGGSSSETVGRDDSTNHLLLWRV
~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~ ~~~~~

sec.str. with predator

> gi|2224615|dbj|BAA20795.1|
              .         .         .         .         .
1    MARRLPRTSALKSSSSELLLTGPGAEEDPLPLIVQDQYVQEARQVFEKIQ   50
     ________________EEEE_____________HHHHHHHHHHHHHHHHH

              .         .         .         .         .
51   RMGAQQDDGSDAPPGSPDWAGDVTRGQRSQEELSGPESSLTDEGIGADPE  100
     HH__________________________HHHHH_________________

              .         .         .         .         .
101  PPVAAFCGLGTTGMWRPLSSSSAQTNHHGPGTEDSLGGWALVSPETPPTP  150
     __EEEEE_______________________________EEEE________

              .         .         .         .         .
151  GALRRRRKVPPSGSGGSEFSNGEAGEAYRSLSDPIPQRHRAATSEEPTGF  200
     _HHHHHH_________EEE_______EEEE______HHHHHHHH______

              .         .         .         .         .
201  SVDSNLLGSLSPKTGLPATSAMDEGLTSGHSDWSVGSEESKGYQEVIQSI  250
     __________________________________________HHHHHHHH

              .         .         .         .         .
251  VQGPGTLGRVVDDRIAGKAPKKKSLSDPSRRGELAGPGFEGPGGEPIREV  300
     H_________HHHHH_______________________________EEE_

              .         .         .         .         .
301  EPMLPPSSSEPILVEQRAEPEEPGATRSRAQSERALPEALPPPATAHRNF  350
     __________HHHHHHH_________HHHHHHHHHH______________

              .         .         .         .         .
351  HLDPKLADILSPRLIRRGSKKRPARSSHQELRRDEGSQDQTGSLSRARPS  400
     ______EEE____EEE__________HHHHHHH_________________

              .         .         .         .         .
401  SRHVRHASVPATFMPIVVPEPPTSVGPPVAVPEPIGFPTRAHPTLQAPSL  450
     _____________EEEEE_________EEEE___________________

              .         .         .         .         .
451  EDVTKQYMLNLHSGEVPAPVPVDMPCLPLAAPPSAEAKPPEAARPADEPT  500
     _HHHHHHHHHH_______________________________________

              .         .         .         .         .
501  PASKCCSKPQVDMRKHVAMTLLDTEQSYVESLRTLMQGYMQPLKQPENSV  550
     _________HHHHHHHHHHHHH_HHHHHHHHHHHHHHHH_________EE

              .         .         .         .         .
551  LCDPSLVDEIFDQIPELLEHHEQFLEQVRHCMQTWHAQQKVGALLVQSFS  600
     E____HHHHHHH_HHHHHHHHHHHHHHHHHHHHHHHHHHHHHHEEE____

              .         .         .         .         .
601  KDVLVNIYSAYIDNFLNAKDAVRVAKEARPAFLKFLEQSMRENKEKQALS  650
     __HHHHHHHHH______HHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHH

              .         .         .         .         .
651  DLMIKPVQRIPRYELLVKDLLKHTPEDHPDHPLLLEAQRNIKQVAERINK  700
     HHH___EEE_HHHHHHHHHHHH__________HHHHHHHHHHHHHHHHH_

              .         .         .         .         .
701  GVRSAEEAERHARVLQEIEAHIEGMEDLQAPLRRFLRQEMVIEVKAIGGK  750
     __HHHHHHHHHHHHHHHHHHHHH___HHHHHHHHHHHHHHHHHHH_____

              .         .         .         .         .
751  KDRSLFLFTDLIVCTTLKRKSGSLRRSSMSLYTAASVIDTASKYKMLWKL  800
     ___EEEEE__EEEE____________________EEEEHHHHHHHHHH__

              .         .         .         .         .
801  PLEDADIIKGASQATNRENIQKAISRLDEDLTTLGQMSKLSESLGFPHQS  850
     _____EEEEHHHHHHHHHHHHHHHHHHH______________________

              .         .         .         .         .
851  LDDALRDLSAAMHRDLSEKQALCYALSFPPTKLELCATRPEGTDSYIFEF  900
     HHHHHHHHHHHHHHH_HHHHHHHHH_______EEEEE________EEEE_

              .         .         .         .         .
901  PHPDARLGFEQAFDEAKRKLASSKSCLDPEFLKAIPIMKTRSGMQFSCAA  950
     _____HHHHHHHHHHHHHHHHHH________EEEEEEEE_______EEE_

              .         .         .         .         .
951  PTLNSCPEPSPEVWVCNSDGYVGQVCLLSLRAEPDVEACIAVCSARILCI 1000
     ____________EEEE____EEEEHHHHHHH_____EEEEHHHHHEEEE_

              .         .         .         .         .
1001 GAVPGLQPRCHREPPPSLRSPPETAPEPAGPELDVEAAADEEAATLAEPG 1050
     ________________________________HHHHHHHHHHHHHH____

              .         .         .         .         .
1051 PQPCLHISIAGSGLEMTPGLGEGDPRPELVPFDSDSDDESSPSPSGTLQS 1100
     ___EEEEEE___________________EEE___________________

              .         .         .         .         .
1101 QASRSTISSSFGNEETPSSKEATAETTSSEEEQEPGFLPLSGSFGPGGPC 1150
     __________________________________________________

              .         .         .         .         .
1151 GTSPMDGRALRRSSHGSFTRGSLEDLLSVDPEAYQSSVWLGTEDGCVHVY 1200
     ________EEEE__________HHHHH_________EEEEE____EEEEE

              .         .         .         .         .
1201 QSSDSIRDRRNSMKLQHAASVTCILYLNNQVFVSLANGELVVYQREAGHF 1250
     E____HHHHHHHHHHHHHH___EEEE____EEEEE____EEEEE______

              .         .         .         .         .
1251 WDPQNFKSVTLGTQGSPITKMVSVGGRLWCGCQNRVLVLSPDTLQLEHMF 1300
     ______EEEE_______EEEEE____EEE____EEEEEE____HHHHHEE

              .         .         .         .         .
1301 YVGQDSSRCVACMVDSSLGVWVTLKGSAHVCLYHPDTFEQLAEVDVTPPV 1350
     EE______EEEEEE______EEEE___EEEEE____HHHHHHHH______

              .         .         .         .         .
1351 HRMLAGSDAIIRQHKAACLRITALLVCEELLWVGTSAGVVLTMPTSPGTV 1400
     EEEE_______HHHHHHHHHHHHHHHHHHHHEEE____EEEEE_______

              .         .         .         .         .
1401 SCPRAPLSPTGLGQGHTGHVRFLAAVQLPDGFNLLCPTPPPPPDTGPEKL 1450
     __________________HHHHHHHH_______EEE______________

              .         .         .         .         .
1451 PSLEHRDSPWHRGPAPARPKMLVISGGDGYEDFRLSSGGGSSSETVGRDD 1500
     ____________________EEEE______EEEE________________

              .
1501 STNHLLLWRV                                         1510
     _HHHHHH___


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


method         :         1
alpha-contents :      37.2 %
beta-contents  :       0.0 %
coil-contents  :      62.8 %
class          :     alpha


method         :         2
alpha-contents :      10.0 %
beta-contents  :       0.0 %
coil-contents  :      90.0 %
class          : irregular


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

GPI: learning from metazoa
-20.28  -2.95  -1.11  -0.32  -4.00   0.00   0.00   0.00  -1.11  -9.32  -2.76 -12.00 -12.00  -4.00 -12.00   0.00  -81.85
 -5.92  -1.83  -1.58  -2.52   0.00   0.00 -16.00   0.00  -0.38  -9.32  -2.76 -12.00 -12.00  -4.00 -12.00   0.00  -80.31
ID: gi|2224615|dbj|BAA20795.1|	AC: xxx Len: 1470 1:I  1457 Sc:  -80.31 Pv: 7.847467e-01 NO_GPI_SITE
GPI: learning from protozoa
-13.38  -2.06  -1.53  -1.49  -4.00   0.00 -16.00   0.00  -0.05  -7.95  -9.85 -12.00 -12.00  -4.00 -12.00   0.00  -96.31
-25.88  -1.89  -0.85  -1.52  -4.00   0.00   0.00   0.00  -0.36  -8.90  -9.85 -12.00 -12.00  -4.00 -12.00   0.00  -93.25
ID: gi|2224615|dbj|BAA20795.1|	AC: xxx Len: 1470 1:I  1452 Sc:  -93.25 Pv: 7.551621e-01 NO_GPI_SITE

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

# SignalP euk predictions
# name       Cmax  pos ?  Ymax  pos ?  Smax  pos ?  Smean ?
gi|2224615|  0.539 488 Y  0.497 1388 Y  0.862 1377 N  0.082 N
# SignalP gram- predictions
# name       Cmax  pos ?  Ymax  pos ?  Smax  pos ?  Smean ?
gi|2224615|  0.676 488 Y  0.395 488 Y  0.979   6 Y  0.128 N
# SignalP gram+ predictions
# name       Cmax  pos ?  Ymax  pos ?  Smax  pos ?  Smean ?
gi|2224615|  0.559 488 Y  0.490 488 Y  0.978   7 Y  0.182 N

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

low complexity regions: SEG 12 2.2 2.5
>gi|2224615|dbj|BAA20795.1| KIAA0337 [Homo sapiens]

                                  1-8    MARRLPRT
                  salksssselll    9-20   
                                 21-56   TGPGAEEDPLPLIVQDQYVQEARQVFEKIQ
                                         RMGAQQ
                  ddgsdappgspd   57-68   
                                 69-161  WAGDVTRGQRSQEELSGPESSLTDEGIGAD
                                         PEPPVAAFCGLGTTGMWRPLSSSSAQTNHH
                                         GPGTEDSLGGWALVSPETPPTPGALRRRRK
                                         VPP
               sgsggsefsngeage  162-176  
                                177-281  AYRSLSDPIPQRHRAATSEEPTGFSVDSNL
                                         LGSLSPKTGLPATSAMDEGLTSGHSDWSVG
                                         SEESKGYQEVIQSIVQGPGTLGRVVDDRIA
                                         GKAPKKKSLSDPSRR
               gelagpgfegpggep  282-296  
                                297-334  IREVEPMLPPSSSEPILVEQRAEPEEPGAT
                                         RSRAQSER
                  alpealpppata  335-346  
                                347-477  HRNFHLDPKLADILSPRLIRRGSKKRPARS
                                         SHQELRRDEGSQDQTGSLSRARPSSRHVRH
                                         ASVPATFMPIVVPEPPTSVGPPVAVPEPIG
                                         FPTRAHPTLQAPSLEDVTKQYMLNLHSGEV
                                         PAPVPVDMPCL
     plaappsaeakppeaarpadeptpa  478-502  
                                503-565  SKCCSKPQVDMRKHVAMTLLDTEQSYVESL
                                         RTLMQGYMQPLKQPENSVLCDPSLVDEIFD
                                         QIP
                  ellehheqfleq  566-577  
                                578-766  VRHCMQTWHAQQKVGALLVQSFSKDVLVNI
                                         YSAYIDNFLNAKDAVRVAKEARPAFLKFLE
                                         QSMRENKEKQALSDLMIKPVQRIPRYELLV
                                         KDLLKHTPEDHPDHPLLLEAQRNIKQVAER
                                         INKGVRSAEEAERHARVLQEIEAHIEGMED
                                         LQAPLRRFLRQEMVIEVKAIGGKKDRSLFL
                                         FTDLIVCTT
               lkrksgslrrssmsl  767-781  
                                782-1011 YTAASVIDTASKYKMLWKLPLEDADIIKGA
                                         SQATNRENIQKAISRLDEDLTTLGQMSKLS
                                         ESLGFPHQSLDDALRDLSAAMHRDLSEKQA
                                         LCYALSFPPTKLELCATRPEGTDSYIFEFP
                                         HPDARLGFEQAFDEAKRKLASSKSCLDPEF
                                         LKAIPIMKTRSGMQFSCAAPTLNSCPEPSP
                                         EVWVCNSDGYVGQVCLLSLRAEPDVEACIA
                                         VCSARILCIGAVPGLQPRCH
repppslrsppetapepagpeldveaaade 1012-1048 
                       eaatlae
                               1049-1082 PGPQPCLHISIAGSGLEMTPGLGEGDPRPE
                                         LVPF
                 dsdsddesspsps 1083-1095 
                               1096-1099 GTLQ
                   sqasrstisss 1100-1110 
                               1111-1113 FGN
eetpsskeataettsseeeqepgflplsgs 1114-1149 
                        fgpggp
                               1150-1436 CGTSPMDGRALRRSSHGSFTRGSLEDLLSV
                                         DPEAYQSSVWLGTEDGCVHVYQSSDSIRDR
                                         RNSMKLQHAASVTCILYLNNQVFVSLANGE
                                         LVVYQREAGHFWDPQNFKSVTLGTQGSPIT
                                         KMVSVGGRLWCGCQNRVLVLSPDTLQLEHM
                                         FYVGQDSSRCVACMVDSSLGVWVTLKGSAH
                                         VCLYHPDTFEQLAEVDVTPPVHRMLAGSDA
                                         IIRQHKAACLRITALLVCEELLWVGTSAGV
                                         VLTMPTSPGTVSCPRAPLSPTGLGQGHTGH
                                         VRFLAAVQLPDGFNLLC
                   ptpppppdtgp 1437-1447 
                               1448-1485 EKLPSLEHRDSPWHRGPAPARPKMLVISGG
                                         DGYEDFRL
                      ssgggsss 1486-1493 
                               1494-1510 ETVGRDDSTNHLLLWRV

low complexity regions: SEG 25 3.0 3.3
>gi|2224615|dbj|BAA20795.1| KIAA0337 [Homo sapiens]

                                  1-1    M
arrlprtsalkssssellltgpgaeedplp    2-32   
                             l
                                 33-135  IVQDQYVQEARQVFEKIQRMGAQQDDGSDA
                                         PPGSPDWAGDVTRGQRSQEELSGPESSLTD
                                         EGIGADPEPPVAAFCGLGTTGMWRPLSSSS
                                         AQTNHHGPGTEDS
lggwalvspetpptpgalrrrrkvppsgsg  136-186  
         gsefsngeageayrslsdpip
                                187-274  QRHRAATSEEPTGFSVDSNLLGSLSPKTGL
                                         PATSAMDEGLTSGHSDWSVGSEESKGYQEV
                                         IQSIVQGPGTLGRVVDDRIAGKAPKKKS
lsdpsrrgelagpgfegpggepirevepml  275-344  
ppsssepilveqraepeepgatrsraqser
                    alpealpppa
                                345-403  TAHRNFHLDPKLADILSPRLIRRGSKKRPA
                                         RSSHQELRRDEGSQDQTGSLSRARPSSRH
vrhasvpatfmpivvpepptsvgppvavpe  404-444  
                   pigfptrahpt
                                445-465  LQAPSLEDVTKQYMLNLHSGE
vpapvpvdmpclplaappsaeakppeaarp  466-509  
                adeptpaskccskp
                                510-663  QVDMRKHVAMTLLDTEQSYVESLRTLMQGY
                                         MQPLKQPENSVLCDPSLVDEIFDQIPELLE
                                         HHEQFLEQVRHCMQTWHAQQKVGALLVQSF
                                         SKDVLVNIYSAYIDNFLNAKDAVRVAKEAR
                                         PAFLKFLEQSMRENKEKQALSDLMIKPVQR
                                         IPRY
       ellvkdllkhtpedhpdhpllle  664-686  
                                687-1001 AQRNIKQVAERINKGVRSAEEAERHARVLQ
                                         EIEAHIEGMEDLQAPLRRFLRQEMVIEVKA
                                         IGGKKDRSLFLFTDLIVCTTLKRKSGSLRR
                                         SSMSLYTAASVIDTASKYKMLWKLPLEDAD
                                         IIKGASQATNRENIQKAISRLDEDLTTLGQ
                                         MSKLSESLGFPHQSLDDALRDLSAAMHRDL
                                         SEKQALCYALSFPPTKLELCATRPEGTDSY
                                         IFEFPHPDARLGFEQAFDEAKRKLASSKSC
                                         LDPEFLKAIPIMKTRSGMQFSCAAPTLNSC
                                         PEPSPEVWVCNSDGYVGQVCLLSLRAEPDV
                                         EACIAVCSARILCIG
avpglqprchrepppslrsppetapepagp 1002-1056 
     eldveaaadeeaatlaepgpqpclh
                               1057-1060 ISIA
gsglemtpglgegdprpelvpfdsdsddes 1061-1173 
spspsgtlqsqasrstisssfgneetpssk
eataettsseeeqepgflplsgsfgpggpc
       gtspmdgralrrsshgsftrgsl
                               1174-1436 EDLLSVDPEAYQSSVWLGTEDGCVHVYQSS
                                         DSIRDRRNSMKLQHAASVTCILYLNNQVFV
                                         SLANGELVVYQREAGHFWDPQNFKSVTLGT
                                         QGSPITKMVSVGGRLWCGCQNRVLVLSPDT
                                         LQLEHMFYVGQDSSRCVACMVDSSLGVWVT
                                         LKGSAHVCLYHPDTFEQLAEVDVTPPVHRM
                                         LAGSDAIIRQHKAACLRITALLVCEELLWV
                                         GTSAGVVLTMPTSPGTVSCPRAPLSPTGLG
                                         QGHTGHVRFLAAVQLPDGFNLLC
ptpppppdtgpeklpslehrdspwhrgpap 1437-1466 
                               1467-1474 ARPKMLVI
   sggdgyedfrlssgggsssetvgrdds 1475-1501 
                               1502-1510 TNHLLLWRV

low complexity regions: SEG 45 3.4 3.75
>gi|2224615|dbj|BAA20795.1| KIAA0337 [Homo sapiens]

                                  1-38   MARRLPRTSALKSSSSELLLTGPGAEEDPL
                                         PLIVQDQY
vqearqvfekiqrmgaqqddgsdappgspd   39-105  
wagdvtrgqrsqeelsgpessltdegigad
                       peppvaa
                                106-227  FCGLGTTGMWRPLSSSSAQTNHHGPGTEDS
                                         LGGWALVSPETPPTPGALRRRRKVPPSGSG
                                         GSEFSNGEAGEAYRSLSDPIPQRHRAATSE
                                         EPTGFSVDSNLLGSLSPKTGLPATSAMDEG
                                         LT
sghsdwsvgseeskgyqeviqsivqgpgtl  228-523  
grvvddriagkapkkkslsdpsrrgelagp
gfegpggepirevepmlppsssepilveqr
aepeepgatrsraqseralpealpppatah
rnfhldpkladilsprlirrgskkrparss
hqelrrdegsqdqtgslsrarpssrhvrha
svpatfmpivvpepptsvgppvavpepigf
ptrahptlqapsledvtkqymlnlhsgevp
apvpvdmpclplaappsaeakppeaarpad
    eptpaskccskpqvdmrkhvamtlld
                                524-946  TEQSYVESLRTLMQGYMQPLKQPENSVLCD
                                         PSLVDEIFDQIPELLEHHEQFLEQVRHCMQ
                                         TWHAQQKVGALLVQSFSKDVLVNIYSAYID
                                         NFLNAKDAVRVAKEARPAFLKFLEQSMREN
                                         KEKQALSDLMIKPVQRIPRYELLVKDLLKH
                                         TPEDHPDHPLLLEAQRNIKQVAERINKGVR
                                         SAEEAERHARVLQEIEAHIEGMEDLQAPLR
                                         RFLRQEMVIEVKAIGGKKDRSLFLFTDLIV
                                         CTTLKRKSGSLRRSSMSLYTAASVIDTASK
                                         YKMLWKLPLEDADIIKGASQATNRENIQKA
                                         ISRLDEDLTTLGQMSKLSESLGFPHQSLDD
                                         ALRDLSAAMHRDLSEKQALCYALSFPPTKL
                                         ELCATRPEGTDSYIFEFPHPDARLGFEQAF
                                         DEAKRKLASSKSCLDPEFLKAIPIMKTRSG
                                         MQF
scaaptlnscpepspevwvcnsdgyvgqvc  947-1182 
llslraepdveaciavcsarilcigavpgl
qprchrepppslrsppetapepagpeldve
aaadeeaatlaepgpqpclhisiagsglem
tpglgegdprpelvpfdsdsddesspspsg
tlqsqasrstisssfgneetpsskeataet
tsseeeqepgflplsgsfgpggpcgtspmd
    gralrrsshgsftrgsledllsvdpe
                               1183-1510 AYQSSVWLGTEDGCVHVYQSSDSIRDRRNS
                                         MKLQHAASVTCILYLNNQVFVSLANGELVV
                                         YQREAGHFWDPQNFKSVTLGTQGSPITKMV
                                         SVGGRLWCGCQNRVLVLSPDTLQLEHMFYV
                                         GQDSSRCVACMVDSSLGVWVTLKGSAHVCL
                                         YHPDTFEQLAEVDVTPPVHRMLAGSDAIIR
                                         QHKAACLRITALLVCEELLWVGTSAGVVLT
                                         MPTSPGTVSCPRAPLSPTGLGQGHTGHVRF
                                         LAAVQLPDGFNLLCPTPPPPPDTGPEKLPS
                                         LEHRDSPWHRGPAPARPKMLVISGGDGYED
                                         FRLSSGGGSSSETVGRDDSTNHLLLWRV


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|2224615|dbj|BAA20795.1| KIAA0337 [Homo sapiens]
MARRLPRTSALKSSSSELLLTGPGAEEDPLPLIVQDQYVQEARQVFEKIQRMGAQQDDGS
DAPPGSPDWAGDVTRGQRSQEELSGPESSLTDEGIGADPEPPVAAFCGLGTTGMWRPLSS
SSAQTNHHGPGTEDSLGGWALVSPETPPTPGALrrrrkvppsgsggsEFSNGEAGEAYRS
LSDPIPQRHRAATSEEPTGFSVDSNLLGSLSPKTGLPATSAMDEGLTSGHSDWSVGSEES
KGYQEVIQSIVQGPGTLGRVVDDRIAGKAPKKKSLSDPSRRGELAGPGFEGPGGEPIREV
EPMLPPSSSEPILVEQRAEPEEPGATRSRAQSERALPEALPPPATAHRNFHLDPKLADIL
SPRLIRRGSKKRPARSSHQELRRDEGSQDQTGSLSRARPSSRHVRHASVPATFMPIVVPE
PPTSVGPPVAVPEPIGFPTRAHPTLQAPSLEDVTKQYMLNLHSGEVPAPVPVDMPCLPLA
APPSAEAKPPEAARPADEPTPASKCCSKPQVDMRKHVAMTLLDTEQSYVESlrtlmqgym
qplkqpENSVLCDPSLVDEIFDQIPELLEHHEQFLEQVRHCMQTWHAQQKVGALLVQSFS
KDVLVNIYSAYIDNFLNAKDAVRVAKEARPAFLKFLEQSMRENKEKQALSDLMIKPVQRI
PRYELLVKDLLKHTPEDHPDHPLLLEAQRNIKQVAERINKGVRSAEEAERHARVLQEIEA
HIEGMEDLQAPLRRFLRQEMVIEVKAIGGKKDRSLFLFTDLIVCTTLKRKSGSLRRSSMS
LYTAASVIDTASKYKMLWKLPLEDADIIKGASQATNRENIQKAISRLDEDLTTLGQMSKL
SESLGFPHQSLDDALRDLSAAMHRDLSEKQALCYALSFPPTKLELCATRPEGTDSYIFEF
PHPDARLGFEQAFDEAKRKLASSKSCLDPEFLKAIPIMKTRSGMQFSCAAPTLNSCPEPS
PEVWVCNSDGYVGQVCLLSLRAEPDVEACIAVCSARILCIGAVPGLQPRCHREPPPSLRS
PPETAPEPAGPELDVEAAADEEAATLAEPGPQPCLHISIAGSGLEMTPGLGEGDPRPELV
PFdsdsddesspspsgtlqsqasrstisssfgneetpsskeataettsseeeqePGFLPL
SGSFGPGGPCGTSPMDGRALRRSSHGSFTRGSLEDLLSVDPEAYQSSVWLGTEDGCVHVY
QSSDSIRDRRNSMKLQHAASVTCILYLNNQVFVSLANGELVVYQREAGHFWDPQNFKSVT
LGTQGSPITKMVSVGGRLWCGCQNRVLVLSPDTLQLEHMFYVGQDSSRCVACMVDSSLGV
WVTLKGSAHVCLYHPDTFEQLAEVDVTPPVHRMLAGSDAIIRQHKAACLRITALLVCEEL
LWVGTSAGVVLTMPTSPGTVSCPRAPLSPTGLGQGHTGHVRFLAAVQLPDGFNLLCptpp
pppDTGPEKLPSLEHRDSPWHRGPAPARPKMLVISGGDGYEDFRLssgggsssetvgrdd
stnhlllWRV
    1 -  153 MARRLPRTSA LKSSSSELLL TGPGAEEDPL PLIVQDQYVQ EARQVFEKIQ RMGAQQDDGS 
             DAPPGSPDWA GDVTRGQRSQ EELSGPESSL TDEGIGADPE PPVAAFCGLG TTGMWRPLSS 
             SSAQTNHHGP GTEDSLGGWA LVSPETPPTP GAL
  154 -  167   rrrrkvp psgsggs
  168 -  531 EFS NGEAGEAYRS LSDPIPQRHR AATSEEPTGF SVDSNLLGSL SPKTGLPATS AMDEGLT
             SGH SDWSVGSEES KGYQEVIQSI VQGPGTLGRV VDDRIAGKAP KKKSLSDPSR RGELAGP
             GFE GPGGEPIREV EPMLPPSSSE PILVEQRAEP EEPGATRSRA QSERALPEAL PPPATAH
             RNF HLDPKLADIL SPRLIRRGSK KRPARSSHQE LRRDEGSQDQ TGSLSRARPS SRHVRHA
             SVP ATFMPIVVPE PPTSVGPPVA VPEPIGFPTR AHPTLQAPSL EDVTKQYMLN LHSGEVP
             APV PVDMPCLPLA APPSAEAKPP EAARPADEPT PASKCCSKPQ VDMRKHVAMT LLDTEQS
             YVE S
  532 -  546   lrtlmqgym qplkqp
  547 - 1082 ENSV LCDPSLVDEI FDQIPELLEH HEQFLEQVRH CMQTWHAQQK VGALLVQSFS KDVLVN
             IYSA YIDNFLNAKD AVRVAKEARP AFLKFLEQSM RENKEKQALS DLMIKPVQRI PRYELL
             VKDL LKHTPEDHPD HPLLLEAQRN IKQVAERINK GVRSAEEAER HARVLQEIEA HIEGME
             DLQA PLRRFLRQEM VIEVKAIGGK KDRSLFLFTD LIVCTTLKRK SGSLRRSSMS LYTAAS
             VIDT ASKYKMLWKL PLEDADIIKG ASQATNRENI QKAISRLDED LTTLGQMSKL SESLGF
             PHQS LDDALRDLSA AMHRDLSEKQ ALCYALSFPP TKLELCATRP EGTDSYIFEF PHPDAR
             LGFE QAFDEAKRKL ASSKSCLDPE FLKAIPIMKT RSGMQFSCAA PTLNSCPEPS PEVWVC
             NSDG YVGQVCLLSL RAEPDVEACI AVCSARILCI GAVPGLQPRC HREPPPSLRS PPETAP
             EPAG PELDVEAAAD EEAATLAEPG PQPCLHISIA GSGLEMTPGL GEGDPRPELV PF
 1083 - 1134   dsdsddes spspsgtlqs qasrstisss fgneetpssk eataettsse eeqe
 1135 - 1436 PGFLPL SGSFGPGGPC GTSPMDGRAL RRSSHGSFTR GSLEDLLSVD PEAYQSSVWL GTED
             GCVHVY QSSDSIRDRR NSMKLQHAAS VTCILYLNNQ VFVSLANGEL VVYQREAGHF WDPQ
             NFKSVT LGTQGSPITK MVSVGGRLWC GCQNRVLVLS PDTLQLEHMF YVGQDSSRCV ACMV
             DSSLGV WVTLKGSAHV CLYHPDTFEQ LAEVDVTPPV HRMLAGSDAI IRQHKAACLR ITAL
             LVCEEL LWVGTSAGVV LTMPTSPGTV SCPRAPLSPT GLGQGHTGHV RFLAAVQLPD GFNL
             LC
 1437 - 1443   ptpp ppp
 1444 - 1485 DTGPEKL PSLEHRDSPW HRGPAPARPK MLVISGGDGY EDFRL
 1486 - 1507   ssggg sssetvgrdd stnhlll
 1508 - 1510 WRV 

low complexity regions: DUST
>gi|2224615|dbj|BAA20795.1| KIAA0337 [Homo sapiens]
MARRLPRTSALKSSSSELLLTGPGAEEDPLPLIVQDQYVQEARQVFEKIQRMGAQQDDGS
DAPPGSPDWAGDVTRGQRSQEELSGPESSLTDEGIGADPEPPVAAFCGLGTTGMWRPLSS
SSAQTNHHGPGTEDSLGGWALVSPETPPTPGALRRRRKVPPSGSGGSEFSNGEAGEAYRS
LSDPIPQRHRAATSEEPTGFSVDSNLLGSLSPKTGLPATSAMDEGLTSGHSDWSVGSEES
KGYQEVIQSIVQGPGTLGRVVDDRIAGKAPKKKSLSDPSRRGELAGPGFEGPGGEPIREV
EPMLPPSSSEPILVEQRAEPEEPGATRSRAQSERALPEALPPPATAHRNFHLDPKLADIL
SPRLIRRGSKKRPARSSHQELRRDEGSQDQTGSLSRARPSSRHVRHASVPATFMPIVVPE
PPTSVGPPVAVPEPIGFPTRAHPTLQAPSLEDVTKQYMLNLHSGEVPAPVPVDMPCLPLA
APPSAEAKPPEAARPADEPTPASKCCSKPQVDMRKHVAMTLLDTEQSYVESLRTLMQGYM
QPLKQPENSVLCDPSLVDEIFDQIPELLEHHEQFLEQVRHCMQTWHAQQKVGALLVQSFS
KDVLVNIYSAYIDNFLNAKDAVRVAKEARPAFLKFLEQSMRENKEKQALSDLMIKPVQRI
PRYELLVKDLLKHTPEDHPDHPLLLEAQRNIKQVAERINKGVRSAEEAERHARVLQEIEA
HIEGMEDLQAPLRRFLRQEMVIEVKAIGGKKDRSLFLFTDLIVCTTLKRKSGSLRRSSMS
LYTAASVIDTASKYKMLWKLPLEDADIIKGASQATNRENIQKAISRLDEDLTTLGQMSKL
SESLGFPHQSLDDALRDLSAAMHRDLSEKQALCYALSFPPTKLELCATRPEGTDSYIFEF
PHPDARLGFEQAFDEAKRKLASSKSCLDPEFLKAIPIMKTRSGMQFSCAAPTLNSCPEPS
PEVWVCNSDGYVGQVCLLSLRAEPDVEACIAVCSARILCIGAVPGLQPRCHREPPPSLRS
PPETAPEPAGPELDVEAAADEEAATLAEPGPQPCLHISIAGSGLEMTPGLGEGDPRPELV
PFDSDSDDESSPSPSGTLQSQASRSTISSSFGNEETPSSKEATAETTSSEEEQEPGFLPL
SGSFGPGGPCGTSPMDGRALRRSSHGSFTRGSLEDLLSVDPEAYQSSVWLGTEDGCVHVY
QSSDSIRDRRNSMKLQHAASVTCILYLNNQVFVSLANGELVVYQREAGHFWDPQNFKSVT
LGTQGSPITKMVSVGGRLWCGCQNRVLVLSPDTLQLEHMFYVGQDSSRCVACMVDSSLGV
WVTLKGSAHVCLYHPDTFEQLAEVDVTPPVHRMLAGSDAIIRQHKAACLRITALLVCEEL
LWVGTSAGVVLTMPTSPGTVSCPRAPLSPTGLGQGHTGHVRFLAAVQLPDGFNLLCPTPP
PPPDTGPEKLPSLEHRDSPWHRGPAPARPKMLVISGGDGYEDFRLSSGGGSSSETVGRDD
STNHLLLWRV

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

coiled coil prediction for gi|2224615|dbj|BAA20795.1|
sequence: 1470 amino acids, 41 residue(s) in coiled coil state
  Coil   1 *  688*   729* RNIKQVAERINKGVRSAEEAERHARVLQEIEAHIEGMEDLQA

    .    |     .    |     .    |     .    |     .    |     .   60
MARRLPRTSA LKSSSSELLL TGPGAEEDPL PLIVQDQYVQ EARQVFEKIQ RMGAQQDDGS
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 M'95 -w border
---------- ---------- ---------- ---------- ---------- ---------- * 21 M'95 -w register
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 M'95 +w polar
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 MTK  -w class
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 28 M'95 -w signif.
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~11 1111111111 11~~~~~~~~ * 14 M'95 -w local

    .    |     .    |     .    |     .    |     .    |     .  120
DAPPGSPDWA GDVTRGQRSQ EELSGPESSL TDEGIGADPE PPVAAFCGLG TTGMWRPLSS
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 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
SSAQTNHHGP GTEDSLGGWA LVSPETPPTP GALRRRRKVP PSGSGGSEFS NGEAGEAYRS
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 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
LSDPIPQRHR AATSEEPTGF SVDSNLLGSL SPKTGLPATS AMDEGLTSGH SDWSVGSEES
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 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
KGYQEVIQSI VQGPGTLGRV VDDRIAGKAP KKKSLSDPSR RGELAGPGFE GPGGEPIREV
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 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
EPMLPPSSSE PILVEQRAEP EEPGATRSRA QSERALPEAL PPPATAHRNF HLDPKLADIL
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 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
SPRLIRRGSK KRPARSSHQE LRRDEGSQDQ TGSLSRARPS SRHVRHASVP ATFMPIVVPE
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 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
PPTSVGPPVA VPEPIGFPTR AHPTLQAPSL EDVTKQYMLN LHSGEVPAPV PVDMPCLPLA
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 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

    .    |     .    |     .    |     .    |     .    |     .  540
APPSAEAKPP EAARPADEPT PASKCCSKPQ VDMRKHVAMT LLDTEQSYVE SLRTLMQGYM
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 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
QPLKQPENSV LCDPSLVDEI FDQIPELLEH HEQFLEQVRH CMQTWHAQQK VGALLVQSFS
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 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
KDVLVNIYSA YIDNFLNAKD AVRVAKEARP AFLKFLEQSM RENKEKQALS DLMIKPVQRI
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 M'95 -w border
---------- ---------- ---------- ---------- ---------- ---------- * 21 M'95 -w register
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 M'95 +w polar
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 MTK  -w class
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 28 M'95 -w signif.
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~66666666 6666662221 ~~~~~~~~~~ * 14 M'95 -w local

    .    |     .    |     .    |     .    |     .    |     .  720
PRYELLVKDL LKHTPEDHPD HPLLLEAQRN IKQVAERINK GVRSAEEAER HARVLQEIEA
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~33366 6666666666 6666699999 9999999999 * 21 M'95 -w border
---------- ---------- -----cdefg abcdefgabc defgafgabc defgabcdef * 21 M'95 -w register
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~266666 6666666666 * 21 M'95 +w polar
~~~~~~~~~~ ~~~~~~~~~~ ~~11144466 6666666666 6666666666 5555555555 * 21 MTK  -w class
~~~~~~~~~~ ~~~~~~~~~~ ~~99999999 9999999999 9999999999 9998888333 * 28 M'95 -w signif.
~~~~~~~~~~ ~~~~~~~~~~ ~~37799999 9999999999 1145588888 8888888885 * 14 M'95 -w local

    .    |     .    |     .    |     .    |     .    |     .  780
HIEGMEDLQA PLRRFLRQEM VIEVKAIGGK KDRSLFLFTD LIVCTTLKRK SGSLRRSSMS
9999999553 ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 M'95 -w border
gabcdefdef ---------- ---------- ---------- ---------- ---------- * 21 M'95 -w register
6666665221 ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 M'95 +w polar
5555554221 ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 MTK  -w class
3333333333 ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 28 M'95 -w signif.
5555555553 ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 14 M'95 -w local

    .    |     .    |     .    |     .    |     .    |     .  840
LYTAASVIDT ASKYKMLWKL PLEDADIIKG ASQATNRENI QKAISRLDED LTTLGQMSKL
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ 3334444444 4444444444 4444441~~~ * 21 M'95 -w border
---------- ---------- ---------- bcdefgabcd efgabcdefg abcdefg--- * 21 M'95 -w register
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ 3338888888 8888888888 888888611~ * 21 M'95 +w polar
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ 1112222222 2222222222 2222221~~~ * 21 MTK  -w class
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 28 M'95 -w signif.
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~1166888 8888888888 8333~~~~~~ * 14 M'95 -w local

    .    |     .    |     .    |     .    |     .    |     .  900
SESLGFPHQS LDDALRDLSA AMHRDLSEKQ ALCYALSFPP TKLELCATRP EGTDSYIFEF
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 M'95 -w border
---------- ---------- ---------- ---------- ---------- ---------- * 21 M'95 -w register
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 M'95 +w polar
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 MTK  -w class
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 28 M'95 -w signif.
~~~~~~~~22 2222222222 2211~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 14 M'95 -w local

    .    |     .    |     .    |     .    |     .    |     .  960
PHPDARLGFE QAFDEAKRKL ASSKSCLDPE FLKAIPIMKT RSGMQFSCAA PTLNSCPEPS
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 M'95 -w border
---------- ---------- ---------- ---------- ---------- ---------- * 21 M'95 -w register
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 M'95 +w polar
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 MTK  -w class
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 28 M'95 -w signif.
~~~~~~~366 6666666666 6652~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 14 M'95 -w local

    .    |     .    |     .    |     .    |     .    |     . 1020
PEVWVCNSDG YVGQVCLLSL RAEPDVEACI AVCSARILCI GAVPGLQPRC HREPPPSLRS
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 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

    .    |     .    |     .    |     .    |     .    |     . 1080
PPETAPEPAG PELDVEAAAD EEAATLAEPG PQPCLHISIA GSGLEMTPGL GEGDPRPELV
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 M'95 -w border
---------- ---------- ---------- ---------- ---------- ---------- * 21 M'95 -w register
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 M'95 +w polar
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 21 MTK  -w class
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 28 M'95 -w signif.
~~~~~~~~~~ ~~~~111111 11111111~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 14 M'95 -w local

    .    |     .    |     .    |     .    |     .    |     . 1140
PFDSDSDDES SPSPSGTLQS QASRSTISSS FGNEETPSSK EATAETTSSE EEQEPGFLPL
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 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
SGSFGPGGPC GTSPMDGRAL RRSSHGSFTR GSLEDLLSVD PEAYQSSVWL GTEDGCVHVY
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 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

    .    |     .    |     .    |     .    |     .    |     . 1260
QSSDSIRDRR NSMKLQHAAS VTCILYLNNQ VFVSLANGEL VVYQREAGHF WDPQNFKSVT
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 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

    .    |     .    |     .    |     .    |     .    |     . 1320
LGTQGSPITK MVSVGGRLWC GCQNRVLVLS PDTLQLEHMF YVGQDSSRCV ACMVDSSLGV
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 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

    .    |     .    |     .    |     .    |     .    |     . 1380
WVTLKGSAHV CLYHPDTFEQ LAEVDVTPPV HRMLAGSDAI IRQHKAACLR ITALLVCEEL
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 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

    .    |     .    |     .    |     .    |     .    |     . 1440
LWVGTSAGVV LTMPTSPGTV SCPRAPLSPT GLGQGHTGHV RFLAAVQLPD GFNLLCPTPP
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ * 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

    .    |     .    |     .    | 
PPPDTGPEKL PSLEHRDSPW HRGPAPARPK 
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ 
---------- ---------- ---------- 
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ 
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ 
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ 
~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ 



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

prediction of transmembrane regions with toppred2

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

tem4.___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: tem4.___inter___

 (1 sequences)
MARRLPRTSALKSSSSELLLTGPGAEEDPLPLIVQDQYVQEARQVFEKIQ
RMGAQQDDGSDAPPGSPDWAGDVTRGQRSQEELSGPESSLTDEGIGADPE
PPVAAFCGLGTTGMWRPLSSSSAQTNHHGPGTEDSLGGWALVSPETPPTP
GALRRRRKVPPSGSGGSEFSNGEAGEAYRSLSDPIPQRHRAATSEEPTGF
SVDSNLLGSLSPKTGLPATSAMDEGLTSGHSDWSVGSEESKGYQEVIQSI
VQGPGTLGRVVDDRIAGKAPKKKSLSDPSRRGELAGPGFEGPGGEPIREV
EPMLPPSSSEPILVEQRAEPEEPGATRSRAQSERALPEALPPPATAHRNF
HLDPKLADILSPRLIRRGSKKRPARSSHQELRRDEGSQDQTGSLSRARPS
SRHVRHASVPATFMPIVVPEPPTSVGPPVAVPEPIGFPTRAHPTLQAPSL
EDVTKQYMLNLHSGEVPAPVPVDMPCLPLAAPPSAEAKPPEAARPADEPT
PASKCCSKPQVDMRKHVAMTLLDTEQSYVESLRTLMQGYMQPLKQPENSV
LCDPSLVDEIFDQIPELLEHHEQFLEQVRHCMQTWHAQQKVGALLVQSFS
KDVLVNIYSAYIDNFLNAKDAVRVAKEARPAFLKFLEQSMRENKEKQALS
DLMIKPVQRIPRYELLVKDLLKHTPEDHPDHPLLLEAQRNIKQVAERINK
GVRSAEEAERHARVLQEIEAHIEGMEDLQAPLRRFLRQEMVIEVKAIGGK
KDRSLFLFTDLIVCTTLKRKSGSLRRSSMSLYTAASVIDTASKYKMLWKL
PLEDADIIKGASQATNRENIQKAISRLDEDLTTLGQMSKLSESLGFPHQS
LDDALRDLSAAMHRDLSEKQALCYALSFPPTKLELCATRPEGTDSYIFEF
PHPDARLGFEQAFDEAKRKLASSKSCLDPEFLKAIPIMKTRSGMQFSCAA
PTLNSCPEPSPEVWVCNSDGYVGQVCLLSLRAEPDVEACIAVCSARILCI
GAVPGLQPRCHREPPPSLRSPPETAPEPAGPELDVEAAADEEAATLAEPG
PQPCLHISIAGSGLEMTPGLGEGDPRPELVPFDSDSDDESSPSPSGTLQS
QASRSTISSSFGNEETPSSKEATAETTSSEEEQEPGFLPLSGSFGPGGPC
GTSPMDGRALRRSSHGSFTRGSLEDLLSVDPEAYQSSVWLGTEDGCVHVY
QSSDSIRDRRNSMKLQHAASVTCILYLNNQVFVSLANGELVVYQREAGHF
WDPQNFKSVTLGTQGSPITKMVSVGGRLWCGCQNRVLVLSPDTLQLEHMF
YVGQDSSRCVACMVDSSLGVWVTLKGSAHVCLYHPDTFEQLAEVDVTPPV
HRMLAGSDAIIRQHKAACLRITALLVCEELLWVGTSAGVVLTMPTSPGTV
SCPRAPLSPTGLGQGHTGHVRFLAAVQLPDGFNLLCPTPPPPPDTGPEKL
PSLEHRDSPWHRGPAPARPKMLVISGGDGYEDFRLSSGGGSSSETVGRDD
STNHLLLWRV


(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 128 structures are to be tested


Candidate membrane-spanning segments:

 Helix Begin   End   Score Certainity
     1    99   119   0.874 Putative
     2   407   427   0.848 Putative
     3   465   485   0.617 Putative
     4   987  1007   0.903 Putative
     5  1052  1072   0.668 Putative
     6  1135  1155   1.224 Certain
     7  1217  1237   0.725 Putative
     8  1308  1328   0.636 Putative
     9  1380  1400   1.539 Certain
Warning! Only printing 30 out of 128 structures


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

Transmembrane segments included in this structure:
     Segment       1     2     3     4     5     6     7     8     9
 Loop length    98   287    37   501    44    62    61    70    51   110
 K+R profile     +        2.00        3.00           +        5.00      
                       +           +           +           +           +      
CYT-EXT prof -0.22           -           -        0.28           -      
                    0.43       -0.50        0.19        1.18        0.57      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: 10.00
Tm probability: 0.00
-> Orientation: N-in

Charge-difference over N-terminal Tm (+-15 residues): -6.00
 (NEG-POS)/(NEG+POS): 0.3571
                 NEG: 19.0000
                 POS: 9.0000
-> Orientation: N-out

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

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

Transmembrane segments included in this structure:
     Segment       2     3     4     5     6     7     8     9
 Loop length   406    37   501    44    62    61    70    51   110
 K+R profile     +           +           +           +           +      
                    2.00        3.00           +        5.00      
CYT-EXT prof  0.32       -0.50        0.19        1.18        0.57      
                       -           -        0.28           -      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: -10.00
Tm probability: 0.00
-> Orientation: N-out

Charge-difference over N-terminal Tm (+-15 residues): 4.00
 (NEG-POS)/(NEG+POS): 0.0769
                 NEG: 56.0000
                 POS: 48.0000
-> Orientation: N-in

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

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

Transmembrane segments included in this structure:
     Segment       1     3     4     5     6     7     8     9
 Loop length    98   345   501    44    62    61    70    51   110
 K+R profile     +           +           +           +           +      
                       +        3.00           +        5.00      
CYT-EXT prof -0.22       -0.50        0.19        1.18        0.57      
                    0.73           -        0.28           -      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: -8.00
Tm probability: 0.00
-> Orientation: N-out

Charge-difference over N-terminal Tm (+-15 residues): -6.00
 (NEG-POS)/(NEG+POS): 0.3571
                 NEG: 19.0000
                 POS: 9.0000
-> Orientation: N-out

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

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

Transmembrane segments included in this structure:
     Segment       3     4     5     6     7     8     9
 Loop length   464   501    44    62    61    70    51   110
 K+R profile     +        3.00           +        5.00      
                       +           +           +           +      
CYT-EXT prof  0.57           -        0.28           -      
                   -0.50        0.19        1.18        0.57      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: 8.00
Tm probability: 0.00
-> Orientation: N-in

Charge-difference over N-terminal Tm (+-15 residues): 0.00
 (NEG-POS)/(NEG+POS): 0.0991
                 NEG: 61.0000
                 POS: 50.0000
-> Orientation: undecided

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

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

Transmembrane segments included in this structure:
     Segment       1     2     4     5     6     7     8     9
 Loop length    98   287   559    44    62    61    70    51   110
 K+R profile     +           +           +           +           +      
                       +        3.00           +        5.00      
CYT-EXT prof -0.22       -0.22        0.19        1.18        0.57      
                    0.43           -        0.28           -      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: -8.00
Tm probability: 0.00
-> Orientation: N-out

Charge-difference over N-terminal Tm (+-15 residues): -6.00
 (NEG-POS)/(NEG+POS): 0.3571
                 NEG: 19.0000
                 POS: 9.0000
-> Orientation: N-out

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

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

Transmembrane segments included in this structure:
     Segment       2     4     5     6     7     8     9
 Loop length   406   559    44    62    61    70    51   110
 K+R profile     +        3.00           +        5.00      
                       +           +           +           +      
CYT-EXT prof  0.32           -        0.28           -      
                   -0.22        0.19        1.18        0.57      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: 8.00
Tm probability: 0.00
-> Orientation: N-in

Charge-difference over N-terminal Tm (+-15 residues): 4.00
 (NEG-POS)/(NEG+POS): 0.0769
                 NEG: 56.0000
                 POS: 48.0000
-> Orientation: N-in

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

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

Transmembrane segments included in this structure:
     Segment       1     4     5     6     7     8     9
 Loop length    98   867    44    62    61    70    51   110
 K+R profile     +        3.00           +        5.00      
                       +           +           +           +      
CYT-EXT prof -0.22           -        0.28           -      
                    0.21        0.19        1.18        0.57      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: 8.00
Tm probability: 0.00
-> Orientation: N-in

Charge-difference over N-terminal Tm (+-15 residues): -6.00
 (NEG-POS)/(NEG+POS): 0.3571
                 NEG: 19.0000
                 POS: 9.0000
-> Orientation: N-out

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

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

Transmembrane segments included in this structure:
     Segment       4     5     6     7     8     9
 Loop length   986    44    62    61    70    51   110
 K+R profile     +           +           +           +      
                    3.00           +        5.00      
CYT-EXT prof  0.17        0.19        1.18        0.57      
                       -        0.28           -      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: -8.00
Tm probability: 0.00
-> Orientation: N-out

Charge-difference over N-terminal Tm (+-15 residues): -4.00
 (NEG-POS)/(NEG+POS): 0.0837
                 NEG: 136.0000
                 POS: 115.0000
-> Orientation: N-out

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

----------------------------------------------------------------------
Structure 9

Transmembrane segments included in this structure:
     Segment       1     2     3     6     7     8     9
 Loop length    98   287    37   649    61    70    51   110
 K+R profile     +        2.00           +        5.00      
                       +           +           +           +      
CYT-EXT prof -0.22           -        0.28           -      
                    0.43       -0.44        1.18        0.57      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: 7.00
Tm probability: 0.00
-> Orientation: N-in

Charge-difference over N-terminal Tm (+-15 residues): -6.00
 (NEG-POS)/(NEG+POS): 0.3571
                 NEG: 19.0000
                 POS: 9.0000
-> Orientation: N-out

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

----------------------------------------------------------------------
Structure 10

Transmembrane segments included in this structure:
     Segment       2     3     6     7     8     9
 Loop length   406    37   649    61    70    51   110
 K+R profile     +           +           +           +      
                    2.00           +        5.00      
CYT-EXT prof  0.32       -0.44        1.18        0.57      
                       -        0.28           -      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: -7.00
Tm probability: 0.00
-> Orientation: N-out

Charge-difference over N-terminal Tm (+-15 residues): 4.00
 (NEG-POS)/(NEG+POS): 0.0769
                 NEG: 56.0000
                 POS: 48.0000
-> Orientation: N-in

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

----------------------------------------------------------------------
Structure 11

Transmembrane segments included in this structure:
     Segment       1     2     3     5     6     8     9
 Loop length    98   287    37   566    62   152    51   110
 K+R profile     +        2.00           +        5.00      
                       +           +           +           +      
CYT-EXT prof -0.22           -        0.19           -      
                    0.43       -0.57        1.05        0.57      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: 7.00
Tm probability: 0.00
-> Orientation: N-in

Charge-difference over N-terminal Tm (+-15 residues): -6.00
 (NEG-POS)/(NEG+POS): 0.3571
                 NEG: 19.0000
                 POS: 9.0000
-> Orientation: N-out

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

----------------------------------------------------------------------
Structure 12

Transmembrane segments included in this structure:
     Segment       2     3     5     6     8     9
 Loop length   406    37   566    62   152    51   110
 K+R profile     +           +           +           +      
                    2.00           +        5.00      
CYT-EXT prof  0.32       -0.57        1.05        0.57      
                       -        0.19           -      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: -7.00
Tm probability: 0.00
-> Orientation: N-out

Charge-difference over N-terminal Tm (+-15 residues): 4.00
 (NEG-POS)/(NEG+POS): 0.0769
                 NEG: 56.0000
                 POS: 48.0000
-> Orientation: N-in

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

----------------------------------------------------------------------
Structure 13

Transmembrane segments included in this structure:
     Segment       1     2     3     4     6     8     9
 Loop length    98   287    37   501   127   152    51   110
 K+R profile     +        2.00           +        5.00      
                       +           +           +           +      
CYT-EXT prof -0.22           -       -0.06           -      
                    0.43       -0.50        1.05        0.57      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: 7.00
Tm probability: 0.00
-> Orientation: N-in

Charge-difference over N-terminal Tm (+-15 residues): -6.00
 (NEG-POS)/(NEG+POS): 0.3571
                 NEG: 19.0000
                 POS: 9.0000
-> Orientation: N-out

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

----------------------------------------------------------------------
Structure 14

Transmembrane segments included in this structure:
     Segment       2     3     4     6     8     9
 Loop length   406    37   501   127   152    51   110
 K+R profile     +           +           +           +      
                    2.00           +        5.00      
CYT-EXT prof  0.32       -0.50        1.05        0.57      
                       -       -0.06           -      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: -7.00
Tm probability: 0.00
-> Orientation: N-out

Charge-difference over N-terminal Tm (+-15 residues): 4.00
 (NEG-POS)/(NEG+POS): 0.0769
                 NEG: 56.0000
                 POS: 48.0000
-> Orientation: N-in

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

----------------------------------------------------------------------
Structure 15

Transmembrane segments included in this structure:
     Segment       1     3     5     6     7     8     9
 Loop length    98   345   566    62    61    70    51   110
 K+R profile     +           +           +        5.00      
                       +           +           +           +      
CYT-EXT prof -0.22       -0.57        0.28           -      
                    0.73        0.19        1.18        0.57      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: 5.00
Tm probability: 0.00
-> Orientation: N-in

Charge-difference over N-terminal Tm (+-15 residues): -6.00
 (NEG-POS)/(NEG+POS): 0.3571
                 NEG: 19.0000
                 POS: 9.0000
-> Orientation: N-out

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

----------------------------------------------------------------------
Structure 16

Transmembrane segments included in this structure:
     Segment       3     5     6     7     8     9
 Loop length   464   566    62    61    70    51   110
 K+R profile     +           +           +           +      
                       +           +        5.00      
CYT-EXT prof  0.57        0.19        1.18        0.57      
                   -0.57        0.28           -      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: -5.00
Tm probability: 0.00
-> Orientation: N-out

Charge-difference over N-terminal Tm (+-15 residues): 0.00
 (NEG-POS)/(NEG+POS): 0.0991
                 NEG: 61.0000
                 POS: 50.0000
-> Orientation: undecided

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

----------------------------------------------------------------------
Structure 17

Transmembrane segments included in this structure:
     Segment       1     2     5     6     7     8     9
 Loop length    98   287   624    62    61    70    51   110
 K+R profile     +           +           +        5.00      
                       +           +           +           +      
CYT-EXT prof -0.22       -0.31        0.28           -      
                    0.43        0.19        1.18        0.57      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: 5.00
Tm probability: 0.00
-> Orientation: N-in

Charge-difference over N-terminal Tm (+-15 residues): -6.00
 (NEG-POS)/(NEG+POS): 0.3571
                 NEG: 19.0000
                 POS: 9.0000
-> Orientation: N-out

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

----------------------------------------------------------------------
Structure 18

Transmembrane segments included in this structure:
     Segment       2     5     6     7     8     9
 Loop length   406   624    62    61    70    51   110
 K+R profile     +           +           +           +      
                       +           +        5.00      
CYT-EXT prof  0.32        0.19        1.18        0.57      
                   -0.31        0.28           -      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: -5.00
Tm probability: 0.00
-> Orientation: N-out

Charge-difference over N-terminal Tm (+-15 residues): 4.00
 (NEG-POS)/(NEG+POS): 0.0769
                 NEG: 56.0000
                 POS: 48.0000
-> Orientation: N-in

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

----------------------------------------------------------------------
Structure 19

Transmembrane segments included in this structure:
     Segment       1     5     6     7     8     9
 Loop length    98   932    62    61    70    51   110
 K+R profile     +           +           +           +      
                       +           +        5.00      
CYT-EXT prof -0.22        0.19        1.18        0.57      
                    0.11        0.28           -      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: -5.00
Tm probability: 0.00
-> Orientation: N-out

Charge-difference over N-terminal Tm (+-15 residues): -6.00
 (NEG-POS)/(NEG+POS): 0.3571
                 NEG: 19.0000
                 POS: 9.0000
-> Orientation: N-out

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

----------------------------------------------------------------------
Structure 20

Transmembrane segments included in this structure:
     Segment       5     6     7     8     9
 Loop length  1051    62    61    70    51   110
 K+R profile     +           +        5.00      
                       +           +           +      
CYT-EXT prof  0.09        0.28           -      
                    0.19        1.18        0.57      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: 5.00
Tm probability: 0.00
-> Orientation: N-in

Charge-difference over N-terminal Tm (+-15 residues): -1.00
 (NEG-POS)/(NEG+POS): 0.1019
                 NEG: 146.0000
                 POS: 119.0000
-> Orientation: N-out

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

----------------------------------------------------------------------
Structure 21

Transmembrane segments included in this structure:
     Segment       1     3     4     6     7     8     9
 Loop length    98   345   501   127    61    70    51   110
 K+R profile     +           +           +        5.00      
                       +           +           +           +      
CYT-EXT prof -0.22       -0.50        0.28           -      
                    0.73       -0.06        1.18        0.57      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: 5.00
Tm probability: 0.00
-> Orientation: N-in

Charge-difference over N-terminal Tm (+-15 residues): -6.00
 (NEG-POS)/(NEG+POS): 0.3571
                 NEG: 19.0000
                 POS: 9.0000
-> Orientation: N-out

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

----------------------------------------------------------------------
Structure 22

Transmembrane segments included in this structure:
     Segment       3     4     6     7     8     9
 Loop length   464   501   127    61    70    51   110
 K+R profile     +           +           +           +      
                       +           +        5.00      
CYT-EXT prof  0.57       -0.06        1.18        0.57      
                   -0.50        0.28           -      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: -5.00
Tm probability: 0.00
-> Orientation: N-out

Charge-difference over N-terminal Tm (+-15 residues): 0.00
 (NEG-POS)/(NEG+POS): 0.0991
                 NEG: 61.0000
                 POS: 50.0000
-> Orientation: undecided

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

----------------------------------------------------------------------
Structure 23

Transmembrane segments included in this structure:
     Segment       1     2     4     6     7     8     9
 Loop length    98   287   559   127    61    70    51   110
 K+R profile     +           +           +        5.00      
                       +           +           +           +      
CYT-EXT prof -0.22       -0.22        0.28           -      
                    0.43       -0.06        1.18        0.57      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: 5.00
Tm probability: 0.01
-> Orientation: N-in

Charge-difference over N-terminal Tm (+-15 residues): -6.00
 (NEG-POS)/(NEG+POS): 0.3571
                 NEG: 19.0000
                 POS: 9.0000
-> Orientation: N-out

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

----------------------------------------------------------------------
Structure 24

Transmembrane segments included in this structure:
     Segment       2     4     6     7     8     9
 Loop length   406   559   127    61    70    51   110
 K+R profile     +           +           +           +      
                       +           +        5.00      
CYT-EXT prof  0.32       -0.06        1.18        0.57      
                   -0.22        0.28           -      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: -5.00
Tm probability: 0.01
-> Orientation: N-out

Charge-difference over N-terminal Tm (+-15 residues): 4.00
 (NEG-POS)/(NEG+POS): 0.0769
                 NEG: 56.0000
                 POS: 48.0000
-> Orientation: N-in

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

----------------------------------------------------------------------
Structure 25

Transmembrane segments included in this structure:
     Segment       1     4     6     7     8     9
 Loop length    98   867   127    61    70    51   110
 K+R profile     +           +           +           +      
                       +           +        5.00      
CYT-EXT prof -0.22       -0.06        1.18        0.57      
                    0.21        0.28           -      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: -5.00
Tm probability: 0.01
-> Orientation: N-out

Charge-difference over N-terminal Tm (+-15 residues): -6.00
 (NEG-POS)/(NEG+POS): 0.3571
                 NEG: 19.0000
                 POS: 9.0000
-> Orientation: N-out

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

----------------------------------------------------------------------
Structure 26

Transmembrane segments included in this structure:
     Segment       4     6     7     8     9
 Loop length   986   127    61    70    51   110
 K+R profile     +           +        5.00      
                       +           +           +      
CYT-EXT prof  0.17        0.28           -      
                   -0.06        1.18        0.57      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: 5.00
Tm probability: 0.02
-> Orientation: N-in

Charge-difference over N-terminal Tm (+-15 residues): -4.00
 (NEG-POS)/(NEG+POS): 0.0837
                 NEG: 136.0000
                 POS: 115.0000
-> Orientation: N-out

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

----------------------------------------------------------------------
Structure 27

Transmembrane segments included in this structure:
     Segment       1     3     6     7     8     9
 Loop length    98   345   649    61    70    51   110
 K+R profile     +           +           +           +      
                       +           +        5.00      
CYT-EXT prof -0.22       -0.44        1.18        0.57      
                    0.73        0.28           -      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: -5.00
Tm probability: 0.00
-> Orientation: N-out

Charge-difference over N-terminal Tm (+-15 residues): -6.00
 (NEG-POS)/(NEG+POS): 0.3571
                 NEG: 19.0000
                 POS: 9.0000
-> Orientation: N-out

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

----------------------------------------------------------------------
Structure 28

Transmembrane segments included in this structure:
     Segment       3     6     7     8     9
 Loop length   464   649    61    70    51   110
 K+R profile     +           +        5.00      
                       +           +           +      
CYT-EXT prof  0.57        0.28           -      
                   -0.44        1.18        0.57      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: 5.00
Tm probability: 0.00
-> Orientation: N-in

Charge-difference over N-terminal Tm (+-15 residues): 0.00
 (NEG-POS)/(NEG+POS): 0.0991
                 NEG: 61.0000
                 POS: 50.0000
-> Orientation: undecided

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

----------------------------------------------------------------------
Structure 29

Transmembrane segments included in this structure:
     Segment       1     2     6     7     8     9
 Loop length    98   287   707    61    70    51   110
 K+R profile     +           +           +           +      
                       +           +        5.00      
CYT-EXT prof -0.22       -0.23        1.18        0.57      
                    0.43        0.28           -      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: -5.00
Tm probability: 0.01
-> Orientation: N-out

Charge-difference over N-terminal Tm (+-15 residues): -6.00
 (NEG-POS)/(NEG+POS): 0.3571
                 NEG: 19.0000
                 POS: 9.0000
-> Orientation: N-out

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

----------------------------------------------------------------------
Structure 30

Transmembrane segments included in this structure:
     Segment       2     6     7     8     9
 Loop length   406   707    61    70    51   110
 K+R profile     +           +        5.00      
                       +           +           +      
CYT-EXT prof  0.32        0.28           -      
                   -0.23        1.18        0.57      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: 5.00
Tm probability: 0.02
-> Orientation: N-in

Charge-difference over N-terminal Tm (+-15 residues): 4.00
 (NEG-POS)/(NEG+POS): 0.0769
                 NEG: 56.0000
                 POS: 48.0000
-> Orientation: N-in

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

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

"tem4" 1510 
 99 119 #f 0.873958
 407 427 #f 0.847917
 465 485 #f 0.616667
 987 1007 #f 0.903125
 1052 1072 #f 0.667708
 1135 1155 #t 1.22396
 1217 1237 #f 0.725
 1308 1328 #f 0.636458
 1380 1400 #t 1.53854



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

tem4.___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: tem4.___inter___

 (1 sequences)
MARRLPRTSALKSSSSELLLTGPGAEEDPLPLIVQDQYVQEARQVFEKIQ
RMGAQQDDGSDAPPGSPDWAGDVTRGQRSQEELSGPESSLTDEGIGADPE
PPVAAFCGLGTTGMWRPLSSSSAQTNHHGPGTEDSLGGWALVSPETPPTP
GALRRRRKVPPSGSGGSEFSNGEAGEAYRSLSDPIPQRHRAATSEEPTGF
SVDSNLLGSLSPKTGLPATSAMDEGLTSGHSDWSVGSEESKGYQEVIQSI
VQGPGTLGRVVDDRIAGKAPKKKSLSDPSRRGELAGPGFEGPGGEPIREV
EPMLPPSSSEPILVEQRAEPEEPGATRSRAQSERALPEALPPPATAHRNF
HLDPKLADILSPRLIRRGSKKRPARSSHQELRRDEGSQDQTGSLSRARPS
SRHVRHASVPATFMPIVVPEPPTSVGPPVAVPEPIGFPTRAHPTLQAPSL
EDVTKQYMLNLHSGEVPAPVPVDMPCLPLAAPPSAEAKPPEAARPADEPT
PASKCCSKPQVDMRKHVAMTLLDTEQSYVESLRTLMQGYMQPLKQPENSV
LCDPSLVDEIFDQIPELLEHHEQFLEQVRHCMQTWHAQQKVGALLVQSFS
KDVLVNIYSAYIDNFLNAKDAVRVAKEARPAFLKFLEQSMRENKEKQALS
DLMIKPVQRIPRYELLVKDLLKHTPEDHPDHPLLLEAQRNIKQVAERINK
GVRSAEEAERHARVLQEIEAHIEGMEDLQAPLRRFLRQEMVIEVKAIGGK
KDRSLFLFTDLIVCTTLKRKSGSLRRSSMSLYTAASVIDTASKYKMLWKL
PLEDADIIKGASQATNRENIQKAISRLDEDLTTLGQMSKLSESLGFPHQS
LDDALRDLSAAMHRDLSEKQALCYALSFPPTKLELCATRPEGTDSYIFEF
PHPDARLGFEQAFDEAKRKLASSKSCLDPEFLKAIPIMKTRSGMQFSCAA
PTLNSCPEPSPEVWVCNSDGYVGQVCLLSLRAEPDVEACIAVCSARILCI
GAVPGLQPRCHREPPPSLRSPPETAPEPAGPELDVEAAADEEAATLAEPG
PQPCLHISIAGSGLEMTPGLGEGDPRPELVPFDSDSDDESSPSPSGTLQS
QASRSTISSSFGNEETPSSKEATAETTSSEEEQEPGFLPLSGSFGPGGPC
GTSPMDGRALRRSSHGSFTRGSLEDLLSVDPEAYQSSVWLGTEDGCVHVY
QSSDSIRDRRNSMKLQHAASVTCILYLNNQVFVSLANGELVVYQREAGHF
WDPQNFKSVTLGTQGSPITKMVSVGGRLWCGCQNRVLVLSPDTLQLEHMF
YVGQDSSRCVACMVDSSLGVWVTLKGSAHVCLYHPDTFEQLAEVDVTPPV
HRMLAGSDAIIRQHKAACLRITALLVCEELLWVGTSAGVVLTMPTSPGTV
SCPRAPLSPTGLGQGHTGHVRFLAAVQLPDGFNLLCPTPPPPPDTGPEKL
PSLEHRDSPWHRGPAPARPKMLVISGGDGYEDFRLSSGGGSSSETVGRDD
STNHLLLWRV


(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 128 structures are to be tested


Candidate membrane-spanning segments:

 Helix Begin   End   Score Certainity
     1    99   119   0.874 Putative
     2   407   427   0.848 Putative
     3   465   485   0.617 Putative
     4   987  1007   0.903 Putative
     5  1052  1072   0.668 Putative
     6  1135  1155   1.224 Certain
     7  1217  1237   0.725 Putative
     8  1308  1328   0.636 Putative
     9  1380  1400   1.539 Certain
Warning! Only printing 30 out of 128 structures


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

Transmembrane segments included in this structure:
     Segment       1     2     3     4     5     6     7     8     9
 Loop length    98   287    37   501    44    62    61    70    51   110
 K+R profile     +        2.00        3.00           +        5.00      
                       +           +           +           +           +      
CYT-EXT prof -0.22           -           -        0.28           -      
                    0.43       -0.50        0.19        1.18        0.57      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: 10.00
Tm probability: 0.00
-> Orientation: N-in

Charge-difference over N-terminal Tm (+-15 residues): -6.00
 (NEG-POS)/(NEG+POS): 0.3571
                 NEG: 19.0000
                 POS: 9.0000
-> Orientation: N-out

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

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

Transmembrane segments included in this structure:
     Segment       2     3     4     5     6     7     8     9
 Loop length   406    37   501    44    62    61    70    51   110
 K+R profile     +           +           +           +           +      
                    2.00        3.00           +        5.00      
CYT-EXT prof  0.32       -0.50        0.19        1.18        0.57      
                       -           -        0.28           -      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: -10.00
Tm probability: 0.00
-> Orientation: N-out

Charge-difference over N-terminal Tm (+-15 residues): 4.00
 (NEG-POS)/(NEG+POS): 0.0769
                 NEG: 56.0000
                 POS: 48.0000
-> Orientation: N-in

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

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

Transmembrane segments included in this structure:
     Segment       1     3     4     5     6     7     8     9
 Loop length    98   345   501    44    62    61    70    51   110
 K+R profile     +           +           +           +           +      
                       +        3.00           +        5.00      
CYT-EXT prof -0.22       -0.50        0.19        1.18        0.57      
                    0.73           -        0.28           -      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: -8.00
Tm probability: 0.00
-> Orientation: N-out

Charge-difference over N-terminal Tm (+-15 residues): -6.00
 (NEG-POS)/(NEG+POS): 0.3571
                 NEG: 19.0000
                 POS: 9.0000
-> Orientation: N-out

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

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

Transmembrane segments included in this structure:
     Segment       3     4     5     6     7     8     9
 Loop length   464   501    44    62    61    70    51   110
 K+R profile     +        3.00           +        5.00      
                       +           +           +           +      
CYT-EXT prof  0.57           -        0.28           -      
                   -0.50        0.19        1.18        0.57      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: 8.00
Tm probability: 0.00
-> Orientation: N-in

Charge-difference over N-terminal Tm (+-15 residues): 0.00
 (NEG-POS)/(NEG+POS): 0.0991
                 NEG: 61.0000
                 POS: 50.0000
-> Orientation: undecided

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

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

Transmembrane segments included in this structure:
     Segment       1     2     4     5     6     7     8     9
 Loop length    98   287   559    44    62    61    70    51   110
 K+R profile     +           +           +           +           +      
                       +        3.00           +        5.00      
CYT-EXT prof -0.22       -0.22        0.19        1.18        0.57      
                    0.43           -        0.28           -      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: -8.00
Tm probability: 0.00
-> Orientation: N-out

Charge-difference over N-terminal Tm (+-15 residues): -6.00
 (NEG-POS)/(NEG+POS): 0.3571
                 NEG: 19.0000
                 POS: 9.0000
-> Orientation: N-out

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

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

Transmembrane segments included in this structure:
     Segment       2     4     5     6     7     8     9
 Loop length   406   559    44    62    61    70    51   110
 K+R profile     +        3.00           +        5.00      
                       +           +           +           +      
CYT-EXT prof  0.32           -        0.28           -      
                   -0.22        0.19        1.18        0.57      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: 8.00
Tm probability: 0.00
-> Orientation: N-in

Charge-difference over N-terminal Tm (+-15 residues): 4.00
 (NEG-POS)/(NEG+POS): 0.0769
                 NEG: 56.0000
                 POS: 48.0000
-> Orientation: N-in

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

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

Transmembrane segments included in this structure:
     Segment       1     4     5     6     7     8     9
 Loop length    98   867    44    62    61    70    51   110
 K+R profile     +        3.00           +        5.00      
                       +           +           +           +      
CYT-EXT prof -0.22           -        0.28           -      
                    0.21        0.19        1.18        0.57      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: 8.00
Tm probability: 0.00
-> Orientation: N-in

Charge-difference over N-terminal Tm (+-15 residues): -6.00
 (NEG-POS)/(NEG+POS): 0.3571
                 NEG: 19.0000
                 POS: 9.0000
-> Orientation: N-out

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

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

Transmembrane segments included in this structure:
     Segment       4     5     6     7     8     9
 Loop length   986    44    62    61    70    51   110
 K+R profile     +           +           +           +      
                    3.00           +        5.00      
CYT-EXT prof  0.17        0.19        1.18        0.57      
                       -        0.28           -      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: -8.00
Tm probability: 0.00
-> Orientation: N-out

Charge-difference over N-terminal Tm (+-15 residues): -4.00
 (NEG-POS)/(NEG+POS): 0.0837
                 NEG: 136.0000
                 POS: 115.0000
-> Orientation: N-out

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

----------------------------------------------------------------------
Structure 9

Transmembrane segments included in this structure:
     Segment       1     2     3     6     7     8     9
 Loop length    98   287    37   649    61    70    51   110
 K+R profile     +        2.00           +        5.00      
                       +           +           +           +      
CYT-EXT prof -0.22           -        0.28           -      
                    0.43       -0.44        1.18        0.57      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: 7.00
Tm probability: 0.00
-> Orientation: N-in

Charge-difference over N-terminal Tm (+-15 residues): -6.00
 (NEG-POS)/(NEG+POS): 0.3571
                 NEG: 19.0000
                 POS: 9.0000
-> Orientation: N-out

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

----------------------------------------------------------------------
Structure 10

Transmembrane segments included in this structure:
     Segment       2     3     6     7     8     9
 Loop length   406    37   649    61    70    51   110
 K+R profile     +           +           +           +      
                    2.00           +        5.00      
CYT-EXT prof  0.32       -0.44        1.18        0.57      
                       -        0.28           -      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: -7.00
Tm probability: 0.00
-> Orientation: N-out

Charge-difference over N-terminal Tm (+-15 residues): 4.00
 (NEG-POS)/(NEG+POS): 0.0769
                 NEG: 56.0000
                 POS: 48.0000
-> Orientation: N-in

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

----------------------------------------------------------------------
Structure 11

Transmembrane segments included in this structure:
     Segment       1     2     3     5     6     8     9
 Loop length    98   287    37   566    62   152    51   110
 K+R profile     +        2.00           +        5.00      
                       +           +           +           +      
CYT-EXT prof -0.22           -        0.19           -      
                    0.43       -0.57        1.05        0.57      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: 7.00
Tm probability: 0.00
-> Orientation: N-in

Charge-difference over N-terminal Tm (+-15 residues): -6.00
 (NEG-POS)/(NEG+POS): 0.3571
                 NEG: 19.0000
                 POS: 9.0000
-> Orientation: N-out

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

----------------------------------------------------------------------
Structure 12

Transmembrane segments included in this structure:
     Segment       2     3     5     6     8     9
 Loop length   406    37   566    62   152    51   110
 K+R profile     +           +           +           +      
                    2.00           +        5.00      
CYT-EXT prof  0.32       -0.57        1.05        0.57      
                       -        0.19           -      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: -7.00
Tm probability: 0.00
-> Orientation: N-out

Charge-difference over N-terminal Tm (+-15 residues): 4.00
 (NEG-POS)/(NEG+POS): 0.0769
                 NEG: 56.0000
                 POS: 48.0000
-> Orientation: N-in

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

----------------------------------------------------------------------
Structure 13

Transmembrane segments included in this structure:
     Segment       1     2     3     4     6     8     9
 Loop length    98   287    37   501   127   152    51   110
 K+R profile     +        2.00           +        5.00      
                       +           +           +           +      
CYT-EXT prof -0.22           -       -0.06           -      
                    0.43       -0.50        1.05        0.57      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: 7.00
Tm probability: 0.00
-> Orientation: N-in

Charge-difference over N-terminal Tm (+-15 residues): -6.00
 (NEG-POS)/(NEG+POS): 0.3571
                 NEG: 19.0000
                 POS: 9.0000
-> Orientation: N-out

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

----------------------------------------------------------------------
Structure 14

Transmembrane segments included in this structure:
     Segment       2     3     4     6     8     9
 Loop length   406    37   501   127   152    51   110
 K+R profile     +           +           +           +      
                    2.00           +        5.00      
CYT-EXT prof  0.32       -0.50        1.05        0.57      
                       -       -0.06           -      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: -7.00
Tm probability: 0.00
-> Orientation: N-out

Charge-difference over N-terminal Tm (+-15 residues): 4.00
 (NEG-POS)/(NEG+POS): 0.0769
                 NEG: 56.0000
                 POS: 48.0000
-> Orientation: N-in

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

----------------------------------------------------------------------
Structure 15

Transmembrane segments included in this structure:
     Segment       1     3     5     6     7     8     9
 Loop length    98   345   566    62    61    70    51   110
 K+R profile     +           +           +        5.00      
                       +           +           +           +      
CYT-EXT prof -0.22       -0.57        0.28           -      
                    0.73        0.19        1.18        0.57      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: 5.00
Tm probability: 0.00
-> Orientation: N-in

Charge-difference over N-terminal Tm (+-15 residues): -6.00
 (NEG-POS)/(NEG+POS): 0.3571
                 NEG: 19.0000
                 POS: 9.0000
-> Orientation: N-out

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

----------------------------------------------------------------------
Structure 16

Transmembrane segments included in this structure:
     Segment       3     5     6     7     8     9
 Loop length   464   566    62    61    70    51   110
 K+R profile     +           +           +           +      
                       +           +        5.00      
CYT-EXT prof  0.57        0.19        1.18        0.57      
                   -0.57        0.28           -      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: -5.00
Tm probability: 0.00
-> Orientation: N-out

Charge-difference over N-terminal Tm (+-15 residues): 0.00
 (NEG-POS)/(NEG+POS): 0.0991
                 NEG: 61.0000
                 POS: 50.0000
-> Orientation: undecided

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

----------------------------------------------------------------------
Structure 17

Transmembrane segments included in this structure:
     Segment       1     2     5     6     7     8     9
 Loop length    98   287   624    62    61    70    51   110
 K+R profile     +           +           +        5.00      
                       +           +           +           +      
CYT-EXT prof -0.22       -0.31        0.28           -      
                    0.43        0.19        1.18        0.57      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: 5.00
Tm probability: 0.00
-> Orientation: N-in

Charge-difference over N-terminal Tm (+-15 residues): -6.00
 (NEG-POS)/(NEG+POS): 0.3571
                 NEG: 19.0000
                 POS: 9.0000
-> Orientation: N-out

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

----------------------------------------------------------------------
Structure 18

Transmembrane segments included in this structure:
     Segment       2     5     6     7     8     9
 Loop length   406   624    62    61    70    51   110
 K+R profile     +           +           +           +      
                       +           +        5.00      
CYT-EXT prof  0.32        0.19        1.18        0.57      
                   -0.31        0.28           -      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: -5.00
Tm probability: 0.00
-> Orientation: N-out

Charge-difference over N-terminal Tm (+-15 residues): 4.00
 (NEG-POS)/(NEG+POS): 0.0769
                 NEG: 56.0000
                 POS: 48.0000
-> Orientation: N-in

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

----------------------------------------------------------------------
Structure 19

Transmembrane segments included in this structure:
     Segment       1     5     6     7     8     9
 Loop length    98   932    62    61    70    51   110
 K+R profile     +           +           +           +      
                       +           +        5.00      
CYT-EXT prof -0.22        0.19        1.18        0.57      
                    0.11        0.28           -      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: -5.00
Tm probability: 0.00
-> Orientation: N-out

Charge-difference over N-terminal Tm (+-15 residues): -6.00
 (NEG-POS)/(NEG+POS): 0.3571
                 NEG: 19.0000
                 POS: 9.0000
-> Orientation: N-out

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

----------------------------------------------------------------------
Structure 20

Transmembrane segments included in this structure:
     Segment       5     6     7     8     9
 Loop length  1051    62    61    70    51   110
 K+R profile     +           +        5.00      
                       +           +           +      
CYT-EXT prof  0.09        0.28           -      
                    0.19        1.18        0.57      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: 5.00
Tm probability: 0.00
-> Orientation: N-in

Charge-difference over N-terminal Tm (+-15 residues): -1.00
 (NEG-POS)/(NEG+POS): 0.1019
                 NEG: 146.0000
                 POS: 119.0000
-> Orientation: N-out

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

----------------------------------------------------------------------
Structure 21

Transmembrane segments included in this structure:
     Segment       1     3     4     6     7     8     9
 Loop length    98   345   501   127    61    70    51   110
 K+R profile     +           +           +        5.00      
                       +           +           +           +      
CYT-EXT prof -0.22       -0.50        0.28           -      
                    0.73       -0.06        1.18        0.57      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: 5.00
Tm probability: 0.00
-> Orientation: N-in

Charge-difference over N-terminal Tm (+-15 residues): -6.00
 (NEG-POS)/(NEG+POS): 0.3571
                 NEG: 19.0000
                 POS: 9.0000
-> Orientation: N-out

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

----------------------------------------------------------------------
Structure 22

Transmembrane segments included in this structure:
     Segment       3     4     6     7     8     9
 Loop length   464   501   127    61    70    51   110
 K+R profile     +           +           +           +      
                       +           +        5.00      
CYT-EXT prof  0.57       -0.06        1.18        0.57      
                   -0.50        0.28           -      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: -5.00
Tm probability: 0.00
-> Orientation: N-out

Charge-difference over N-terminal Tm (+-15 residues): 0.00
 (NEG-POS)/(NEG+POS): 0.0991
                 NEG: 61.0000
                 POS: 50.0000
-> Orientation: undecided

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

----------------------------------------------------------------------
Structure 23

Transmembrane segments included in this structure:
     Segment       1     2     4     6     7     8     9
 Loop length    98   287   559   127    61    70    51   110
 K+R profile     +           +           +        5.00      
                       +           +           +           +      
CYT-EXT prof -0.22       -0.22        0.28           -      
                    0.43       -0.06        1.18        0.57      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: 5.00
Tm probability: 0.01
-> Orientation: N-in

Charge-difference over N-terminal Tm (+-15 residues): -6.00
 (NEG-POS)/(NEG+POS): 0.3571
                 NEG: 19.0000
                 POS: 9.0000
-> Orientation: N-out

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

----------------------------------------------------------------------
Structure 24

Transmembrane segments included in this structure:
     Segment       2     4     6     7     8     9
 Loop length   406   559   127    61    70    51   110
 K+R profile     +           +           +           +      
                       +           +        5.00      
CYT-EXT prof  0.32       -0.06        1.18        0.57      
                   -0.22        0.28           -      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: -5.00
Tm probability: 0.01
-> Orientation: N-out

Charge-difference over N-terminal Tm (+-15 residues): 4.00
 (NEG-POS)/(NEG+POS): 0.0769
                 NEG: 56.0000
                 POS: 48.0000
-> Orientation: N-in

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

----------------------------------------------------------------------
Structure 25

Transmembrane segments included in this structure:
     Segment       1     4     6     7     8     9
 Loop length    98   867   127    61    70    51   110
 K+R profile     +           +           +           +      
                       +           +        5.00      
CYT-EXT prof -0.22       -0.06        1.18        0.57      
                    0.21        0.28           -      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: -5.00
Tm probability: 0.01
-> Orientation: N-out

Charge-difference over N-terminal Tm (+-15 residues): -6.00
 (NEG-POS)/(NEG+POS): 0.3571
                 NEG: 19.0000
                 POS: 9.0000
-> Orientation: N-out

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

----------------------------------------------------------------------
Structure 26

Transmembrane segments included in this structure:
     Segment       4     6     7     8     9
 Loop length   986   127    61    70    51   110
 K+R profile     +           +        5.00      
                       +           +           +      
CYT-EXT prof  0.17        0.28           -      
                   -0.06        1.18        0.57      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: 5.00
Tm probability: 0.02
-> Orientation: N-in

Charge-difference over N-terminal Tm (+-15 residues): -4.00
 (NEG-POS)/(NEG+POS): 0.0837
                 NEG: 136.0000
                 POS: 115.0000
-> Orientation: N-out

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

----------------------------------------------------------------------
Structure 27

Transmembrane segments included in this structure:
     Segment       1     3     6     7     8     9
 Loop length    98   345   649    61    70    51   110
 K+R profile     +           +           +           +      
                       +           +        5.00      
CYT-EXT prof -0.22       -0.44        1.18        0.57      
                    0.73        0.28           -      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: -5.00
Tm probability: 0.00
-> Orientation: N-out

Charge-difference over N-terminal Tm (+-15 residues): -6.00
 (NEG-POS)/(NEG+POS): 0.3571
                 NEG: 19.0000
                 POS: 9.0000
-> Orientation: N-out

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

----------------------------------------------------------------------
Structure 28

Transmembrane segments included in this structure:
     Segment       3     6     7     8     9
 Loop length   464   649    61    70    51   110
 K+R profile     +           +        5.00      
                       +           +           +      
CYT-EXT prof  0.57        0.28           -      
                   -0.44        1.18        0.57      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: 5.00
Tm probability: 0.00
-> Orientation: N-in

Charge-difference over N-terminal Tm (+-15 residues): 0.00
 (NEG-POS)/(NEG+POS): 0.0991
                 NEG: 61.0000
                 POS: 50.0000
-> Orientation: undecided

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

----------------------------------------------------------------------
Structure 29

Transmembrane segments included in this structure:
     Segment       1     2     6     7     8     9
 Loop length    98   287   707    61    70    51   110
 K+R profile     +           +           +           +      
                       +           +        5.00      
CYT-EXT prof -0.22       -0.23        1.18        0.57      
                    0.43        0.28           -      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: -5.00
Tm probability: 0.01
-> Orientation: N-out

Charge-difference over N-terminal Tm (+-15 residues): -6.00
 (NEG-POS)/(NEG+POS): 0.3571
                 NEG: 19.0000
                 POS: 9.0000
-> Orientation: N-out

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

----------------------------------------------------------------------
Structure 30

Transmembrane segments included in this structure:
     Segment       2     6     7     8     9
 Loop length   406   707    61    70    51   110
 K+R profile     +           +        5.00      
                       +           +           +      
CYT-EXT prof  0.32        0.28           -      
                   -0.23        1.18        0.57      
For CYT-EXT profile neg. values indicate cytoplasmic preference.


K+R difference: 5.00
Tm probability: 0.02
-> Orientation: N-in

Charge-difference over N-terminal Tm (+-15 residues): 4.00
 (NEG-POS)/(NEG+POS): 0.0769
                 NEG: 56.0000
                 POS: 48.0000
-> Orientation: N-in

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

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

"tem4" 1510 
 99 119 #f 0.873958
 407 427 #f 0.847917
 465 485 #f 0.616667
 987 1007 #f 0.903125
 1052 1072 #f 0.667708
 1135 1155 #t 1.22396
 1217 1237 #f 0.725
 1308 1328 #f 0.636458
 1380 1400 #t 1.53854



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

SAPS.  Version of April 11, 1996.
Date run: Wed Oct 11 20:38:41 2000

File: /people/maria/tem4.___saps___
ID   gi|2224615|dbj|BAA20795.1|
DE   KIAA0337 [Homo sapiens]

number of residues: 1510;   molecular weight: 163.6 kdal
 
         1  MARRLPRTSA LKSSSSELLL TGPGAEEDPL PLIVQDQYVQ EARQVFEKIQ RMGAQQDDGS 
        61  DAPPGSPDWA GDVTRGQRSQ EELSGPESSL TDEGIGADPE PPVAAFCGLG TTGMWRPLSS 
       121  SSAQTNHHGP GTEDSLGGWA LVSPETPPTP GALRRRRKVP PSGSGGSEFS NGEAGEAYRS 
       181  LSDPIPQRHR AATSEEPTGF SVDSNLLGSL SPKTGLPATS AMDEGLTSGH SDWSVGSEES 
       241  KGYQEVIQSI VQGPGTLGRV VDDRIAGKAP KKKSLSDPSR RGELAGPGFE GPGGEPIREV 
       301  EPMLPPSSSE PILVEQRAEP EEPGATRSRA QSERALPEAL PPPATAHRNF HLDPKLADIL 
       361  SPRLIRRGSK KRPARSSHQE LRRDEGSQDQ TGSLSRARPS SRHVRHASVP ATFMPIVVPE 
       421  PPTSVGPPVA VPEPIGFPTR AHPTLQAPSL EDVTKQYMLN LHSGEVPAPV PVDMPCLPLA 
       481  APPSAEAKPP EAARPADEPT PASKCCSKPQ VDMRKHVAMT LLDTEQSYVE SLRTLMQGYM 
       541  QPLKQPENSV LCDPSLVDEI FDQIPELLEH HEQFLEQVRH CMQTWHAQQK VGALLVQSFS 
       601  KDVLVNIYSA YIDNFLNAKD AVRVAKEARP AFLKFLEQSM RENKEKQALS DLMIKPVQRI 
       661  PRYELLVKDL LKHTPEDHPD HPLLLEAQRN IKQVAERINK GVRSAEEAER HARVLQEIEA 
       721  HIEGMEDLQA PLRRFLRQEM VIEVKAIGGK KDRSLFLFTD LIVCTTLKRK SGSLRRSSMS 
       781  LYTAASVIDT ASKYKMLWKL PLEDADIIKG ASQATNRENI QKAISRLDED LTTLGQMSKL 
       841  SESLGFPHQS LDDALRDLSA AMHRDLSEKQ ALCYALSFPP TKLELCATRP EGTDSYIFEF 
       901  PHPDARLGFE QAFDEAKRKL ASSKSCLDPE FLKAIPIMKT RSGMQFSCAA PTLNSCPEPS 
       961  PEVWVCNSDG YVGQVCLLSL RAEPDVEACI AVCSARILCI GAVPGLQPRC HREPPPSLRS 
      1021  PPETAPEPAG PELDVEAAAD EEAATLAEPG PQPCLHISIA GSGLEMTPGL GEGDPRPELV 
      1081  PFDSDSDDES SPSPSGTLQS QASRSTISSS FGNEETPSSK EATAETTSSE EEQEPGFLPL 
      1141  SGSFGPGGPC GTSPMDGRAL RRSSHGSFTR GSLEDLLSVD PEAYQSSVWL GTEDGCVHVY 
      1201  QSSDSIRDRR NSMKLQHAAS VTCILYLNNQ VFVSLANGEL VVYQREAGHF WDPQNFKSVT 
      1261  LGTQGSPITK MVSVGGRLWC GCQNRVLVLS PDTLQLEHMF YVGQDSSRCV ACMVDSSLGV 
      1321  WVTLKGSAHV CLYHPDTFEQ LAEVDVTPPV HRMLAGSDAI IRQHKAACLR ITALLVCEEL 
      1381  LWVGTSAGVV LTMPTSPGTV SCPRAPLSPT GLGQGHTGHV RFLAAVQLPD GFNLLCPTPP 
      1441  PPPDTGPEKL PSLEHRDSPW HRGPAPARPK MLVISGGDGY EDFRLSSGGG SSSETVGRDD 
      1501  STNHLLLWRV

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

A  :121( 8.0%); C  : 31( 2.1%); D  : 81( 5.4%); E  :116( 7.7%); F  : 38( 2.5%)
G  :112( 7.4%); H  : 40( 2.6%); I  : 47( 3.1%); K  : 56( 3.7%); L  :150( 9.9%)
M  : 32( 2.1%); N- : 25( 1.7%); P+ :144( 9.5%); Q  : 68( 4.5%); R  : 95( 6.3%)
S  :156(10.3%); T  : 74( 4.9%); V  : 89( 5.9%); W  : 14( 0.9%); Y  : 21( 1.4%)

KR      :  151 ( 10.0%);   ED      :  197 ( 13.0%);   AGP     :  377 ( 25.0%);
KRED    :  348 ( 23.0%);   KR-ED   :  -46 ( -3.0%);   FIKMNY- :  219 ( 14.5%);
LVIFM   :  356 ( 23.6%);   ST      :  230 ( 15.2%).

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

A. CHARGE CLUSTERS.


Positive charge clusters (cmin = 11/30 or 14/45 or 17/60):  none


Negative charge clusters (cmin = 13/30 or 17/45 or 20/60):  none


Mixed charge clusters (cmin = 18/30 or 24/45 or 30/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.
There are no high scoring uncharged segments.


C. CHARGE RUNS AND PATTERNS.

pattern  (+)|  (-)|  (*)|  (0)| (+0)| (-0)| (*0)|(+00)|(-00)|(*00)| (H.)|(H..)|
lmin0     5 |   6 |   8 |  40 |  10 |  11 |  14 |  12 |  13 |  17 |   7 |  10 | 
lmin1     7 |   7 |  10 |  48 |  12 |  13 |  17 |  15 |  16 |  20 |   9 |  12 | 
lmin2     8 |   9 |  11 |  53 |  14 |  15 |  19 |  17 |  18 |  23 |  10 |  13 | 
 (Significance level: 0.010000; Minimal displayed length:  6)
There are no charge runs or patterns exceeding the given minimal lengths.

Run count statistics:

  +  runs >=   3:   5, at  154;  271;  370;  768;  917;
  -  runs >=   4:   0
  *  runs >=   5:   1, at  154;
  0  runs >=  27:   0

--------------------------------------------------------------------------------
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-106-C-368-C-28-C-C-45-C-28-C-182-C-108-C-12-C-39-C-21-C-7-C-9-C-9-C-12-C-3-C-5-C-10-C-43-C-95-C-45-C-26-C-56-C-1-C-26-C-2-C-18-C-36-C-8-C-24-C-33-C-74-COOH


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


H2N-106-C-19-H-H-60-H-40-H-116-H-3-H-26-H-24-H-2-H-35-H-19-H-13-C-28-C-C-9-H-35-C-17-H-H-8-H-C-4-H-86-H-4-H-2-H-29-H-9-H-42-C-83-H-14-H-9-C-12-C-15-H-23-C-21-C-7-C-9-C-9-C-12-C-3-C-5-C-10-C-H-42-C-1-H-93-C-14-H-30-C-1-H-18-H-5-C-25-H-30-C-1-C-15-H-10-C-2-C-16-H-1-C-2-H-16-H-12-H-3-C-8-C-24-C-13-H-2-H-16-C-18-H-5-H-42-H-6-COOH

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

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

Aligned matching blocks:


[ 180- 184]   SLSDP
[ 274- 278]   SLSDP

______________________________

[ 645- 649]   EKQAL
[ 868- 872]   EKQAL

______________________________

[ 774- 778]   LRRSS
[1160-1164]   LRRSS


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: 111  (Expected range:  65--122)

2. Histogram of spacings between consecutive amino acid multiplets:
   (1-5) 40   (6-10) 26   (11-20) 24   (>=21) 22

3. Clusters of amino acid multiplets (cmin = 12/30 or 16/45 or 19/60):  none

4. Long amino acid multiplets (>= 5; Letter/Length/Position):
    P/5/1439


B. CHARGE ALPHABET.

1. Total number of charge multiplets:  45  (Expected range:  18-- 54)
   15 +plets (f+: 10.0%), 30 -plets (f-: 13.0%)
   Total number of charge altplets: 25 (Critical number: 59)

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

3. Long charge multiplets (>= 5; Letter/Length/Position):
    +/5/154

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

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

Location	Period	Element		Copies	Core	Errors
  13-  16	 1	S         	 4	 4  	 0
  53-  76	 6	G.....    	 4	 4  	 0
 119- 122	 1	S         	 4	 4  	 0
 154- 157	 1	R         	 4	 4  	 0
 228- 242	 3	S..       	 5	 5 !	 0
 278- 313	 9	P....E... 	 4	 4  	/0/././././1/./././
 352- 367	 4	L...      	 4	 4  	 0
 483- 506	 6	P.....    	 4	 4  	 0
 676- 715	10	E.........	 4	 4  	 0
1016-1035	 5	P....     	 4	 4  	 0
1090-1114	 5	S....     	 5	 5 !	 0
1142-1153	 3	G..       	 4	 4  	 0
1439-1443	 1	P         	 5	 5 !	 0
1439-1454	 4	P...      	 4	 4  	 0


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

Location	Period	Element		Copies	Core	Errors
 154- 158	 1	+         	 5	 5  	 0
 522- 570	 7	i.0..00   	 7	 7  	/0/./2/././2/2/
 703- 762	10	*...*.....	 6	 6  	/0/./././2/./././././
1179-1250	 9	i...00... 	 8	 8  	/0/./././2/1/./././
1278-1325	 8	i0..00.0  	 6	 6  	/0/1/././2/2/./2/
1314-1361	 8	i000.0..  	 6	 6  	/0/2/2/2/./1/././


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

Location (Quartile) Spacing     Rank       P-value   Interpretation

 107- 476  (1.)     C( 369)C     1 of  32   0.0049   large maximal spacing
 787- 963  (3.)     V( 176)V     1 of  90   0.0010   large  1. maximal spacing
1080-1179  (3.)     V(  99)V     2 of  90   0.0094   large  2. 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:            tem4
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query:  gi|2224615|dbj|BAA20795.1|  KIAA0337 [Homo sapiens]

Scores for sequence family classification (score includes all domains):
Model      Description                                  Score    E-value  N 
--------   -----------                                  -----    ------- ---
RhoGEF     RhoGEF domain                                155.8    7.5e-43   1
SH2        Src homology domain 2                          4.3        3.3   1
NGF        Nerve growth factor family                     2.1         19   1
Gastrin    Gastrin/cholecystokinin family                 1.1         56   1
COLFI      Fibrillar collagen C-terminal domain          -1.7         78   1
STAT       STAT protein                                  -2.8         79   1
SpoA       Surface presentation of antigens (SPOA) pr   -25.4         87   1
DB         DB module                                    -49.0         64   1
DUF61      Protein of unknown function DUF61            -64.4         52   1
Collagen   Collagen triple helix repeat (20 copies)     -86.3         83   1
Keratin_B2 Keratin, high sulfur B2 protein              -91.0         28   1
Translin   Translin family                              -94.7         73   1
DUF47      Protein of unknown function DUF47           -106.1         62   1

Parsed for domains:
Model      Domain  seq-f seq-t    hmm-f hmm-t      score  E-value
--------   ------- ----- -----    ----- -----      -----  -------
NGF          1/1     276   285 ..     1    10 [.     2.1       19
Collagen     1/1     253   311 ..     1    60 []   -86.3       83
COLFI        1/1     359   376 ..     1    15 [.    -1.7       78
STAT         1/1     452   461 ..   616   625 .]    -2.8       79
SH2          1/1     561   571 ..    69    79 .]     4.3      3.3
RhoGEF       1/1     517   700 ..     1   207 []   155.8  7.5e-43
DUF47        1/1     558   706 ..     1   224 []  -106.1       62
SpoA         1/1     645   711 ..     1    80 []   -25.4       87
Gastrin      1/1     703   712 ..   117   126 .]     1.1       56
DUF61        1/1     686   794 ..     1   140 []   -64.4       52
Translin     1/1     684   839 ..     1   224 []   -94.7       73
DB           1/1     939  1012 ..     1   124 []   -49.0       64
Keratin_B2   1/1     918  1029 ..     1   177 []   -91.0       28

Alignments of top-scoring domains:
NGF: domain 1 of 1, from 276 to 285: score 2.1, E = 19
                   *->epvsrRGElS<-*
                      +++srRGEl+   
  gi|2224615   276    SDPSRRGELA    285  

Collagen: domain 1 of 1, from 253 to 311: score -86.3, E = 83
                   *->GppGppGppGppGppGppGppGpaGapGppGppGepGpPGppGppGp
                        pG  G+  +   +G++  +     p  +G    pG+ Gp G+p  
  gi|2224615   253    -GPGTLGRVVDDRIAGKAPKKKSLSDPSRRGELAGPGFEGPGGEPIR 298  

                   pGppGapGapGpp<-*
                     +p  p +  +p   
  gi|2224615   299 EVEPMLPPSSSEP    311  

COLFI: domain 1 of 1, from 359 to 376: score -1.7, E = 78
                   *->lksP....eGksrknPARt<-*
                      ++sP+  + G s+k PAR+   
  gi|2224615   359    ILSPrlirRG-SKKRPARS    376  

STAT: domain 1 of 1, from 452 to 461: score -2.8, E = 79
                   *->DLiKkHllpl<-*
                      D++K ++l l   
  gi|2224615   452    DVTKQYMLNL    461  

SH2: domain 1 of 1, from 561 to 571: score 4.3, E = 3.3
                   *->FssLpeLvehY<-*
                      F+ +peL eh+   
  gi|2224615   561    FDQIPELLEHH    571  

RhoGEF: domain 1 of 1, from 517 to 700: score 155.8, E = 7.5e-43
                   *->vlkElleTEkkYvrdLeildnvymkpLreaaisskpvltpddietiF
                      v+  ll TE  Yv++L++l ++ym+pL++++ s   + +p  ++ iF
  gi|2224615   517    VAMTLLDTEQSYVESLRTLMQGYMQPLKQPENSV--LCDPSLVDEIF 561  

                   sNiediyefhreFLkssLearisssqfedlDekkiepsaprlGdlFlkl.
                     i++++e h +FL      +++++             ++ +G+l+ ++ 
  gi|2224615   562 DQIPELLEHHEQFLE-QVRHCMQTW-----------HAQQKVGALLVQSf 599  

                   .kepflqvYgeYcsnkpyAqelleklrqaasnpqFaefldeveassntgA
                   +k+   ++Y +Y++n+ +A++ ++    ++++p+F +fl++ + + n++ 
  gi|2224615   600 sKDVLVNIYSAYIDNFLNAKDAVRV--AKEARPAFLKFLEQSMRE-NKE- 645  

                   kddavkltLqsLLlkPvqRilrYpLLLkeLLkhtpegedqpdredLkkal
                        k  L++L++kPvqRi+rY+LL+k+LLkhtpe  d+pd+++L +a+
  gi|2224615   646 -----KQALSDLMIKPVQRIPRYELLVKDLLKHTPE--DHPDHPLLLEAQ 688  

                   dllqdlaksiNe<-*
                   ++++ +a++iN+   
  gi|2224615   689 RNIKQVAERINK    700  

DUF47: domain 1 of 1, from 558 to 706: score -106.1, E = 62
                   *->klflFekspfkplieHlekvvetvlkleelfealirGdyervEellk
                             + f+ + e+le +++ ++  +  ++++ +    + ++ + 
  gi|2224615   558    ------DEIFDQIPELLEHHEQFLEQVRHCMQTWHA----QQKVGAL 594  

                   eIselEreADeikreirlkltsGaFLPvdRtDllelvelqDkiaDaaEdi
                    + +         ++++ +++s                 +D+ ++ a+d+
  gi|2224615   595 LVQSF-------SKDVLVNIYS---------------AYIDNFLN-AKDA 621  

                   AkllllkknLvsiPeelkdeflamlkktvkaaellaeaieeLekllesaf
                      ++++k       e ++ fl++l  ++   +  + a+  L+       
  gi|2224615   622 ---VRVAK-------EARPAFLKFLEQSMRE-NKEKQALSDLM------I 654  

                   sgkEvkkardivkeieelEdeiDqlqrkllkklFsaYredainpkdflil
                   +          v+ i + E         l k+l+    +  ++  d  +l
  gi|2224615   655 K---------PVQRIPRYE--------LLVKDLLK---HTPEDHPDHPLL 684  

                   leiieklgnvSDraEDaadrisiillk<-*
                   le  +++  v+ r+ +      + +++   
  gi|2224615   685 LEAQRNIKQVAERINK-----GVRSAE    706  

SpoA: domain 1 of 1, from 645 to 711: score -25.4, E = 87
                   *->llalLmdipvdlvvelGrteltlreLlrLqpGsVieLdklagdpgvt
                      + + L+d+ ++ v ++ r+el +++Ll+ +p +    +  ++     
  gi|2224615   645    EKQALSDLMIKPVQRIPRYELLVKDLLKHTPED----HPDHPLL--- 684  

                   VdilvngrligkGelvvvgdnlGVrIteiiesk<-*
                   ++   n + ++     ++++  GVr  e  e++   
  gi|2224615   685 LEAQRNIKQVA----ERINK--GVRSAEEAERH    711  

Gastrin: domain 1 of 1, from 703 to 712: score 1.1, E = 56
                   *->RSaEEydyas<-*
                      RSaEE++ ++   
  gi|2224615   703    RSAEEAERHA    712  

DUF61: domain 1 of 1, from 686 to 794: score -64.4, E = 52
                   *->kvEKiikkeIealNkHLPrFkRrsLkE......LLkEkePhVtlkDG
                         + ik   e++N    +   rs +E +++ + L+E e+h+ +   
  gi|2224615   686    EAQRNIKQVAERIN----K-GVRSAEEaerharVLQEIEAHIEG--- 724  

                   keHyiKrrELEfLaelldEDeisrLrlPIvlEiesiGeeraIkIrGkvEV
                             E+L   l       Lr   v+E+  iG++           
  gi|2224615   725 ---------MEDLQAPLRRF----LRQEMVIEVKAIGGK----------- 750  

                   KVIkkvLGkdeegylEKdelylYkPElklAEiRr...kLPTtTQviFrvs
                           kd+ ++l  d ++      k   +Rr++  L T+  vi   s
  gi|2224615   751 --------KDRSLFLFTDLIVCTTLKRKSGSLRRssmSLYTAASVIDTAS 792  

                   ll<-*
                    +   
  gi|2224615   793 KY    794  

Translin: domain 1 of 1, from 684 to 839: score -94.7, E = 73
                   *->ddkREeilKvsReItieSkklIfllqGVHrgsskdeae...erLeka
                            +l   R I + + ++ + +    r++  +eae++  +L ++
  gi|2224615   684    ------LLEAQRNIKQVAERINKGV----RSA--EEAErhaRVLQEI 718  

                   rekleevKiHqhlaslkkelegedyekFsracstglQEyVEAvaFykyle
                   + ++e+         l  +        ++r ++   QE V  v     ++
  gi|2224615   719 EAHIEGME------DLQAP--------LRRFLR---QEMVIEVKAIGGKK 751  

                   tetLpsfdelnKqLiFtTeilgienKtPssDApDkelgtWRLrinpeDYl
                   +++L+ f +l    i  T++ +++                          
  gi|2224615   752 DRSLFLFTDL----IVCTTLKRKS-------------------------- 771  

                   LGvaDltGELrRmaltsvsnGdldeaeriskFmeqlyeeLmeleylPnkv
                          G LrR +    ++ ++      +k + +l  e   +     k+
  gi|2224615   772 -------GSLRRSSMSLYTAASVIDTASKYKMLWKLPLEDADI----IKG 810  

                   vseLrkKldvlrqsleKveedcyalkvrgl<-*
                    s    +  + +  +++  ed+++l +  +   
  gi|2224615   811 ASQATNRENIQK-AISRLDEDLTTLGQMSK    839  

DB: domain 1 of 1, from 939 to 1012: score -49.0, E = 64
                   *->CCesdrkvpeeClskkCsfntysksalnmvvqkFkqlmffkqddsCp
                        ++  ++  +C             +l                 sCp
  gi|2224615   939    --KTRSGMQFSCAA----------PTL----------------NSCP 957  

                   lkgkalpemlfCAAqgrDHtdCClrn....GVttPqTlaGkkCLtfCdqs
                   +     pe++ C ++g     C l+ + +++V +        C + C+++
  gi|2224615   958 E---PSPEVWVCNSDGYVGQVCLLSLraepDVEA--------CIAVCSAR 996  

                   pgPRRFWfvtqldksYvpCldrfesIk.qCFy<-*
                              ++     C++ +    ++C++   
  gi|2224615   997 -----------IL-----CIGAVPGLQpRCHR    1012 

Keratin_B2: domain 1 of 1, from 918 to 1029: score -91.0, E = 28
                   *->qtScCGfptCStlgtrPsCGsscCQPsCCeSCCQpsCcqpSCCqPtc
                           +    S+ +        c+ P            +     P  
  gi|2224615   918    -----RK-LASSKS--------CLDP---------EFLKAI---PIM 938  

                   sqtscCqPtcfqssCCrPsCcqTSCCqPtCcqsssCqtgCgigGsiGyGQ
                    + s  q+     sC  P+     C  P+  + ++   g         GQ
  gi|2224615   939 KTRSGMQF-----SCAAPTLNS--CPEPSPEVWVCNSDGY-------VGQ 974  

                   eGsSGAvScrirWCRPdCrvegtClPpC..CvvsCtaPTCCqpvsaQasC
                   +                C  +    P ++ C+  C a       +  a  
  gi|2224615   975 V----------------CLLSLRAEPDVeaCIAVCSAR----ILCIGAV- 1003 

                   CRPsCqPyCgqsCCRPaCccsvtCtrTccePc<-*
                     P  qP C    + P+ + ++    T  eP    
  gi|2224615  1004 --PGLQPRCH-REPPPSLRSPP---ETAPEPA    1029 

//

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:            tem4
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query:  gi|2224615|dbj|BAA20795.1|  KIAA0337 [Homo sapiens]

Scores for sequence family classification (score includes all domains):
Model          Description                              Score    E-value  N 
--------       -----------                              -----    ------- ---
RhoGEF         RhoGEF domain                            154.0    5.5e-44   1
SH2            Src homology domain 2                      4.3        3.3   1
WD40           WD domain, G-beta repeat                   2.5         85   1
UreD           UreD urease accessory protein              2.2         19   1
NGF            Nerve growth factor family                 2.1         19   1
PPR            PPR repeat                                 1.9         87   1
Gastrin        Gastrin/cholecystokinin family             1.1         56   1
DX             DX module                                  0.5         49   1
Ran_BP1        RanBP1 domain.                             0.4         73   1
Ribosomal_S28e Ribosomal protein S28e                     0.4         77   1
Lambda_CIII    Lambda Phage CIII                          0.4         92   1
HECT           HECT-domain (ubiquitin-transferase).       0.0         53   1
Sua5_yciO_yrdC SUA5/yciO/yrdC family                     -0.1      1e+02   1
tRNA-synt_1d   tRNA synthetases class I (R)              -0.3         22   1
lig_chan       Ligand-gated ion channel                  -0.5         99   1
COLFI          Fibrillar collagen C-terminal domain      -1.7         78   1
STAT           STAT protein                              -2.8         79   1

Parsed for domains:
Model          Domain  seq-f seq-t    hmm-f hmm-t      score  E-value
--------       ------- ----- -----    ----- -----      -----  -------
PPR              1/1      39    58 ..    16    35 .]     1.9       87
NGF              1/1     276   285 ..     1    10 [.     2.1       19
Sua5_yciO_yrdC   1/1     358   368 ..   179   191 .]    -0.1    1e+02
COLFI            1/1     359   376 ..     1    15 [.    -1.7       78
STAT             1/1     452   461 ..   616   625 .]    -2.8       79
SH2              1/1     561   571 ..    69    79 .]     4.3      3.3
UreD             1/1     569   583 ..   171   185 .]     2.2       19
tRNA-synt_1d     1/1     646   660 ..   595   609 .]    -0.3       22
RhoGEF           1/1     517   700 ..     1   207 []   154.0  5.5e-44
Gastrin          1/1     703   712 ..   117   126 .]     1.1       56
Ribosomal_S28e   1/1     707   716 ..    62    71 .]     0.4       77
Ran_BP1          1/1     909   919 ..   117   127 .]     0.4       73
lig_chan         1/1     960   966 ..     1     7 [.    -0.5       99
Lambda_CIII      1/1    1055  1063 ..     1     9 [.     0.4       92
HECT             1/1    1172  1190 ..   104   122 ..     0.0       53
WD40             1/1    1217  1244 ..     7    36 .]     2.5       85
DX               1/1    1329  1335 ..    95   101 .]     0.5       49

Alignments of top-scoring domains:
PPR: domain 1 of 1, from 39 to 58: score 1.9, E = 87
                   *->leeAlelfeeMkekGikPdv<-*
                      ++eA ++fe+++++G   d+   
  gi|2224615    39    VQEARQVFEKIQRMGAQQDD    58   

NGF: domain 1 of 1, from 276 to 285: score 2.1, E = 19
                   *->epvsrRGElS<-*
                      +++srRGEl+   
  gi|2224615   276    SDPSRRGELA    285  

Sua5_yciO_yrdC: domain 1 of 1, from 358 to 368: score -0.1, E = 1e+02
                   *->dltddtpvvlRrG<-*
                      d++   p+++RrG   
  gi|2224615   358    DIL--SPRLIRRG    368  

COLFI: domain 1 of 1, from 359 to 376: score -1.7, E = 78
                   *->lksP....eGksrknPARt<-*
                      ++sP+  + G s+k PAR+   
  gi|2224615   359    ILSPrlirRG-SKKRPARS    376  

STAT: domain 1 of 1, from 452 to 461: score -2.8, E = 79
                   *->DLiKkHllpl<-*
                      D++K ++l l   
  gi|2224615   452    DVTKQYMLNL    461  

SH2: domain 1 of 1, from 561 to 571: score 4.3, E = 3.3
                   *->FssLpeLvehY<-*
                      F+ +peL eh+   
  gi|2224615   561    FDQIPELLEHH    571  

UreD: domain 1 of 1, from 569 to 583: score 2.2, E = 19
                   *->vadeallkqlrelle<-*
                      +++e +l+q+r++++   
  gi|2224615   569    EHHEQFLEQVRHCMQ    583  

tRNA-synt_1d: domain 1 of 1, from 646 to 660: score -0.3, E = 22
                   *->angLdlLGistverM<-*
                      +++L+ L i++v+r    
  gi|2224615   646    KQALSDLMIKPVQRI    660  

RhoGEF: domain 1 of 1, from 517 to 700: score 154.0, E = 5.5e-44
                   *->vlkElleTEkkYvrdLeildnvymkpLreaaisskpvltpddietiF
                      v+  ll TE  Yv++L++l ++ym+pL++++ s   + +p  ++ iF
  gi|2224615   517    VAMTLLDTEQSYVESLRTLMQGYMQPLKQPENSV--LCDPSLVDEIF 561  

                   sNiediyefhreFLkssLearisssqfedlDekkiepsaprlGdlFlkl.
                     i++++e h +FL      +++++             ++ +G+l+ ++ 
  gi|2224615   562 DQIPELLEHHEQFLE-QVRHCMQTW-----------HAQQKVGALLVQSf 599  

                   .kepflqvYgeYcsnkpyAqelleklrqaasnpqFaefldeveassntgA
                   +k+   ++Y +Y++n+ +A++ ++    ++++p+F +fl++ + + n++ 
  gi|2224615   600 sKDVLVNIYSAYIDNFLNAKDAVRV--AKEARPAFLKFLEQSMRE-NKE- 645  

                   kddavkltLqsLLlkPvqRilrYpLLLkeLLkhtpegedqpdredLkkal
                        k  L++L++kPvqRi+rY+LL+k+LLkhtpe  d+pd+++L +a+
  gi|2224615   646 -----KQALSDLMIKPVQRIPRYELLVKDLLKHTPE--DHPDHPLLLEAQ 688  

                   dllqdlaksiNe<-*
                   ++++ +a++iN+   
  gi|2224615   689 RNIKQVAERINK    700  

Gastrin: domain 1 of 1, from 703 to 712: score 1.1, E = 56
                   *->RSaEEydyas<-*
                      RSaEE++ ++   
  gi|2224615   703    RSAEEAERHA    712  

Ribosomal_S28e: domain 1 of 1, from 707 to 716: score 0.4, E = 77
                   *->EsEREARrLr<-*
                      E ER AR L    
  gi|2224615   707    EAERHARVLQ    716  

Ran_BP1: domain 1 of 1, from 909 to 919: score 0.4, E = 73
                   *->FKkkFEeaqke<-*
                      F ++F+ea+++   
  gi|2224615   909    FEQAFDEAKRK    919  

lig_chan: domain 1 of 1, from 960 to 966: score -0.5, E = 99
                   *->SseiWmm<-*
                      S+e+W++   
  gi|2224615   960    SPEVWVC    966  

Lambda_CIII: domain 1 of 1, from 1055 to 1063: score 0.4, E = 92
                   *->MhfaiAGsg<-*
                      +h+ iAGsg   
  gi|2224615  1055    LHISIAGSG    1063 

HECT: domain 1 of 1, from 1172 to 1190: score 0.0, E = 53
                   *->sfhDlesvDpelYnsLiki<-*
                      s++Dl svDpe Y+s +++   
  gi|2224615  1172    SLEDLLSVDPEAYQSSVWL    1190 

WD40: domain 1 of 1, from 1217 to 1244: score 2.5, E = 85
                   *->hsssvtslafdpdggllatgsaDgtvriwd<-*
                      h +svt++ + ++  +++++ a+g + +++   
  gi|2224615  1217    HAASVTCILYLNN--QVFVSLANGELVVYQ    1244 

DX: domain 1 of 1, from 1329 to 1335: score 0.5, E = 49
                   *->aeCyedP<-*
                      ++C+++P   
  gi|2224615  1329    HVCLYHP    1335 

//

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:            tem4
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query:  gi|2224615|dbj|BAA20795.1|  KIAA0337 [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 Oct 11 20:41:38 2000

Sequence file: tem4

----------------------------------------
Sequence gi|2224615|dbj|BAA20795.1| (1510 residues):

Matching pattern PS00004 CAMP_PHOSPHO_SITE:
  271: KKKS
  366: RRGS
  768: KRKS
  775: RRSS
 1161: RRSS
 1209: RRNS
Total matches: 6

Matching pattern PS00005 PKC_PHOSPHO_SITE:
  211: SPK
  279: SRR
  332: SER
  361: SPR
  369: SKK
  400: SSR
  531: SLR
  639: SMR
  766: TLK
  773: SLR
  815: TNR
  867: SEK
  922: SSK
  979: SLR
  994: SAR
 1017: SLR
 1118: SSK
 1205: SIR
 1212: SMK
 1306: SSR
 1323: TLK
Total matches: 21

Matching pattern PS00006 CK2_PHOSPHO_SITE:
   14: SSSE
   79: SQEE
   84: SGPE
   89: SLTD
  170: SNGE
  180: SLSD
  193: TSEE
  220: SAMD
  274: SLSD
  307: SSSE
  377: SHQE
  449: SLED
  520: TLLD
  527: SYVE
  555: SLVD
  639: SMRE
  674: TPED
  704: SAEE
  786: SVID
  815: TNRE
  825: SRLD
  850: SLDD
  881: TKLE
  888: TRPE
  925: SCLD
  955: SCPE
 1020: SPPE
 1024: TAPE
 1045: TLAE
 1062: SGLE
 1084: SDSD
 1086: SDDE
 1118: SSKE
 1127: TSSE
 1128: SSEE
 1129: SEEE
 1153: SPMD
 1172: SLED
 1205: SIRD
 1445: TGPE
 1475: SGGD
 1491: SSSE
Total matches: 42

Matching pattern PS00008 MYRISTYL:
  108: GLGTTG
  131: GTEDSL
  163: GSGGSE
  166: GSEFSN
  215: GLPATS
  225: GLTSGH
  236: GSEESK
  324: GATRSR
  701: GVRSAE
  810: GASQAT
  943: GMQFSC
 1001: GAVPGL
 1069: GLGEGD
 1096: GTLQSQ
 1147: GGPCGT
 1191: GTEDGC
 1265: GSPITK
 1303: GQDSSR
 1319: GVWVTL
 1384: GTSAGV
 1388: GVVLTM
 1411: GLGQGH
 1413: GQGHTG
 1488: GGGSSS
 1489: GGSSSE
Total matches: 25

Matching pattern PS00009 AMIDATION:
  748: GGKK
Total matches: 1

Total no of hits in this sequence: 95

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

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

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

Start with Profile Search

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

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, 1510 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:            tem4
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query:  gi|2224615|dbj|BAA20795.1|  KIAA0337 [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:            tem4
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query:  gi|2224615|dbj|BAA20795.1|  KIAA0337 [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]
//

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

L. Aravind's signalling DB
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|2224615|dbj|BAA20795.1| KIAA0337 [Homo sapiens]
         (1510 letters)

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


                                                                   Score     E
Sequences producing significant alignments:                        (bits)  Value

RHOGAP RHO-type GTPase GTP hydrolysis activating protein          172  2e-45
RHOGEF RHO-type GTPase GDP exchange factor                        135  3e-34
SNARE Alpha helical domains which are involved in vesicle fu...    25  0.79
HSF Heat Shock Factor Domain (Winged helix type DNA binding ...    25  0.97
POZ  Pox virus Zinc finger domain (Also called BTB domain; a...    24  1.2
ARR Arrestin domain                                                24  1.4
KELCH Kelch repeat- beta propeller like domain                     23  1.9
PAS domain (Per-Arnt-Sim Domain (A ligand binding and dimeri...    23  2.4
SEC7 Sec7 like GDP exchange factor for ARF like GTPases            23  2.8
PH  Pleckstrin homology domain (lipid and protein interactio...    23  3.3
PUM Pumilio repeat RNA binding domain                              23  3.3
S1  S1 RNA binding domain                                          22  5.0
CNMP Cyclic nucleotide binding domain                              22  6.4
HECT  A ubiquitin conjugating enzyme domain                        22  7.5
MIZFIN  MIZ type Cysteine zinc DNA binding domain                  21  9.2

>RHOGAP RHO-type GTPase GTP hydrolysis activating protein 
          Length = 194

 Score =  172 bits (434), Expect = 2e-45
 Identities = 32/196 (16%), Positives = 72/196 (36%), Gaps = 16/196 (8%)

Query: 517 VAMTLLDTEQSYVESLRTLMQGYMQPLKQPENSVLCDPSLVDEIF--DQIPELLEHHEQF 574
           V   L+ TE  +V +L+ + + +   + +       +P +V  +F       +   H   
Sbjct: 1   VMYELIQTELHHVRTLKIMTRLFRTGMLEEL---HLEPGVVQGLFPLRGRASVTSIHASS 57

Query: 575 LEQVR-------HCMQTWHAQQKVGALLVQSFSK---DVLVNIYSAYIDNFLNAKDAVRV 624
                                 ++G LL+  FS    + +   YS +      A    + 
Sbjct: 58  ASYQNADARPCALGSTGNFVIPRLGDLLISQFSGGSAEQMCKTYSEFCSRHSKALKFYKE 117

Query: 625 AKEARPAFLKFLEQSMRE-NKEKQALSDLMIKPVQRIPRYELLVKDLLKHTPEDHPDHPL 683
                  F +F+ +  R    ++  + + ++   QRI +Y LL+  +L+H+     +   
Sbjct: 118 LYARDKRFQQFIRKVTRPAVLKRHGVQECILLVTQRITKYPLLISRILQHSHGIEEERQD 177

Query: 684 LLEAQRNIKQVAERIN 699
           L  A   +K++   ++
Sbjct: 178 LTTALGLVKELLSNVD 193


>RHOGEF RHO-type GTPase GDP exchange factor 
          Length = 176

 Score =  135 bits (339), Expect = 3e-34
 Identities = 45/162 (27%), Positives = 81/162 (49%), Gaps = 15/162 (9%)

Query: 517 VAMTLLDTEQSYVESLRTLMQGYMQPLKQPENSVLCDPSLVDE---IFDQIPELLE-HHE 572
           V   L+ TE+ YV  L T++ GY   +  PE   L  P L ++   +F  + E+ E H++
Sbjct: 1   VLNELIQTERVYVRELYTVLLGYRAEMDNPEMFDLMPPLLRNKKDILFGNMAEIYEFHND 60

Query: 573 QFLEQVRHCMQTWHAQQKVGALLVQSFSKDVLVNIYSAYIDNFLNAKDAVRVAKEARPAF 632
            FL  + +C       ++VG   ++         +Y+ Y  N   ++   R     + + 
Sbjct: 61  IFLSSLENCAHAP---ERVGPCFLERKD---DFQMYAKYCQNKPRSETIWR-----KYSE 109

Query: 633 LKFLEQSMRENKEKQALSDLMIKPVQRIPRYELLVKDLLKHT 674
             F ++  R+ K +  L   ++KPVQRI +Y+LL+K+LLK++
Sbjct: 110 CAFFQECQRKLKHRLRLDSYLLKPVQRITKYQLLLKELLKYS 151


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

 Score = 24.9 bits (54), Expect = 0.79
 Identities = 12/88 (13%), Positives = 20/88 (22%), Gaps = 15/88 (17%)

Query: 789 DTASKYKMLWKLPLEDADIIKGASQATNRENIQKAISRLDEDLTTLGQMSKLSESLGFPH 848
           +   + K  +  P   +  I       N  ++Q  I                  S     
Sbjct: 117 NFIQRTKQRYNNPRSLSTKI-------NLSDMQMEIKL-RPPYQIPMCE---LGSA---- 161

Query: 849 QSLDDALRDLSAAMHRDLSEKQALCYAL 876
             +  A         +  S  Q L  A 
Sbjct: 162 NGVTSAFSVDCKGAGKISSAHQRLEPAT 189


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

 Score = 24.5 bits (53), Expect = 0.97
 Identities = 6/27 (22%), Positives = 10/27 (36%)

Query: 899 EFPHPDARLGFEQAFDEAKRKLASSKS 925
           EF +     G +      KR+    +S
Sbjct: 73  EFANEGFLGGQKHLLKTIKRRRNVGQS 99


>POZ  Pox virus Zinc finger domain (Also called BTB domain; a protein-protein interaction domain) 
          Length = 229

 Score = 24.4 bits (52), Expect = 1.2
 Identities = 18/150 (12%), Positives = 47/150 (31%), Gaps = 11/150 (7%)

Query: 507 SKPQVDMRKHVAMTLLDTEQSYVESLRTLMQGYMQPLKQPENSVLCDPSLVDEI-FDQIP 565
           S     +       L+D    +  ++   ++     L Q    +     L+ E  +  I 
Sbjct: 58  SGRMEVLTDKEGWILIDRCGKHFGTILNYLRDDTITLPQNRQEI---KELMAEAKYYLIQ 114

Query: 566 ELLEHHEQFLEQVRHCMQT------WHAQQKVGALLVQSFSKDVLVNIYSAYIDNFLNAK 619
            L+   +  L+  +   Q         + ++   L+  S  K V+  +Y+   + +    
Sbjct: 115 GLVNMCQSALQDKKDSYQPVCNIPIITSLKEEERLIESST-KPVVKLLYNRSNNKYSYTS 173

Query: 620 DAVRVAKEARPAFLKFLEQSMRENKEKQAL 649
           ++     +    F K   +        + +
Sbjct: 174 NSDDHLLKNIELFDKLSLRFNGRVLFIKDV 203


>ARR Arrestin domain 
          Length = 454

 Score = 24.1 bits (52), Expect = 1.4
 Identities = 7/30 (23%), Positives = 15/30 (49%)

Query: 875 ALSFPPTKLELCATRPEGTDSYIFEFPHPD 904
           ++S P ++L+    R  G +++ F F    
Sbjct: 117 SISRPLSRLQERLKRKLGANAFPFWFEVAP 146


 Score = 21.4 bits (45), Expect = 8.7
 Identities = 6/30 (20%), Positives = 9/30 (30%), Gaps = 1/30 (3%)

Query: 997  ILCIGAVPGLQPRCHREPPPSLRSPPETAP 1026
             L     P   P       PS+ +   + P
Sbjct: 376  TLTHSK-PPESPERTDRGLPSIEATNGSEP 404


>KELCH Kelch repeat- beta propeller like domain 
          Length = 319

 Score = 23.4 bits (50), Expect = 1.9
 Identities = 13/121 (10%), Positives = 36/121 (29%), Gaps = 14/121 (11%)

Query: 1188 VWLGTEDGCVHVYQSSDSIRDRRNSMKLQHAASVTCILYLNNQVFVSLANGELVVYQREA 1247
                 +   V  Y   ++      +       +   +  L NQ+F+   + E   +   +
Sbjct: 183  SVSHDKLPKVQCYDQCENRWTVPATCPQPWRYTAAAV--LGNQIFIMGGDTE---FSACS 237

Query: 1248 GHFWDPQNFKSVTLGTQGSPIT--KMVSVGGRL--WCGCQ-----NRVLVLSPDTLQLEH 1298
             + ++ + ++   +G   +       V+ G +L    G         +    P       
Sbjct: 238  AYKFNSETYQWTKVGDVTAKRMSCHAVASGNKLYVVGGYFGIQRCKTLDCYDPTLDVWNS 297

Query: 1299 M 1299
            +
Sbjct: 298  I 298


>PAS domain (Per-Arnt-Sim Domain (A ligand binding and dimerization domain)) 
          Length = 117

 Score = 23.4 bits (50), Expect = 2.4
 Identities = 7/88 (7%), Positives = 22/88 (24%), Gaps = 9/88 (10%)

Query: 695 AERINKGVRSAEEAERHARVLQEIEAHIEGMEDLQAPLRRFLR--QEMVIEVKAI----G 748
            E   +        +       ++   ++  E     ++   +      +   A+     
Sbjct: 30  QELQGQPHNMVRHPDMPKAAFADMWFTLKKGEPWSGIVKNRRKNGDHYWVRANAVPMVRE 89

Query: 749 GKKDRSLFLFTDLIVCTTLKRKSGSLRR 776
           GK    + + T     T  +  +     
Sbjct: 90  GKISGYMSIRTRA---TDEEIAAVEPLY 114


>SEC7 Sec7 like GDP exchange factor for ARF like GTPases 
          Length = 192

 Score = 22.9 bits (49), Expect = 2.8
 Identities = 22/108 (20%), Positives = 43/108 (39%), Gaps = 14/108 (12%)

Query: 655 KPVQRIPRYELLVKDLLKHTPEDHPDHPLLLEAQRNIKQVAERINKGVRSAEEAERHARV 714
           KP + +    L+ +  +  TP     H LL     + + + E +        + + +  V
Sbjct: 18  KPEKGVQY--LIERGFVPDTPVGVA-HFLLQRKGLSRQMIGEFL-----GNRQKQFNRDV 69

Query: 715 LQEIEAHIE--GMEDLQAPLRRFLRQEMVIEVKAIGGKKDRSLFLFTD 760
           L  +   ++   ME  +A LR+F      I V+    K +R +  F+ 
Sbjct: 70  LDCVVDEMDFSTMELDEA-LRKFQ---AHIRVQGEAQKVERLIEAFSQ 113


 Score = 21.3 bits (45), Expect = 7.9
 Identities = 5/21 (23%), Positives = 10/21 (46%)

Query: 1186 SSVWLGTEDGCVHVYQSSDSI 1206
            S  +     G V  +++ D+I
Sbjct: 112  SQRYCICNPGVVRQFRNPDTI 132


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

 Score = 22.7 bits (48), Expect = 3.3
 Identities = 10/103 (9%), Positives = 28/103 (26%), Gaps = 4/103 (3%)

Query: 722 IEGMEDLQAPLRRFLRQEMVIEVKAIGGKKDRSLFLFTDLIVCTTLKRKSGSLRRSSMSL 781
            +G +             ++     +   K R   L  + +     +  +    R  + L
Sbjct: 3   DDGNDLTHTFFNPDREGWLLKLGGRVKTWKRRWFILTDNCLYYF--EYTTDKEPRGIIPL 60

Query: 782 YTAASVIDTASKYKMLWKL--PLEDADIIKGASQATNRENIQK 822
              +       +    ++L  P     +IK      +   ++ 
Sbjct: 61  ENLSIREVDDPRKPNCFELYIPNNKGQLIKACKTEADGRVVEG 103


>PUM Pumilio repeat RNA binding domain 
          Length = 337

 Score = 22.9 bits (49), Expect = 3.3
 Identities = 11/34 (32%), Positives = 15/34 (43%)

Query: 797 LWKLPLEDADIIKGASQATNRENIQKAISRLDED 830
           +  LPL D D IK A+       +QK +    E 
Sbjct: 4   ISALPLRDLDYIKLATDQFGCRFLQKKLETPSES 37


>S1  S1 RNA binding domain 
          Length = 305

 Score = 22.2 bits (47), Expect = 5.0
 Identities = 8/50 (16%), Positives = 21/50 (42%), Gaps = 3/50 (6%)

Query: 703 RSAEEAERHARVLQEIEAHIEGMEDLQAPLRRFLRQEMVIEVKAIGGKKD 752
            +A   +  A+ L + + H E  + +   + +   + + ++    GGK  
Sbjct: 7   SAAFSLDEFAKALDKHDYHAEKGQTVHGKICQHANEGVYVD---FGGKSP 53


 Score = 21.4 bits (45), Expect = 9.5
 Identities = 10/50 (20%), Positives = 17/50 (34%), Gaps = 5/50 (10%)

Query: 660 IPRYELLVKDLLKHTPEDHPDHPLLLEAQRNIKQVAERINKGVRSAEEAE 709
           IPR  L+ KD +             +       Q   ++    R  ++AE
Sbjct: 142 IPRSHLMHKDNMDALVGQ--VLKAHILE---ANQDNNKLVLTQRRIQQAE 186


>CNMP Cyclic nucleotide binding domain 
          Length = 118

 Score = 21.7 bits (46), Expect = 6.4
 Identities = 10/59 (16%), Positives = 16/59 (26%), Gaps = 2/59 (3%)

Query: 1126 TTSSEEEQEPGFLPLSGSFGPGG--PCGTSPMDGRALRRSSHGSFTRGSLEDLLSVDPE 1182
             +    EQ    L      G  G       P   R  R S          ++++  D +
Sbjct: 50   ASKGVSEQVLAKLGPGSMAGEIGVVFNIPQPFTVRTRRLSQVIRIGHHKFKEMVQSDND 108


>HECT  A ubiquitin conjugating enzyme domain 
          Length = 255

 Score = 21.7 bits (46), Expect = 7.5
 Identities = 7/19 (36%), Positives = 10/19 (51%)

Query: 1172 SLEDLLSVDPEAYQSSVWL 1190
            +LEDL  VD    Q+   +
Sbjct: 116  TLEDLEEVDLLYVQTLNSI 134


>MIZFIN  MIZ type Cysteine zinc DNA binding domain 
          Length = 172

 Score = 21.1 bits (44), Expect = 9.2
 Identities = 14/60 (23%), Positives = 20/60 (33%), Gaps = 5/60 (8%)

Query: 1081 PFDSDSDDESSPSPSGTLQSQASRSTISSSFGNEETPSSKEA-----TAETTSSEEEQEP 1135
            P     +     S   T    +S  +   S       S K+      T E++S EEE  P
Sbjct: 91   PMRPKKEAMKVTSQPCTKVESSSVFSKPCSVTVASDASKKKIDVIDLTIESSSDEEEDPP 150


Underlying Matrix: BLOSUM62
Number of sequences tested against query: 105
Number of sequences better than 10.0: 15 
Number of calls to ALIGN: 18 
Length of query: 1510 
Total length of test sequences: 20182  
Effective length of test sequences: 16136.0
Effective search space size: 23740223.4
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|2224615|dbj|BAA20795.1| KIAA0337 [Homo sapiens]
         (1510 letters)

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


                                                                   Score     E
Sequences producing significant alignments:                        (bits)  Value

gi|2746866 [80..273] 7-bladed beta-propeller                       51  1e-07
gi|1556399 [647..916] 7-bladed beta-propeller                      48  1e-06
gi|2832331 [157..459] 7-bladed beta-propeller                      47  2e-06
gi|2137729 [37..266] 7-bladed beta-propeller                       46  5e-06
gi|2622601 [25..254] 7-bladed beta-propeller                       43  3e-05
gi|2102675 [28..258] 7-bladed beta-propeller                       43  4e-05
gi|1001440 [1028..1666] 7-bladed beta-propeller                    43  4e-05
gi|549752 [8..246] 7-bladed beta-propeller                         42  8e-05
gi|2815132 [909..1205] 7-bladed beta-propeller                     42  8e-05
gi|1181604 [175..454] 7-bladed beta-propeller                      40  2e-04
gi|1001939 [228..429] 7-bladed beta-propeller                      40  3e-04
gi|1086900 [26..235] 7-bladed beta-propeller                       37  0.002
gi|586529 [222..473] 7-bladed beta-propeller                       37  0.002
gi|1363708 [373..683] 7-bladed beta-propeller                      36  0.004
gi|2104446 [1048..1249] 7-bladed beta-propeller                    36  0.004
gi|452523 [41..225] 7-bladed beta-propeller                        35  0.007
gi|1495018 [84..244] 7-bladed beta-propeller                       35  0.008
gi|2131405 [8..233] 7-bladed beta-propeller                        35  0.011
gi|2736365 [143..399] 7-bladed beta-propeller                      35  0.012
gi|2494897 [186..459] 7-bladed beta-propeller                      32  0.059
gi|2662028 [241..421] 7-bladed beta-propeller                      32  0.059
gi|2129405 [8..201] S-adenosyl-L-methionine-dependent methyl...    32  0.075
gi|1066956 [784..993] 7-bladed beta-propeller                      31  0.15
gi|549660 [99..278] 7-bladed beta-propeller                        29  0.68
gi|4006 [276..439] P-loop containing nucleotide triphosphate...    28  0.90
gi|135122 [4..325] Class II aaRS and biotin synthetases            28  1.0
gi|123030 [1..120] Four-helical up-and-down bundle                 28  1.2
gi|121663 [123..314] ATP-grasp                                     28  1.3
gi|433311 [68..521] DNA/RNA polymerases                            28  1.5
gi|417569 [3..336] Aspartate/ornithine carbamoyltransferase,...    27  2.1
gi|1360139 [98..407] 7-bladed beta-propeller                       27  2.6
gi|1568595 [1..151] Ribosomal protein L9                           26  3.1
gi|2117634 [147..843] Lipoxigenase, C-terminal domain              26  4.0
gi|125988 [12..184] EF Hand-like                                   26  4.9
gi|2811252 [82..382] Periplasmic binding protein-like II           25  6.2
gi|153548 [125..455] beta-Lactamase/D-ala carboxypeptidase         26  6.4
gi|2497466 [46..370] Purine nucleoside hydrolase                   25  7.8
gi|2463100 [37..265] Trypsin-like serine proteases                 25  7.9
gi|1938424 [576..818] Acyl-CoA dehydrogenase (flavoprotein),...    25  8.2
gi|478310 [306..563] beta/alpha (TIM)-barrel                       25  8.6
gi|1420859 [36..611] Periplasmic binding protein-like II           25  9.0
gi|68234 [34..414] Citrate synthase                                25  9.4
gi|544020 [86..429] Nitrogenase iron-molybdenum protein, alp...    25  9.7

>gi|2746866 [80..273] 7-bladed beta-propeller 
          Length = 194

 Score = 51.2 bits (121), Expect = 1e-07
 Identities = 28/190 (14%), Positives = 56/190 (28%), Gaps = 25/190 (13%)

Query: 1191 GTEDGCVHVYQSSDSIRDRRNSMKLQHAASVTCILYLNNQVFVSLANGELVVYQREAGHF 1250
               +G +  +  +    D  N +   H   VT ++ +N   F+S +    V         
Sbjct: 10   AGRNGTLTHWMFNHEPGDLPNVVFDGHKFGVTGVVSINENQFMSGSRDCSV-------KI 62

Query: 1251 WDPQNFKSV-TLGTQGSPITKMVSVGGRLWCGCQNR-VLVLSPDTLQLEHMFYVGQDSSR 1308
            WD ++ K V T     + +T M   G  +    +++ V +  P    +   F   +    
Sbjct: 63   WDLESAKCVRTQTINRNLVTHMAYNGNLVAQTSEDKSVRLWDPRNASVVTEFPRKRHIQM 122

Query: 1309 CVA--------------CMVDSSLGVWVTLKGS-AHVCLYHPDTFEQLAEVDVTPPVHRM 1353
                                   +  +       +     H      L+ V +      +
Sbjct: 123  FCEFDGANRFFSCSNGFNNDGCEITSYDVRNPRQSREARGHEGNVTSLSIVQMDQTKRLL 182

Query: 1354 LAGS-DAIIR 1362
             + S D  IR
Sbjct: 183  ASVSADRQIR 192


>gi|1556399 [647..916] 7-bladed beta-propeller 
          Length = 270

 Score = 48.1 bits (113), Expect = 1e-06
 Identities = 21/163 (12%), Positives = 48/163 (28%), Gaps = 10/163 (6%)

Query: 1188 VWLGTEDGCVHVYQSSDSIRDRRNSMKLQHAASVTCILYLNNQVFVSLANGELVVYQREA 1247
            V   + D  V+ Y  +   R         H  +V+ I + +N+++ +  +  + V+    
Sbjct: 41   VITSSWDNNVYFYSIAFGRRQDTLMG---HDDAVSKICWHDNRLYSASWDSTVKVWSGVP 97

Query: 1248 GHFWDPQNFKSVTLGTQG--SPITKM-VSVGGRL--WCGCQNRVLVLSPDTLQLEHMFYV 1302
                  +      L        +  + ++    L      +  V +    T  L H    
Sbjct: 98   AEMPGTKRHHFDLLAELEHDVSVDTISLNAASTLLVSGTKEGTVNIWDLTTATLMHQIPC 157

Query: 1303 GQDSSRCVACMVDSSLGVWVTLKGSAHVCLYHPDTFEQLAEVD 1345
                    A   DS   +     G   + +    T   ++ + 
Sbjct: 158  HSGIVCDTAFSPDSRHVLSTGTDGC--LNVIDVQTGMLISSMT 198


>gi|2832331 [157..459] 7-bladed beta-propeller 
          Length = 303

 Score = 47.4 bits (111), Expect = 2e-06
 Identities = 32/254 (12%), Positives = 77/254 (29%), Gaps = 34/254 (13%)

Query: 1188 VWLGTEDGCVHVYQSSDSIRDRRNSMKLQHAASVTCILYL--NNQVFVSLANGELVVYQR 1245
            +  G     + ++++         +    H  +V+ + +    +Q++ +  +  + V   
Sbjct: 58   LASGDRSKLILIWEAQSCQHLYTFTG---HRDAVSGLAFRRGTHQLYSTSHDRSVKV--- 111

Query: 1246 EAGHFWDPQNFKSV-TLGTQGSPIT--KMVSVGGRLWCGCQNRVLVLSPDTLQLEHMFYV 1302
                 W+      V TL      +     +S    +  G ++  + +     + + +FY 
Sbjct: 112  -----WNVAENSYVETLFGHQDAVAALDALSRECCVTAGGRDGTVRVWKIPEESQLVFYG 166

Query: 1303 GQDSSRCVACMVDSSLGVWVTLKGSAHVCLYHPDTFEQLAEVDVTPPVHRMLAGSDAIIR 1362
             Q S  C+  + +  +   V+      V L+       LA       +            
Sbjct: 167  HQGSIDCIHLINEEHM---VSGADDGSVALWGLSKKRPLALQREAHGLRGE------PGL 217

Query: 1363 QHKAACLRITALLVCEELLWVGTSAGVVLTMPTSPGTVSCPRAPLSPTGLGQGHTGHVRF 1422
            +       + A L+  +L+  G+ +  V       G     +  L          G +  
Sbjct: 218  EQPFWISSV-AALLNTDLVATGSHSSCVRLWQCGEGF---RQLDLLCD---IPLVGFINS 270

Query: 1423 LAAVQLPDGFNLLC 1436
            L       G  L+ 
Sbjct: 271  LKF--SSSGDFLVA 282


>gi|2137729 [37..266] 7-bladed beta-propeller 
          Length = 230

 Score = 45.7 bits (107), Expect = 5e-06
 Identities = 21/195 (10%), Positives = 50/195 (24%), Gaps = 31/195 (15%)

Query: 1188 VWLGTEDGCVHVYQSSDSIRDRRNSMKLQHAASVTCILYL--NNQVFVSLANGELVVYQR 1245
            +   +  G   V  S          ++L      TC+     +  +     N  + +   
Sbjct: 13   LVTASSSGDKIVVSSCKC--KPVPLLELAEGQKQTCVDLNSTSMYLASGGLNNTVNI--- 67

Query: 1246 EAGHFWDPQNFKSV-TLGTQGSPITKM--VSVGGRLWCGCQNRVL-------------VL 1289
                 WD ++ +   +L      +T +        +  G  +  +               
Sbjct: 68   -----WDLKSKRLHRSLKDHKCEVTCVAYNWNDCYIASGSLSGEIILHSVTTNTSSTPFG 122

Query: 1290 SPDTLQLEHMFYVGQDSSRCVACMVDSSLGVWVTLKGSAHVCLYHPDTFEQLAEVDVTPP 1349
                  + H+ Y     S   +   +  + +W     S++              +  +P 
Sbjct: 123  HGSKQPIRHIKYSLFRKSLLGSVSDNGVVTLWDVNSQSSYHTFDSTHKAPASG-ICFSPV 181

Query: 1350 VHRMLA--GSDAIIR 1362
               +    G D  I 
Sbjct: 182  NELLFVTIGLDKRII 196


>gi|2622601 [25..254] 7-bladed beta-propeller 
          Length = 230

 Score = 43.3 bits (101), Expect = 3e-05
 Identities = 16/157 (10%), Positives = 42/157 (26%), Gaps = 19/157 (12%)

Query: 1188 VWLGTEDGCVHVYQSSDSIRDRRNSMKLQ-HAASVTCILYLNNQVFVSLANGELVVYQRE 1246
            V++G+ DG ++               K +      +      + V++   +  +      
Sbjct: 87   VFVGSWDGYLYAI----DTDTGDLEWKFKTGNRIESSPAVSGDTVYIGSDDCRVYA---- 138

Query: 1247 AGHFWDPQNFKSV-TLGTQGSPITKMVSVGGRLWCGCQ-NRVLVLSPDTLQLEHMFYVGQ 1304
                 D  +        T  +  +  + V G L+ G    +V  LS         +    
Sbjct: 139  ----VDRDDGSKKWEFYTGDAVKSSPLPVNGTLYVGSFNGKVYALSESDGNEVWSYT--S 192

Query: 1305 DSSRCVACMVDSSLGVWVTLKGSAHVCLYHPDTFEQL 1341
              +   +    +      +  G+  +          +
Sbjct: 193  GDAIRSSPAFWNGTVYVGSDDGN--IYALSESDGNVI 227


>gi|2102675 [28..258] 7-bladed beta-propeller 
          Length = 231

 Score = 43.0 bits (100), Expect = 4e-05
 Identities = 16/165 (9%), Positives = 42/165 (24%), Gaps = 20/165 (12%)

Query: 1191 GTEDGCVHVYQSSDSIRDRRNSMKLQ-HAASVTCILYLNNQVFVSLANGELVVYQREAGH 1249
               DG + ++       +  +    +  +       +      V  A     +       
Sbjct: 4    CGSDGDIRIW-------ESLDDDDPKSISIGEKAYSFALKNGKVVTAASNNAIQ------ 50

Query: 1250 FWDPQNFKSV-TLGTQGSPITKMV--SVGGRLWCGCQNRVL-VLSPDTLQLEHMFYVGQD 1305
                 + +    L    +    +V  + G R+  G  + ++ VL  +    +        
Sbjct: 51   LHTFPDGEPDGILTRFTTNANHVVFNTDGTRIAAGSGDFLVKVLQVEDSTQQKTLRGHSA 110

Query: 1306 SSRCVACMVDSSLGVWVTLKGSAHVCLYHPDTFEQLAEVDVTPPV 1350
                V+           +  GS  V ++        A + +    
Sbjct: 111  PVLSVSFDPKDIYLASASCDGS--VRIWKISDQTCEAVLPLLEKC 153


>gi|1001440 [1028..1666] 7-bladed beta-propeller 
          Length = 639

 Score = 42.5 bits (98), Expect = 4e-05
 Identities = 26/205 (12%), Positives = 46/205 (21%), Gaps = 18/205 (8%)

Query: 1143 SFGP-GGPCGTSPMDGRALRRSSHG---SFTRGSLEDLLSVDPEAYQSSVWLGTEDGCVH 1198
             F   G    T+  D  A      G      RG    +      A    +   + D    
Sbjct: 445  EFSADGQFIVTASDDNTAGIWEIVGREVGICRGHEGPVYFAQFSADSRYILTASVDNTAR 504

Query: 1199 VYQSSDSIRDRRNSMKLQHAASVTCILYLNNQVFVSLANGELVVYQREAGHFWDPQNFKS 1258
            ++        R       H + V    +      ++  + +           WD      
Sbjct: 505  IW----DFLGRPLLTLAGHQSIVYQARFSPEGNLIATVSADHTA------RLWDRSGKTV 554

Query: 1259 VTLGTQGSPITKMV--SVGGRLWCGCQNRVLVLSPDTLQLEHMFYVGQDSSRCVACMVDS 1316
              L      +  +     G  L     +    L   + +         +  R      D 
Sbjct: 555  AVLYGHQGLVGTVDWSPDGQMLVTASNDGTARLWDLSGRELLTLEGHGNWVRSAEFSPDG 614

Query: 1317 SLGVWVTLKGSAHVCLYHPDTFEQL 1341
                 +T        L+   T  QL
Sbjct: 615  RW--VLTSSADGTAKLWPVKTLPQL 637


 Score = 41.7 bits (96), Expect = 8e-05
 Identities = 20/234 (8%), Positives = 55/234 (22%), Gaps = 12/234 (5%)

Query: 1171 GSLEDLLSVDPEAYQSSVWLGTEDGCVHVYQSSDSIRDRRNSMKLQHAASVTCILYLNNQ 1230
            G  + + + +   +   +   + DG   ++     +  R   +   H + V    +  + 
Sbjct: 108  GHADWVRNAEFNCHGQILLTASRDGTARLW----DLEGREIGLCQGHTSWVRNAQFSPDG 163

Query: 1231 VFVSLANGELVVYQREAGHFWDPQNFKSVTLGTQGSPITKMV--SVGGRLWCGCQNRVLV 1288
             ++   + +           WD  +     L    + +   +    G  +     +    
Sbjct: 164  QWIVTCSADGTA------RLWDLSSQCFAVLKGHQNWVNNALWSPDGQHIITSSSDGTAR 217

Query: 1289 LSPDTLQLEHMFYVGQDSSRCVACMVDSSLGVWVTLKGSAHVCLYHPDTFEQLAEVDVTP 1348
            +     +          +       +D    V  +   +A +          L       
Sbjct: 218  VWSRHGKCLGTLRGHDHNIHGARFSLDGQKIVTYSTDNTARLWTKEGTLLTILRGHQKEV 277

Query: 1349 PVHRMLAGSDAIIRQHKAACLRITALLVCEELLWVGTSAGVVLTMPTSPGTVSC 1402
                  A    +         R   +   + +   G S  V        G    
Sbjct: 278  YDADFSADGRFVFTVSADQTARQWDISQKDTITLTGHSHWVRNAHFNPKGDRLL 331


 Score = 29.3 bits (64), Expect = 0.43
 Identities = 11/103 (10%), Positives = 31/103 (29%), Gaps = 16/103 (15%)

Query: 1191 GTEDGCVHVYQSSDSIRDRRNSMKLQHAASVTCILYLNN--QVFVSLANGELVVYQREAG 1248
             T +G  +++    S      ++   H  ++T + +  +      + A+  + +      
Sbjct: 5    ATTNGIAYLW----SFHGELINVLRGHQEAITALDWSADGQYFATASADHTVKL------ 54

Query: 1249 HFWDPQNFKSVTLGTQGSPITKMV--SVGGRLWCGCQNRVLVL 1289
              W     +  TL      +  +        L    Q+    +
Sbjct: 55   --WQRHGEEVATLRGHEDWVRSVHFSPHHQFLVTSGQDNTARI 95


>gi|549752 [8..246] 7-bladed beta-propeller 
          Length = 239

 Score = 41.9 bits (97), Expect = 8e-05
 Identities = 18/188 (9%), Positives = 49/188 (25%), Gaps = 25/188 (13%)

Query: 1188 VWLGTEDGCVHVYQSSDSIRDRRNSMKLQHAASVTCILYLNNQVFVSLANGELVVYQREA 1247
            +  G +D  ++      +  +      + H  +V  + + +  V     +    V     
Sbjct: 62   LLFGGKDTMINGVPLFATSGEDPLYTLIGHQGNVCSLSFQDGVVISGSWDKTAKV----- 116

Query: 1248 GHFWDPQNFKSVTLGTQGSPITKMV---SVGGRLWCGCQNRVLVLSPDTLQLEHMFYVGQ 1304
               W   +            +           +      ++ + L     ++   F    
Sbjct: 117  ---WKEGSLVYNLQAHNA-SVWDAKVVSFSENKFLTASADKTIKLW-QNDKVIKTFSGIH 171

Query: 1305 DSSRCVACMVDSSLGVWVTLKGSAHVCLYHPDTFEQLAE----------VDVTPPVHRML 1354
            +    V  +     G +++      + L    T + L            + + P    + 
Sbjct: 172  NDV--VRHLAVVDDGHFISCSNDGLIKLVDMHTGDVLRTYEGHESFVYCIKLLPNGDIVS 229

Query: 1355 AGSDAIIR 1362
             G D  +R
Sbjct: 230  CGEDRTVR 237


>gi|2815132 [909..1205] 7-bladed beta-propeller 
          Length = 297

 Score = 41.6 bits (96), Expect = 8e-05
 Identities = 23/190 (12%), Positives = 52/190 (27%), Gaps = 18/190 (9%)

Query: 1152 TSPMDGRA----LRRSSHGSFTRGSLEDLLSVDPEAYQSSVWLGTEDGCVHVYQSSDSIR 1207
            T   D       +      S   G +  + ++   +  + V +  ED  + V+    S  
Sbjct: 117  TGSGDFVVQMWDVTNGKCISRMGGLMAPVSTLAITSNDAFVVVACEDETLKVF----STV 172

Query: 1208 DRRNSMKLQ-HAASVTCILYL-NNQVFVSLANGELVVYQREAGHFWDPQNFKSVTLGTQG 1265
              +   +L  H   V  ++   ++    +    ++  Y    G          V    Q 
Sbjct: 173  GGQELHELMGHEGKVNSLVCAQDDCQLFAATKSKVFCYDIHNGQ------MIDVLDTAQP 226

Query: 1266 SPITKMVSVGGRLW--CGCQNRVLVLSPDTLQLEHMFYVGQDSSRCVACMVDSSLGVWVT 1323
             PI  +       +    C  +V + +      +             A  + +       
Sbjct: 227  FPICSLKISSDNYFLISPCGPKVTIWNVTKRNHDAHDVHADKEGFLTAVALSNDDKYAAC 286

Query: 1324 LKGSAHVCLY 1333
               +  V L+
Sbjct: 287  GTNNGIVALW 296


 Score = 28.4 bits (62), Expect = 0.99
 Identities = 14/122 (11%), Positives = 30/122 (24%), Gaps = 16/122 (13%)

Query: 1251 WDPQNFKSV-TLGTQGSPITKM-VSVGGRL--WCGCQNRVLVLSPDTLQLEHMFYVGQDS 1306
            +   + + + T     + +T +  S    L         V V                  
Sbjct: 1    YHIASEQLIGTFKGHTAAVTCLCSSNDSSLFVSTSFDKTVNVWVFSQSTPTMSLTHHTAK 60

Query: 1307 SRCVACMVDSSLGVWVTLKGSAHVCLYHPDTFEQLAE----------VDVTPPVHRMLAG 1356
              C     D    +  +   S    +   +T E +            + +T     ++ G
Sbjct: 61   VTCAILTSDDQYLITASADSS--AKMIKLETGEVMRSFNDHTGSVVSLQLTSNNQFLITG 118

Query: 1357 SD 1358
            S 
Sbjct: 119  SG 120


 Score = 26.0 bits (56), Expect = 3.8
 Identities = 13/63 (20%), Positives = 23/63 (35%), Gaps = 5/63 (7%)

Query: 1188 VWLGTEDGCVHVYQSSDSIRDRRNSMKLQHAASVTCILYL--NNQVFVSLANGELVVYQR 1245
                + D  V+V+  S S           H A VTC +    +  +  + A+    + + 
Sbjct: 31   FVSTSFDKTVNVWVFSQSTPTMS---LTHHTAKVTCAILTSDDQYLITASADSSAKMIKL 87

Query: 1246 EAG 1248
            E G
Sbjct: 88   ETG 90


>gi|1181604 [175..454] 7-bladed beta-propeller 
          Length = 280

 Score = 40.4 bits (93), Expect = 2e-04
 Identities = 17/197 (8%), Positives = 51/197 (25%), Gaps = 27/197 (13%)

Query: 1178 SVDPEAYQSSVWLGTEDGCVHVYQSSDSIRDRRNSMKLQHAASVTCILYL--NNQVFVSL 1235
            S+ P A+         +   +V    D +    N     ++  V  I +     ++    
Sbjct: 1    SLTPIAHIVKFGSAATESHANVTDIHDGLDFSSNDDG-GYSFGVFSIKFSTDGREIVAGT 59

Query: 1236 ANGELVVYQREAGHFWDPQNFKSV-TLGTQGSPITKMV---SVGGRLWCGCQNRVLVL-- 1289
            ++  + V        +D +  +    +    S +  +      G  ++ G  + +  +  
Sbjct: 60   SDESICV--------YDLEADRLSLRISAHESDVNSVCFADESGHLIYSGSDDNLCKVWD 111

Query: 1290 ---SPDTLQLEHMFYVGQDSSRCVACMVDSSLGVWVTLKGSAHVCLYHPDTFEQLAE--- 1343
                    +   +     +    +    D    +      +  + L+        A    
Sbjct: 112  RRCFNAKGKPAGILMGHLEGITFIDSRGDGRYFISNGKDQT--IKLWDIRKMSSNAGGTI 169

Query: 1344 --VDVTPPVHRMLAGSD 1358
               +       M    +
Sbjct: 170  QSRNSEWDYRWMEYPQE 186


>gi|1001939 [228..429] 7-bladed beta-propeller 
          Length = 202

 Score = 39.9 bits (92), Expect = 3e-04
 Identities = 19/186 (10%), Positives = 48/186 (25%), Gaps = 18/186 (9%)

Query: 1191 GTEDGCVHVY--QSSDSIRDRRNSMKLQHAASVTCILYLNNQVFVSLA-NGELVVYQREA 1247
             + +  +H++   + +     R    L    +   + +  +   +    N  + V+    
Sbjct: 9    SSRENPIHIWDAFTGELRASFRAYNHLDELTAAHSLCFSPDGSQLFCGFNRTVRVFSTSR 68

Query: 1248 GHFWDPQNFKSVTLGTQGSPITKM--VSVGGRLWCGCQ-NRVLVLSPDTLQLEHMFYVGQ 1304
                            Q   I+ +          CG     + + + D      +    Q
Sbjct: 69   PGRDCEVRATFAKKQGQSGIISCIAFSPSQPLYACGSYGRTLGLYAWDDGSPLALLGGHQ 128

Query: 1305 DSSRCVA-----------CMVDSSLGVWVTLKGSAHVCLY-HPDTFEQLAEVDVTPPVHR 1352
                 +               D+ L  W   +    +       T  Q    D+ P    
Sbjct: 129  GGITHLCFHPDGNLFFSGARKDAELLCWDLRQPGHLLWSLSREVTTNQRIYFDLDPSGQF 188

Query: 1353 MLAGSD 1358
            +++G+ 
Sbjct: 189  LVSGNT 194


>gi|1086900 [26..235] 7-bladed beta-propeller 
          Length = 210

 Score = 37.1 bits (85), Expect = 0.002
 Identities = 18/163 (11%), Positives = 48/163 (29%), Gaps = 17/163 (10%)

Query: 1176 LLSVDPEAYQSSVWLGTEDGCVHVYQSSDSIRDRRNSMKLQ-HAASVTCILY--LNNQVF 1232
            + +V      S V+ G+ D  +  +      R R        + ++   +      +++ 
Sbjct: 59   VNAVAFNEESSVVFSGSMDCTMQAF----DCRSRSEKPIQIFNESTDGILSIDVNGHEIV 114

Query: 1233 VSLANGELVVYQREAGHFWDPQNFKSVTLGTQGSPITKMVSVGGRLW--CGCQNRVLVLS 1290
               A+G            +  ++          S  +   +              V ++ 
Sbjct: 115  AGSADGNY--------RVYSIRDGNMTVDYMGDSVNSVSFTPDSNCLLAGVMGGIVRLID 166

Query: 1291 PDTLQLEHMFYVGQDSSRCVACMVDSSLGVWVTLKGSAHVCLY 1333
              + +L   +   Q++   + C V  S+    +      V +Y
Sbjct: 167  KSSGKLLASYKGHQNTEYKLDCRVLQSIEHVASGSEDGFVYVY 209


>gi|586529 [222..473] 7-bladed beta-propeller 
          Length = 252

 Score = 36.9 bits (84), Expect = 0.002
 Identities = 21/146 (14%), Positives = 40/146 (27%), Gaps = 18/146 (12%)

Query: 1191 GTEDGCVHVYQSSDSIRDRRNSMKLQ-HAASVTCILYL--NNQVFVSLANGELVVYQREA 1247
                G + VY     I ++  + KL  H   ++ + +   N  +  +  +G L +     
Sbjct: 114  PGPKGAIFVY----QITEKTPTGKLIGHHGPISVLEFNDTNKLLLSASDDGTLRI----- 164

Query: 1248 GHFWDPQNFKSV-TLGTQGSPITKM-VSVGGRLWCGCQ-NRVLVLSPDTLQLEHMFYVGQ 1304
               W   N  S  +       I         ++        V + S     L  +  V  
Sbjct: 165  ---WHGGNGNSQNSFYGHSQSIVSASWVGDDKVISCSMDGSVRLWSLKQNTLLALSIVDG 221

Query: 1305 DSSRCVACMVDSSLGVWVTLKGSAHV 1330
                      D        + G  +V
Sbjct: 222  VPIFAGRISQDGQKYAVAFMDGQVNV 247


 Score = 28.3 bits (62), Expect = 0.94
 Identities = 17/107 (15%), Positives = 36/107 (32%), Gaps = 13/107 (12%)

Query: 1221 VTCILYL--NNQVFVSLANGELVVYQREAGHFWDPQNFKSVTLGTQGSPITKM-VSVGGR 1277
            VTC+ +    N +   + NGEL +        W+        L    +PI  +  +  G 
Sbjct: 2    VTCLAWSHDGNSIVTGVENGELRL--------WNKTGALLNVLNFHRAPIVSVKWNKDGT 53

Query: 1278 --LWCGCQNRVLVLSPDTLQLEHMFYVGQDSSRCVACMVDSSLGVWV 1322
              +    +N  ++ +  +  +   F + +     +     S  G   
Sbjct: 54   HIISMDVENVTILWNVISGTVMQHFELKETGGSSINAENHSGDGSLG 100


>gi|1363708 [373..683] 7-bladed beta-propeller 
          Length = 311

 Score = 36.1 bits (82), Expect = 0.004
 Identities = 17/151 (11%), Positives = 35/151 (22%), Gaps = 21/151 (13%)

Query: 1188 VWLGTEDGCVHVYQSSDSIRDRRNSMKLQ----HAASVTCILYLNNQVFVSLANGELVVY 1243
            + +   D  V V+          +         H   V  I    +   +  ++ +  + 
Sbjct: 174  LAISLLDNTVKVF-------FLDSMKFYLSLYGHKLPVLSIDISFDSKMIITSSADKNIK 226

Query: 1244 QREAGHFWDPQNFKSV-TLGTQGSPITKM-VSVGGRLW--CGCQNRVLVLSPDTLQLEHM 1299
                   W         +L      I  +        +  C     V     +  +    
Sbjct: 227  ------IWGLDFGDCHKSLFAHQDSIMNVKFLPQSHNFFSCSKDAVVKYWDGEKFECIQK 280

Query: 1300 FYVGQDSSRCVACMVDSSLGVWVTLKGSAHV 1330
             Y  Q     +A   D    V  +   S  +
Sbjct: 281  LYAHQSEVWALAVATDGGFVVSSSHDHSIRI 311


 Score = 36.1 bits (82), Expect = 0.005
 Identities = 28/209 (13%), Positives = 54/209 (25%), Gaps = 27/209 (12%)

Query: 1217 HAASVTCILYLNNQVFVSLANGELVVYQREAGHFWDPQNFKSVTLGTQGSPITKMV--SV 1274
                V  I   ++   ++ A+   +         W+ +  K +    +            
Sbjct: 18   QRTDVRSIDISDDNKLLATASNGSL-------KIWNIKTHKCIRTF-ECGYALTCKFLPG 69

Query: 1275 GGRLWCGCQNRVLVL-------SPDTLQLEHMFYVG-----QDSSRCVACMVDSSLGVWV 1322
            G  +  G +N  L L         DT++  H   +       D  R V    D ++  W 
Sbjct: 70   GLLVILGTRNGELQLFDLASSSLLDTIEDAHDAAIWSLDLTSDGKRLVTGSADKTVKFWD 129

Query: 1323 TLKGSAHVCLYHPDTFEQLAEVDVTPPVHRMLAGSDAIIRQHKAACLRITALLVCEELLW 1382
                ++ V          L     T             +R          +LL     ++
Sbjct: 130  FKVENSLVPGTKNKFLPVLKLHHDTTLELT---DDILCVRVSPDDRYLAISLLDNTVKVF 186

Query: 1383 VGTSAGVVLTMPTSPGTVSCPRAPLSPTG 1411
               S    L++      V      +S   
Sbjct: 187  FLDSMKFYLSLYGHKLPVLSID--ISFDS 213


>gi|2104446 [1048..1249] 7-bladed beta-propeller 
          Length = 202

 Score = 36.0 bits (82), Expect = 0.004
 Identities = 13/151 (8%), Positives = 34/151 (21%), Gaps = 15/151 (9%)

Query: 1188 VWLGTEDGCVHVYQSSDSIRDRRNSMKLQ--HAASVTCILYL---NNQVFVSLANGELVV 1242
            + +  +   + ++              L    + S+T +       N +    ++G L V
Sbjct: 57   LLVAGDVRVIRIW----DASKEICYANLPVRSSNSITSLTSDLVGCNIIVAGFSDGVLRV 112

Query: 1243 YQREAGHFWDPQNFKSVTLGTQGSPITKMVSVGGRLW---CGCQNRVLVLSPDTLQLEHM 1299
            Y +             V        +   +   G            V +           
Sbjct: 113  YDKRLPARDSL---TDVWKEHSSEIVNVEMQSSGMRELISASSDGEVKLWDIRMNHSLQT 169

Query: 1300 FYVGQDSSRCVACMVDSSLGVWVTLKGSAHV 1330
            F         +     + +    +   S  +
Sbjct: 170  FSTDNSGLTSLTVHSHAPVYATGSSNQSIKI 200


>gi|452523 [41..225] 7-bladed beta-propeller 
          Length = 185

 Score = 35.2 bits (80), Expect = 0.007
 Identities = 19/145 (13%), Positives = 46/145 (31%), Gaps = 11/145 (7%)

Query: 1170 RGSLEDLLSVDPEAYQSSVWLGTEDGCVHVYQSSDSIRDRRNSMKLQHAASVTCILYL-N 1228
            +       +V        + LGT  G + +Y +       +      H   V CI +  +
Sbjct: 48   QRKKRKSEAVGMSNQTDLLALGTAVGSILLYSTVKGELHSKLISG-GHDNRVNCIQWHQD 106

Query: 1229 NQVFVSLANGELVVYQREAGHFWDPQNFKSVTLGT-QGSPITKM-VSVGGRLWCGCQNRV 1286
            +    S ++ + +         W+ Q  K         S ++ + +S  G++       +
Sbjct: 107  SGCLYSCSDDKHI-------VEWNVQTCKVKCKWKGDNSSVSSLCISPDGKMLLSAGRTI 159

Query: 1287 LVLSPDTLQLEHMFYVGQDSSRCVA 1311
             +   +T ++   F         + 
Sbjct: 160  KLWVLETKEVYRHFTGHATPVSSLM 184


>gi|1495018 [84..244] 7-bladed beta-propeller 
          Length = 161

 Score = 35.1 bits (80), Expect = 0.008
 Identities = 20/139 (14%), Positives = 48/139 (34%), Gaps = 15/139 (10%)

Query: 1191 GTEDGCVHVYQSSDSIRDRRNSMKLQHAASVTCILYL--NNQVFVSLANGELVVYQREAG 1248
             + D  +  +  + S  DR++ +   H  SV  I +   +  VFV+      V       
Sbjct: 29   ISGDSTIRCWDLNQSTLDRKSQVFFGHEGSVRSICFAPDDPNVFVTGGRDFQV------- 81

Query: 1249 HFWDPQNFKSVTLGTQG--SPITKMV--SVGGRLW--CGCQNRVLVLSPDTLQLEHMFYV 1302
              WD +      +      + IT         ++      +++    + +  ++  + ++
Sbjct: 82   KIWDMRVSTVKKMEEDCRMATITYKTAHPKPSKVLTSGTPKSKAKAKTIEGYKVTSVLFL 141

Query: 1303 GQDSSRCVACMVDSSLGVW 1321
             +      +   DS + VW
Sbjct: 142  DEHHVASASENADSGIRVW 160


>gi|2131405 [8..233] 7-bladed beta-propeller 
          Length = 226

 Score = 34.8 bits (79), Expect = 0.011
 Identities = 15/202 (7%), Positives = 51/202 (24%), Gaps = 31/202 (15%)

Query: 1188 VWLGTEDGCVHVYQSSDSIRDRRNSMKLQ-----HAASVTCILYL--NNQVFVSLANGEL 1240
            +  G+ D  + +     S++    ++        H  ++  + +    + +     +  +
Sbjct: 20   LATGSTDRKIKLV----SVKYDDFTLIDVLDETAHKKAIRSVAWRPHTSLLAAGSFDSTV 75

Query: 1241 VVYQREAGHFWDPQNFKSVTLGTQGSPITKM--VSVGGRLWCGCQ-NRVLVLSPDTLQLE 1297
             ++ +E       +      +    + +  +   + G  L    +   V +   D    E
Sbjct: 76   SIWAKEESADRTFEMDLLAIIEGHENEVKGVAWSNDGYYLATCSRDKSVWIWETDESGEE 135

Query: 1298 HMFYVG--------------QDSSRCVACMVDSSLGVWVTLKGSAHVCLYHPDTFEQLAE 1343
            +                      +   +   D ++ +W                   +  
Sbjct: 136  YECISVLQEHSQDVKHVIWHPSEALLASSSYDDTVRIWKDYDDDWECVAVLNGHEGTVWS 195

Query: 1344 VDVTPPVHRML---AGSDAIIR 1362
             D               D+ +R
Sbjct: 196  SDFDKTEGVFRLCSGSDDSTVR 217


>gi|2736365 [143..399] 7-bladed beta-propeller 
          Length = 257

 Score = 34.5 bits (78), Expect = 0.012
 Identities = 21/254 (8%), Positives = 50/254 (19%), Gaps = 40/254 (15%)

Query: 1188 VWLGTEDGCVHVYQSSDSIRDRRNSMKLQ-HAASVTCILYL-------NNQVFVSLANGE 1239
            +  G+ D  + +   ++        M+L  H  +V  ++++          V     N  
Sbjct: 13   LATGSNDKTIRLMAFNEDSCRIGAEMELTCHDGTVRDLIFIESSINRSTILVSGGAGNCH 72

Query: 1240 LVVYQREAGH-FWDPQNFKSVTLGTQGSPITKMVSVGG-----RLWC---GCQNRVLVLS 1290
            L V     G      +      +         +          +                
Sbjct: 73   LNVTDCNTGQLIQSMKGHTGSCMIVHLMLKISVTKFHILNCTVKTMIYHKTFSTLFHSFR 132

Query: 1291 PDTLQLEHMFYVGQDSSRCVACMVDSSLGVWVTLKGSAHVCLYHPDTFEQLAEVDVTPPV 1350
             D    +                   +  V  +      +  +       +         
Sbjct: 133  EDEKVQKSGSLGLPPILGLYTWSQAGNQFVSCSQDK--TIRFWDLRQQTAVN-------- 182

Query: 1351 HRMLAGSDAIIRQHKAACLRITALLVCEELLWVGTSAGVVLTMPTSPGTVSCPRAPLSPT 1410
                     I   + +    +  +    +LL  G     V         +          
Sbjct: 183  --------VISPGNASPVTSV-CVDPNGKLLVSGHEDASVALFDIGGNRILQTF----RY 229

Query: 1411 GLGQGHTGHVRFLA 1424
                 H   VR + 
Sbjct: 230  LRSLPHGDEVRTVR 243


>gi|2494897 [186..459] 7-bladed beta-propeller 
          Length = 274

 Score = 32.1 bits (72), Expect = 0.059
 Identities = 27/274 (9%), Positives = 57/274 (19%), Gaps = 59/274 (21%)

Query: 1188 VWLGTEDGCVHVYQ--------------SSDSIRDRRNSMKLQ----HAASVTCILY--L 1227
            + +G     + V+               S  S + ++   K      H  +V  + +  L
Sbjct: 23   IAVGNMTPVIEVWDLDIVDSLEPVFTLGSKLSKKKKKKGKKSSSAEGHTDAVLDLSWNKL 82

Query: 1228 NNQVFVSLANGELVVYQREAGHFWDPQNFKSV-TLGTQGSPITKMV--SVGGRLW--CGC 1282
               V  S +    V         WD    K   +L      +  +       +       
Sbjct: 83   IRNVLASASADNTV-------ILWDMSLGKPAASLAVHTDKVQTLQFHPFEAQTLISGSY 135

Query: 1283 QNRVLVLSPDTLQLEHMFYVGQDSSRCVACMVDSSLGVWVTLKGSAHVCLYHPDTFEQLA 1342
               V +    +    H  +        V     S     +       V      + + + 
Sbjct: 136  DKSVALYDCRSPDESHRMWRFSGQIERVTWNHFSPCHF-LASTDDGFVYNLDARSDKPIF 194

Query: 1343 EVDVTPPVHRMLAGSDAIIRQHKAACLRITALLVCEELLWVGTSAGVVLTMPTSPGTVSC 1402
                              +  H      +      +  L   ++   V            
Sbjct: 195  T-----------------LNAHNDEISGLDLSSQIKGCLVTASADKYVKIWDILGDR--- 234

Query: 1403 PRAPLSPTGLGQGHTGHVRFLAAVQLPDGFNLLC 1436
                 S         G +   +    PD   +  
Sbjct: 235  ----PSLVHSRDMKMGVLFCSSC--CPDLPFIYA 262


>gi|2662028 [241..421] 7-bladed beta-propeller 
          Length = 181

 Score = 32.4 bits (73), Expect = 0.059
 Identities = 12/150 (8%), Positives = 32/150 (21%), Gaps = 25/150 (16%)

Query: 1186 SSVWLGTEDGCVHVYQSSDSIRDRRNSMKLQ-HAASVTCILYLNNQVFVSLANGELVVYQ 1244
            +                   +       +   H A +  + +  N   +S +  +     
Sbjct: 40   TCSSDPNNPKLNLKA----PVFSEAPIREYAGHTADILDLSWSRNNFLLSSSMDKTA--- 92

Query: 1245 REAGHFWDPQNFKSVTLGTQGSPITKMV--SVGGRLW--CGCQNRVLVLSPDTLQL---- 1296
                  W P     +        +T +       R +       ++ + S     +    
Sbjct: 93   ----RLWHPVRKDCLCCFEHSDFVTSIAFHPKDDRFFLSGSLDCKLRLWSIKEKAVSFWN 148

Query: 1297 EHMFYV-----GQDSSRCVACMVDSSLGVW 1321
            E    +       D    +A         +
Sbjct: 149  ELPELITAVAFSPDGGLAIAGTFVGLCLFY 178


>gi|2129405 [8..201] S-adenosyl-L-methionine-dependent methyltransferases 
          Length = 194

 Score = 31.8 bits (72), Expect = 0.075
 Identities = 14/58 (24%), Positives = 20/58 (34%), Gaps = 3/58 (5%)

Query: 661 PRYELL--VKDLLKHTPEDHPDHPLLL-EAQRNIKQVAERINKGVRSAEEAERHARVL 715
           P + LL  VK+     P D   +P              E +   V   E+ E +  VL
Sbjct: 114 PEHVLLENVKEFRSWGPIDEDGNPTKDGSIFERWISTLEALGYSVVYDEKDEDYGVVL 171


>gi|1066956 [784..993] 7-bladed beta-propeller 
          Length = 210

 Score = 31.0 bits (69), Expect = 0.15
 Identities = 15/145 (10%), Positives = 39/145 (26%), Gaps = 18/145 (12%)

Query: 1188 VWLGTEDGCVHVYQSSDSIRDRRNSMKLQ-HAASVTCILYL-NNQVFVSLANGELVVYQR 1245
            +W           + +           L  H AS++ +     + + VS     ++    
Sbjct: 73   LWHWNGTTGFIAGEYNQPGEVPSPRSILTGHEASISALCVSAEHGLVVSGCEDGVI---- 128

Query: 1246 EAGHFWDPQNFKSVTLGTQGSPITKMVSVGGRLWCGCQNRVLVL----SPDTLQLEHMFY 1301
                     +     +   G      +S    L     ++ +V     +    ++     
Sbjct: 129  ---LIHTTSSDLLRRIRGHGIVTQLSMSRECILLSLFDSKRMVTYSATAKKLDEVLVDDK 185

Query: 1302 VG-----QDSSRCVACMVDSSLGVW 1321
            +      +D    V   V+  + +W
Sbjct: 186  IECVTVTRDGEFAVTGAVNGRINIW 210


>gi|549660 [99..278] 7-bladed beta-propeller 
          Length = 180

 Score = 28.6 bits (63), Expect = 0.68
 Identities = 12/107 (11%), Positives = 31/107 (28%), Gaps = 12/107 (11%)

Query: 1191 GTEDGCVHVYQSSDSIRDRRNSMK--LQHAASVTCILYL-NNQVFVSLANGELVVYQREA 1247
              E+  + +Y            +K      +    + +  N +  +  ++          
Sbjct: 74   NPENFEIGLYNLKKIQEGPFLIIKINDATFSQWNKLEFSNNGKYLLVGSSIGKH------ 127

Query: 1248 GHFWDPQNFKSVTL--GTQGSPITKMVSVGGRLWCGCQNRVLVLSPD 1292
               +D    + +    GT+  P+ + +  G   +      VL    D
Sbjct: 128  -LIFDAFTGQQLFELIGTRAFPMREFLDSGSACFTPDGEFVLGTDYD 173


>gi|4006 [276..439] P-loop containing nucleotide triphosphate hydrolases 
          Length = 164

 Score = 28.4 bits (62), Expect = 0.90
 Identities = 6/77 (7%), Positives = 20/77 (25%), Gaps = 5/77 (6%)

Query: 1184 YQSSVWLGTEDGCVHVYQSSDSI---RDRRNSMKLQHAASVTCILYLNNQVFVSLANGEL 1240
             ++       D C+ +   +D      +   +           I+ + N+    L + + 
Sbjct: 71   DRAKKKSVQSDLCLFIVDPTDLSKLLPEDILAHLSSKTFGNKRIIIVVNK--SDLVSDDE 128

Query: 1241 VVYQREAGHFWDPQNFK 1257
            +              + 
Sbjct: 129  MTKVLNKLQTRLGSKYP 145


>gi|135122 [4..325] Class II aaRS and biotin synthetases 
          Length = 322

 Score = 28.2 bits (62), Expect = 1.0
 Identities = 25/218 (11%), Positives = 52/218 (23%), Gaps = 39/218 (17%)

Query: 698 INKGVRSAEEAERHARVLQEIEAHI------EGMEDLQAPLRRFLRQ-----EMVIEVKA 746
           I    R     +   R   ++   +      +   +L     R+ R       + +E+ +
Sbjct: 105 IGPMFRHERPQKGRYRQFHQLGCEVFGLQGPDIDAELIMLTARWWRALGISEHVTLELNS 164

Query: 747 IGGKKDRSLFLFTDLIVCTTLKRKSGSLRRSSMSLYTAASVIDTASKYKMLWKLPLEDAD 806
           IG  + R                    L +    L          +  ++L     E   
Sbjct: 165 IGSLEAR----------ANYRDALVAFLEQHKEKLDEDCKRRMYTNPLRVLDSKNPEVQA 214

Query: 807 IIKGASQATNRENIQKAISRLDEDLTTLGQMSKLSESLGFPHQSLDDALRDLSAAMHRDL 866
           ++  A    +           +E       + KL ES G  +        +       D 
Sbjct: 215 LLNDAPALGDYL--------DEESREHFAGLCKLLESAGIAYTV------NQRLVRGLDY 260

Query: 867 -SEKQALCYALSFPPTKLELCATRPEGTDSYIFEFPHP 903
            +         S           R +G    + +    
Sbjct: 261 YNRTVFEWVTNSLGSQGTVCAGGRYDG---LVEQLGGR 295


>gi|123030 [1..120] Four-helical up-and-down bundle 
          Length = 120

 Score = 27.7 bits (61), Expect = 1.2
 Identities = 14/47 (29%), Positives = 17/47 (35%), Gaps = 10/47 (21%)

Query: 563 QIPELLEH---HEQFLEQVRHCMQTWHAQQKVGALLVQSFSKDVLVN 606
              E   H   HE FL  +R       A      L    ++KD LVN
Sbjct: 68  YKSEHPGHKKKHEDFLAVIRG----LSAPVPNDKL---LYAKDWLVN 107


>gi|121663 [123..314] ATP-grasp 
          Length = 192

 Score = 27.7 bits (61), Expect = 1.3
 Identities = 16/144 (11%), Positives = 37/144 (25%), Gaps = 13/144 (9%)

Query: 617 NAKDAVRVAKEARPAFLKFLEQSMRENKEK-QALSDLMIKPV-------QRIPRYELLVK 668
           N K       +  P  L    ++  + K   +  SD+++KP+           +      
Sbjct: 1   NEKLFTAWFSDLTPETLVT--RNKAQLKAFWEKHSDIILKPLDGMGGASIFRVKEGDPNL 58

Query: 669 DLLKHTPEDHPDHPLLLEAQRNIKQ---VAERINKGVRSAEEAERHARVLQEIEAHIEGM 725
            ++  T  +H     + +      +       +  G        R  +  +       G 
Sbjct: 59  GVIAETLTEHGTRYCMAQNYLPAIKDGDKRVLVVDGEPVPYCLARIPQGGETRGNLAAGG 118

Query: 726 EDLQAPLRRFLRQEMVIEVKAIGG 749
                PL     +        +  
Sbjct: 119 RGEPRPLTESDWKIARQIGPTLKE 142


>gi|433311 [68..521] DNA/RNA polymerases 
          Length = 454

 Score = 27.7 bits (61), Expect = 1.5
 Identities = 13/99 (13%), Positives = 36/99 (36%), Gaps = 6/99 (6%)

Query: 488 KPPEAARPADEPTPASKCCSKP-QVDMRKHVAMTLLDTEQSYVESLRTLMQGYMQPLKQP 546
                 +P    T  +K    P +  +   +    + +   Y  +L  ++Q + +     
Sbjct: 114 PLDPDYQPYTAFTLPNKNNQGPGRRYVWCSLPQGWVLSPLIYQSTLDNILQPWRK----- 168

Query: 547 ENSVLCDPSLVDEIFDQIPELLEHHEQFLEQVRHCMQTW 585
           +   +     +D+I+         HE+ ++++R  +  W
Sbjct: 169 KYPNIDVYQYMDDIYIGSDFSRLEHEKIIQELRDLLIFW 207


>gi|417569 [3..336] Aspartate/ornithine carbamoyltransferase, Catalytic-chain 
          Length = 334

 Score = 27.2 bits (60), Expect = 2.1
 Identities = 12/95 (12%), Positives = 25/95 (25%), Gaps = 9/95 (9%)

Query: 1256 FKSVTLGTQGSPITKMVSVGGRLWCGCQNRVLVLSPDTLQLEHMFYVGQDSSRCVACMVD 1315
            F   +  T+ S                  ++  L+ ++           D+   +A   +
Sbjct: 57   FLEPSTRTKESFRNAAKFHQV--------KLSDLAAESSSFNKGESY-ADTFNTLAGYQN 107

Query: 1316 SSLGVWVTLKGSAHVCLYHPDTFEQLAEVDVTPPV 1350
            S   V   ++G           F Q   +   P  
Sbjct: 108  SIFIVRSEVEGVCRWLEDEAQAFYQRNNLKRKPAF 142


>gi|1360139 [98..407] 7-bladed beta-propeller 
          Length = 310

 Score = 26.9 bits (59), Expect = 2.6
 Identities = 18/173 (10%), Positives = 37/173 (20%), Gaps = 20/173 (11%)

Query: 1197 VHVYQSSDSIRDRRNSMKLQHAASVTCILYLNNQVFV-SLANGELVVYQREAGHFWDP-- 1253
              ++       D    M              N  V V S       ++       +    
Sbjct: 60   AIIFNPKSGAWDNVPDM-NDKRWYPNNTTLANGDVLVLSGETDGEGLFNELP-QRYVAAT 117

Query: 1254 QNFKSVTLGTQGSPI--TKMVSVGGRLWC-GCQNRVLVLSPDTLQLEH-----MFYVGQD 1305
             +++++T   +  P      ++   +L+  G       L PD            F     
Sbjct: 118  NSWQNLTTAQRKIPYYPHMFLAPNNKLFFSGPWRSSQWLDPDGTGTWFEAPYSHFGGR-- 175

Query: 1306 SSRCVACMVDSSLGVWVTLKG---SAHVCLYHPDTFEQLAEVDVTPPVHRMLA 1355
             S       D  +           +  +   +       A         R   
Sbjct: 176  -SYGGHVYFDGKVLPVGGGNPPTETVELIDLNLPL-PTWAYQTPMSVARRQHN 226


>gi|1568595 [1..151] Ribosomal protein L9 
          Length = 151

 Score = 26.5 bits (58), Expect = 3.1
 Identities = 12/61 (19%), Positives = 24/61 (38%), Gaps = 5/61 (8%)

Query: 348 RNFHLDPKLADILSPRLIRRGSKKRPARSSHQELRRDEGSQ-----DQTGSLSRARPSSR 402
           RNF L   LA + S    ++  + R AR +      +  ++     +  G ++    +S 
Sbjct: 27  RNFLLPRGLAIVASRGAQKQADEIRRARETKSVRDLEHANEIKAAIEALGPIALPVKTSA 86

Query: 403 H 403
            
Sbjct: 87  D 87


 Score = 25.3 bits (55), Expect = 6.6
 Identities = 14/85 (16%), Positives = 27/85 (31%), Gaps = 5/85 (5%)

Query: 687 AQRNIKQVAERIN-KGVRSAEEAERHARVLQE---IEAHIEGMEDLQAPLRRFLRQEMVI 742
           AQ+   ++      K VR  E A      ++    I   ++   D           ++V 
Sbjct: 43  AQKQADEIRRARETKSVRDLEHANEIKAAIEALGPIALPVKTSADSGKLFGSVTAADVVA 102

Query: 743 EVKAIGGKK-DRSLFLFTDLIVCTT 766
            +K  GG   D+ +       +   
Sbjct: 103 AIKKAGGPNLDKRIVRLPKTHIKAV 127


>gi|2117634 [147..843] Lipoxigenase, C-terminal domain 
          Length = 697

 Score = 26.0 bits (57), Expect = 4.0
 Identities = 11/46 (23%), Positives = 16/46 (33%), Gaps = 4/46 (8%)

Query: 1015 PPSLRSPPETAPEPAGPELDVEAAADEEAATLAEPGP-QPCLHISI 1059
            P   R+     PEP  PE +      ++          Q  L IS+
Sbjct: 586  PTLSRNF---MPEPGSPEYEELKTNPDKVFLKTITPQLQTLLGISL 628


>gi|125988 [12..184] EF Hand-like 
          Length = 173

 Score = 25.9 bits (56), Expect = 4.9
 Identities = 17/118 (14%), Positives = 33/118 (27%), Gaps = 5/118 (4%)

Query: 715 LQEIEAHIEGMEDLQAPLRRFLRQEMVIEVKAIGGKKDRSLFL--FTDLIVCTTLKRKSG 772
            + I   I  +  L A      RQE +     +G      + +         + LK K+ 
Sbjct: 24  YELIAVRIAKIAKLSAEKAEETRQEFLRVADQLGLAPGVRISVEEAAVNATDSLLKMKAE 83

Query: 773 SLRRSSMSLYTAASVIDTASKYKMLWKLPLEDADIIKGASQATNRENIQKAISRLDED 830
               + +        IDT     +      E    ++        +      + LD +
Sbjct: 84  EKAMAVIQSLIMYDCIDTDKDGYVSLP---EFKAFLQAVGPDITDDKAITCFNTLDFN 138


>gi|2811252 [82..382] Periplasmic binding protein-like II 
          Length = 301

 Score = 25.3 bits (55), Expect = 6.2
 Identities = 7/60 (11%), Positives = 15/60 (24%), Gaps = 7/60 (11%)

Query: 765 TTLKRKSGSLRRSSMSLYTAASVIDTASKYKMLWKLPLEDADIIKGASQATNRENIQKAI 824
             +K K G    +          +            P    ++ +  +   NR  I +  
Sbjct: 188 DPVKDKPGFSYTNYYRHGGGTIALHCNFG-------PTSLTEVRQAFAYVLNRPKIVELF 240


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

 Score = 25.7 bits (56), Expect = 6.4
 Identities = 12/69 (17%), Positives = 26/69 (37%), Gaps = 5/69 (7%)

Query: 1125 ETTSSEEEQEPGFLPLSGSFGPGGPCGTSPMDGRALRRSSHGSFTRGSLEDLLSVDPEAY 1184
            +T   +     G  P++  +             +AL  SS+    +      L +  + Y
Sbjct: 88   QTLPDQSIVFQGSAPINSWYPAFSRPTPITGV-QALEYSSNAYMVQ----TALGLMGQTY 142

Query: 1185 QSSVWLGTE 1193
            Q ++++GT 
Sbjct: 143  QPNMFVGTS 151


>gi|2497466 [46..370] Purine nucleoside hydrolase 
          Length = 325

 Score = 25.3 bits (55), Expect = 7.8
 Identities = 6/65 (9%), Positives = 17/65 (25%), Gaps = 4/65 (6%)

Query: 1293 TLQLEHMFYVGQDSSRCVACMVDSSLGVWVTLKG--SAHVCLYHPDTFEQLAEVDVTPPV 1350
             L  E   ++  D++  +    D  +     +      H  +      E +        +
Sbjct: 162  NLSAEFNVWIDPDAANYIFR--DPDVKDKCIVVPLNLTHKAIATYKVNEMIYNEKNNSKL 219

Query: 1351 HRMLA 1355
              +  
Sbjct: 220  RELFL 224


>gi|2463100 [37..265] Trypsin-like serine proteases 
          Length = 229

 Score = 25.1 bits (54), Expect = 7.9
 Identities = 18/149 (12%), Positives = 28/149 (18%), Gaps = 34/149 (22%)

Query: 1265 GSPITK------MVSV--GGRLWCGCQNRVLVLSPDTLQLEHMFYVGQDSSRCVACMVDS 1316
            G P+T       M +V    R  C                     V  +           
Sbjct: 7    GKPLTSIVKYPFMATVWQNDRKLCTA-----------------SIVSPNYILSAGHCFVK 49

Query: 1317 SLGVWVTLKGSAHVCLYHPDTFEQLAEVDVTPPVHRMLAGSD-AIIRQHKAACLRITALL 1375
                   +               Q  +V+        + G D AI++   +      A  
Sbjct: 50   MSEENYIILVGTV-NAKLEKGNGQQFKVEKAHVYSETVFGQDIAIVKLKNSIDFSDNATQ 108

Query: 1376 VC----EELLWVGTS---AGVVLTMPTSP 1397
                              AG       S 
Sbjct: 109  PITLSRRSNFTKTDLAFIAGWGRITDWSS 137


>gi|1938424 [576..818] Acyl-CoA dehydrogenase (flavoprotein), N-terminal and middle domains 
          Length = 243

 Score = 25.2 bits (55), Expect = 8.2
 Identities = 4/32 (12%), Positives = 9/32 (27%)

Query: 1371 ITALLVCEELLWVGTSAGVVLTMPTSPGTVSC 1402
            +    +CE +     +  V        G +  
Sbjct: 74   VEYAHICELMGRSIFAPEVFNCQAPDTGNMEV 105


>gi|478310 [306..563] beta/alpha (TIM)-barrel 
          Length = 258

 Score = 24.9 bits (54), Expect = 8.6
 Identities = 9/51 (17%), Positives = 14/51 (26%), Gaps = 3/51 (5%)

Query: 1232 FVSLANGELVVYQREAGHF---WDPQNFKSVTLGTQGSPITKMVSVGGRLW 1279
            F         V           +    FK+     Q +    ++SVGG   
Sbjct: 8    FADATTTPGAVTFNLDSAGLGGYTVDQFKADVRAKQAAGKKVIISVGGEKG 58


>gi|1420859 [36..611] Periplasmic binding protein-like II 
          Length = 576

 Score = 24.9 bits (54), Expect = 9.0
 Identities = 6/26 (23%), Positives = 10/26 (38%)

Query: 799 KLPLEDADIIKGASQATNRENIQKAI 824
           K  L + D  +    A +R + Q   
Sbjct: 294 KKALNNKDFRQAIQFAFDRASFQAQT 319


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

 Score = 24.8 bits (54), Expect = 9.4
 Identities = 5/24 (20%), Positives = 11/24 (45%)

Query: 661 PRYELLVKDLLKHTPEDHPDHPLL 684
           PR +++ +   +   E   + P L
Sbjct: 281 PRAKVMKQTCDEVLQELGINDPQL 304


>gi|544020 [86..429] Nitrogenase iron-molybdenum protein, alpha and beta chains 
          Length = 344

 Score = 24.8 bits (54), Expect = 9.7
 Identities = 12/69 (17%), Positives = 23/69 (32%), Gaps = 2/69 (2%)

Query: 682 PLLLEAQRNIKQVAERINKGVRSAEEAER--HARVLQEIEAHIEGMEDLQAPLRRFLRQE 739
            +    + + +   +   + V  A    R    + L   E  + G     A + +   +E
Sbjct: 183 HISSSKEVDYEPYIDGQTRFVSQAAWFSRSIDCQNLTGKETVVFGDATHAASITKIPARE 242

Query: 740 MVIEVKAIG 748
           M I V   G
Sbjct: 243 MGIRVSYTG 251


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