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: saccharomyces
Comp Matrix: blosum62.bla
Gap open: 12
Gap extension: 2
Start/End local
Target Sequence HS307871
Strand: forward
Start/End (protein) local
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 268.00 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
saccharomyces 4 FPAPKNDLILRAAKGEKVERPPCWIMRQAGRYLP
FP KND +LRAA GE+ + P W MRQAGRYLP
FPELKNDTFLRAAWGEETDYTPVWCMRQAGRYLP
HS307871 1758 tcgcaagatccggtgggagtacgttaccggcttc
tcataaacttgccggaacaacctggtgacggatc
tggggtcacgaacgagaacctctgcgcgactcaa
saccharomyces 38 YHEVKNNRDFFQTCRDAEIASEI
+ E + +DFF TCR E E+
E:E[gag] FRETRAAQDFFSTCRSPEACCEL
HS307871 1860 GGTAAGAG Intron 1 CAGAGtagacggcgttaatctcggttgc
<1-----[1861 : 1975]-1> tgacgccaattgcggccacggat
tgacgtcgctccgtcttgcctag
saccharomyces 62 TIQ PVRRYRGLIDAAIIFSDILVIPQ
T+Q P+RR+ L+DAAIIFSDILV+PQ
TLQ PLRRF--LLDAAIIFSDILVVPQ
HS307871 2047 accGTGAGGG Intron 2 CAGcccct ccgggaattgacggcc
cta<0-----[2056 : 2131]-0>ctggt ttacctttcattttca
tgg agtcc tgttcctccccttacg
saccharomyces 88 AMGMRVEMLEGKGPHFPEPLRNPEDL
A+GM V M+ GKGP FPEPLR +DL
ALGMEVTMVPGKGPSFPEPLREEQDL
HS307871 2195 GTACCCA Intron 3 CAGgcgaggaagcgagcatcgctaggcgc
<0-----[2195 : 2433]-0>ctgtatcttcgagcgtcactgaaaat
agcgggcgatcaacccagaaaaggca
saccharomyces 114 QTVLDYKVDVLKELDWAFKAITMTRIKLDGEVPLFGFCGGP
+ + D +V V EL + F+AIT+TR +L G VPL+GF G P
ERLRDPEV-VASELGYVFQAITLTRQRLAGRVPLIGFAGAP
HS307871 2512 gcccgcgg ggtgcgtgtcgaacaccccggcgccagtgggc
agtgacat tccatgattactctcgagtcggtcttgtcgcc
acagtaag actgactgcaccctcaaagtatgggtctttca
saccharomyces 155 WTLMVYMTEGGGSRLFRFAKQWINMY
WTLM YM EGGGS + AK+W+ +
WTLMTYMVEGGGSSTMAQAKRWLYQR
HS307871 2632 GTAATGT Intron 4 TAGtacaatagggggtaaagcgactctca
<0-----[2632 : 2748]-0>gcttcattagggcgctcacaggtaag
gcggacgtgttcaccgtgcgcgctga
saccharomyces 181 PELSHKLLQKITDVAVEFLSQQVVAGAQ
P+ SH+LL+ +TD V +L QVVAGAQ
PQASHQLLRILTDALVPYLVGQVVAGAQ
HS307871 2827 ccgacccccacaggcgctcggcgggggcGTGAGTC Intron 5
cacgaattgttcacttcattgattcgca<0-----[2911 : 3278]
tgttcggtccctttgcatgaaaggttcg
saccharomyces 209 ILQVFESWGGELSSVDFDEFSLPYLRQIAERVPKRLQELGIMEQIP
LQ+FES G L F++F+LPY+R +A++V RL+E G+ +P
ALQLFESHAGHLGPQLFNKFALPYIRDVAKQVKARLREAGLA-PVP
HS307871 3276 TAGgtcctgtcggccgccctaatgcctacgggacgagatcgggcg cgc
-0>ctattacacgatgcattaatctcatgatcaatacgtgacgtc ctc
agggtgctagttcagcccgtagtccttgcgaggcggggacga agc
saccharomyces 255 M IVFAKGSWYALDKLCCSGFDVVSLD
M I+FAK +AL++L +G++VV LD
M IIFAKDGHFALEELAQAGYEVVGLD
HS307871 3414 aGTGAGGA Intron 6 TAGaatgaggctgcggcgcggtggggcg
t<0-----[3417 : 3575]-0>tttcaagatctaatcacgaattgta
g ccttgtgttcggggcatctggtgtc
saccharomyces 281 WSWDPREA KINKNRVTLQGNLDP
W+ P++A + VTLQGNLDP
WTVAPKKA V:R[cgg] ECVGKTVTLQGNLDP
HS307871 3651 taggcaagCGGTAAGCC Intron 7 CAGGgtggaagatcgacgc
gctccaac <2-----[3677 : 3779]-2> agtgactctagatac
gagcagac gtgggggaggccgcc
saccharomyces 305 GVMYGSK EVITKKVKQMIEAFGGGKS
+Y S+ E I + VKQM++ FG
CALYASE EEIGQLVKQMLDDFG--PH
HS307871 3826 tgttgtgGTAACAG Intron 8 CAGggagctgacacggtg cc
gctacca<0-----[3847 : 4178]-0>aatgattaattaatg ca
tcgtatg ggcggggggggtcta at
saccharomyces 331 RYIVNFGHGTHPFMDPDVIKFFLEECHR
RYI N+GHG +P MDP+ + F++ H+
RYIANLGHGLYPDMDPEHVGAFVDAVHK
HS307871 4230 ctagatgcgctcgagcgcgggtggggca
gatcatgagtacatacaatgcttactaa
cctccgctgtttcgcaatgcctgttgta
//
Gene 1
Gene 1758 4313
Exon 1758 1860 phase 0
Exon 1976 2055 phase 1
Exon 2132 2194 phase 0
Exon 2434 2631 phase 0
Exon 2749 2910 phase 0
Exon 3279 3416 phase 0
Exon 3576 3676 phase 0
Exon 3780 3846 phase 2
Exon 4179 4313 phase 0
//
|
