# /hgtech/tools/solaris8/bin/fasta34_t -T 8 -b50 -d10 -E0.01 -H -Oeh00728.fasta.nr -Q eh00728.ptfa /cdna2/lib/nr/nr 2 FASTA searches a protein or DNA sequence data bank version 34.26.5 April 26, 2007 Please cite: W.R. Pearson & D.J. Lipman PNAS (1988) 85:2444-2448 eh00728, 164 aa vs /cdna2/lib/nr/nr library 2355815254 residues in 6825066 sequences statistics sampled from 60000 to 6824814 sequences Expectation_n fit: rho(ln(x))= 5.1112+/-0.000182; mu= 7.9857+/- 0.010 mean_var=67.5887+/-13.323, 0's: 45 Z-trim: 46 B-trim: 2907 in 1/66 Lambda= 0.156005 FASTA (3.5 Sept 2006) function [optimized, BL50 matrix (15:-5)] ktup: 2 join: 36, opt: 24, open/ext: -10/-2, width: 16 The best scores are: opt bits E(6825066) gi|62089396|dbj|BAD93142.1| mitochondrial ribosoma ( 164) 1090 253.2 1.1e-65 gi|13633907|sp|P82914|RT15_HUMAN 28S ribosomal pro ( 257) 838 196.6 1.9e-48 gi|114555553|ref|XP_524665.2| PREDICTED: mitochond ( 257) 826 193.9 1.2e-47 gi|90083306|dbj|BAE90735.1| unnamed protein produc ( 132) 735 173.2 1e-41 gi|109002021|ref|XP_001111055.1| PREDICTED: mitoch ( 253) 621 147.8 9.3e-34 gi|149694025|ref|XP_001503663.1| PREDICTED: simila ( 258) 621 147.8 9.4e-34 gi|76614526|ref|XP_585783.2| PREDICTED: similar to ( 259) 569 136.1 3.1e-30 gi|12005982|gb|AAG44697.1|AF265439_1 DC37 [Homo sa ( 283) 524 126.0 3.8e-27 gi|148698367|gb|EDL30314.1| mitochondrial ribosoma ( 136) 513 123.3 1.2e-26 gi|149023938|gb|EDL80435.1| mitochondrial ribosoma ( 135) 503 121.0 5.5e-26 gi|12842922|dbj|BAB25785.1| unnamed protein produc ( 235) 483 116.7 2e-24 gi|26341146|dbj|BAC34235.1| unnamed protein produc ( 258) 483 116.7 2.1e-24 gi|12845693|dbj|BAB26857.1| unnamed protein produc ( 258) 483 116.7 2.1e-24 gi|33585821|gb|AAH55861.1| Mitochondrial ribosomal ( 258) 483 116.7 2.1e-24 gi|13959567|sp|Q9DC71|RT15_MOUSE 28S ribosomal pro ( 258) 483 116.7 2.1e-24 gi|12833594|dbj|BAB22586.1| unnamed protein produc ( 258) 483 116.7 2.1e-24 gi|53734559|gb|AAH83856.1| Mitochondrial ribosomal ( 257) 473 114.5 1e-23 gi|126330382|ref|XP_001380812.1| PREDICTED: simila ( 362) 418 102.2 7e-20 gi|73977047|ref|XP_532559.2| PREDICTED: similar to ( 221) 350 86.7 1.9e-15 gi|149638508|ref|XP_001511018.1| PREDICTED: simila ( 191) 279 70.7 1.1e-10 gi|62858843|ref|NP_001017066.1| mitochondrial ribo ( 226) 222 57.9 9.1e-07 gi|89266810|emb|CAJ81537.1| mitochondrial ribosoma ( 248) 222 58.0 9.8e-07 gi|161611760|gb|AAI55947.1| Unknown (protein for I ( 253) 213 55.9 4.1e-06 gi|119627766|gb|EAX07361.1| mitochondrial ribosoma ( 161) 199 52.7 2.5e-05 gi|47215059|emb|CAG03494.1| unnamed protein produc ( 254) 190 50.8 0.00015 gi|49900813|gb|AAH76324.1| Zgc:92862 [Danio rerio] ( 263) 188 50.3 0.00021 gi|24212640|sp|P82913.2|RT15_BOVIN Mitochondrial 2 ( 96) 177 47.6 0.00052 >>gi|62089396|dbj|BAD93142.1| mitochondrial ribosomal pr (164 aa) initn: 1090 init1: 1090 opt: 1090 Z-score: 1334.4 bits: 253.2 E(): 1.1e-65 Smith-Waterman score: 1090; 100.000% identity (100.000% similar) in 164 aa overlap (1-164:1-164) 10 20 30 40 50 60 eh0072 PVRAAKRTAAAARAPRGLEVTMLRVAWRTLSLIRTRAVTQVLVPGLPGGGSAKFPFNQWG :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: gi|620 PVRAAKRTAAAARAPRGLEVTMLRVAWRTLSLIRTRAVTQVLVPGLPGGGSAKFPFNQWG 10 20 30 40 50 60 70 80 90 100 110 120 eh0072 LQPRSLLLQAARGYVVRKPAQSRLDDDPPPSTLLKDYQNVPGIEKVDDVVKRLLSLEMAN :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: gi|620 LQPRSLLLQAARGYVVRKPAQSRLDDDPPPSTLLKDYQNVPGIEKVDDVVKRLLSLEMAN 70 80 90 100 110 120 130 140 150 160 eh0072 KKEMLKIKQEQFMKKIVANPEDTRSLEARSWSGTLGEWGCREAC :::::::::::::::::::::::::::::::::::::::::::: gi|620 KKEMLKIKQEQFMKKIVANPEDTRSLEARSWSGTLGEWGCREAC 130 140 150 160 >>gi|13633907|sp|P82914|RT15_HUMAN 28S ribosomal protein (257 aa) initn: 838 init1: 838 opt: 838 Z-score: 1025.1 bits: 196.6 E(): 1.9e-48 Smith-Waterman score: 838; 100.000% identity (100.000% similar) in 128 aa overlap (22-149:1-128) 10 20 30 40 50 60 eh0072 PVRAAKRTAAAARAPRGLEVTMLRVAWRTLSLIRTRAVTQVLVPGLPGGGSAKFPFNQWG ::::::::::::::::::::::::::::::::::::::: gi|136 MLRVAWRTLSLIRTRAVTQVLVPGLPGGGSAKFPFNQWG 10 20 30 70 80 90 100 110 120 eh0072 LQPRSLLLQAARGYVVRKPAQSRLDDDPPPSTLLKDYQNVPGIEKVDDVVKRLLSLEMAN :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: gi|136 LQPRSLLLQAARGYVVRKPAQSRLDDDPPPSTLLKDYQNVPGIEKVDDVVKRLLSLEMAN 40 50 60 70 80 90 130 140 150 160 eh0072 KKEMLKIKQEQFMKKIVANPEDTRSLEARSWSGTLGEWGCREAC ::::::::::::::::::::::::::::: gi|136 KKEMLKIKQEQFMKKIVANPEDTRSLEARIIALSVKIRSYEEHLEKHRKDKAHKRYLLMS 100 110 120 130 140 150 >>gi|114555553|ref|XP_524665.2| PREDICTED: mitochondrial (257 aa) initn: 826 init1: 826 opt: 826 Z-score: 1010.5 bits: 193.9 E(): 1.2e-47 Smith-Waterman score: 826; 97.656% identity (100.000% similar) in 128 aa overlap (22-149:1-128) 10 20 30 40 50 60 eh0072 PVRAAKRTAAAARAPRGLEVTMLRVAWRTLSLIRTRAVTQVLVPGLPGGGSAKFPFNQWG ::::::::::::::.::::.::::::::::::::::::: gi|114 MLRVAWRTLSLIRTQAVTQALVPGLPGGGSAKFPFNQWG 10 20 30 70 80 90 100 110 120 eh0072 LQPRSLLLQAARGYVVRKPAQSRLDDDPPPSTLLKDYQNVPGIEKVDDVVKRLLSLEMAN :::::::::::::::.:::::::::::::::::::::::::::::::::::::::::::: gi|114 LQPRSLLLQAARGYVIRKPAQSRLDDDPPPSTLLKDYQNVPGIEKVDDVVKRLLSLEMAN 40 50 60 70 80 90 130 140 150 160 eh0072 KKEMLKIKQEQFMKKIVANPEDTRSLEARSWSGTLGEWGCREAC ::::::::::::::::::::::::::::: gi|114 KKEMLKIKQEQFMKKIVANPEDTRSLEARIIALSVKIRSYEEHLEKHRKDKAHKRYLLMS 100 110 120 130 140 150 >>gi|90083306|dbj|BAE90735.1| unnamed protein product [M (132 aa) initn: 674 init1: 674 opt: 735 Z-score: 904.0 bits: 173.2 E(): 1e-41 Smith-Waterman score: 735; 85.294% identity (92.647% similar) in 136 aa overlap (22-157:1-132) 10 20 30 40 50 60 eh0072 PVRAAKRTAAAARAPRGLEVTMLRVAWRTLSLIRTRAVTQVLVPGLPGGGSAKFPFNQWG :::::::.::::::.: :: ::::::::.::::: : gi|900 MLRVAWRALSLIRTQA--PVL--GLPGGGSARFPFNQRG 10 20 30 70 80 90 100 110 120 eh0072 LQPRSLLLQAARGYVVRKPAQSRLDDDPPPSTLLKDYQNVPGIEKVDDVVKRLLSLEMAN :::::::::::::::...::: .::::::::::::::::::::::::::::::::::: gi|900 LQPRSLLLQAARGYVIQNPAQPGQNDDPPPSTLLKDYQNVPGIEKVDDVVKRLLSLEMAN 40 50 60 70 80 90 130 140 150 160 eh0072 KKEMLKIKQEQFMKKIVANPEDTRSLEARSWSGTLGEWGCREAC ::::::.:.::.:::: :::::::::::::::::::: gi|900 KKEMLKMKREQLMKKIDANPEDTRSLEARSWSGTLGE 100 110 120 130 >>gi|109002021|ref|XP_001111055.1| PREDICTED: mitochondr (253 aa) initn: 560 init1: 560 opt: 621 Z-score: 761.2 bits: 147.8 E(): 9.3e-34 Smith-Waterman score: 621; 78.125% identity (85.938% similar) in 128 aa overlap (22-149:1-124) 10 20 30 40 50 60 eh0072 PVRAAKRTAAAARAPRGLEVTMLRVAWRTLSLIRTRAVTQVLVPGLPGGGSAKFPFNQWG :::::::.::::::.: :: ::::::::.::::: : gi|109 MLRVAWRALSLIRTQA--PVL--GLPGGGSARFPFNQRG 10 20 30 70 80 90 100 110 120 eh0072 LQPRSLLLQAARGYVVRKPAQSRLDDDPPPSTLLKDYQNVPGIEKVDDVVKRLLSLEMAN :::::::::::::::...::: .:::::::::::::::::::: ::::::: gi|109 LQPRSLLLQAARGYVIQNPAQPGQNDDPPPSTLLKDYQNVPGIEKXXXXXXXXLSLEMAN 40 50 60 70 80 90 130 140 150 160 eh0072 KKEMLKIKQEQFMKKIVANPEDTRSLEARSWSGTLGEWGCREAC ::::::.:.::.:::: :::::::::::: gi|109 KKEMLKMKREQLMKKIDANPEDTRSLEARIVALSVKIRSYEEHMQKHRKDKAHKRYLLMS 100 110 120 130 140 150 >>gi|149694025|ref|XP_001503663.1| PREDICTED: similar to (258 aa) initn: 622 init1: 465 opt: 621 Z-score: 761.1 bits: 147.8 E(): 9.4e-34 Smith-Waterman score: 621; 72.308% identity (89.231% similar) in 130 aa overlap (22-149:1-130) 10 20 30 40 50 60 eh0072 PVRAAKRTAAAARAPRGLEVTMLRVAWRTLSLIRTRAVTQVLVPGLPGGGSAKFPFNQWG :::.:::.:: .::::::.. : :::: :.:..: . :: gi|149 MLRAAWRALSSVRTRAVTHAPVLGLPGRGGARLPCDPWG 10 20 30 70 80 90 100 110 eh0072 L--QPRSLLLQAARGYVVRKPAQSRLDDDPPPSTLLKDYQNVPGIEKVDDVVKRLLSLEM : ::.:.:::::::.::::.. .::::::::::::::.:::::::::::::::::: gi|149 LLSPPRGLILQAARGYAVRKPVRPSQEDDPPPSTLLKDYQNIPGIEKVDDVVKRLLSLEM 40 50 60 70 80 90 120 130 140 150 160 eh0072 ANKKEMLKIKQEQFMKKIVANPEDTRSLEARSWSGTLGEWGCREAC :..:: ::.::::.:.:.::::::: ::::: gi|149 ASRKEKLKVKQEQLMNKVVANPEDTSSLEARIVALTVKIRNYEEHMRKHRKDKAHKRYLL 100 110 120 130 140 150 gi|149 MSIDQRNKMLKNLRKTNYVVFEKACKELGIEYTIPPFYKRKAHRRWVTKKALCIRVFQEA 160 170 180 190 200 210 >>gi|76614526|ref|XP_585783.2| PREDICTED: similar to mit (259 aa) initn: 589 init1: 426 opt: 569 Z-score: 697.8 bits: 136.1 E(): 3.1e-30 Smith-Waterman score: 569; 70.769% identity (86.154% similar) in 130 aa overlap (22-149:1-130) 10 20 30 40 50 60 eh0072 PVRAAKRTAAAARAPRGLEVTMLRVAWRTLSLIRTRAVTQVLVPGLPGGGSAKFPFNQWG :::.:::.:: :::.:::: : :::::: ::. : gi|766 MLRAAWRALSSIRTQAVTQPPVLGLPGGGCAKLLSVQRD 10 20 30 70 80 90 100 110 eh0072 L--QPRSLLLQAARGYVVRKPAQSRLDDDPPPSTLLKDYQNVPGIEKVDDVVKRLLSLEM : .: .:.:::::::...::.:. .:::::: :::::::::::::::::::::.:::: gi|766 LPSSPGGLILQAARGYATQKPVQQSQEDDPPPSMLLKDYQNVPGIEKVDDVVKRLISLEM 40 50 60 70 80 90 120 130 140 150 160 eh0072 ANKKEMLKIKQEQFMKKIVANPEDTRSLEARSWSGTLGEWGCREAC ::::::::::.::.:.:.: ::.:: ::::. gi|766 ANKKEMLKIKKEQLMSKVVENPKDTSSLEAQIVALTVKIRSYEEHMQKHQKDKAHKRYLL 100 110 120 130 140 150 gi|766 MSIDKRQKMLKNLRKTNYPVFEKTCKELGIEYTFPPPYHRKIHRRLATKKALCIRVFQEV 160 170 180 190 200 210 >>gi|12005982|gb|AAG44697.1|AF265439_1 DC37 [Homo sapien (283 aa) initn: 496 init1: 496 opt: 524 Z-score: 642.6 bits: 126.0 E(): 3.8e-27 Smith-Waterman score: 525; 66.013% identity (75.817% similar) in 153 aa overlap (6-149:11-154) 10 20 30 40 50 eh0072 PVRAAKRTAAAARAPRGLEVTMLRVAWRTLSLIRTRAVTQVLVPGLPGGGSAKFP .:::::::::::::::::::::::::::::::::::::::::. :. .: gi|120 MGAHVTRSVPQRTAAAARAPRGLEVTMLRVAWRTLSLIRTRAVTQVLVPGLPAVGAPSF- 10 20 30 40 50 60 70 80 90 100 eh0072 FNQWGLQPRSLLLQAARGYVVRKPAQSRLDDDPPP-----STLL----KDYQNVPGIEKV .. : :: ....: : . . . : . :: : . :::::: gi|120 LSTSGAC--SLEVSSSR------PRADMSSGNQPSLGWMMTHLLLRCSKTTRMSPGIEKV 60 70 80 90 100 110 110 120 130 140 150 160 eh0072 DDVVKRLLSLEMANKKEMLKIKQEQFMKKIVANPEDTRSLEARSWSGTLGEWGCREAC ::::::::::::::::::::::::::::::::::::::::::: gi|120 DDVVKRLLSLEMANKKEMLKIKQEQFMKKIVANPEDTRSLEARIIALSVKIRSYEEHLEK 120 130 140 150 160 170 gi|120 HRKDKAHKRYLLMSIDQRKKMLKNLRNTNYDVFEKICWGLGIEYTFPPLYYRRAHRRFVT 180 190 200 210 220 230 >>gi|148698367|gb|EDL30314.1| mitochondrial ribosomal pr (136 aa) initn: 552 init1: 273 opt: 513 Z-score: 633.7 bits: 123.3 E(): 1.2e-26 Smith-Waterman score: 513; 60.432% identity (85.612% similar) in 139 aa overlap (22-158:1-136) 10 20 30 40 50 60 eh0072 PVRAAKRTAAAARAPRGLEVTMLRVAWRTLSLIRTRAVTQVLVPGLPGGGSAKFPFNQWG :::.:::.:: .:..:::.. ::.: ::.::.. . : gi|148 MLRAAWRALSSVRAQAVTRAPVPALRGGSSASLLSARCG 10 20 30 70 80 90 100 110 eh0072 LQPRSLLLQAARGYV-VRKPAQSRLDDDPPPSTLLKDYQNV-PGIEKVDDVVKRLLSLEM ::: ::: .:::.:. :.::.: . ::.:: :...:.:... :.::::::::::.::::: gi|148 LQPPSLL-RAARAYAAVQKPVQPKQDDEPPSSAFIKEYKDIIPNIEKVDDVVKRILSLEM 40 50 60 70 80 90 120 130 140 150 160 eh0072 ANKKEMLKIKQEQFMKKIVANPEDTRSLEARSWSGTLGEWGCREAC :..:: :::::::.:.::: ::::.:.:::.::: : : gi|148 ASRKEKLKIKQEQLMNKIVENPEDSRTLEAQSWSR--GPW 100 110 120 130 >>gi|149023938|gb|EDL80435.1| mitochondrial ribosomal pr (135 aa) initn: 387 init1: 387 opt: 503 Z-score: 621.6 bits: 121.0 E(): 5.5e-26 Smith-Waterman score: 503; 59.420% identity (84.783% similar) in 138 aa overlap (22-158:1-135) 10 20 30 40 50 60 eh0072 PVRAAKRTAAAARAPRGLEVTMLRVAWRTLSLIRTRAVTQVLVPGLPGGGSAKFPFNQWG :::.:::.:: .:..::::. ::.: . .::..: . : gi|149 MLRAAWRALSSVRVQAVTQAPVPALRARSSASLPSARCG 10 20 30 70 80 90 100 110 eh0072 LQPRSLLLQAARGYVVRKPAQSRLDDDPPPSTLLKDYQNV-PGIEKVDDVVKRLLSLEMA :: ::: .:::.:.:.::.:.. ::.: ::..:.:.:. :..:::::::::.:::::: gi|149 LQTPSLL-NAARAYAVQKPVQAKQDDEPASSTFIKEYKNIIPNMEKVDDVVKRILSLEMA 40 50 60 70 80 90 120 130 140 150 160 eh0072 NKKEMLKIKQEQFMKKIVANPEDTRSLEARSWSGTLGEWGCREAC ..:: ::::.::.:.::. :::: :.::::::: : : gi|149 SRKEKLKIKREQLMNKIAENPEDYRTLEARSWSR--GSW 100 110 120 130 164 residues in 1 query sequences 2355815254 residues in 6825066 library sequences Tcomplib [34.26] (8 proc) start: Fri Aug 8 17:13:22 2008 done: Fri Aug 8 17:16:29 2008 Total Scan time: 600.990 Total Display time: 0.020 Function used was FASTA [version 34.26.5 April 26, 2007]