FASTA searches a protein or DNA sequence data bank 36.3.4 Apr, 2011 Please cite: W.R. Pearson & D.J. Lipman PNAS (1988) 85:2444-2448 Query: pF1KE2419, 690 aa 1>>>pF1KE2419 690 - 690 aa - 690 aa Library: human.CCDS.faa 18511270 residues in 32554 sequences Statistics: Expectation_n fit: rho(ln(x))= 5.3789+/-0.00108; mu= 20.9321+/- 0.064 mean_var=98.0598+/-19.605, 0's: 0 Z-trim(106.3): 69 B-trim: 0 in 0/52 Lambda= 0.129518 statistics sampled from 8825 (8891) to 8825 sequences Algorithm: FASTA (3.7 Nov 2010) [optimized] Parameters: BL50 matrix (15:-5), open/ext: -10/-2 ktup: 2, E-join: 1 (0.617), E-opt: 0.2 (0.273), width: 16 Scan time: 2.330 The best scores are: opt bits E(32554) CCDS3435.1 SLC34A2 gene_id:10568|Hs108|chr4 ( 690) 4602 871.0 0 CCDS54750.1 SLC34A2 gene_id:10568|Hs108|chr4 ( 689) 4579 866.7 0 CCDS4418.1 SLC34A1 gene_id:6569|Hs108|chr5 ( 639) 1257 245.9 1.3e-64 CCDS7038.1 SLC34A3 gene_id:142680|Hs108|chr9 ( 599) 1102 216.9 6.7e-56 CCDS54953.1 SLC34A1 gene_id:6569|Hs108|chr5 ( 340) 958 189.8 5.7e-48 >>CCDS3435.1 SLC34A2 gene_id:10568|Hs108|chr4 (690 aa) initn: 4602 init1: 4602 opt: 4602 Z-score: 4650.8 bits: 871.0 E(32554): 0 Smith-Waterman score: 4602; 99.9% identity (99.9% similar) in 690 aa overlap (1-690:1-690) 10 20 30 40 50 60 pF1KE2 MAPWPELGDAQPNPDKYLEGAAGQQPTAPDKSKETNKTDNTEAPVTKIELLPSYSTATLI :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: CCDS34 MAPWPELGDAQPNPDKYLEGAAGQQPTAPDKSKETNKTDNTEAPVTKIELLPSYSTATLI 10 20 30 40 50 60 70 80 90 100 110 120 pF1KE2 DEPTEVDDPWNLPTLQDSGIKWSERDTKGKILCFFQGIGRLILLLGFLYFFVCSLDILSS :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: CCDS34 DEPTEVDDPWNLPTLQDSGIKWSERDTKGKILCFFQGIGRLILLLGFLYFFVCSLDILSS 70 80 90 100 110 120 130 140 150 160 170 180 pF1KE2 AFQLVGGKMAGQFFSNSSIMSNPLLGLVIGVLVTVLVQSSSTSTSIVVSMVSSSLLTVRA :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: CCDS34 AFQLVGGKMAGQFFSNSSIMSNPLLGLVIGVLVTVLVQSSSTSTSIVVSMVSSSLLTVRA 130 140 150 160 170 180 190 200 210 220 230 240 pF1KE2 AIPIIMGANIGTSITNTIVALMQVGDRSEFRRAFAGATVHDFFNWLSVLVLLPVEVATHY :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: CCDS34 AIPIIMGANIGTSITNTIVALMQVGDRSEFRRAFAGATVHDFFNWLSVLVLLPVEVATHY 190 200 210 220 230 240 250 260 270 280 290 300 pF1KE2 LEIITQLIVESFHFKNGEDAPDLLKVITKPFTKLIVQLDKKVISQIAMNDEKAKNKSLVK :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: CCDS34 LEIITQLIVESFHFKNGEDAPDLLKVITKPFTKLIVQLDKKVISQIAMNDEKAKNKSLVK 250 260 270 280 290 300 310 320 330 340 350 360 pF1KE2 IWCKTFTNKTQINVTVPSTANCTSPSLCWTDGIQNWTMKNVTYKENIAKCQHIFVNFHLP :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: CCDS34 IWCKTFTNKTQINVTVPSTANCTSPSLCWTDGIQNWTMKNVTYKENIAKCQHIFVNFHLP 310 320 330 340 350 360 370 380 390 400 410 420 pF1KE2 DLAVGTILLILSLLVLCGCLIMIVKILGSVLKGQVATVIKKTINTDFPFPFAWLTGYLAI :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: CCDS34 DLAVGTILLILSLLVLCGCLIMIVKILGSVLKGQVATVIKKTINTDFPFPFAWLTGYLAI 370 380 390 400 410 420 430 440 450 460 470 480 pF1KE2 LVGAGMTFIVQSSSVFTSALTPLIGIGVITIERAYPLTLGSNIGTTTTAILAALASPGNA :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: CCDS34 LVGAGMTFIVQSSSVFTSALTPLIGIGVITIERAYPLTLGSNIGTTTTAILAALASPGNA 430 440 450 460 470 480 490 500 510 520 530 540 pF1KE2 LRSSLQIALCHFFFNISGILLWYPIPFTRLPIRMAKGLGNISAKYRWFAVFYLIIFFFLI :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: CCDS34 LRSSLQIALCHFFFNISGILLWYPIPFTRLPIRMAKGLGNISAKYRWFAVFYLIIFFFLI 490 500 510 520 530 540 550 560 570 580 590 600 pF1KE2 PLTVFGLSLAGWRVLVGVGVPVVFIIILVLCLRLLQSRCPRVLPKKLQNWNFLPLWMRSL :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: CCDS34 PLTVFGLSLAGWRVLVGVGVPVVFIIILVLCLRLLQSRCPRVLPKKLQNWNFLPLWMRSL 550 560 570 580 590 600 610 620 630 640 650 660 pF1KE2 KPWDAVVSKFTGCFQMRCCCCCRVCCRACCLLCGCPKCCRCSKCCEDLEEAQEGQDVPVK ::::::::::::::::::::::::::::::::: :::::::::::::::::::::::::: CCDS34 KPWDAVVSKFTGCFQMRCCCCCRVCCRACCLLCDCPKCCRCSKCCEDLEEAQEGQDVPVK 610 620 630 640 650 660 670 680 690 pF1KE2 APETFDNITISREAQGEVPASDSKTECTAL :::::::::::::::::::::::::::::: CCDS34 APETFDNITISREAQGEVPASDSKTECTAL 670 680 690 >>CCDS54750.1 SLC34A2 gene_id:10568|Hs108|chr4 (689 aa) initn: 4326 init1: 4326 opt: 4579 Z-score: 4627.6 bits: 866.7 E(32554): 0 Smith-Waterman score: 4579; 99.6% identity (99.7% similar) in 690 aa overlap (1-690:1-689) 10 20 30 40 50 60 pF1KE2 MAPWPELGDAQPNPDKYLEGAAGQQPTAPDKSKETNKTDNTEAPVTKIELLPSYSTATLI ::::::::::::::::::::::::::::::::::::: .::::::::::::::::::::: CCDS54 MAPWPELGDAQPNPDKYLEGAAGQQPTAPDKSKETNK-NNTEAPVTKIELLPSYSTATLI 10 20 30 40 50 70 80 90 100 110 120 pF1KE2 DEPTEVDDPWNLPTLQDSGIKWSERDTKGKILCFFQGIGRLILLLGFLYFFVCSLDILSS :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: CCDS54 DEPTEVDDPWNLPTLQDSGIKWSERDTKGKILCFFQGIGRLILLLGFLYFFVCSLDILSS 60 70 80 90 100 110 130 140 150 160 170 180 pF1KE2 AFQLVGGKMAGQFFSNSSIMSNPLLGLVIGVLVTVLVQSSSTSTSIVVSMVSSSLLTVRA :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: CCDS54 AFQLVGGKMAGQFFSNSSIMSNPLLGLVIGVLVTVLVQSSSTSTSIVVSMVSSSLLTVRA 120 130 140 150 160 170 190 200 210 220 230 240 pF1KE2 AIPIIMGANIGTSITNTIVALMQVGDRSEFRRAFAGATVHDFFNWLSVLVLLPVEVATHY :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: CCDS54 AIPIIMGANIGTSITNTIVALMQVGDRSEFRRAFAGATVHDFFNWLSVLVLLPVEVATHY 180 190 200 210 220 230 250 260 270 280 290 300 pF1KE2 LEIITQLIVESFHFKNGEDAPDLLKVITKPFTKLIVQLDKKVISQIAMNDEKAKNKSLVK :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: CCDS54 LEIITQLIVESFHFKNGEDAPDLLKVITKPFTKLIVQLDKKVISQIAMNDEKAKNKSLVK 240 250 260 270 280 290 310 320 330 340 350 360 pF1KE2 IWCKTFTNKTQINVTVPSTANCTSPSLCWTDGIQNWTMKNVTYKENIAKCQHIFVNFHLP :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: CCDS54 IWCKTFTNKTQINVTVPSTANCTSPSLCWTDGIQNWTMKNVTYKENIAKCQHIFVNFHLP 300 310 320 330 340 350 370 380 390 400 410 420 pF1KE2 DLAVGTILLILSLLVLCGCLIMIVKILGSVLKGQVATVIKKTINTDFPFPFAWLTGYLAI :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: CCDS54 DLAVGTILLILSLLVLCGCLIMIVKILGSVLKGQVATVIKKTINTDFPFPFAWLTGYLAI 360 370 380 390 400 410 430 440 450 460 470 480 pF1KE2 LVGAGMTFIVQSSSVFTSALTPLIGIGVITIERAYPLTLGSNIGTTTTAILAALASPGNA :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: CCDS54 LVGAGMTFIVQSSSVFTSALTPLIGIGVITIERAYPLTLGSNIGTTTTAILAALASPGNA 420 430 440 450 460 470 490 500 510 520 530 540 pF1KE2 LRSSLQIALCHFFFNISGILLWYPIPFTRLPIRMAKGLGNISAKYRWFAVFYLIIFFFLI :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: CCDS54 LRSSLQIALCHFFFNISGILLWYPIPFTRLPIRMAKGLGNISAKYRWFAVFYLIIFFFLI 480 490 500 510 520 530 550 560 570 580 590 600 pF1KE2 PLTVFGLSLAGWRVLVGVGVPVVFIIILVLCLRLLQSRCPRVLPKKLQNWNFLPLWMRSL :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: CCDS54 PLTVFGLSLAGWRVLVGVGVPVVFIIILVLCLRLLQSRCPRVLPKKLQNWNFLPLWMRSL 540 550 560 570 580 590 610 620 630 640 650 660 pF1KE2 KPWDAVVSKFTGCFQMRCCCCCRVCCRACCLLCGCPKCCRCSKCCEDLEEAQEGQDVPVK ::::::::::::::::::::::::::::::::: :::::::::::::::::::::::::: CCDS54 KPWDAVVSKFTGCFQMRCCCCCRVCCRACCLLCDCPKCCRCSKCCEDLEEAQEGQDVPVK 600 610 620 630 640 650 670 680 690 pF1KE2 APETFDNITISREAQGEVPASDSKTECTAL :::::::::::::::::::::::::::::: CCDS54 APETFDNITISREAQGEVPASDSKTECTAL 660 670 680 >>CCDS4418.1 SLC34A1 gene_id:6569|Hs108|chr5 (639 aa) initn: 2194 init1: 1250 opt: 1257 Z-score: 1273.3 bits: 245.9 E(32554): 1.3e-64 Smith-Waterman score: 2176; 64.8% identity (86.8% similar) in 514 aa overlap (100-613:103-599) 70 80 90 100 110 120 pF1KE2 WNLPTLQDSGIKWSERDTKGKILCFFQGIGRLILLLGFLYFFVCSLDILSSAFQLVGGKM .. :.: :::.::::::.:::::::.:::. CCDS44 PLPAKLALEEEQKPESRLVPKLRQAGAMLLKVPLMLTFLYLFVCSLDMLSSAFQLAGGKV 80 90 100 110 120 130 130 140 150 160 170 180 pF1KE2 AGQFFSNSSIMSNPLLGLVIGVLVTVLVQSSSTSTSIVVSMVSSSLLTVRAAIPIIMGAN ::..:....:.:::. :::.:.:::::::::::::::.::::::.:: : .:::::::.: CCDS44 AGDIFKDNAILSNPVAGLVVGILVTVLVQSSSTSTSIIVSMVSSGLLEVSSAIPIIMGSN 140 150 160 170 180 190 190 200 210 220 230 240 pF1KE2 IGTSITNTIVALMQVGDRSEFRRAFAGATVHDFFNWLSVLVLLPVEVATHYLEIITQLIV ::::.:::::::::.:::..:::::::::::: :::::::::::.:.:: ::. ::.:.: CCDS44 IGTSVTNTIVALMQAGDRTDFRRAFAGATVHDCFNWLSVLVLLPLEAATGYLHHITRLVV 200 210 220 230 240 250 250 260 270 280 290 300 pF1KE2 ESFHFKNGEDAPDLLKVITKPFTKLIVQLDKKVISQIAMNDEKAKNKSLVKIWCKTFTNK ::....:.:::::::.::.::::::.:::..::..:: .::. .:.::..:::. CCDS44 ASFNIHGGRDAPDLLKIITEPFTKLIIQLDESVITSIATGDESLRNHSLIQIWCH----- 260 270 280 290 300 310 320 330 340 350 360 pF1KE2 TQINVTVPSTANCTSPSLCWTDGIQNWTMKNVTYKENIAKCQHIFVNFHLPDLAVGTILL :.. . . :. ... .. :. :.: . ::.::::. ::::::: ::: CCDS44 -------PDSLQAPT-SMSRAEANSSQTLGNAT----MEKCNHIFVDTGLPDLAVGLILL 310 320 330 340 350 370 380 390 400 410 420 pF1KE2 ILSLLVLCGCLIMIVKILGSVLKGQVATVIKKTINTDFPFPFAWLTGYLAILVGAGMTFI ::..:: :::..::.:.:.:::::: ::.:.:::::: ::.:.:::.:..:::.:::. CCDS44 AGSLVLLCTCLILLVKMLNSLLKGQVAKVIQKVINTDFPAPFTWVTGYFAMVVGASMTFV 360 370 380 390 400 410 430 440 450 460 470 480 pF1KE2 VQSSSVFTSALTPLIGIGVITIERAYPLTLGSNIGTTTTAILAALASPGNALRSSLQIAL ::::::::::.:::::.:::.::::::::::::::::::::::::::: . : :..:::: CCDS44 VQSSSVFTSAITPLIGLGVISIERAYPLTLGSNIGTTTTAILAALASPREKLSSAFQIAL 420 430 440 450 460 470 490 500 510 520 530 540 pF1KE2 CHFFFNISGILLWYPIPFTRLPIRMAKGLGNISAKYRWFAVFYLIIFFFLIPLTVFGLSL :::::::::::::::.: :::::::::.::. .::::::::.::.. :.:.: :::.:. CCDS44 CHFFFNISGILLWYPVPCTRLPIRMAKALGKRTAKYRWFAVLYLLVCFLLLPSLVFGISM 480 490 500 510 520 530 550 560 570 580 590 600 pF1KE2 AGWRVLVGVGVPVVFIIILVLCLRLLQSRCPRVLPKKLQNWNFLPLWMRSLKPWDAVVSK :::.:.::::.: .. .:. . .:::: : ::: ::.:.::: ::.:::: : .... CCDS44 AGWQVMVGVGTPFGALLAFVVLINVLQSRSPGHLPKWLQTWDFLPRWMHSLKPLDHLITR 540 550 560 570 580 590 610 620 630 640 650 660 pF1KE2 FTGCFQMRCCCCCRVCCRACCLLCGCPKCCRCSKCCEDLEEAQEGQDVPVKAPETFDNIT : : CCDS44 ATLCCARPEPRSPPLPPRVFLEELPPATPSPRLALPAHHNATRL 600 610 620 630 >>CCDS7038.1 SLC34A3 gene_id:142680|Hs108|chr9 (599 aa) initn: 1979 init1: 1086 opt: 1102 Z-score: 1117.1 bits: 216.9 E(32554): 6.7e-56 Smith-Waterman score: 1945; 52.9% identity (78.1% similar) in 571 aa overlap (55-625:32-575) 30 40 50 60 70 80 pF1KE2 QPTAPDKSKETNKTDNTEAPVTKIELLPSYSTATLIDEPTEVDDPWNLPTLQDSGIKWSE :.: ...: :::.:: :.:.. :.: CCDS70 PSSLPGSQVPHPTLDAVDLVEKTLRNEGTSSSAPVLEEGD--TDPWTLPQLKDTSQPWKE 10 20 30 40 50 90 100 110 120 130 140 pF1KE2 RDTKGKILCFFQGIGRLILLLGFLYFFVCSLDILSSAFQLVGGKMAGQFFSNSSIMSNPL . :.. .. . ::: ::::.::::.:::::::.:.:.::..:... ..:::. CCDS70 LRVAGRLRRVAGSVLKACGLLGSLYFFICSLDVLSSAFQLLGSKVAGDIFKDNVVLSNPV 60 70 80 90 100 110 150 160 170 180 190 200 pF1KE2 LGLVIGVLVTVLVQSSSTSTSIVVSMVSSSLLTVRAAIPIIMGANIGTSITNTIVALMQV :::::::::.::::::::.:::::::...:::::...:::::.:.:::::.:.:.. : CCDS70 AGLVIGVLVTALVQSSSTSSSIVVSMVAAKLLTVRVSVPIIMGVNVGTSITSTLVSMAQS 120 130 140 150 160 170 210 220 230 240 250 260 pF1KE2 GDRSEFRRAFAGATVHDFFNWLSVLVLLPVEVATHYLEIITQLIVESFHFKNGEDAPDLL :::.::.:::.:..:: .::::.::::::.: :: :: ...: . . . .:::.: CCDS70 GDRDEFQRAFSGSAVHGIFNWLTVLVLLPLESATALLERLSELALGAASLTPRAQAPDIL 180 190 200 210 220 230 270 280 290 300 310 320 pF1KE2 KVITKPFTKLIVQLDKKVISQIAMNDEKAKNKSLVKIWCKTFTNKTQINVTVPSTANCTS ::.:::.:.::::::. .: . : .. : :.::.: :: : . :: : ..: . CCDS70 KVLTKPLTHLIVQLDSDMIMSSATGN--ATNSSLIKHWCGTTGQPTQEN------SSCGA 240 250 260 270 280 290 330 340 350 360 370 380 pF1KE2 PSLCWTDGIQNWTMKNVTYKENIAKCQHIFVNFHLPDLAVGTILLILSLLVLCGCLIMIV . : : :: : . :.:.:.. .: ::::: ::: :::::::::..:: CCDS70 FGPC--------TEKNSTAPADRLPCRHLFAGTELTDLAVGCILLAGSLLVLCGCLVLIV 300 310 320 330 340 390 400 410 420 430 440 pF1KE2 KILGSVLKGQVATVIKKTINTDFPFPFAWLTGYLAILVGAGMTFIVQSSSVFTSALTPLI :.:.:::.:.:: :.. .::.:::::..:: ::::.:.:::.:: .:::::::.:..::. CCDS70 KLLNSVLRGRVAQVVRTVINADFPFPLGWLGGYLAVLAGAGLTFALQSSSVFTAAVVPLM 350 360 370 380 390 400 450 460 470 480 490 500 pF1KE2 GIGVITIERAYPLTLGSNIGTTTTAILAALASPGNALRSSLQIALCHFFFNISGILLWYP :.:::...::::: :::::::::::.:::::::.. . :.::.:: :::::..:::::: CCDS70 GVGVISLDRAYPLLLGSNIGTTTTALLAALASPADRMLSALQVALIHFFFNLAGILLWYL 410 420 430 440 450 460 510 520 530 540 550 560 pF1KE2 IPFTRLPIRMAKGLGNISAKYRWFAVFYLIIFFFLIPLTVFGLSLAGWRVLVGVGVPVVF .: :::: .:. .: ..:.::: : ::.. :.:.::..::::::: :..:: :.: CCDS70 VPALRLPIPLARHFGVVTARYRWVAGVYLLLGFLLLPLAAFGLSLAGGMELAAVGGPLVG 470 480 490 500 510 520 570 580 590 600 610 620 pF1KE2 IIILVLCLRLLQSRCPRVLPKKLQNWNFLPLWMRSLKPWDAVVSKFTGCFQMRCCCCCRV ...::. . .:: : : :: .:..: .::.:..::.::: .:.. :: : : CCDS70 LVLLVILVTVLQRRRPAWLPVRLRSWAWLPVWLHSLEPWDRLVTR---------CCPCNV 530 540 550 560 570 630 640 650 660 670 680 pF1KE2 CCRACCLLCGCPKCCRCSKCCEDLEEAQEGQDVPVKAPETFDNITISREAQGEVPASDSK : CCDS70 CSPPKATTKEAYCYENPEILASQQL 580 590 >>CCDS54953.1 SLC34A1 gene_id:6569|Hs108|chr5 (340 aa) initn: 958 init1: 958 opt: 958 Z-score: 974.8 bits: 189.8 E(32554): 5.7e-48 Smith-Waterman score: 958; 71.1% identity (94.6% similar) in 204 aa overlap (100-303:103-306) 70 80 90 100 110 120 pF1KE2 WNLPTLQDSGIKWSERDTKGKILCFFQGIGRLILLLGFLYFFVCSLDILSSAFQLVGGKM .. :.: :::.::::::.:::::::.:::. CCDS54 PLPAKLALEEEQKPESRLVPKLRQAGAMLLKVPLMLTFLYLFVCSLDMLSSAFQLAGGKV 80 90 100 110 120 130 130 140 150 160 170 180 pF1KE2 AGQFFSNSSIMSNPLLGLVIGVLVTVLVQSSSTSTSIVVSMVSSSLLTVRAAIPIIMGAN ::..:....:.:::. :::.:.:::::::::::::::.::::::.:: : .:::::::.: CCDS54 AGDIFKDNAILSNPVAGLVVGILVTVLVQSSSTSTSIIVSMVSSGLLEVSSAIPIIMGSN 140 150 160 170 180 190 190 200 210 220 230 240 pF1KE2 IGTSITNTIVALMQVGDRSEFRRAFAGATVHDFFNWLSVLVLLPVEVATHYLEIITQLIV ::::.:::::::::.:::..:::::::::::: :::::::::::.:.:: ::. ::.:.: CCDS54 IGTSVTNTIVALMQAGDRTDFRRAFAGATVHDCFNWLSVLVLLPLEAATGYLHHITRLVV 200 210 220 230 240 250 250 260 270 280 290 300 pF1KE2 ESFHFKNGEDAPDLLKVITKPFTKLIVQLDKKVISQIAMNDEKAKNKSLVKIWCKTFTNK ::....:.:::::::.::.::::::.:::..::..:: .::. .:.::..::: CCDS54 ASFNIHGGRDAPDLLKIITEPFTKLIIQLDESVITSIATGDESLRNHSLIQIWCHPDSLQ 260 270 280 290 300 310 310 320 330 340 350 360 pF1KE2 TQINVTVPSTANCTSPSLCWTDGIQNWTMKNVTYKENIAKCQHIFVNFHLPDLAVGTILL CCDS54 QNLEGREITHFDLRKKQAMEDSSVPHCP 320 330 340 690 residues in 1 query sequences 18511270 residues in 32554 library sequences Tcomplib [36.3.4 Apr, 2011] (8 proc) start: Mon Nov 7 19:37:34 2016 done: Mon Nov 7 19:37:35 2016 Total Scan time: 2.330 Total Display time: 0.030 Function used was FASTA [36.3.4 Apr, 2011]