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:       Sel2
Comp Matrix:         blosum62.bla
Gap open:            12
Gap extension:       2
Start/End            global
Target Sequence      subseq(1128000,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 137.65 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

Sel2               1 MDTNENMKNVKDMMSFTTKNIIYIFIGISLLIFIYKILKKNKKDAEVKR 
                     MD   N    K +  ++ KN+I + +G++  I+IYK++K  KK +E +R 
                     MDREGNNAPTKGVDLLSVKNVILLVLGVTSFILIYKFMKYRKKVSEFER 
subseq(1128000, 3131 agaggaagcaagggtcagaagaccgtggattacatataatcaagtgtga 
                     tagagaacccagtattgtaattttttgtccttttaattaagaatcatag 
                     gcggccttcagagtagcgatgcgaggctattatttgcgatagatagtgg 


Sel2              50 NNEISIKMKLSREKQLQELDKEMMINKEKMKEQNIKKNEEKKKDADQAK 
                      N+I  KMK+SREK+LQEL+KEM++NKEKM+EQ  K N++K K  D  K 
                     KNKIDEKMKMSREKRLQELEKEMIVNKEKMREQMNKGNDKKDKALDSDK 
subseq(1128000, 3278 aaaaggaaaatagaaccgtgagaagaagaaagcaaagagaagagcgtga 
                     aaataaatatcgaagtaataaatttaaaatgaataagaaaaaactacaa 
                     acaataggagagaaggaagaaagcttaaagaaagcgatcaataaacaca 


Sel2              99 PKLGSKDNSSFNHLNDYSNYYRPSLKNRFYNNRKSR              
                      K  SKDNSSF+H  D SNYYRPS++NR+      +              
                     AKPDSKDNSSFSHFRDLSNYYRPSIRNRLCE----K              
subseq(1128000, 3425 gacgtagatttactagttattactaaaattg    a              
                     cacacaaacctgatgatcaaagcctgagtga    a              
                     aaatcgttggccttgtaatccggaaatgaca    g              


//
Gene 1
Gene 3131 3520 
  Exon 3131 3520 phase 0
//
>subseq(1128000,5000).[3131:3520].sp.tr
MDREGNNAPTKGVDLLSVKNVILLVLGVTSFILIYKFMKYRKKVSEFERKNKIDEKMKMS
REKRLQELEKEMIVNKEKMREQMNKGNDKKDKALDSDKAKPDSKDNSSFSHFRDLSNYYR
PSIRNRLCEK
//
>subseq(1128000,5000).[3131:3520].sp
ATGGACAGGGAGGGCAACAATGCTCCCACAAAGGGAGTGGATTTACTGAGCGTGAAAAAT
GTGATCCTGCTAGTGTTGGGCGTTACATCTTTTATACTTATTTATAAGTTCATGAAATAT
CGAAAGAAAGTTTCAGAGTTTGAGAGGAAAAACAAAATAGATGAAAAGATGAAAATGTCA
AGGGAAAAAAGGCTGCAAGAATTGGAAAAAGAAATGATCGTTAATAAAGAAAAAATGAGA
GAACAAATGAACAAGGGAAATGACAAAAAAGATAAAGCACTAGACTCAGACAAAGCAAAA
CCAGATTCCAAGGATAATTCGTCGTTCAGCCATTTTAGGGATTTATCAAATTACTACAGG
CCGTCAATAAGAAATAGGTTATGCGAAAAG
//