FASTA searches a protein or DNA sequence data bank version 3.4t11 Apr 17, 2002 Please cite: W.R. Pearson & D.J. Lipman PNAS (1988) 85:2444-2448 ../query/mFLJ00017.ptfa, 408 aa vs ./tmplib.26680 library 1768609 residues in 2168 sequences Expectation_n fit: rho(ln(x))= 4.2909+/-0.00454; mu= 16.8664+/- 0.305 mean_var=104.4461+/-23.822, 0's: 0 Z-trim: 17 B-trim: 0 in 0/36 Lambda= 0.1255 FASTA (3.45 Mar 2002) function [optimized, BL50 matrix (15:-5)] ktup: 2 join: 37, opt: 25, open/ext: -10/-2, width: 16 The best scores are: opt bits E(2168) mKIAA4241 ( 453 res) mfj17221 ( 453) 1957 365 5.9e-102 mKIAA4240 ( 416 res) mia02001 ( 416) 1893 353 1.7e-98 mKIAA0763 ( 530 res) mbp04234 ( 530) 501 101 1.5e-22 mKIAA1110 ( 911 res) mbg05269 ( 911) 487 99 1.1e-21 mKIAA0248 ( 1803 res) mbg12530 (1803) 471 97 1.2e-20 mKIAA2011 ( 771 res) mfj31262 ( 771) 322 69 1e-12 mKIAA0522 ( 646 res) mbp09103 ( 646) 280 61 1.8e-10 mKIAA0969 ( 1106 res) mic02056 (1106) 190 46 1.9e-05 mKIAA1200 ( 1447 res) mbg09531 (1447) 165 41 0.00052 mKIAA1686 ( 1205 res) mbg21605 (1205) 147 38 0.0045 mKIAA0902 ( 1058 res) mbg09814 (1058) 146 38 0.0047 mKIAA0041 ( 807 res) mbg03660 ( 807) 144 37 0.0053 >>mKIAA4241 ( 453 res) mfj17221 (453 aa) initn: 1971 init1: 1328 opt: 1957 Z-score: 1922.7 bits: 364.8 E(): 5.9e-102 Smith-Waterman score: 1957; 71.835% identity (72.021% ungapped) in 387 aa overlap (23-408:66-452) 10 20 30 40 50 mFLJ00 AGKASSWATFPLGFRMDVCHTDPAELSSGEAKELQQIKWHRKQLLEDIQKLK : .:: : .:: .:. ..:.:..::..:: mKIAA4 PSRSRRLRARRLAGSGLKMDEGGGGEGGSVPEDLSLEEREELLDIRRRKKELIDDIERLK 40 50 60 70 80 90 60 70 80 90 100 110 mFLJ00 DEIADVFAQIDCFESTEESRMAQKEKEMCIGRKKFNMDPNKGIQYLIEHKLLTSDVQDIA :::.:...:: . :.:::. .:..:.. .:::::::::.::::.:::. :: :. .:.: mKIAA4 YEIAEVMTEIDNLTSVEESKTTQRNKQIAMGRKKFNMDPKKGIQFLIENDLLQSSPEDVA 100 110 120 130 140 150 120 130 140 150 160 170 mFLJ00 QFLYKGDGLNKTAIGTYLGEKDPINLQVLQAFVDCHEFANLNLVQALRQFLWSFRLPGEA ::::::.:::::.:: ::::.: .:..::::::. ::::.:::::::::::::::::::: mKIAA4 QFLYKGEGLNKTVIGDYLGERDDFNIKVLQAFVELHEFADLNLVQALRQFLWSFRLPGEA 160 170 180 190 200 210 180 190 200 210 220 230 mFLJ00 QKIDRMMEAFAARYCLCNPGVFRSTDTCYVLSFSVIMLNTGLHNPNVRDRPPFERFVTMN :::::::::::.::::::::::.::::::::::..:::::.::: ::::.: :::.::: mKIAA4 QKIDRMMEAFASRYCLCNPGVFQSTDTCYVLSFAIIMLNTSLHNHNVRDKPTAERFITMN 220 230 240 250 260 270 240 250 260 270 280 290 mFLJ00 RGINSGSDLPEEQLRNLFDSIKSEPFSIPEDDGGDLTHTFFNPDREGWLLKLGG-RVKTW :::: :.::::: ::::..:::.:::.::::::.:::::::::::::::::::: ::::: mKIAA4 RGINEGGDLPEELLRNLYESIKNEPFKIPEDDGNDLTHTFFNPDREGWLLKLGGGRVKTW 280 290 300 310 320 330 300 310 320 330 340 350 mFLJ00 KRRWFILTDNCLYYFEFTTDKEPRGIIPLENLSVQKVEDPKKPFCLELYNPSCRGQKIKA ::::::::::::::::.::::::::::::::::...::::.:: :.:::::: .:: ::: mKIAA4 KRRWFILTDNCLYYFEYTTDKEPRGIIPLENLSIREVEDPRKPNCFELYNPSHKGQVIKA 340 350 360 370 380 390 360 370 380 390 400 mFLJ00 CKTDGDGKVVEGKHESYRISAANAEERDQWIEAIRASITRVPFYDLLSARKKKIVGK :::..::.::::.: ::::: . ::...:...:.:::.: ::::.:..::..:..: mKIAA4 CKTEADGRVVEGNHVVYRISAPSPEEKEEWMKSIKASISRDPFYDMLATRKRRIANKK 400 410 420 430 440 450 >>mKIAA4240 ( 416 res) mia02001 (416 aa) initn: 1908 init1: 1302 opt: 1893 Z-score: 1860.4 bits: 353.2 E(): 1.7e-98 Smith-Waterman score: 1893; 70.312% identity (70.496% ungapped) in 384 aa overlap (23-405:27-410) 10 20 30 40 50 mFLJ00 AGKASSWATFPLGFRMDVCHTDPAELSSGEAKELQQIKWHRKQLLEDIQKLKDEIA :..:.. : .::..:. ....:: :::.::.::: mKIAA4 RGLREVSGGRGPEPVASTMEDDDSYVPSDLTAEERQELENIRRRKQELLADIQRLKEEIA 10 20 30 40 50 60 60 70 80 90 100 110 mFLJ00 DVFAQIDCFESTEESRMAQKEKEMCIGRKKFNMDPNKGIQYLIEHKLLTSDVQDIAQFLY .: .:. . :::: . :..:.. .:::::::::.::::.:::. :: . .::::::: mKIAA4 EVANEIESLGSTEERKNMQRNKQVAMGRKKFNMDPKKGIQFLIENGLLKNTCEDIAQFLY 70 80 90 100 110 120 120 130 140 150 160 170 mFLJ00 KGDGLNKTAIGTYLGEKDPINLQVLQAFVDCHEFANLNLVQALRQFLWSFRLPGEAQKID ::.:::::::: ::::.: ...::: :::. :::..:::::::::::::::::::::::: mKIAA4 KGEGLNKTAIGDYLGERDEFSIQVLYAFVELHEFTDLNLVQALRQFLWSFRLPGEAQKID 130 140 150 160 170 180 180 190 200 210 220 230 mFLJ00 RMMEAFAARYCLCNPGVFRSTDTCYVLSFSVIMLNTGLHNPNVRDRPPFERFVTMNRGIN ::::::: ::: :: :::.::::::::::..:::::.::::::.:.: :::..:::::: mKIAA4 RMMEAFAQRYCQCNTGVFQSTDTCYVLSFAIIMLNTSLHNPNVKDKPTVERFIAMNRGIN 190 200 210 220 230 240 240 250 260 270 280 290 mFLJ00 SGSDLPEEQLRNLFDSIKSEPFSIPEDDGGDLTHTFFNPDREGWLLKLGG-RVKTWKRRW .:.::::: ::::..:::.:::.::::::.:::::::::::::::::::: ::::::::: mKIAA4 DGGDLPEELLRNLYESIKNEPFKIPEDDGNDLTHTFFNPDREGWLLKLGGGRVKTWKRRW 250 260 270 280 290 300 300 310 320 330 340 350 mFLJ00 FILTDNCLYYFEFTTDKEPRGIIPLENLSVQKVEDPKKPFCLELYNPSCRGQKIKACKTD ::::::::::::.::::::::::::::::...::: ::: :.::: :. . : ::::::. mKIAA4 FILTDNCLYYFEYTTDKEPRGIIPLENLSIREVEDSKKPNCFELYIPDNKDQVIKACKTE 310 320 330 340 350 360 360 370 380 390 400 mFLJ00 GDGKVVEGKHESYRISAANAEERDQWIEAIRASITRVPFYDLLSARKKKIVGK .::.::::.: ::::: . ::...::. :.:.:.: :::..:.:::::. mKIAA4 ADGRVVEGNHTVYRISAPTPEEKEDWIKCIKAAISRDPFYEMLAARKKKVSSTKRH 370 380 390 400 410 >>mKIAA0763 ( 530 res) mbp04234 (530 aa) initn: 458 init1: 212 opt: 501 Z-score: 497.4 bits: 101.3 E(): 1.5e-22 Smith-Waterman score: 501; 38.164% identity (39.303% ungapped) in 207 aa overlap (63-263:78-284) 40 50 60 70 80 90 mFLJ00 ELQQIKWHRKQLLEDIQKLKDEIADVFAQIDCFESTEESRMAQKEKEMCIGRKKFNMDPN . ..: : . ..... :: . :: :. mKIAA0 DTINCSSESSSRDSLREQTLSKQTYHKETRNSWDSPAFSNDVIRKRHYRIGLNLFNKKPE 50 60 70 80 90 100 100 110 120 130 140 150 mFLJ00 KGIQYLIEHKLLTSDVQDIAQFLYKGDGLNKTAIGTYLGEKDP-INLQVLQAFVDCHEFA ::::::::. .. . .:.:: . ::.. :: .::... .: .::. :: .:. mKIAA0 KGIQYLIERGFVPDTPVGVAHFLLQRKGLSRQMIGEFLGNRQKQFNRDVLDCVVDEMDFS 110 120 130 140 150 160 160 170 180 190 200 mFLJ00 NLNLVQALRQFLWSFRLPGEAQKIDRMMEAFAARYCLCNPGV---FRSTDTCYVLSFSVI ..: .:::.: .:. :::::..:..:::. :::.::::: ::. :: ..:.:..: mKIAA0 AMELDEALRKFQAHIRVQGEAQKVERLIEAFSQRYCVCNPGVVRQFRNPDTIFILAFAII 170 180 190 200 210 220 210 220 230 240 250 260 mFLJ00 MLNTGLHNPNVRD--RPPFERFVTMNRGINSGSDLPEEQLRNLFDSIKSEPFSIPEDDGG .::: ...:::. . .: :: ::...: :.:.: : .... :... .. :: mKIAA0 LLNTDMYSPNVKPERKMKLEDFVKNLRGVDDGEDIPRETLIGIYERIRKRELKTNEDHVS 230 240 250 260 270 280 270 280 290 300 310 320 mFLJ00 DLTHTFFNPDREGWLLKLGGRVKTWKRRWFILTDNCLYYFEFTTDKEPRGIIPLENLSVQ mKIAA0 QVQKVEKLIVGKKPIGSLHHGLGCVLSLPHRRLVCYCRLFQVPDPNKPQKLGLHQREIFL 290 300 310 320 330 340 >>mKIAA1110 ( 911 res) mbg05269 (911 aa) initn: 398 init1: 201 opt: 487 Z-score: 481.4 bits: 99.2 E(): 1.1e-21 Smith-Waterman score: 487; 39.362% identity (40.659% ungapped) in 188 aa overlap (82-263:377-564) 60 70 80 90 100 110 mFLJ00 KDEIADVFAQIDCFESTEESRMAQKEKEMCIGRKKFNMDPNKGIQYLIEHKLLTSDVQDI :: . ::..:.::::.:: . .. . . mKIAA1 LSLPRYHCENPASCRSPTLSTDTLRKRLYRIGLNLFNINPDKGIQFLISRGFIPDTPIGV 350 360 370 380 390 400 120 130 140 150 160 170 mFLJ00 AQFLYKGDGLNKTAIGTYLGE-KDPINLQVLQAFVDCHEFANLNLVQALRQFLWSFRLPG :.:: . ::.. :: .::. : .: .::. :: .:.:..: .:::.: .:. : mKIAA1 AHFLLQRKGLSRQMIGEFLGNSKKQFNRDVLDCVVDEMDFSNMELDEALRKFQAHIRVQG 410 420 430 440 450 460 180 190 200 210 220 mFLJ00 EAQKIDRMMEAFAARYCLCNPGV---FRSTDTCYVLSFSVIMLNTGLHNPNVR-DRPP-F ::::..:..:::. :::.::: : :.. :: ..:.:..:.::: ...::.. :: . mKIAA1 EAQKVERLIEAFSQRYCMCNPEVVQQFHNPDTIFILAFAIILLNTDMYSPNIKPDRKMML 470 480 490 500 510 520 230 240 250 260 270 280 mFLJ00 ERFVTMNRGINSGSDLPEEQLRNLFDSIKSEPFSIPEDDGGDLTHTFFNPDREGWLLKLG : :. ::...:.:.:.: . .... :... .. :: mKIAA1 EDFIRNLRGVDDGADIPRELVVGIYERIQQKELKSNEDHVTYVTKVEKSIVGMKTVLSMP 530 540 550 560 570 580 290 300 310 320 330 340 mFLJ00 GRVKTWKRRWFILTDNCLYYFEFTTDKEPRGIIPLENLSVQKVEDPKKPFCLELYNPSCR mKIAA1 HRRLVCCSRLFEVTDVNKLQKQAAHQREVFLFNDLLVILKLCPKKKSSFTYTFCKAVGLL 590 600 610 620 630 640 >>mKIAA0248 ( 1803 res) mbg12530 (1803 aa) initn: 315 init1: 315 opt: 471 Z-score: 462.9 bits: 96.7 E(): 1.2e-20 Smith-Waterman score: 471; 35.565% identity (37.778% ungapped) in 239 aa overlap (40-266:602-838) 10 20 30 40 50 60 mFLJ00 FPLGFRMDVCHTDPAELSSGEAKELQQIKWHRKQLLEDIQKLKDEIADV-FAQ----IDC : : .:... : :: :.. ..: mKIAA0 ERATSDGKSTGVALDARGLHFPSGGWLSTEHGKPGCRDLEEAGDSGADKKFTRKPPRFSC 580 590 600 610 620 630 70 80 90 100 110 120 mFLJ00 F--ESTEESRMAQKEKEMCIGRKKFNMDPNKGIQYLIEHKLLT--SDVQDIAQFLYKGDG . . : .. .:.: . : ..::. :.::::.: :. ::: : ..::.: .. mKIAA0 LLPDPRELIEIKNKKKLLITGTEQFNQKPKKGIQFLQEKGLLTIPMDNTEVAQWLRENPR 640 650 660 670 680 690 130 140 150 160 170 180 mFLJ00 LNKTAIGTYLGEKDPINLQVLQAFVDCHEFANLNLVQALRQFLWSFRLPGEAQKIDRMME :.: :: ..... :...:..::. : .: : .::: .: .::::::: : :..: mKIAA0 LDKKMIGEFVSDRK--NMDLLESFVSTFSFQGLRLDEALRLYLEAFRLPGEAPVIHRLLE 700 710 720 730 740 190 200 210 220 230 mFLJ00 AFAARYCLCNPGVFRSTDTCYVLSFSVIMLNTGLHNPNVRDR--P-PFERFVTMNRGINS .:. .. :: . : ..:.:..:...:::::: :: ::: . : .:.: .:.:. mKIAA0 VFTEHWRSCNGSPFANSDACFALAYAVIMLNTDQHNHNVRKQNAPMTLEEFRKNLKGVNG 750 760 770 780 790 800 240 250 260 270 280 290 mFLJ00 GSDLPEEQLRNLFDSIKSEPFSIPEDDGGDLTHTFFNPDREGWLLKLGGRVKTWKRRWFI :.:. .. :.... .::.: . .::.. : mKIAA0 GKDFEQDILEDMYHAIKNEEIVMPEEQTGLVRENYVWSVLLHRGASPEGVFLRVPPGSYD 810 820 830 840 850 860 >>mKIAA2011 ( 771 res) mfj31262 (771 aa) initn: 309 init1: 309 opt: 322 Z-score: 320.7 bits: 69.2 E(): 1e-12 Smith-Waterman score: 322; 33.974% identity (33.974% ungapped) in 156 aa overlap (109-264:306-461) 80 90 100 110 120 130 mFLJ00 EMCIGRKKFNMDPNKGIQYLIEHKLLTSDVQDIAQFLYKGDGLNKTAIGTYLGEKDPINL : .:. ::. ::. :. .. .::... .. mKIAA2 SDSELDSTERLALGSTDTLSNGQKADLEAAQRLAKRLYRLDGFRKADVARHLGKNNDFSK 280 290 300 310 320 330 140 150 160 170 180 190 mFLJ00 QVLQAFVDCHEFANLNLVQALRQFLWSFRLPGEAQKIDRMMEAFAARYCLCNPGVFRSTD : .. :....: :::: :: . : ::.:. .:.. :. :: ::: .. : : mKIAA2 LVAGEYLKFFVFTGMTLDQALRVFLKELALMGETQERERVLAHFSQRYFQCNPEALSSED 340 350 360 370 380 390 200 210 220 230 240 250 mFLJ00 TCYVLSFSVIMLNTGLHNPNVRDRPPFERFVTMNRGINSGSDLPEEQLRNLFDSIKSEPF ..:. ....::: ::. :. : :. .:.:.:.:.:.: :. :..:::.: . mKIAA2 GAHTLTCALMLLNTDLHGHNIGKRMTCGDFIGNLEGLNDGGDFPRELLKALYSSIKNEKL 400 410 420 430 440 450 260 270 280 290 300 310 mFLJ00 SIPEDDGGDLTHTFFNPDREGWLLKLGGRVKTWKRRWFILTDNCLYYFEFTTDKEPRGII . :. mKIAA2 QWAIDEEELRRSLSELADPNPKVIKRVSGGSGSSSSPFLDLTPEPGAAVYKHGALVRKVH 460 470 480 490 500 510 >>mKIAA0522 ( 646 res) mbp09103 (646 aa) initn: 223 init1: 135 opt: 280 Z-score: 280.3 bits: 61.5 E(): 1.8e-10 Smith-Waterman score: 280; 36.275% identity (38.144% ungapped) in 102 aa overlap (167-263:3-104) 140 150 160 170 180 190 mFLJ00 NLQVLQAFVDCHEFANLNLVQALRQFLWSFRLPGEAQKIDRMMEAFAARYCLCNPGV--- :. :::::..:..:::. :::.:::.. mKIAA0 LIRVQGEAQKVERLIEAFSQRYCVCNPALVRQ 10 20 30 200 210 220 230 240 250 mFLJ00 FRSTDTCYVLSFSVIMLNTGLHNPNVRD--RPPFERFVTMNRGINSGSDLPEEQLRNLFD ::. :: ..:.:..:.::: ...:.:. . .. :. ::...: :.:.. : .... mKIAA0 FRNPDTIFILAFAIILLNTDMYSPSVKAERKMKLDDFIKNLRGVDNGEDIPRDLLVGIYQ 40 50 60 70 80 90 260 270 280 290 300 310 mFLJ00 SIKSEPFSIPEDDGGDLTHTFFNPDREGWLLKLGGRVKTWKRRWFILTDNCLYYFEFTTD :... . .: mKIAA0 RIQGRELRTNDDHVSQVQAVERMIVGKKPVLSLPHRRLVCCCQLYEVPDPNRPQRLGLHQ 100 110 120 130 140 150 >>mKIAA0969 ( 1106 res) mic02056 (1106 aa) initn: 167 init1: 129 opt: 190 Z-score: 190.0 bits: 45.5 E(): 1.9e-05 Smith-Waterman score: 190; 27.568% identity (32.278% ungapped) in 185 aa overlap (226-393:45-219) 200 210 220 230 240 250 mFLJ00 STDTCYVLSFSVIMLNTGLHNPNVRDRPPFERFVTMNRGI---NSGSDLPEEQ--LRNLF .: .:.: : : :..: :. .: mKIAA0 RHPVTGQIPEENKKFDLKTSALDMSNKTGGKRSATINSDIANHNMVSEVPPERPNIRATR 20 30 40 50 60 70 260 270 280 290 300 mFLJ00 DSIKSEPFSIPEDDGGDLTHTF-FNPD----REGWLLKLGGR-VKTWKRRWFILTDNCLY : :. : : .:.. ::. . ::: : .. :: :..:::.:.: ::. mKIAA0 TSRKAIAF-------GKRAHSMKRNPNAPVTKAGWLYKQASSGVKQWNKRWFVLVDRCLF 80 90 100 110 120 310 320 330 340 350 mFLJ00 YFEFTTDKEPRGIIPLENLSVQKVEDP-----KKPFCLELYNPSCRGQKIKACKTDGDGK :.. .. : ::: .. : :. :. : . .. . : . : . . mKIAA0 YYKDEKQESILGSIPLLSFRVAAVQPSDNISRKHTFKVTVHWVDEAGASSTHCLSP---Q 130 140 150 160 170 180 360 370 380 390 400 mFLJ00 VVEGKHESYRISAANAEERDQWIEAI-RASITRVPFYDLLSARKKKIVGK . .. ..: .:: . ::.. ::.:. .:. ...: mKIAA0 AEHAGVRTYFFSAESPEEQEAWIQAMGEAARVQIPPAQKSVPQPVRHSLEKPDSENIPPS 190 200 210 220 230 240 mKIAA0 KHHQQPPHNNLTKLEPEAKTRGEGDGRGCEKAERRPERPEVKKETLVKANGLPSGPETAS 250 260 270 280 290 300 >>mKIAA1200 ( 1447 res) mbg09531 (1447 aa) initn: 147 init1: 147 opt: 165 Z-score: 164.4 bits: 41.2 E(): 0.00052 Smith-Waterman score: 165; 41.538% identity (44.262% ungapped) in 65 aa overlap (276-337:665-728) 250 260 270 280 290 300 mFLJ00 LRNLFDSIKSEPFSIPEDDGGDLTHTFFNPDREGWLLKLGGRVKTWKRRWFILTDNCLYY .. :.:::.:.::::::::::.: .. . : mKIAA1 PDACSLDSDYSEPEHKLQRTSSYSTDGEALEKSGYLLKMGSRVKTWKRRWFVLRQGQILY 640 650 660 670 680 690 310 320 330 340 350 360 mFLJ00 FEFTTD--KEPRGIIPLENLSVQKV-EDPKKPFCLELYNPSCRGQKIKACKTDGDGKVVE .. .: ..:.: . : : : : :. . : : mKIAA1 YKSPSDVIRKPQGQVDL-NSHCQIVREEEAQTFQLISGNKTYYLTAESPSLLEEWIRVLQ 700 710 720 730 740 750 370 380 390 400 mFLJ00 GKHESYRISAANAEERDQWIEAIRASITRVPFYDLLSARKKKIVGK mKIAA1 SLLKVQVTGPPALHQGGTKPTVKGWLTKVKHGHSKLVWCALVGKTFYYYRSHEDKRPLGC 760 770 780 790 800 810 >>mKIAA1686 ( 1205 res) mbg21605 (1205 aa) initn: 112 init1: 112 opt: 147 Z-score: 147.5 bits: 37.8 E(): 0.0045 Smith-Waterman score: 147; 31.325% identity (37.681% ungapped) in 83 aa overlap (274-345:102-181) 250 260 270 280 290 mFLJ00 EQLRNLFDSIKSEPFSIPEDDGGDLTHTFFNPD----REGWLLKLGGR-VKTWKRRWFIL ::. :.::: : . .: ::.:::.: mKIAA1 AVHPMSPVGRTSRASKKVHNFGKRSNSIKRNPNAPVVRRGWLYKQDSTGMKLWKKRWFVL 80 90 100 110 120 130 300 310 320 330 340 350 mFLJ00 TDNCLYYFEFTTDKEPRGII-----PLENLSVQKVEDP-KKPFCLELYNPSCRGQKIKAC .: ::.:.. :.. .::. : .... .:: .. . .. .:. : mKIAA1 SDLCLFYYR---DEKEEGILGSILLPSFQIAMLTAEDHINRKYAFKAAHPNMRTYYFCTD 140 150 160 170 180 360 370 380 390 400 mFLJ00 KTDGDGKVVEGKHESYRISAANAEERDQWIEAIRASITRVPFYDLLSARKKKIVGK mKIAA1 TGKEMELWMKAMLDAALVQTEPVKRITFNFRVDKITTDSASTKETNNIPNHRVLIRPEVQ 190 200 210 220 230 240 408 residues in 1 query sequences 1768609 residues in 2168 library sequences Scomplib [34t11] start: Mon Mar 27 10:06:03 2006 done: Mon Mar 27 10:06:04 2006 Scan time: 0.600 Display time: 0.130 Function used was FASTA [version 3.4t11 Apr 17, 2002]