Result of FASTA (ccds) for
pF1KA1343
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: pF1KA1343, 495 aa
1>>>pF1KA1343 495 - 495 aa - 495 aa
Library: human.CCDS.faa
18511270 residues in 32554 sequences
Statistics: Expectation_n fit: rho(ln(x))= 13.3113+/-0.00111; mu= -17.4843+/- 0.067
mean_var=431.7403+/-87.916, 0's: 0 Z-trim(115.1): 16 B-trim: 0 in 0/53
Lambda= 0.061725
statistics sampled from 15613 (15623) to 15613 sequences
Algorithm: FASTA (3.7 Nov 2010) [optimized]
Parameters: BL50 matrix (15:-5), open/ext: -10/-2
ktup: 2, E-join: 1 (0.776), E-opt: 0.2 (0.48), width: 16
Scan time: 3.820
The best scores are: opt bits E(32554)
CCDS31016.1 RCOR3 gene_id:55758|Hs108|chr1 ( 495) 3362 313.4 3.7e-85
CCDS44312.1 RCOR3 gene_id:55758|Hs108|chr1 ( 553) 3362 313.4 4e-85
CCDS44314.1 RCOR3 gene_id:55758|Hs108|chr1 ( 449) 2557 241.7 1.3e-63
CCDS44313.1 RCOR3 gene_id:55758|Hs108|chr1 ( 436) 2047 196.3 5.9e-50
CCDS8052.1 RCOR2 gene_id:283248|Hs108|chr11 ( 523) 1674 163.1 6.8e-40
CCDS9974.2 RCOR1 gene_id:23186|Hs108|chr14 ( 485) 909 94.9 2.1e-19
>>CCDS31016.1 RCOR3 gene_id:55758|Hs108|chr1 (495 aa)
initn: 3362 init1: 3362 opt: 3362 Z-score: 1642.4 bits: 313.4 E(32554): 3.7e-85
Smith-Waterman score: 3362; 100.0% identity (100.0% similar) in 495 aa overlap (1-495:1-495)
10 20 30 40 50 60
pF1KA1 MRVGAEYQARIPEFDPGATKYTDKDNGGMLVWSPYHSIPDAKLDEYIAIAKEKHGYNVEQ
::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
CCDS31 MRVGAEYQARIPEFDPGATKYTDKDNGGMLVWSPYHSIPDAKLDEYIAIAKEKHGYNVEQ
10 20 30 40 50 60
70 80 90 100 110 120
pF1KA1 ALGMLFWHKHNIEKSLADLPNFTPFPDEWTVEDKVLFEQAFSFHGKSFHRIQQMLPDKTI
::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
CCDS31 ALGMLFWHKHNIEKSLADLPNFTPFPDEWTVEDKVLFEQAFSFHGKSFHRIQQMLPDKTI
70 80 90 100 110 120
130 140 150 160 170 180
pF1KA1 ASLVKYYYSWKKTRSRTSLMDRQARKLANRHNQGDSDDDVEETHPMDGNDSDYDPKKEAK
::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
CCDS31 ASLVKYYYSWKKTRSRTSLMDRQARKLANRHNQGDSDDDVEETHPMDGNDSDYDPKKEAK
130 140 150 160 170 180
190 200 210 220 230 240
pF1KA1 KEGNTEQPVQTSKIGLGRREYQSLQHRHHSQRSKCRPPKGMYLTQEDVVAVSCSPNAANT
::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
CCDS31 KEGNTEQPVQTSKIGLGRREYQSLQHRHHSQRSKCRPPKGMYLTQEDVVAVSCSPNAANT
190 200 210 220 230 240
250 260 270 280 290 300
pF1KA1 ILRQLDMELISLKRQVQNAKQVNSALKQKMEGGIEEFKPPESNQKINARWTTEEQLLAVQ
::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
CCDS31 ILRQLDMELISLKRQVQNAKQVNSALKQKMEGGIEEFKPPESNQKINARWTTEEQLLAVQ
250 260 270 280 290 300
310 320 330 340 350 360
pF1KA1 GVRKYGKDFQAIADVIGNKTVGQVKNFFVNYRRRFNLEEVLQEWEAEQGTQASNGDASTL
::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
CCDS31 GVRKYGKDFQAIADVIGNKTVGQVKNFFVNYRRRFNLEEVLQEWEAEQGTQASNGDASTL
310 320 330 340 350 360
370 380 390 400 410 420
pF1KA1 GEETKSASNVPSGKSTDEEEEAQTPQAPRTLGPSPPAPSSTPTPTAPIATLNQPPPLLRP
::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
CCDS31 GEETKSASNVPSGKSTDEEEEAQTPQAPRTLGPSPPAPSSTPTPTAPIATLNQPPPLLRP
370 380 390 400 410 420
430 440 450 460 470 480
pF1KA1 TLPAAPALHRQPPPLQQQARFIQPRPTLNQPPPPLIRPANSMPPRLNPRPVLSTVGGQQP
::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
CCDS31 TLPAAPALHRQPPPLQQQARFIQPRPTLNQPPPPLIRPANSMPPRLNPRPVLSTVGGQQP
430 440 450 460 470 480
490
pF1KA1 PSLIGIQTDSQSSLH
:::::::::::::::
CCDS31 PSLIGIQTDSQSSLH
490
>>CCDS44312.1 RCOR3 gene_id:55758|Hs108|chr1 (553 aa)
initn: 3362 init1: 3362 opt: 3362 Z-score: 1641.7 bits: 313.4 E(32554): 4e-85
Smith-Waterman score: 3362; 100.0% identity (100.0% similar) in 495 aa overlap (1-495:59-553)
10 20 30
pF1KA1 MRVGAEYQARIPEFDPGATKYTDKDNGGML
::::::::::::::::::::::::::::::
CCDS44 GGSGASSTNGGLHYSEPESGCSSDDEHDVGMRVGAEYQARIPEFDPGATKYTDKDNGGML
30 40 50 60 70 80
40 50 60 70 80 90
pF1KA1 VWSPYHSIPDAKLDEYIAIAKEKHGYNVEQALGMLFWHKHNIEKSLADLPNFTPFPDEWT
::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
CCDS44 VWSPYHSIPDAKLDEYIAIAKEKHGYNVEQALGMLFWHKHNIEKSLADLPNFTPFPDEWT
90 100 110 120 130 140
100 110 120 130 140 150
pF1KA1 VEDKVLFEQAFSFHGKSFHRIQQMLPDKTIASLVKYYYSWKKTRSRTSLMDRQARKLANR
::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
CCDS44 VEDKVLFEQAFSFHGKSFHRIQQMLPDKTIASLVKYYYSWKKTRSRTSLMDRQARKLANR
150 160 170 180 190 200
160 170 180 190 200 210
pF1KA1 HNQGDSDDDVEETHPMDGNDSDYDPKKEAKKEGNTEQPVQTSKIGLGRREYQSLQHRHHS
::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
CCDS44 HNQGDSDDDVEETHPMDGNDSDYDPKKEAKKEGNTEQPVQTSKIGLGRREYQSLQHRHHS
210 220 230 240 250 260
220 230 240 250 260 270
pF1KA1 QRSKCRPPKGMYLTQEDVVAVSCSPNAANTILRQLDMELISLKRQVQNAKQVNSALKQKM
::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
CCDS44 QRSKCRPPKGMYLTQEDVVAVSCSPNAANTILRQLDMELISLKRQVQNAKQVNSALKQKM
270 280 290 300 310 320
280 290 300 310 320 330
pF1KA1 EGGIEEFKPPESNQKINARWTTEEQLLAVQGVRKYGKDFQAIADVIGNKTVGQVKNFFVN
::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
CCDS44 EGGIEEFKPPESNQKINARWTTEEQLLAVQGVRKYGKDFQAIADVIGNKTVGQVKNFFVN
330 340 350 360 370 380
340 350 360 370 380 390
pF1KA1 YRRRFNLEEVLQEWEAEQGTQASNGDASTLGEETKSASNVPSGKSTDEEEEAQTPQAPRT
::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
CCDS44 YRRRFNLEEVLQEWEAEQGTQASNGDASTLGEETKSASNVPSGKSTDEEEEAQTPQAPRT
390 400 410 420 430 440
400 410 420 430 440 450
pF1KA1 LGPSPPAPSSTPTPTAPIATLNQPPPLLRPTLPAAPALHRQPPPLQQQARFIQPRPTLNQ
::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
CCDS44 LGPSPPAPSSTPTPTAPIATLNQPPPLLRPTLPAAPALHRQPPPLQQQARFIQPRPTLNQ
450 460 470 480 490 500
460 470 480 490
pF1KA1 PPPPLIRPANSMPPRLNPRPVLSTVGGQQPPSLIGIQTDSQSSLH
:::::::::::::::::::::::::::::::::::::::::::::
CCDS44 PPPPLIRPANSMPPRLNPRPVLSTVGGQQPPSLIGIQTDSQSSLH
510 520 530 540 550
>>CCDS44314.1 RCOR3 gene_id:55758|Hs108|chr1 (449 aa)
initn: 2557 init1: 2557 opt: 2557 Z-score: 1255.6 bits: 241.7 E(32554): 1.3e-63
Smith-Waterman score: 2557; 100.0% identity (100.0% similar) in 381 aa overlap (1-381:59-439)
10 20 30
pF1KA1 MRVGAEYQARIPEFDPGATKYTDKDNGGML
::::::::::::::::::::::::::::::
CCDS44 GGSGASSTNGGLHYSEPESGCSSDDEHDVGMRVGAEYQARIPEFDPGATKYTDKDNGGML
30 40 50 60 70 80
40 50 60 70 80 90
pF1KA1 VWSPYHSIPDAKLDEYIAIAKEKHGYNVEQALGMLFWHKHNIEKSLADLPNFTPFPDEWT
::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
CCDS44 VWSPYHSIPDAKLDEYIAIAKEKHGYNVEQALGMLFWHKHNIEKSLADLPNFTPFPDEWT
90 100 110 120 130 140
100 110 120 130 140 150
pF1KA1 VEDKVLFEQAFSFHGKSFHRIQQMLPDKTIASLVKYYYSWKKTRSRTSLMDRQARKLANR
::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
CCDS44 VEDKVLFEQAFSFHGKSFHRIQQMLPDKTIASLVKYYYSWKKTRSRTSLMDRQARKLANR
150 160 170 180 190 200
160 170 180 190 200 210
pF1KA1 HNQGDSDDDVEETHPMDGNDSDYDPKKEAKKEGNTEQPVQTSKIGLGRREYQSLQHRHHS
::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
CCDS44 HNQGDSDDDVEETHPMDGNDSDYDPKKEAKKEGNTEQPVQTSKIGLGRREYQSLQHRHHS
210 220 230 240 250 260
220 230 240 250 260 270
pF1KA1 QRSKCRPPKGMYLTQEDVVAVSCSPNAANTILRQLDMELISLKRQVQNAKQVNSALKQKM
::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
CCDS44 QRSKCRPPKGMYLTQEDVVAVSCSPNAANTILRQLDMELISLKRQVQNAKQVNSALKQKM
270 280 290 300 310 320
280 290 300 310 320 330
pF1KA1 EGGIEEFKPPESNQKINARWTTEEQLLAVQGVRKYGKDFQAIADVIGNKTVGQVKNFFVN
::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
CCDS44 EGGIEEFKPPESNQKINARWTTEEQLLAVQGVRKYGKDFQAIADVIGNKTVGQVKNFFVN
330 340 350 360 370 380
340 350 360 370 380 390
pF1KA1 YRRRFNLEEVLQEWEAEQGTQASNGDASTLGEETKSASNVPSGKSTDEEEEAQTPQAPRT
:::::::::::::::::::::::::::::::::::::::::::::::::::
CCDS44 YRRRFNLEEVLQEWEAEQGTQASNGDASTLGEETKSASNVPSGKSTDEEEEVCLCMEFEL
390 400 410 420 430 440
400 410 420 430 440 450
pF1KA1 LGPSPPAPSSTPTPTAPIATLNQPPPLLRPTLPAAPALHRQPPPLQQQARFIQPRPTLNQ
CCDS44 I
>>CCDS44313.1 RCOR3 gene_id:55758|Hs108|chr1 (436 aa)
initn: 2073 init1: 2044 opt: 2047 Z-score: 1010.4 bits: 196.3 E(32554): 5.9e-50
Smith-Waterman score: 2047; 98.7% identity (99.0% similar) in 306 aa overlap (1-306:59-364)
10 20 30
pF1KA1 MRVGAEYQARIPEFDPGATKYTDKDNGGML
::::::::::::::::::::::::::::::
CCDS44 GGSGASSTNGGLHYSEPESGCSSDDEHDVGMRVGAEYQARIPEFDPGATKYTDKDNGGML
30 40 50 60 70 80
40 50 60 70 80 90
pF1KA1 VWSPYHSIPDAKLDEYIAIAKEKHGYNVEQALGMLFWHKHNIEKSLADLPNFTPFPDEWT
::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
CCDS44 VWSPYHSIPDAKLDEYIAIAKEKHGYNVEQALGMLFWHKHNIEKSLADLPNFTPFPDEWT
90 100 110 120 130 140
100 110 120 130 140 150
pF1KA1 VEDKVLFEQAFSFHGKSFHRIQQMLPDKTIASLVKYYYSWKKTRSRTSLMDRQARKLANR
::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
CCDS44 VEDKVLFEQAFSFHGKSFHRIQQMLPDKTIASLVKYYYSWKKTRSRTSLMDRQARKLANR
150 160 170 180 190 200
160 170 180 190 200 210
pF1KA1 HNQGDSDDDVEETHPMDGNDSDYDPKKEAKKEGNTEQPVQTSKIGLGRREYQSLQHRHHS
::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
CCDS44 HNQGDSDDDVEETHPMDGNDSDYDPKKEAKKEGNTEQPVQTSKIGLGRREYQSLQHRHHS
210 220 230 240 250 260
220 230 240 250 260 270
pF1KA1 QRSKCRPPKGMYLTQEDVVAVSCSPNAANTILRQLDMELISLKRQVQNAKQVNSALKQKM
::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
CCDS44 QRSKCRPPKGMYLTQEDVVAVSCSPNAANTILRQLDMELISLKRQVQNAKQVNSALKQKM
270 280 290 300 310 320
280 290 300 310 320 330
pF1KA1 EGGIEEFKPPESNQKINARWTTEEQLLAVQGVRKYGKDFQAIADVIGNKTVGQVKNFFVN
:::::::::::::::::::::::::::::::. :
CCDS44 EGGIEEFKPPESNQKINARWTTEEQLLAVQGTDPTGSSDTGSITSCPIIHSNTNSPYCHS
330 340 350 360 370 380
340 350 360 370 380 390
pF1KA1 YRRRFNLEEVLQEWEAEQGTQASNGDASTLGEETKSASNVPSGKSTDEEEEAQTPQAPRT
CCDS44 EPASTTSSSNTACCPGSSPAASSTPAAGSVHPAPANFKSASTTSYSPC
390 400 410 420 430
>>CCDS8052.1 RCOR2 gene_id:283248|Hs108|chr11 (523 aa)
initn: 1570 init1: 756 opt: 1674 Z-score: 829.7 bits: 163.1 E(32554): 6.8e-40
Smith-Waterman score: 1898; 60.3% identity (80.2% similar) in 491 aa overlap (1-485:46-513)
10 20
pF1KA1 MRVGAEYQARIPEFDP-GATKYTDKDNGGM
.:::..::: ::: : . ..:..:. ::
CCDS80 SRSRAKTVPNGGQPHSEDDSSEEEHSHDSMIRVGTNYQAVIPECKPESPARYSNKELKGM
20 30 40 50 60 70
30 40 50 60 70 80
pF1KA1 LVWSPYHSIPDAKLDEYIAIAKEKHGYNVEQALGMLFWHKHNIEKSLADLPNFTPFPDEW
::::: : . :::::.:::.::::::::.:::::::.::::..::::::: :::::::::
CCDS80 LVWSPNHCVSDAKLDKYIAMAKEKHGYNIEQALGMLLWHKHDVEKSLADLANFTPFPDEW
80 90 100 110 120 130
90 100 110 120 130 140
pF1KA1 TVEDKVLFEQAFSFHGKSFHRIQQMLPDKTIASLVKYYYSWKKTRSRTSLMDRQARKLAN
::::::::::::.:::: :.::::::::: : :::::::::::::::::.::::::.:..
CCDS80 TVEDKVLFEQAFGFHGKCFQRIQQMLPDKLIPSLVKYYYSWKKTRSRTSVMDRQARRLGG
140 150 160 170 180 190
150 160 170 180 190 200
pF1KA1 RHNQGDSDDDVEETHPMDGNDSDYDPKK-EAKKEGNTEQPVQTSKIGLGRREYQSLQHRH
:... :::. .:: . :. :. .: . :.: .:.. .. : :..: : :.::
CCDS80 RKDKEDSDE-LEEGR---GGVSEGEPDPADPKREPLPSRPLN-ARPGPGKKEVQVSQYRH
200 210 220 230 240 250
210 220 230 240 250 260
pF1KA1 HSQRSKCRPPKGMYLTQEDVVAVSCSPNAANTILRQLDMELISLKRQVQNAKQVNSALKQ
: :.. ::::::::. : ..::: ::. :: :: :: .:::::::::. ::.::.:.:
CCDS80 HPLRTRRRPPKGMYLSPEGLTAVSGSPDLANLTLRGLDSQLISLKRQVQSMKQTNSSLRQ
260 270 280 290 300 310
270 280 290 300 310 320
pF1KA1 KMEGGIEEFKPPESNQKINARWTTEEQLLAVQGVRKYGKDFQAIADVIGNKTVGQVKNFF
.::::. ..:::.: :.:.::::.:::::::..:.::::: :::.::::::. :::.::
CCDS80 ALEGGIDPLRPPEANTKFNSRWTTDEQLLAVQAIRRYGKDFGAIAEVIGNKTLTQVKTFF
320 330 340 350 360 370
330 340 350 360 370 380
pF1KA1 VNYRRRFNLEEVLQEWEAEQGTQASNGDASTLGEETKSASNVPSGKSTDEEEEAQ----T
:.:::::::::::::::::: ... : : . ::.. .. .:. . .:..:.: .
CCDS80 VSYRRRFNLEEVLQEWEAEQ--DGAPG-APVPMEEARRGAPLPA-PALEEDDEVQITSVS
380 390 400 410 420
390 400 410 420 430 440
pF1KA1 PQAPRTLGPSPPAPSSTPTPTAPIATLNQPPPLLRPTLPAAPALHRQPPPLQQQARFIQP
..::.. :.:: : : :: .:.:::::::: ::.::.: :::::::: .::.::
CCDS80 TSVPRSVPPAPPPP---PPPT----SLSQPPPLLRPPLPTAPTLLRQPPPLQQ-GRFLQP
430 440 450 460 470
450 460 470 480 490
pF1KA1 RPTLNQPPPPLIRPANSMPPRLNPRPVLSTVGGQQPPSLIGIQTDSQSSLH
: . ::::::::::: . : : . :: : : ::.:::
CCDS80 RLAPNQPPPPLIRPALAAP-RHSARP-----GPQPPPTLIGTPLEPPAPSL
480 490 500 510 520
>>CCDS9974.2 RCOR1 gene_id:23186|Hs108|chr14 (485 aa)
initn: 1567 init1: 759 opt: 909 Z-score: 462.0 bits: 94.9 E(32554): 2.1e-19
Smith-Waterman score: 1616; 64.6% identity (84.4% similar) in 384 aa overlap (1-382:105-474)
10 20
pF1KA1 MRVGAEYQARIPEFDPG--ATKYTDKDNGG
:::: .::: .:.:::. : . ..:: :
CCDS99 AAAPNGNSSSNSWEEGSSGSSSDEEHGGGGMRVGPQYQAVVPDFDPAKLARRSQERDNLG
80 90 100 110 120 130
30 40 50 60 70 80
pF1KA1 MLVWSPYHSIPDAKLDEYIAIAKEKHGYNVEQALGMLFWHKHNIEKSLADLPNFTPFPDE
:::::: ... .:::::::::::::::::.::::::::::::::::::::::::::::::
CCDS99 MLVWSPNQNLSEAKLDEYIAIAKEKHGYNMEQALGMLFWHKHNIEKSLADLPNFTPFPDE
140 150 160 170 180 190
90 100 110 120 130 140
pF1KA1 WTVEDKVLFEQAFSFHGKSFHRIQQMLPDKTIASLVKYYYSWKKTRSRTSLMDRQARKLA
::::::::::::::::::.:::::::::::.::::::.::::::::..::.:::.:::
CCDS99 WTVEDKVLFEQAFSFHGKTFHRIQQMLPDKSIASLVKFYYSWKKTRTKTSVMDRHARK--
200 210 220 230 240 250
150 160 170 180 190 200
pF1KA1 NRHNQGDSDDDVEETHPMDGNDSDYDPKKEAKKEGNTEQPVQTSKIGLGRREYQSLQHRH
..... .:.:..::.. . : . : .::.::: :..: . ..: .: : ..
CCDS99 QKREREESEDELEEANGNNPIDIEVDQNKESKKE---VPPTET--VPQVKKEKHSTQAKN
260 270 280 290 300
210 220 230 240 250 260
pF1KA1 HSQRSKCRPPKGMYLTQEDVVAVSCSPNAANTILRQLDMELISLKRQVQNAKQVNSALKQ
:.: .:::::.:.:::: ::: . .::.:.::::::::.:.:::.:: ::.:::::.
CCDS99 ---RAKRKPPKGMFLSQEDVEAVSANATAATTVLRQLDMELVSVKRQIQNIKQTNSALKE
310 320 330 340 350 360
270 280 290 300 310 320
pF1KA1 KMEGGIEEFKPPESNQKINARWTTEEQLLAVQGVRKYGKDFQAIADVIGNKTVGQVKNFF
:..:::: .. :: :: ::::::::::::::..::::.:::::.::::::.: ::::::
CCDS99 KLDGGIEPYRLPEVIQKCNARWTTEEQLLAVQAIRKYGRDFQAISDVIGNKSVVQVKNFF
370 380 390 400 410 420
330 340 350 360 370 380
pF1KA1 VNYRRRFNLEEVLQEWEAEQGTQASNGDASTLGEETKSASNVPSGKSTDEEEEAQTPQAP
::::::::..:::::::::.: . .:: .. . .:: .: : : .::.::
CCDS99 VNYRRRFNIDEVLQEWEAEHGKEETNGPSNQ--KPVKSPDN--SIKMPEEEDEAPVLDVR
430 440 450 460 470 480
390 400 410 420 430 440
pF1KA1 RTLGPSPPAPSSTPTPTAPIATLNQPPPLLRPTLPAAPALHRQPPPLQQQARFIQPRPTL
CCDS99 YASAS
495 residues in 1 query sequences
18511270 residues in 32554 library sequences
Tcomplib [36.3.4 Apr, 2011] (8 proc)
start: Fri Nov 4 01:29:21 2016 done: Fri Nov 4 01:29:21 2016
Total Scan time: 3.820 Total Display time: 0.040
Function used was FASTA [36.3.4 Apr, 2011]