genewise $Name: wise2-2-0 $ (unreleased release)
This program is freely distributed under a GPL. See source directory
Copyright (c) GRL limited: portions of the code are from separate copyright

Query protein:       Sel1
Comp Matrix:         blosum62.bla
Gap open:            12
Gap extension:       2
Start/End            global
Target Sequence      subseq(2083000,5000)
Strand:              forward
Start/End (protein)  global
Gene Paras:          human.gf
Codon Table:         codon.table
Subs error:          1e-05
Indel error:         1e-05
Model splice?        model
Model codon bias?    flat
Model intron bias?   tied
Null model           syn
Algorithm            623

genewise output
Score 105.18 bits over entire alignment
Scores as bits over a synchronous coding model

Warning: The bits scores is not probablistically correct for single seqs
See WWW help for more info

Sel1               1 MDDRKENKKN-ADYEIFLRDRKLKYEKIRKRNALHNKFILPIINFIKGI 
                     M+ +KEN++N  +YE  LRDRKLKYE +RK++ LHN+F+ PI+N+I    
                     MNRKKENQENEGEYENVLRDRKLKYENLRKKSTLHNRFMSPIFNLINSA 
subseq(2083000, 2944 aacaagacgagggtgagtagcacatgaccaaaaccaatatcatacaaag 
                     tagaaaaaaaagaaaattgagataaaatgaagctaagttccttattagc 
                     gttaaacaatggatgctaattaggtgtggaatcctcgtggcatcaatta 


Sel1              49 INFIKTVINI                       LAIFFKTLLGQSNTSQ 
                     INF+KT++NI                       LAIFFKTL G SN S  
                     INFLKTLFNI                       LAIFFKTLFGLSNGSN 
subseq(2083000, 3091 aatcaactaaGTAATGC  Intron 1   CAGcgattaattgttagta 
                     tattacttat<0-----[3121 : 3321]-0>tctttacttgtcagca 
                     accgacttca                       attttaaatcgttcat 


Sel1              75                        SSNGKNNDDDDGFFKKKRT-GGLPKS 
                                            +  G+N D DD FF KKR  G L KS 
                                            TLGGRNGDGDDDFFNKKRKMGTLQKS 
subseq(2083000, 3370 GTAACCC  Intron 2   CAGacggcaggggggttaaaaaagatcaa 
                     <0-----[3370 : 3531]-0>ctgggagagaaattaaagatgctaag 
                                            cgcaccccatttcccggaagatgaac 


Sel1             100 RIMELKNLGTCLGST                                   
                     R+MELKNLGTCLGS+                                   
                     RMMELKNLGTCLGSS                                   
subseq(2083000, 3610 aaagcaatgattgat                                   
                     gttataatgcgtggc                                   
                     aggaggtgcctgacc                                   


//
Gene 1
Gene 2944 3654 
  Exon 2944 3120 phase 0
  Exon 3322 3369 phase 0
  Exon 3532 3654 phase 0
//
>subseq(2083000,5000).[2944:3654].sp.tr
MNRKKENQENEGEYENVLRDRKLKYENLRKKSTLHNRFMSPIFNLINSAINFLKTLFNIL
AIFFKTLFGLSNGSNTLGGRNGDGDDDFFNKKRKMGTLQKSRMMELKNLGTCLGSS
//
>subseq(2083000,5000).[2944:3654].sp
ATGAATCGTAAAAAAGAAAACCAAGAAAATGAGGGGGAATATGAGAACGTTTTAAGAGAT
CGTAAACTGAAGTATGAGAATCTGCGGAAAAAAAGTACCCTCCATAACAGGTTTATGTCG
CCCATATTTAACCTAATAAATAGTGCAATAAACTTCCTGAAAACCCTTTTTAACATACTA
GCTATTTTTTTTAAAACATTATTTGGCTTGTCTAATGGCTCAAATACCCTGGGCGGACGC
AACGGCGACGGAGATGATGATTTCTTCAACAAGAAGAGAAAAATGGGAACTTTGCAAAAA
AGCAGAATGATGGAACTGAAGAATTTGGGCACCTGTTTGGGAAGCTCC TGAGGGTAA
//