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: pF1KB3998, 191 aa 1>>>pF1KB3998 191 - 191 aa - 191 aa Library: /omim/omim.rfq.tfa 60827320 residues in 85289 sequences Statistics: Expectation_n fit: rho(ln(x))= 4.4082+/-0.000333; mu= 18.7524+/- 0.021 mean_var=78.1262+/-15.654, 0's: 0 Z-trim(116.9): 375 B-trim: 711 in 1/52 Lambda= 0.145103 statistics sampled from 28069 (28500) to 28069 sequences Algorithm: FASTA (3.7 Nov 2010) [optimized] Parameters: BL50 matrix (15:-5), open/ext: -10/-2 ktup: 2, E-join: 1 (0.729), E-opt: 0.2 (0.334), width: 16 Scan time: 5.300 The best scores are: opt bits E(85289) NP_426359 (OMIM: 116952,616737) cell division cont ( 191) 1283 277.4 9.9e-75 NP_001034891 (OMIM: 116952,616737) cell division c ( 191) 1224 265.0 5.1e-71 NP_001782 (OMIM: 116952,616737) cell division cont ( 191) 1224 265.0 5.1e-71 NP_005043 (OMIM: 602050) ras-related C3 botulinum ( 192) 946 206.8 1.7e-53 NP_008839 (OMIM: 602048) ras-related C3 botulinum ( 192) 941 205.8 3.5e-53 NP_002863 (OMIM: 602049,608203) ras-related C3 bot ( 192) 925 202.4 3.6e-52 NP_036381 (OMIM: 605857) rho-related GTP-binding p ( 205) 880 193.1 2.6e-49 NP_065714 (OMIM: 607653) rho-related GTP-binding p ( 214) 821 180.7 1.4e-45 NP_001656 (OMIM: 179505) rho-related GTP-binding p ( 191) 814 179.2 3.5e-45 XP_005252973 (OMIM: 179505) PREDICTED: rho-related ( 191) 814 179.2 3.5e-45 XP_016873208 (OMIM: 179505) PREDICTED: rho-related ( 191) 814 179.2 3.5e-45 NP_067028 (OMIM: 606366) rho-related GTP-binding p ( 258) 769 169.9 2.9e-42 NP_001303236 (OMIM: 602050) ras-related C3 botulin ( 170) 764 168.7 4.6e-42 XP_006724349 (OMIM: 602049,608203) PREDICTED: ras- ( 161) 754 166.5 1.9e-41 XP_011531028 (OMIM: 605857) PREDICTED: rho-related ( 192) 743 164.3 1.1e-40 XP_016859205 (OMIM: 605857) PREDICTED: rho-related ( 172) 733 162.2 4.2e-40 NP_001300870 (OMIM: 165390) transforming protein R ( 193) 696 154.5 9.7e-38 NP_001655 (OMIM: 165390) transforming protein RhoA ( 193) 696 154.5 9.7e-38 XP_006720277 (OMIM: 607653) PREDICTED: rho-related ( 176) 687 152.6 3.4e-37 NP_001036143 (OMIM: 165380) rho-related GTP-bindin ( 193) 684 152.0 5.5e-37 NP_001036144 (OMIM: 165380) rho-related GTP-bindin ( 193) 684 152.0 5.5e-37 NP_786886 (OMIM: 165380) rho-related GTP-binding p ( 193) 684 152.0 5.5e-37 NP_004031 (OMIM: 165370) rho-related GTP-binding p ( 196) 660 147.0 1.8e-35 NP_055393 (OMIM: 605781) rho-related GTP-binding p ( 210) 595 133.4 2.4e-31 XP_011531030 (OMIM: 605857) PREDICTED: rho-related ( 157) 583 130.7 1.1e-30 NP_001300872 (OMIM: 165390) transforming protein R ( 187) 575 129.2 4e-30 XP_011511994 (OMIM: 602037) PREDICTED: rho-related ( 191) 524 118.5 6.7e-27 NP_001265298 (OMIM: 602037) rho-related GTP-bindin ( 191) 524 118.5 6.7e-27 NP_001265293 (OMIM: 602037) rho-related GTP-bindin ( 191) 524 118.5 6.7e-27 NP_001265296 (OMIM: 602037) rho-related GTP-bindin ( 191) 524 118.5 6.7e-27 XP_016863678 (OMIM: 602037) PREDICTED: rho-related ( 191) 524 118.5 6.7e-27 NP_001265288 (OMIM: 602037) rho-related GTP-bindin ( 191) 524 118.5 6.7e-27 NP_001265294 (OMIM: 602037) rho-related GTP-bindin ( 191) 524 118.5 6.7e-27 NP_001265289 (OMIM: 602037) rho-related GTP-bindin ( 191) 524 118.5 6.7e-27 NP_001265292 (OMIM: 602037) rho-related GTP-bindin ( 191) 524 118.5 6.7e-27 NP_001265297 (OMIM: 602037) rho-related GTP-bindin ( 191) 524 118.5 6.7e-27 NP_001265295 (OMIM: 602037) rho-related GTP-bindin ( 191) 524 118.5 6.7e-27 NP_004301 (OMIM: 602037) rho-related GTP-binding p ( 191) 524 118.5 6.7e-27 NP_001265290 (OMIM: 602037) rho-related GTP-bindin ( 191) 524 118.5 6.7e-27 XP_016863677 (OMIM: 602037) PREDICTED: rho-related ( 191) 524 118.5 6.7e-27 NP_001265291 (OMIM: 602037) rho-related GTP-bindin ( 191) 524 118.5 6.7e-27 NP_061485 (OMIM: 602048) ras-related C3 botulinum ( 211) 523 118.3 8.2e-27 NP_005159 (OMIM: 602924) rho-related GTP-binding p ( 244) 513 116.3 3.8e-26 NP_001241667 (OMIM: 602924) rho-related GTP-bindin ( 244) 513 116.3 3.8e-26 NP_005431 (OMIM: 601555) rho-related GTP-binding p ( 227) 483 110.0 2.9e-24 XP_011523618 (OMIM: 601555) PREDICTED: rho-related ( 234) 483 110.0 2.9e-24 NP_001300873 (OMIM: 165390) transforming protein R ( 173) 470 107.1 1.6e-23 NP_055285 (OMIM: 609038) rho-related GTP-binding p ( 232) 435 100.0 3.1e-21 XP_011523619 (OMIM: 601555) PREDICTED: rho-related ( 206) 407 94.0 1.7e-19 NP_001300874 (OMIM: 165390) transforming protein R ( 112) 360 83.9 1e-16 >>NP_426359 (OMIM: 116952,616737) cell division control (191 aa) initn: 1283 init1: 1283 opt: 1283 Z-score: 1462.7 bits: 277.4 E(85289): 9.9e-75 Smith-Waterman score: 1283; 100.0% identity (100.0% similar) in 191 aa overlap (1-191:1-191) 10 20 30 40 50 60 pF1KB3 MQTIKCVVVGDGAVGKTCLLISYTTNKFPSEYVPTVFDNYAVTVMIGGEPYTLGLFDTAG :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: NP_426 MQTIKCVVVGDGAVGKTCLLISYTTNKFPSEYVPTVFDNYAVTVMIGGEPYTLGLFDTAG 10 20 30 40 50 60 70 80 90 100 110 120 pF1KB3 QEDYDRLRPLSYPQTDVFLVCFSVVSPSSFENVKEKWVPEITHHCPKTPFLLVGTQIDLR :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: NP_426 QEDYDRLRPLSYPQTDVFLVCFSVVSPSSFENVKEKWVPEITHHCPKTPFLLVGTQIDLR 70 80 90 100 110 120 130 140 150 160 170 180 pF1KB3 DDPSTIEKLAKNKQKPITPETAEKLARDLKAVKYVECSALTQRGLKNVFDEAILAALEPP :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: NP_426 DDPSTIEKLAKNKQKPITPETAEKLARDLKAVKYVECSALTQRGLKNVFDEAILAALEPP 130 140 150 160 170 180 190 pF1KB3 ETQPKRKCCIF ::::::::::: NP_426 ETQPKRKCCIF 190 >>NP_001034891 (OMIM: 116952,616737) cell division contr (191 aa) initn: 1271 init1: 1224 opt: 1224 Z-score: 1395.9 bits: 265.0 E(85289): 5.1e-71 Smith-Waterman score: 1224; 95.3% identity (98.4% similar) in 191 aa overlap (1-191:1-191) 10 20 30 40 50 60 pF1KB3 MQTIKCVVVGDGAVGKTCLLISYTTNKFPSEYVPTVFDNYAVTVMIGGEPYTLGLFDTAG :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: NP_001 MQTIKCVVVGDGAVGKTCLLISYTTNKFPSEYVPTVFDNYAVTVMIGGEPYTLGLFDTAG 10 20 30 40 50 60 70 80 90 100 110 120 pF1KB3 QEDYDRLRPLSYPQTDVFLVCFSVVSPSSFENVKEKWVPEITHHCPKTPFLLVGTQIDLR :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: NP_001 QEDYDRLRPLSYPQTDVFLVCFSVVSPSSFENVKEKWVPEITHHCPKTPFLLVGTQIDLR 70 80 90 100 110 120 130 140 150 160 170 180 pF1KB3 DDPSTIEKLAKNKQKPITPETAEKLARDLKAVKYVECSALTQRGLKNVFDEAILAALEPP ::::::::::::::::::::::::::::::::::::::::::.::::::::::::::::: NP_001 DDPSTIEKLAKNKQKPITPETAEKLARDLKAVKYVECSALTQKGLKNVFDEAILAALEPP 130 140 150 160 170 180 190 pF1KB3 ETQPKRKCCIF : . .:.: .. NP_001 EPKKSRRCVLL 190 >>NP_001782 (OMIM: 116952,616737) cell division control (191 aa) initn: 1271 init1: 1224 opt: 1224 Z-score: 1395.9 bits: 265.0 E(85289): 5.1e-71 Smith-Waterman score: 1224; 95.3% identity (98.4% similar) in 191 aa overlap (1-191:1-191) 10 20 30 40 50 60 pF1KB3 MQTIKCVVVGDGAVGKTCLLISYTTNKFPSEYVPTVFDNYAVTVMIGGEPYTLGLFDTAG :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: NP_001 MQTIKCVVVGDGAVGKTCLLISYTTNKFPSEYVPTVFDNYAVTVMIGGEPYTLGLFDTAG 10 20 30 40 50 60 70 80 90 100 110 120 pF1KB3 QEDYDRLRPLSYPQTDVFLVCFSVVSPSSFENVKEKWVPEITHHCPKTPFLLVGTQIDLR :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: NP_001 QEDYDRLRPLSYPQTDVFLVCFSVVSPSSFENVKEKWVPEITHHCPKTPFLLVGTQIDLR 70 80 90 100 110 120 130 140 150 160 170 180 pF1KB3 DDPSTIEKLAKNKQKPITPETAEKLARDLKAVKYVECSALTQRGLKNVFDEAILAALEPP ::::::::::::::::::::::::::::::::::::::::::.::::::::::::::::: NP_001 DDPSTIEKLAKNKQKPITPETAEKLARDLKAVKYVECSALTQKGLKNVFDEAILAALEPP 130 140 150 160 170 180 190 pF1KB3 ETQPKRKCCIF : . .:.: .. NP_001 EPKKSRRCVLL 190 >>NP_005043 (OMIM: 602050) ras-related C3 botulinum toxi (192 aa) initn: 994 init1: 917 opt: 946 Z-score: 1081.4 bits: 206.8 E(85289): 1.7e-53 Smith-Waterman score: 946; 70.3% identity (88.0% similar) in 192 aa overlap (1-191:1-192) 10 20 30 40 50 60 pF1KB3 MQTIKCVVVGDGAVGKTCLLISYTTNKFPSEYVPTVFDNYAVTVMIGGEPYTLGLFDTAG ::.::::::::::::::::::::::: ::.::.:::::::...::. :.: .:::.:::: NP_005 MQAIKCVVVGDGAVGKTCLLISYTTNAFPGEYIPTVFDNYSANVMVDGKPVNLGLWDTAG 10 20 30 40 50 60 70 80 90 100 110 120 pF1KB3 QEDYDRLRPLSYPQTDVFLVCFSVVSPSSFENVKEKWVPEITHHCPKTPFLLVGTQIDLR :::::::::::::::::::.:::.:::.:::::. :: ::. ::::.::.:::::..::: NP_005 QEDYDRLRPLSYPQTDVFLICFSLVSPASFENVRAKWYPEVRHHCPHTPILLVGTKLDLR 70 80 90 100 110 120 130 140 150 160 170 180 pF1KB3 DDPSTIEKLAKNKQKPITPETAEKLARDLKAVKYVECSALTQRGLKNVFDEAILAALEPP :: .:::.: .: ::: . .::.. .:::.:::::::::::.:::::: :.: :: NP_005 DDKDTIERLRDKKLAPITYPQGLAMAREIGSVKYLECSALTQRGLKTVFDEAIRAVLCPP 130 140 150 160 170 180 190 pF1KB3 ETQ-PKRKCCIF .. : .:: .: NP_005 PVKKPGKKCTVF 190 >>NP_008839 (OMIM: 602048) ras-related C3 botulinum toxi (192 aa) initn: 1011 init1: 941 opt: 941 Z-score: 1075.7 bits: 205.8 E(85289): 3.5e-53 Smith-Waterman score: 941; 69.1% identity (89.0% similar) in 191 aa overlap (1-191:1-191) 10 20 30 40 50 60 pF1KB3 MQTIKCVVVGDGAVGKTCLLISYTTNKFPSEYVPTVFDNYAVTVMIGGEPYTLGLFDTAG ::.::::::::::::::::::::::: ::.::.:::::::...::. :.: .:::.:::: NP_008 MQAIKCVVVGDGAVGKTCLLISYTTNAFPGEYIPTVFDNYSANVMVDGKPVNLGLWDTAG 10 20 30 40 50 60 70 80 90 100 110 120 pF1KB3 QEDYDRLRPLSYPQTDVFLVCFSVVSPSSFENVKEKWVPEITHHCPKTPFLLVGTQIDLR :::::::::::::::::::.:::.:::.:::::. :: ::. ::::.::..::::..::: NP_008 QEDYDRLRPLSYPQTDVFLICFSLVSPASFENVRAKWYPEVRHHCPNTPIILVGTKLDLR 70 80 90 100 110 120 130 140 150 160 170 180 pF1KB3 DDPSTIEKLAKNKQKPITPETAEKLARDLKAVKYVECSALTQRGLKNVFDEAILAALEPP :: .::::: ..: ::: . .:... ::::.:::::::::::.:::::: :.: :: NP_008 DDKDTIEKLKEKKLTPITYPQGLAMAKEIGAVKYLECSALTQRGLKTVFDEAIRAVLCPP 130 140 150 160 170 180 190 pF1KB3 ETQPKRKCCIF .. ... :.. NP_008 PVKKRKRKCLLL 190 >>NP_002863 (OMIM: 602049,608203) ras-related C3 botulin (192 aa) initn: 954 init1: 919 opt: 925 Z-score: 1057.6 bits: 202.4 E(85289): 3.6e-52 Smith-Waterman score: 925; 69.3% identity (87.5% similar) in 192 aa overlap (1-191:1-192) 10 20 30 40 50 60 pF1KB3 MQTIKCVVVGDGAVGKTCLLISYTTNKFPSEYVPTVFDNYAVTVMIGGEPYTLGLFDTAG ::.::::::::::::::::::::::: ::.::.:::::::...::. ..: .:::.:::: NP_002 MQAIKCVVVGDGAVGKTCLLISYTTNAFPGEYIPTVFDNYSANVMVDSKPVNLGLWDTAG 10 20 30 40 50 60 70 80 90 100 110 120 pF1KB3 QEDYDRLRPLSYPQTDVFLVCFSVVSPSSFENVKEKWVPEITHHCPKTPFLLVGTQIDLR :::::::::::::::::::.:::.:::.:.:::. :: ::. ::::.::..::::..::: NP_002 QEDYDRLRPLSYPQTDVFLICFSLVSPASYENVRAKWFPEVRHHCPSTPIILVGTKLDLR 70 80 90 100 110 120 130 140 150 160 170 pF1KB3 DDPSTIEKLAKNKQKPITPETAEKLARDLKAVKYVECSALTQRGLKNVFDEAILAALEP- :: .::::: ..: ::: . ::... .:::.:::::::::::.:::::: :.: : NP_002 DDKDTIEKLKEKKLAPITYPQGLALAKEIDSVKYLECSALTQRGLKTVFDEAIRAVLCPQ 130 140 150 160 170 180 180 190 pF1KB3 PETQPKRKCCIF : : :: : .. NP_002 PTRQQKRACSLL 190 >>NP_036381 (OMIM: 605857) rho-related GTP-binding prote (205 aa) initn: 888 init1: 868 opt: 880 Z-score: 1006.4 bits: 193.1 E(85289): 2.6e-49 Smith-Waterman score: 880; 64.9% identity (84.5% similar) in 194 aa overlap (4-190:10-203) 10 20 30 40 50 pF1KB3 MQTIKCVVVGDGAVGKTCLLISYTTNKFPSEYVPTVFDNYAVTVMIGGEPYTLG .::::::::::::::::.::... :: ::::::::.:::.: .::. : :: NP_036 MAHGPGALMLKCVVVGDGAVGKTCLLMSYANDAFPEEYVPTVFDHYAVSVTVGGKQYLLG 10 20 30 40 50 60 60 70 80 90 100 110 pF1KB3 LFDTAGQEDYDRLRPLSYPQTDVFLVCFSVVSPSSFENVKEKWVPEITHHCPKTPFLLVG :.:::::::::::::::::.:::::.:::::.:.::.::::.::::. .. :..::::.: NP_036 LYDTAGQEDYDRLRPLSYPMTDVFLICFSVVNPASFQNVKEEWVPELKEYAPNVPFLLIG 70 80 90 100 110 120 120 130 140 150 160 170 pF1KB3 TQIDLRDDPSTIEKLAKNKQKPITPETAEKLARDLKAVKYVECSALTQRGLKNVFDEAIL :::::::::.:. .: :.::: : ..:::... : :::::::::.:::.::::::. NP_036 TQIDLRDDPKTLARLNDMKEKPICVEQGQKLAKEIGACCYVECSALTQKGLKTVFDEAII 130 140 150 160 170 180 180 190 pF1KB3 AALEPPETQPKRK-------CCIF : : : . :.. ::. NP_036 AILTPKKHTVKKRIGSRCINCCLIT 190 200 >>NP_065714 (OMIM: 607653) rho-related GTP-binding prote (214 aa) initn: 843 init1: 810 opt: 821 Z-score: 939.4 bits: 180.7 E(85289): 1.4e-45 Smith-Waterman score: 821; 62.6% identity (85.6% similar) in 187 aa overlap (2-188:20-204) 10 20 30 40 pF1KB3 MQTIKCVVVGDGAVGKTCLLISYTTNKFPSEYVPTVFDNYAV . .::::::::::::::::.::... :: ::::::::.::: NP_065 MNCKEGTDSSCGCRGNDEKKMLKCVVVGDGAVGKTCLLMSYANDAFPEEYVPTVFDHYAV 10 20 30 40 50 60 50 60 70 80 90 100 pF1KB3 TVMIGGEPYTLGLFDTAGQEDYDRLRPLSYPQTDVFLVCFSVVSPSSFENVKEKWVPEIT :: .::. . :::.::::::::..:::::::.:::::.:::::.:.:..::.:.::::. NP_065 TVTVGGKQHLLGLYDTAGQEDYNQLRPLSYPNTDVFLICFSVVNPASYHNVQEEWVPELK 70 80 90 100 110 120 110 120 130 140 150 160 pF1KB3 HHCPKTPFLLVGTQIDLRDDPSTIEKLAKNKQKPITPETAEKLARDLKAVKYVECSALTQ :..:..:.::::::::::.:. .: :.::.: : . :::. . : :.::::::: NP_065 DCMPHVPYVLIGTQIDLRDDPKTLARLLYMKEKPLTYEHGVKLAKAIGAQCYLECSALTQ 130 140 150 160 170 180 170 180 190 pF1KB3 RGLKNVFDEAILAALEPPETQPKRKCCIF .::: :::::::. ..: . :..: NP_065 KGLKAVFDEAILTIFHPK--KKKKRCSEGHSCCSII 190 200 210 >>NP_001656 (OMIM: 179505) rho-related GTP-binding prote (191 aa) initn: 814 init1: 799 opt: 814 Z-score: 932.1 bits: 179.2 E(85289): 3.5e-45 Smith-Waterman score: 814; 59.2% identity (84.3% similar) in 191 aa overlap (1-191:1-191) 10 20 30 40 50 60 pF1KB3 MQTIKCVVVGDGAVGKTCLLISYTTNKFPSEYVPTVFDNYAVTVMIGGEPYTLGLFDTAG ::.:::::::::::::::::: :::: ::.::.:::::::.. . :. .:.:.:::: NP_001 MQSIKCVVVGDGAVGKTCLLICYTTNAFPKEYIPTVFDNYSAQSAVDGRTVNLNLWDTAG 10 20 30 40 50 60 70 80 90 100 110 120 pF1KB3 QEDYDRLRPLSYPQTDVFLVCFSVVSPSSFENVKEKWVPEITHHCPKTPFLLVGTQIDLR ::.::::: ::::::.::..:::..:: :.:::..:: ::. :::: .:.:::::. ::: NP_001 QEEYDRLRTLSYPQTNVFVICFSIASPPSYENVRHKWHPEVCHHCPDVPILLVGTKKDLR 70 80 90 100 110 120 130 140 150 160 170 180 pF1KB3 DDPSTIEKLAKNKQKPITPETAEKLARDLKAVKYVECSALTQRGLKNVFDEAILAALEPP .:.:...: .. : ::::. .. ::....::.:.::::: : :.:.:: ::. :.:.: NP_001 AQPDTLRRLKEQGQAPITPQQGQALAKQIHAVRYLECSALQQDGVKEVFAEAVRAVLNPT 130 140 150 160 170 180 190 pF1KB3 ETQPKRKCCIF . :.: .. NP_001 PIKRGRSCILL 190 >>XP_005252973 (OMIM: 179505) PREDICTED: rho-related GTP (191 aa) initn: 814 init1: 799 opt: 814 Z-score: 932.1 bits: 179.2 E(85289): 3.5e-45 Smith-Waterman score: 814; 59.2% identity (84.3% similar) in 191 aa overlap (1-191:1-191) 10 20 30 40 50 60 pF1KB3 MQTIKCVVVGDGAVGKTCLLISYTTNKFPSEYVPTVFDNYAVTVMIGGEPYTLGLFDTAG ::.:::::::::::::::::: :::: ::.::.:::::::.. . :. .:.:.:::: XP_005 MQSIKCVVVGDGAVGKTCLLICYTTNAFPKEYIPTVFDNYSAQSAVDGRTVNLNLWDTAG 10 20 30 40 50 60 70 80 90 100 110 120 pF1KB3 QEDYDRLRPLSYPQTDVFLVCFSVVSPSSFENVKEKWVPEITHHCPKTPFLLVGTQIDLR ::.::::: ::::::.::..:::..:: :.:::..:: ::. :::: .:.:::::. ::: XP_005 QEEYDRLRTLSYPQTNVFVICFSIASPPSYENVRHKWHPEVCHHCPDVPILLVGTKKDLR 70 80 90 100 110 120 130 140 150 160 170 180 pF1KB3 DDPSTIEKLAKNKQKPITPETAEKLARDLKAVKYVECSALTQRGLKNVFDEAILAALEPP .:.:...: .. : ::::. .. ::....::.:.::::: : :.:.:: ::. :.:.: XP_005 AQPDTLRRLKEQGQAPITPQQGQALAKQIHAVRYLECSALQQDGVKEVFAEAVRAVLNPT 130 140 150 160 170 180 190 pF1KB3 ETQPKRKCCIF . :.: .. XP_005 PIKRGRSCILL 190 191 residues in 1 query sequences 60827320 residues in 85289 library sequences Tcomplib [36.3.4 Apr, 2011] (8 proc) start: Sat Nov 5 05:31:21 2016 done: Sat Nov 5 05:31:22 2016 Total Scan time: 5.300 Total Display time: -0.010 Function used was FASTA [36.3.4 Apr, 2011]