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:       NP_649556.2
Comp Matrix:         blosum62.bla
Gap open:            12
Gap extension:       2
Start/End            global
Target Sequence      subseq(143000,4000)
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 222.05 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

NP_649556.2        1 MNFEQVNIPEKLVPDNYLQLGLAAQRSKQRSFKELLEKRKLPKNGWSDE 
                      N       +  V      L L    +    +  +L   ++P  G ++  
                     HNEHSHKCMQLRVIRRGNGLPLLRHATLPPPLLSILNYGRVPNEGLNEI 
subseq(143000,4  531 cagctcatactcgacagagccctccgaccccccaacatgcgcaggcaga 
                     aaaacaagtatgttgggagtcttgacctcccttgttaaggtcaagtaat 
                     ctgtctacgagcccaatttcatgttaacctatcccgccgtcatgcacaa 


NP_649556.2       50 RIEELVHMLASLD                          NNYPHKVGLG 
                      I  ++H+++S +                           N    V +G 
                     TIMSILHQISSQN           S:L[ttg]       CNSEKNVKIG 
subseq(143000,4  678 aaaaatccattcaTTGTGCGTA  Intron 1   CAGGtaagaagaag 
                     cttgttaatccaa  <2-----[719  :  823]-2> gagaaatatg 
                     ccgctgcgattat                          cccggccacg 


NP_649556.2       74 EREARIACKLVARRHYNFGHGIG                           
                     ERE RI   LV  ++  FGHGIG                           
                     ERENRIYSALVRNKYIGFGHGIG           R:R[aga]        
subseq(143000,4  855 gcgacatagcgcaatagtgcgagAGGTGAGGA  Intron 2   CAGA 
                     agaagtagcttgaaatgtgagtg  <2-----[926  : 1031]-2>  
                     gggtcaccgcggcgccgcgcccc                           


NP_649556.2       98 SGDLLEAQPKAAGSTLLARLTNALILDLIRGIGLPSCAGCFLVPMCTGM 
                     SG+L + QPK+AG+++LA+ T + + DLI+ +G+  C   +++P  TGM 
                     SGNLDDVQPKSAGNSVLAKATTSFVKDLIKSFGIKGCEDVYILPYATGM 
subseq(143000,4 1033 tgacgggccaaggaagcgagaaatgagcaaatgaagtgggtacctgaga 
                     cgataatacagcgagttcacccgttaattagtgtaggaatattcaccgt 
                     ggccccggggcggccggggccgctagcgcactaaaccatgcagactgcg 


NP_649556.2      147 TLTLCLQSLRKRRPGARYVLWSRIDQKSCFKAITATGLVPVVIPCLIKG 
                      L+ C+   +K R  + YV+ SRID K+C+K I    L  +V+  + +  
                     CLSTCILYTKKEREKSEYVIVSRIDHKTCYKCIDFCALKYLVVDMVYRD 
subseq(143000,4 1180 tcaatattaaagagaagtgagtaagcaattatagttgtattgggagtag 
                     gtgcgttacaaagaagaatttcgtaaacgaagtatgctaatttattaga 
                     tattctgtaagaaaatataagcaaccgcccgcctttcaataagcgatgc 


NP_649556.2      196 ESLNTNVDLFREKIKSLGVDSILCLYTTTSCFAPRNSDDIAEVSKLSKQ 
                     E L+TN+  + + I++ G   I C+ + TS +APRNSDDI ++  + ++ 
                     EELHTNLSEIEKLIQTYG-EKICCVMSVTSSYAPRNSDDIVKIGHMCRR 
subseq(143000,4 1327 ggtcaacagagacacatg gaattgatgatatgcaatggagaagcataa 
                     aatacatgataattacag aatggttctccgaccgacaattatgatggg 
                     aaatatgcaaaataaatt agcttagactattccacatttcaagtgtag 


NP_649556.2      245 WQIPHLVNNAYGLQAKEIVNQLE-CANRVGRIDYFVQSSDKNLLVPVGS 
                     + IPH++NNA+GLQ   +  +++ C    GR+D+ VQS DKN LVPV   
                     YNIPHIINNAFGLQCNYLCKEIQKCFETKGRVDFVVQSCDKNFLVPVNG 
subseq(143000,4 1471 taaccaaaagtgtctatttagacattgaagaggtggcatgaatcgcgag 
                     aatcattaactgtagaatgaataagtacaggtatttaggaaatttctag 
                     ttattttttccgaacttacagaaattgagtaaccacatttgttccatta 


NP_649556.2      293 AIVASFNESVLHDVASTYAGRASGSQSLDVLMTLLSLGRNGFRLLFDQR 
                      IV S ++  + ++   Y GR      LD+ +TLL LG+  +  L  +R 
                     GIVFSSDKKKMKELKKHYPGRTPVHAYLDLFITLLELGKRKILNLRKER 
subseq(143000,4 1618 gagtatgaaaaagtaactcgaacgcgttgctaatcgtgaaaacacaaga 
                     gtttgcaaaataataaaacggcctacatatttcttatgagattatgaag 
                     attccctggggaagggccttactgtatataccaaaaagagatacaaaag 


NP_649556.2      342 GENFNYLRENLRKFAEPRGEIVIDSRFNSISLAITLATLAG----DQMK 
                      ENF +L+  +  +       +I +  N IS+AI L  L         + 
                     EENFAWLKNKVSTLCSKYNLSLIKASKNKISMAINLNELYKIYHVENPR 
subseq(143000,4 1765 ggatgttaaagtacttatatataagaaaaatagaatagctaatcggaca 
                     aaatcgtaaatcctgcaaatgttacgaaatctctataataataataacg 
                     gacccgaatatttaccgctatacatcataacggttgtgtcatctcactg 


NP_649556.2      387 SITKLGSMLHMRGVSGARVIVPGQNKTIDGHEFLGKFRLPPQQPASSVP 
                     SIT LGS+L  R V+G RVI         G +  G       +P  +   
                     SITLLGSLLFYRNVTGHRVICSPLLIRNGGVDPNGVDTTEQGEPKPN-- 
subseq(143000,4 1912 aaaccgttcttaagagcagattcctacaggggcagggaagcggcaca   
                     gtcttgctttagatcgagttgcctttgaggtacagtaccaagacaca   
                     ttcggagacttatctgtgcaccgggaacggatctgtcttggtacact   


NP_649556.2      436 DCSGRNWNNPRGD                                     
                        G ++ NPR                                       
                     ---GVSFPNPRNS                                     
subseq(143000,4 2053    ggttcacaaa                                     
                        gtctcacgag                                     
                        ggtctttatc                                     


//
Gene 1
Gene 531 2082 
  Exon 531 718 phase 0
  Exon 824 925 phase 2
  Exon 1032 2082 phase 2
//
>subseq(143000,4000).[531:2082].sp.tr
HNEHSHKCMQLRVIRRGNGLPLLRHATLPPPLLSILNYGRVPNEGLNEITIMSILHQISS
QNLCNSEKNVKIGERENRIYSALVRNKYIGFGHGIGRSGNLDDVQPKSAGNSVLAKATTS
FVKDLIKSFGIKGCEDVYILPYATGMCLSTCILYTKKEREKSEYVIVSRIDHKTCYKCID
FCALKYLVVDMVYRDEELHTNLSEIEKLIQTYGEKICCVMSVTSSYAPRNSDDIVKIGHM
CRRYNIPHIINNAFGLQCNYLCKEIQKCFETKGRVDFVVQSCDKNFLVPVNGGIVFSSDK
KKMKELKKHYPGRTPVHAYLDLFITLLELGKRKILNLRKEREENFAWLKNKVSTLCSKYN
LSLIKASKNKISMAINLNELYKIYHVENPRSITLLGSLLFYRNVTGHRVICSPLLIRNGG
VDPNGVDTTEQGEPKPNGVSFPNPRNS
//
>subseq(143000,4000).[531:2082].sp
CACAATGAGCATTCCCATAAATGCATGCAATTGCGCGTCATCCGAAGAGGTAATGGTCTC
CCACTTTTGCGTCATGCAACACTCCCCCCTCCACTTCTCAGCATCCTGAACTACGGGCGT
GTCCCAAATGAGGGCCTAAACGAAATAACCATCATGAGCATTTTGCACCAGATATCTTCT
CAAAATTTGTGCAACAGCGAGAAGAACGTCAAAATCGGGGAGCGGGAGAATCGCATATAC
AGCGCGCTCGTGCGGAACAAGTACATCGGGTTCGGGCACGGCATCGGCAGATCGGGGAAC
CTCGACGACGTGCAGCCGAAGAGCGCGGGGAACAGCGTGCTGGCGAAGGCCACCACGAGC
TTTGTAAAGGACCTGATCAAAAGCTTTGGAATAAAAGGCTGCGAAGATGTGTACATACTG
CCATACGCTACGGGCATGTGTCTAAGTACTTGCATTTTGTATACAAAAAAGGAAAGAGAA
AAAAGTGAATATGTAATAGTGTCCAGAATAGACCACAAGACCTGCTACAAGTGCATCGAT
TTTTGTGCCTTAAAATATTTAGTAGTGGACATGGTATATAGGGACGAAGAATTACATACA
AATCTGAGCGAAATAGAAAAACTTATACAAACATATGGTGAAAAGATCTGTTGTGTAATG
TCAGTCACTTCAAGTTATGCCCCCAGAAACTCAGATGATATTGTCAAAATAGGGCATATG
TGTAGAAGGTATAATATACCTCATATTATTAATAATGCCTTCGGGTTACAATGCAATTAT
TTATGCAAAGAGATACAAAAATGTTTTGAGACAAAGGGTAGAGTAGACTTCGTAGTCCAA
AGTTGTGATAAGAATTTTCTCGTCCCAGTTAATGGAGGAATTGTTTTCAGCTCCGATAAG
AAGAAGATGAAAGAATTGAAGAAGCACTACCCTGGTAGAACCCCTGTGCATGCATATTTA
GATCTATTCATCACATTACTAGAATTAGGGAAAAGGAAAATTCTAAACCTAAGAAAAGAA
AGGGAGGAAAACTTCGCCTGGTTAAAAAATAAAGTTTCTACTCTATGCTCCAAGTACAAT
TTAAGTTTAATCAAAGCTAGCAAAAATAAAATATCCATGGCGATTAATTTGAATGAGCTT
TACAAAATTTACCATGTCGAAAACCCTAGGAGTATTACCCTGCTGGGATCGTTACTCTTT
TATAGAAATGTCACTGGGCATAGGGTCATATGCTCCCCGCTGTTGATACGAAACGGGGGG
GTAGATCCCAATGGGGTTGACACTACTGAGCAGGGTGAACCCAAACCCAATGGGGTGTCT
TTCCCTAATCCTAGAAATAGC
//