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: sla/lp_dsp
Comp Matrix: blosum62.bla
Gap open: 12
Gap extension: 2
Start/End global
Target Sequence subseq(150000,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 215.56 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
sla/lp_dsp 1 MNFEQVNIPEKLVPDNYLQLGLAAQRSKQRSFKELLEKRKLPKNGWSDE
+ + + + N L L Q + +L ++P G ++
NEYSHKYMHVRGIQKNN-DLTLLCQTTINLPPYSILNYGRVPNEGLNEI
subseq(150000,5 2914 agttcatacgcgacaaa gcacttcaaaaccctaacatgcgcaggtaga
aaacaaatatggtaaaa atcttgacctatccagttaaggtcaagtaat
tgtttacgcgcaaaatt tcatgtaaacctctcccatcctcatgaacaa
sla/lp_dsp 50 RIEELVHMLASLD NNYPHKVGLG
I ++H+++S + N V +G
TIMSILHQISSQN S:L[tta] CNSEKNVKIG
subseq(150000,5 3058 aaaaatccattcaTTGTACGTA Intron 1 CAGAtaagaagaag
cttgttaatccaa <2-----[3099 : 3219]-2> gagaaatatg
ccgctgtgattat cccggctaag
sla/lp_dsp 74 EREARIACKLVARRHYNFGHGIG
ERE RI LV ++ FGHGIG
ERENRIYSALVRNKYIGFGHGIG R:R[agg]
subseq(150000,5 3251 gcgacatagcgcaatagtgcgagAGGTAAGCA Intron 2 CAGG
agaagtagcttgaaatgtgagtg <2-----[3322 : 3429]-2>
aggtaatcgcgtcgcaatgcctc
sla/lp_dsp 98 SGDLLEAQPKAAGSTLLARLTNALILDLIRGIGLPSCAGCFLVPMCTGM
SG+L + QPK+AG+++LA+ T + + DLIR +G+ C +++P TGM
SGNLDDVQPKSAGNSILAKATTSFVKDLIRSFGIKGCEDVYILPYATGM
subseq(150000,5 3431 tgatgggccaaggaaacgagaaatgagtaaatgaagtgggtacctgaga
cgataatacagcgagttcacccgttaattggtgtaggaatattcaccgt
agcaccggtgctactcgagccaccagtacactaaaccgtgcaaactagg
sla/lp_dsp 147 TLTLCLQSLRKRRPGARYVLWSRIDQKSCFKAITATGLVPVVIPCLIKG
L+ CL +K R + YV+ SRID K+C+K I L +V+ + K
CLSTCLLYTKKEREKSEYVIVSRIDHKTCYKCIDFCALKYLVVDMVYKD
subseq(150000,5 3578 ttaatcttaaagagaagtgagtaagcaattatagttgtattgggagtag
gtgcgttacaaagaagaatttcgtaaacgaagtatgctaatttattaaa
tattcggcaagaaaatataagtaattgcccgtctttcggtaagtgtcgt
sla/lp_dsp 196 ESLNTNVDLFREKIKSLGVDSILCLYTTTSCFAPRNSDDIAEVSKLSKQ
E L TN++ + I++ G I C+ + TS +APRNSDDI +++ + K+
EELFTNLNEIERLIQTYG-EKICCVMSVTSSYAPRNSDDIVKIAHICKR
subseq(150000,5 3725 ggttaacagagacacatg gaattgatgatatgcaatggagaagcataa
aattcataatagttacag aatggttctccgaccgacaattatcatgag
aaatacgcaagagaaatt gaattcgagatttaagtatctaaattttgg
sla/lp_dsp 245 WQIPHLVNNAYGLQAKEIVNQLE-CANRVGRIDYFVQSSDKNLLVPVGS
+ +PH++NNA+GLQ + +++ C GR+D+ VQS DKN LVPV
YNVPHIINNAFGLQCNYLCKEIQKCFESKGRVDFVVQSCDKNFLVPVNG
subseq(150000,5 3869 tagccaaaagtgtctatttagacattgtagcggtggcatgaatcgcgag
aatcattaactgtagaatgaataagtacaggtatttaggaaatttctag
ttaacatttctaagcttgtaaagattagatagttgtgtttgtctcagcc
sla/lp_dsp 293 AIVASFNESVLHDVASTYAGRASGSQSLDVLMTLLSLGRNGFRLLFDQR
IV S ++ + ++ Y GR LD+ +TLL LG+ + L +R
GIVFSSDKKKMKELKKHYPGRTPVHAYLDLFITLLELGKKKILNLRKER
subseq(150000,5 4016 gagtatgaaaaagtaactcgaacgcgttgctaatcgtgaaaacacaaga
gtttgcaaaataataaaacggcctacatatttcttatgaaattatgaag
gcactctgaggaaaagccttaccgtatattctaaaaagaaatacaaagg
sla/lp_dsp 342 GENFNYLRENLRKFAEPRGEIVIDSRFNSISLAITLATLAGD----QMK
ENF++L+E + + +I + N IS+AI L L K
EENFSWLQEKVSTLCSKYNLTLIKASKNKISMAINLNELYNICHILNPK
subseq(150000,5 4163 ggattttcgagtacttatatataagaaaaatagaatagctaatcacaca
aaatcgtaaatcctgcaaatcttacgaaatctctataataatgattaca
gattcgggaatttccagctatgcatcataacgcttatatcttcttacgg
sla/lp_dsp 387 SITKLGSMLHMRGVSGARVIVPGQNKTIDGHEFLGKFRLPPQQPASSVP
SIT LGS+L R V+G RVI + + + + R QQ ++V
SITLLGSLLFYRNVTGHRVIC---SPLLIRNGGVAQRR--DQQIQTNVE
subseq(150000,5 4310 aaaccgtccttaagagcagat tcttacaggggcaa gccacaagg
gtcttgctttagatcgagttg cctttgaggtcagg aaatacata
ttcacttatttattcgtgtac cgggactgagagaa cgatgttaa
sla/lp_dsp 436 DCSGRNWNNPRGD
+ +
PAFPIHHGAVDPN
subseq(150000,5 4442 cgtcaccggggca
cctctaagctaca
tatcttcgaatct
//
Gene 1
Gene 2914 4480
Exon 2914 3098 phase 0
Exon 3220 3321 phase 2
Exon 3430 4480 phase 2
//
>subseq(150000,5000).[2914:4480].sp.tr
NEYSHKYMHVRGIQKNNDLTLLCQTTINLPPYSILNYGRVPNEGLNEITIMSILHQISSQ
NLCNSEKNVKIGERENRIYSALVRNKYIGFGHGIGRSGNLDDVQPKSAGNSILAKATTSF
VKDLIRSFGIKGCEDVYILPYATGMCLSTCLLYTKKEREKSEYVIVSRIDHKTCYKCIDF
CALKYLVVDMVYKDEELFTNLNEIERLIQTYGEKICCVMSVTSSYAPRNSDDIVKIAHIC
KRYNVPHIINNAFGLQCNYLCKEIQKCFESKGRVDFVVQSCDKNFLVPVNGGIVFSSDKK
KMKELKKHYPGRTPVHAYLDLFITLLELGKKKILNLRKEREENFSWLQEKVSTLCSKYNL
TLIKASKNKISMAINLNELYNICHILNPKSITLLGSLLFYRNVTGHRVICSPLLIRNGGV
AQRRDQQIQTNVEPAFPIHHGAVDPN
//
>subseq(150000,5000).[2914:4480].sp
AATGAGTATTCTCATAAATACATGCACGTGCGCGGAATACAAAAAAATAATGATCTCACA
CTTTTGTGTCAAACAACAATCAACCTTCCCCCTTACAGCATCCTAAATTACGGCCGTGTC
CCAAATGAGGGATTAAACGAAATAACCATCATGAGCATTTTGCATCAGATATCTTCTCAA
AATTTATGCAACAGCGAGAAGAACGTTAAAATAGGGGAACGGGAGAATCGAATATATAGC
GCGCTCGTGCGTAACAAGTACATAGGATTTGGGCACGGCATTGGCAGGTCAGGGAACTTA
GACGACGTGCAGCCTAAGAGCGCTGGAAACAGTATCCTGGCAAAGGCCACCACAAGCTTC
GTAAAGGATTTAATCAGAAGCTTTGGAATAAAAGGCTGCGAGGATGTGTACATACTACCA
TACGCTACAGGGATGTGTTTAAGTACTTGCCTGTTGTACACAAAAAAGGAAAGAGAAAAA
AGTGAATATGTAATAGTGTCTAGAATAGATCATAAGACCTGCTACAAGTGTATCGATTTT
TGTGCCTTGAAGTATTTAGTAGTGGATATGGTTTACAAGGATGAAGAATTATTTACAAAC
CTGAACGAAATAGAGAGACTGATACAAACATATGGTGAGAAAATATGTTGTGTCATGTCA
GTGACATCTAGTTATGCACCAAGGAATTCAGATGACATTGTAAAAATAGCTCATATTTGT
AAGAGGTATAATGTACCACACATAATTAATAATGCCTTTGGATTACAGTGCAATTATTTG
TGTAAAGAAATACAGAAATGTTTTGAATCGAAAGGTCGAGTGGATTTTGTGGTTCAGAGT
TGTGATAAGAATTTCCTTGTCCCAGTGAACGGCGGGATCGTATTCAGTTCCGATAAGAAA
AAGATGAAAGAATTAAAAAAGCACTACCCTGGTAGAACCCCCGTGCATGCATATTTAGAT
CTTTTCATTACATTACTAGAATTAGGGAAAAAAAAAATTCTAAACCTAAGAAAAGAGAGG
GAGGAAAATTTTTCCTGGTTGCAGGAAAAAGTTTCTACTCTCTGCTCAAAGTACAATTTA
ACTTTGATCAAAGCTAGCAAAAATAAAATATCCATGGCCATTAATTTAAATGAACTTTAC
AATATTTGCCATATTCTAAACCCGAAGAGTATTACCCTACTCGGTTCTCTACTTTTTTAT
AGAAATGTTACCGGGCATAGGGTTATATGCTCCCCGTTGTTGATACGCAATGGGGGAGTG
GCACAGAGAAGAGACCAGCAAATTCAGACTAATGTAGAACCTGCATTTCCCATTCATCAC
GGGGCAGTAGATCCCAAT
//
