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: pF1KE6715, 748 aa 1>>>pF1KE6715 748 - 748 aa - 748 aa Library: human.CCDS.faa 18511270 residues in 32554 sequences Statistics: Expectation_n fit: rho(ln(x))= 7.5833+/-0.000885; mu= 10.8846+/- 0.054 mean_var=133.2374+/-26.518, 0's: 0 Z-trim(110.9): 23 B-trim: 8 in 1/52 Lambda= 0.111112 statistics sampled from 11955 (11974) to 11955 sequences Algorithm: FASTA (3.7 Nov 2010) [optimized] Parameters: BL50 matrix (15:-5), open/ext: -10/-2 ktup: 2, E-join: 1 (0.704), E-opt: 0.2 (0.368), width: 16 Scan time: 2.990 The best scores are: opt bits E(32554) CCDS10962.1 JPH3 gene_id:57338|Hs108|chr16 ( 748) 4999 812.9 0 CCDS6217.1 JPH1 gene_id:56704|Hs108|chr8 ( 661) 2187 362.1 1.6e-99 CCDS13325.1 JPH2 gene_id:57158|Hs108|chr20 ( 696) 1956 325.1 2.3e-88 CCDS9603.1 JPH4 gene_id:84502|Hs108|chr14 ( 628) 1543 258.9 1.8e-68 CCDS13326.1 JPH2 gene_id:57158|Hs108|chr20 ( 129) 824 143.3 2.3e-34 >>CCDS10962.1 JPH3 gene_id:57338|Hs108|chr16 (748 aa) initn: 4999 init1: 4999 opt: 4999 Z-score: 4335.7 bits: 812.9 E(32554): 0 Smith-Waterman score: 4999; 99.6% identity (99.6% similar) in 748 aa overlap (1-748:1-748) 10 20 30 40 50 60 pF1KE6 MSSGGRFNFDDGGSYCGGWEDGKAHGHGVCTGPKGQGEYTGSWSHGFEVLGVYTWPSGNT :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: CCDS10 MSSGGRFNFDDGGSYCGGWEDGKAHGHGVCTGPKGQGEYTGSWSHGFEVLGVYTWPSGNT 10 20 30 40 50 60 70 80 90 100 110 120 pF1KE6 YQGTWAQGKRHGIGLESKGKWVYKGEWTHGFKGRYGVRECAGNGAKYEGTWSNGLQDGYG :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: CCDS10 YQGTWAQGKRHGIGLESKGKWVYKGEWTHGFKGRYGVRECAGNGAKYEGTWSNGLQDGYG 70 80 90 100 110 120 130 140 150 160 170 180 pF1KE6 TETYSDGGTYQGQWVGGMRQGYGVRQSVPYGMAAVIRSPLRTSINSLRSEHTNGTALHPD :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: CCDS10 TETYSDGGTYQGQWVGGMRQGYGVRQSVPYGMAAVIRSPLRTSINSLRSEHTNGTALHPD 130 140 150 160 170 180 190 200 210 220 230 240 pF1KE6 ASPAVAGSPAVSRGGFVLVAHSDSEILKSKKKGLFRRSLLSGLKLRKSESKSSLASQRSK :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: CCDS10 ASPAVAGSPAVSRGGFVLVAHSDSEILKSKKKGLFRRSLLSGLKLRKSESKSSLASQRSK 190 200 210 220 230 240 250 260 270 280 290 300 pF1KE6 QSSFRSEAGMSTVSSTASDIHSTISLGEAEAELAVIEDDIDATTTETYVGEWKNDKRSGF :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: CCDS10 QSSFRSEAGMSTVSSTASDIHSTISLGEAEAELAVIEDDIDATTTETYVGEWKNDKRSGF 250 260 270 280 290 300 310 320 330 340 350 360 pF1KE6 GVSQRSDGLKYEGEWASNRRHGYGCMTFPDGTKEEGKYKQNILVGGKRKNLIPLRASKIR :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: CCDS10 GVSQRSDGLKYEGEWASNRRHGYGCMTFPDGTKEEGKYKQNILVGGKRKNLIPLRASKIR 310 320 330 340 350 360 370 380 390 400 410 420 pF1KE6 EKVDRAVEAAERAATTAKQKAEIAASRTSHSRAKAEAALTAAQKAQEEARIARITAKEFS ::::::::::::::: :::::::::::::::::::::::::::::::::::::::::::: CCDS10 EKVDRAVEAAERAATIAKQKAEIAASRTSHSRAKAEAALTAAQKAQEEARIARITAKEFS 370 380 390 400 410 420 430 440 450 460 470 480 pF1KE6 PSFQHRENGLEYQRPKRQTSCDDIEVLSTGTPLQQESPELYRKGTTPSDLTTDDSPLQSF ::::::::::::::::::::::::::::::::::::::::::::::::::: :::::::: CCDS10 PSFQHRENGLEYQRPKRQTSCDDIEVLSTGTPLQQESPELYRKGTTPSDLTPDDSPLQSF 430 440 450 460 470 480 490 500 510 520 530 540 pF1KE6 PTSPAATPPPAPAARNKVAHFSRQVSVDEERGGDIQMLLEGRAGDCARSSWGEEQAGGSR :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: CCDS10 PTSPAATPPPAPAARNKVAHFSRQVSVDEERGGDIQMLLEGRAGDCARSSWGEEQAGGSR 490 500 510 520 530 540 550 560 570 580 590 600 pF1KE6 GVRSGALRGGLLVDDFRTRGSGRKQPGNPKPRERRTESPPVFTWTSHHRASNHSPGGSRL :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: CCDS10 GVRSGALRGGLLVDDFRTRGSGRKQPGNPKPRERRTESPPVFTWTSHHRASNHSPGGSRL 550 560 570 580 590 600 610 620 630 640 650 660 pF1KE6 LELQEEKLSNYRMEMKPLLRMETHPQKRRYSKGGACRGLGDDHRLEDRGFGVQRLRSKAQ :::::::::::::::::::::::::::::::::::::::::::: ::::::::::::::: CCDS10 LELQEEKLSNYRMEMKPLLRMETHPQKRRYSKGGACRGLGDDHRPEDRGFGVQRLRSKAQ 610 620 630 640 650 660 670 680 690 700 710 720 pF1KE6 NKENFRPASSAEPAVQKLASLRLGGAEPRLLRWDLTFSPPQKSLPVALESDEENGDELKS :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: CCDS10 NKENFRPASSAEPAVQKLASLRLGGAEPRLLRWDLTFSPPQKSLPVALESDEENGDELKS 670 680 690 700 710 720 730 740 pF1KE6 STGSAPILVVMVILLNIGVAILFINFFI :::::::::::::::::::::::::::: CCDS10 STGSAPILVVMVILLNIGVAILFINFFI 730 740 >>CCDS6217.1 JPH1 gene_id:56704|Hs108|chr8 (661 aa) initn: 2168 init1: 1120 opt: 2187 Z-score: 1900.4 bits: 362.1 E(32554): 1.6e-99 Smith-Waterman score: 2187; 63.4% identity (83.2% similar) in 517 aa overlap (3-510:2-505) 10 20 30 40 50 60 pF1KE6 MSSGGRFNFDDGGSYCGGWEDGKAHGHGVCTGPKGQGEYTGSWSHGFEVLGVYTWPSGNT .::::.:::::.::::::.:::::::.::::::::::.:::::::::.: :::::::: CCDS62 MTGGRFDFDDGGTYCGGWEEGKAHGHGICTGPKGQGEYSGSWSHGFEVVGGYTWPSGNT 10 20 30 40 50 70 80 90 100 110 120 pF1KE6 YQGTWAQGKRHGIGLESKGKWVYKGEWTHGFKGRYGVRECAGNGAKYEGTWSNGLQDGYG ::: ::::::::.:.:.::::.:.:::.::::::::::. . :.:::::::::::::: CCDS62 YQGYWAQGKRHGLGVETKGKWMYRGEWSHGFKGRYGVRQSLCTPARYEGTWSNGLQDGYG 60 70 80 90 100 110 130 140 150 160 170 180 pF1KE6 TETYSDGGTYQGQWVGGMRQGYGVRQSVPYGMAAVIRSPLRTSINSLRSEHTNGTALHPD .:::.:::::::::.::::.:::::::::::::.:::::::::. :::::..::..:: : CCDS62 VETYGDGGTYQGQWAGGMRHGYGVRQSVPYGMATVIRSPLRTSLASLRSEQSNGSVLH-D 120 130 140 150 160 170 190 200 210 220 230 pF1KE6 ASPAVAGSPAVSRGGFVLVAHSDSEILKSKKKGLFRR-SLLSGLKLRKSESKSSLASQRS :. :.: ::: .:::::: :.:.:. .:: ::::: :::...::::::::::..:.:: CCDS62 AA-AAADSPAGTRGGFVLNFHADAELAGKKKGGLFRRGSLLGSMKLRKSESKSSISSKRS 180 190 200 210 220 230 240 250 260 270 280 290 pF1KE6 KQSSFRSEAGMSTVSSTASDIHSTISLGEAEAELAVIEDDIDATTTETYVGEWKNDKRSG : ::.:.:: .:: :: .::::.:... .. .:: .::::::::.::::::::.: CCDS62 ---SVRSDAAMSRISS--SDANSTISFGDVDCDFCPVEDHVDATTTETYMGEWKNDKRNG 240 250 260 270 280 290 300 310 320 330 340 350 pF1KE6 FGVSQRSDGLKYEGEWASNRRHGYGCMTFPDGTKEEGKYKQNILVGGKRKNLIPLRASKI ::::.::.:.:::::::.:.:::::: .::::.:::::::.:::: : ::.:::.: .: CCDS62 FGVSERSNGMKYEGEWANNKRHGYGCTVFPDGSKEEGKYKNNILVRGIRKQLIPIRHTKT 300 310 320 330 340 350 360 370 380 390 400 410 pF1KE6 REKVDRAVEAAERAATTAKQKAEIAASRTSHSRAKAEAALTAAQKAQEEARIARITAKEF :::::::.:.:.:::. :. :.::: :::.:.::::.:: :: :..: ::: .:.:. CCDS62 REKVDRAIEGAQRAAAMARTKVEIANSRTAHARAKADAADQAALAARQECDIARAVAREL 360 370 380 390 400 410 420 430 440 450 460 470 pF1KE6 SPSFQHRENGLEYQRPKRQTSCDDIEV----LSTGTPLQQESPELYRKGTTPSDLTTDDS ::.: . : .: . . : . : : . : .:::..::::::: . . : CCDS62 SPDFY--QPGPDYVKQRFQEGVDAKENPEEKVPEKPPTPKESPHFYRKGTTPPR-SPEAS 420 430 440 450 460 480 490 500 510 520 530 pF1KE6 PLQSFPTSPAATPPP----APAARNKVAHFSRQVSVDEERGGDIQMLLEGRAGDCARSSW : .: :::..: : :.. .. . .:.:. ::. CCDS62 PKHSH--SPASSPKPLKKQNPSSGARLNQDKRSVA-DEQVTAIVNKPLMSKAPTKEAGAV 470 480 490 500 510 520 540 550 560 570 580 590 pF1KE6 GEEQAGGSRGVRSGALRGGLLVDDFRTRGSGRKQPGNPKPRERRTESPPVFTWTSHHRAS CCDS62 VPQSKYSGRHHIPNPSNGELHSQYHGYYVKLNAPQHPPVDVEDGDGSSQSSSALVHKPSA 530 540 550 560 570 580 >>CCDS13325.1 JPH2 gene_id:57158|Hs108|chr20 (696 aa) initn: 1820 init1: 1099 opt: 1956 Z-score: 1699.9 bits: 325.1 E(32554): 2.3e-88 Smith-Waterman score: 1989; 46.8% identity (66.7% similar) in 771 aa overlap (3-747:2-695) 10 20 30 40 50 60 pF1KE6 MSSGGRFNFDDGGSYCGGWEDGKAHGHGVCTGPKGQGEYTGSWSHGFEVLGVYTWPSGNT :::::.:::::.:::::: :::::::.::::::::::.:::. :::: :::::::::: CCDS13 MSGGRFDFDDGGAYCGGWEGGKAHGHGLCTGPKGQGEYSGSWNFGFEVAGVYTWPSGNT 10 20 30 40 50 70 80 90 100 110 120 pF1KE6 YQGTWAQGKRHGIGLESKGKWVYKGEWTHGFKGRYGVRECAGNGAKYEGTWSNGLQDGYG ..: :.::::::.:.:.::.:.::::::::::::::.:. ...::::::::.:::::::: CCDS13 FEGYWSQGKRHGLGIETKGRWLYKGEWTHGFKGRYGIRQSSSSGAKYEGTWNNGLQDGYG 60 70 80 90 100 110 130 140 150 160 170 180 pF1KE6 TETYSDGGTYQGQWVGGMRQGYGVRQSVPYGMAAVIRSPLRTSINSLRSEHTNGTALHPD ::::.::::::::...:::.:::::::::::::.:.:::::::..::::::.::: . :: CCDS13 TETYADGGTYQGQFTNGMRHGYGVRQSVPYGMAVVVRSPLRTSLSSLRSEHSNGT-VAPD 120 130 140 150 160 170 190 200 210 220 230 pF1KE6 --ASPAVAG----SPAVSRGGFVL--VAHSDSEILKSKKKGLFRRSLLSGLKLRKSESKS :::: : :::. ::::.: .:.... : :::.:. : : :::..::.. CCDS13 SPASPASDGPALPSPAIPRGGFALSLLANAEAAARAPKGGGLFQRGALLG-KLRRAESRT 180 190 200 210 220 230 240 250 260 270 280 pF1KE6 SLASQRSKQSSFRSEAGMSTVSSTASDIHSTISLGEAEA---ELAVIEDDIDATTTETYV :..::::. : ..:. .:: ::: :: ::::: : : .: ::::::::::. CCDS13 SVGSQRSRVSFLKSD-----LSSGASDAASTASLGEAAEGADEAAPFEADIDATTTETYM 240 250 260 270 280 290 290 300 310 320 330 340 pF1KE6 GEWKNDKRSGFGVSQRSDGLKYEGEWASNRRHGYGCMTFPDGTKEEGKYKQNILVGGKRK ::::::::::::::.::.::.::::: .: :::::: :.::: .:::::..:.:: .. CCDS13 GEWKNDKRSGFGVSERSSGLRYEGEWLDNLRHGYGCTTLPDGHREEGKYRHNVLVKDTKR 300 310 320 330 340 350 350 360 370 380 390 400 pF1KE6 NLIPLRASKIREKVDRAVEAAERAATTAKQKAEIAASRTSHSRAKAEAALTAAQKAQEEA .. :...:.:.::...::.:.:::. :.::::::::::::..:::::: :: :..:. CCDS13 RMLQLKSNKVRQKVEHSVEGAQRAAAIARQKAEIAASRTSHAKAKAEAAEQAALAANQES 360 370 380 390 400 410 410 420 430 440 450 pF1KE6 RIARITAKEFSPSFQHRENGLEYQRPKR-QTSCDDIEVL---------STGTPLQ-QESP ::: :.:..:.: . : :::. . : .. : : ..: : .::: CCDS13 NIARTLARELAPDF--YQPGPEYQKRRLLQEILENSESLLEPPDRGAGAAGLPQPPRESP 420 430 440 450 460 470 460 470 480 490 500 510 pF1KE6 ELYRKGTTPSDLTTDDSPLQSFPTSPAATPPPAPAARNKVAHFSRQVSVDEERGGDIQML .:... :: : . : :::.::: : :. . : .: CCDS13 QLHER-ETPR-------PEGGSP-SPAGTPPQPKRPRPGVS-----------KDG---LL 480 490 500 520 530 540 550 560 570 pF1KE6 LEGRAGDCARSSWGEEQAG-GSRGVRSGALRGGLLVDDFRTRGSGRKQPGNPKPRERRTE : .:. : .: :::.: ..:.::..:. : .. : CCDS13 SPG--------AWNGEPSGEGSRSVTP-------------SEGAGRRSPARPATERMAIE 510 520 530 540 580 590 600 610 620 630 pF1KE6 SPPVFTWTSHHRASNHSPGGSRLLELQ-EEKLSNYRMEMKPLLRMETHPQKRRYSKGGAC : . :. :: :. . .: .. : . : . .:. CCDS13 ------------ALQAPPAPSREPEVALYQGYHSYAVRTTPPEPPPFEDQPEPEVSGSES 550 560 570 580 590 640 650 660 670 680 690 pF1KE6 RGLGDDHRLEDRGFGVQRLRSKAQNKENFRPAS-SAEPAVQKLAS-LRLGGAEPRLLRWD . . . ::. .:. ::. .: . : . .: : : CCDS13 APSSP----ATAPLQAPTLRGPEPARET--PAKLEPKPIIPKAEPRAKARKTEARGLTKA 600 610 620 630 640 700 710 720 730 740 pF1KE6 LTFSPPQKSLPVALESDEENGDELKSSTGSAPILVVMVILLNIGVAILFINFFI . . .: .: :.. : .:. .. ::. ::::::::.::::.... CCDS13 GAKKKARKEAALAAEAEVEV-EEVPNT-----ILICMVILLNIGLAILFVHLLT 650 660 670 680 690 >>CCDS9603.1 JPH4 gene_id:84502|Hs108|chr14 (628 aa) initn: 1417 init1: 525 opt: 1543 Z-score: 1342.8 bits: 258.9 E(32554): 1.8e-68 Smith-Waterman score: 1620; 47.0% identity (69.4% similar) in 627 aa overlap (1-605:1-591) 10 20 30 40 50 60 pF1KE6 MSSGGRFNFDDGGSYCGGWEDGKAHGHGVCTGPKGQGEYTGSWSHGFEVLGVYTWPSGNT :: ::.:.::::: : :::: :.:::.:::::: .::::.: :.:::: :::.: :.:.. CCDS96 MSPGGKFDFDDGGCYVGGWEAGRAHGYGVCTGPGAQGEYSGCWAHGFESLGVFTGPGGHS 10 20 30 40 50 60 70 80 90 100 110 120 pF1KE6 YQGTWAQGKRHGIGLESKGKWVYKGEWTHGFKGRYGVRECAGNGAKYEGTWSNGLQDGYG ::: : ::::.:.:.: :..:.:.::: :.::: :: : . .: .: : :..:.::::: CCDS96 YQGHWQQGKREGLGVERKSRWTYRGEWLGGLKGRSGVWESV-SGLRYAGLWKDGFQDGYG 70 80 90 100 110 130 140 150 160 170 180 pF1KE6 TETYSDGGTYQGQWVGGMRQGYGVRQSVPYGMAAVIRSPLRTSINSLRSEHTNGTALHPD :::::::::::::: .: :.:::::::::: .::..::: :::..: .:. : : CCDS96 TETYSDGGTYQGQWQAGKRHGYGVRQSVPYHQAALLRSPRRTSLDSGHSDPP--TPPPPL 120 130 140 150 160 170 190 200 210 220 230 pF1KE6 ASPA-VAGSPAV-SRGGFVLVAHSDSEILKSKKK-----GLFRRSLL-SGLKLRKSESKS :. .:::: :::::::.. .:.. .:.:. :.:::::: ::: : . .: CCDS96 PLPGDEGGSPASGSRGGFVLAGPGDADGASSRKRTPAAGGFFRRSLLLSGL--RAGGRRS 180 190 200 210 220 230 240 250 260 270 280 290 pF1KE6 SLASQRSKQSSFRSEAGMSTVSSTASDIHSTISLGEAEAELAVIEDDIDATTTETYVGEW ::.:.:. :.:::.. : :.::. : : : : :.:: ...::.:.::: CCDS96 SLGSKRG---SLRSEVS-SEVGSTGPP-GSEASGPPAAAPPALIE----GSATEVYAGEW 240 250 260 270 280 300 310 320 330 340 350 pF1KE6 KNDKRSGFGVSQRSDGLKYEGEWASNRRHGYGCMTFPDGTKEEGKYKQNILV-GGKRKNL . :.::::::::::.::.::::: .::::::: : :::..::::::.: :: ::. ..: CCDS96 RADRRSGFGVSQRSNGLRYEGEWLGNRRHGYGRTTRPDGSREEGKYKRNRLVHGGRVRSL 290 300 310 320 330 340 360 370 380 390 400 pF1KE6 IPL--RASKIREKVDRAVEAAERAATTAKQKAEIAASRTSHSRAKAEAALTAAQKAQEEA .:: : .:..::::::::.:.::...:.:. ::::.:.. . :: :: ..:.:: : : CCDS96 LPLALRRGKVKEKVDRAVEGARRAVSAARQRQEIAAARAADALLKAVAASSVAEKAVEAA 350 360 370 380 390 400 410 420 430 440 450 460 pF1KE6 RIARITAKEFSPSFQHRENGLEYQRPKRQTSCDDIEVLSTGTPLQQESPELYRKGTTPSD :.:.. :....: .. : .::.... .: : ::...:: .:..: :::. CCDS96 RMAKLIAQDLQPMLE--APG---RRPRQDSEGSDTE------PLDEDSPGVYENGLTPSE 410 420 430 440 450 470 480 490 500 510 520 pF1KE6 LTTD--DSPLQSF-PTSPAATPPPAPAARNKVAHFSRQVSVDEERGGDIQMLLEGRAGDC . . .:: .: : : : : : . .. . :: . :: :: : .. CCDS96 GSPELPSSPASSRQPWRPPACRSPLPPGGDQ-GPFSSPKAWPEEWGG------AGAQAEE 460 470 480 490 500 530 540 550 560 570 pF1KE6 ARSSWGEEQAG-GSRGVRSGA-LRGGLLVDDFRTRGSGRKQPGN---PKPRERR---TES . .:..:: . : :.:. : :: .: . :: :.. :. : : : :: CCDS96 LAGYEAEDEAGMQGPGPRDGSPLLGG--CSD--SSGSLREEEGEDEEPLPPLRAPAGTEP 510 520 530 540 550 560 580 590 600 610 620 630 pF1KE6 PPVFTWTSHHRASNHSPGGSRLLELQEEKLSNYRMEMKPLLRMETHPQKRRYSKGGACRG :. . . .: .: :: : CCDS96 EPIAMLVLRGSSSRGPDAGCLTEELGEPAATERPAQPGAANPLVVGAVALLDLSLAFLFS 570 580 590 600 610 620 >>CCDS13326.1 JPH2 gene_id:57158|Hs108|chr20 (129 aa) initn: 824 init1: 824 opt: 824 Z-score: 730.3 bits: 143.3 E(32554): 2.3e-34 Smith-Waterman score: 824; 80.2% identity (96.0% similar) in 126 aa overlap (3-128:2-127) 10 20 30 40 50 60 pF1KE6 MSSGGRFNFDDGGSYCGGWEDGKAHGHGVCTGPKGQGEYTGSWSHGFEVLGVYTWPSGNT :::::.:::::.:::::: :::::::.::::::::::.:::. :::: :::::::::: CCDS13 MSGGRFDFDDGGAYCGGWEGGKAHGHGLCTGPKGQGEYSGSWNFGFEVAGVYTWPSGNT 10 20 30 40 50 70 80 90 100 110 120 pF1KE6 YQGTWAQGKRHGIGLESKGKWVYKGEWTHGFKGRYGVRECAGNGAKYEGTWSNGLQDGYG ..: :.::::::.:.:.::.:.::::::::::::::.:. ...::::::::.:::::::: CCDS13 FEGYWSQGKRHGLGIETKGRWLYKGEWTHGFKGRYGIRQSSSSGAKYEGTWNNGLQDGYG 60 70 80 90 100 110 130 140 150 160 170 180 pF1KE6 TETYSDGGTYQGQWVGGMRQGYGVRQSVPYGMAAVIRSPLRTSINSLRSEHTNGTALHPD ::::.::: CCDS13 TETYADGGMC 120 748 residues in 1 query sequences 18511270 residues in 32554 library sequences Tcomplib [36.3.4 Apr, 2011] (8 proc) start: Tue Nov 8 15:38:52 2016 done: Tue Nov 8 15:38:52 2016 Total Scan time: 2.990 Total Display time: 0.060 Function used was FASTA [36.3.4 Apr, 2011]