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: pF1KE5208, 244 aa 1>>>pF1KE5208 244 - 244 aa - 244 aa Library: /omim/omim.rfq.tfa 60827320 residues in 85289 sequences Statistics: Expectation_n fit: rho(ln(x))= 6.0015+/-0.000355; mu= 11.3651+/- 0.022 mean_var=82.1044+/-15.813, 0's: 0 Z-trim(114.4): 49 B-trim: 0 in 0/52 Lambda= 0.141544 statistics sampled from 24125 (24174) to 24125 sequences Algorithm: FASTA (3.7 Nov 2010) [optimized] Parameters: BL50 matrix (15:-5), open/ext: -10/-2 ktup: 2, E-join: 1 (0.659), E-opt: 0.2 (0.283), width: 16 Scan time: 4.310 The best scores are: opt bits E(85289) NP_001272 (OMIM: 601303) tubulin-folding cofactor ( 244) 1610 338.2 7.8e-93 NP_001287900 (OMIM: 601303) tubulin-folding cofact ( 193) 1294 273.6 1.7e-73 NP_001128512 (OMIM: 105400,168605,601143,607641) d (1253) 241 59.0 4.4e-08 NP_001177766 (OMIM: 105400,168605,601143,607641) d (1271) 241 59.0 4.5e-08 NP_004073 (OMIM: 105400,168605,601143,607641) dyna (1278) 241 59.0 4.5e-08 NP_001177765 (OMIM: 105400,168605,601143,607641) d (1236) 239 58.6 5.8e-08 XP_016875281 (OMIM: 179838) PREDICTED: CAP-Gly dom (1355) 239 58.6 6.3e-08 XP_016875280 (OMIM: 179838) PREDICTED: CAP-Gly dom (1390) 239 58.6 6.4e-08 NP_937883 (OMIM: 179838) CAP-Gly domain-containing (1392) 239 58.6 6.5e-08 XP_016875279 (OMIM: 179838) PREDICTED: CAP-Gly dom (1392) 239 58.6 6.5e-08 XP_016875277 (OMIM: 179838) PREDICTED: CAP-Gly dom (1394) 239 58.6 6.5e-08 XP_016875278 (OMIM: 179838) PREDICTED: CAP-Gly dom (1394) 239 58.6 6.5e-08 XP_016875276 (OMIM: 179838) PREDICTED: CAP-Gly dom (1401) 239 58.6 6.5e-08 XP_016875275 (OMIM: 179838) PREDICTED: CAP-Gly dom (1427) 239 58.6 6.6e-08 NP_002947 (OMIM: 179838) CAP-Gly domain-containing (1427) 239 58.6 6.6e-08 XP_006719615 (OMIM: 179838) PREDICTED: CAP-Gly dom (1429) 239 58.6 6.6e-08 NP_001234926 (OMIM: 179838) CAP-Gly domain-contain (1438) 239 58.6 6.6e-08 XP_005253650 (OMIM: 179838) PREDICTED: CAP-Gly dom (1440) 239 58.6 6.6e-08 XP_011542761 (OMIM: 607350) PREDICTED: kinesin-lik (1783) 240 58.9 6.9e-08 XP_011542762 (OMIM: 607350) PREDICTED: kinesin-lik (1783) 240 58.9 6.9e-08 XP_011542763 (OMIM: 607350) PREDICTED: kinesin-lik (1783) 240 58.9 6.9e-08 XP_005273515 (OMIM: 607350) PREDICTED: kinesin-lik (1805) 240 58.9 7e-08 NP_056069 (OMIM: 607350) kinesin-like protein KIF1 (1826) 240 58.9 7.1e-08 XP_011542760 (OMIM: 607350) PREDICTED: kinesin-lik (1826) 240 58.9 7.1e-08 XP_011542759 (OMIM: 607350) PREDICTED: kinesin-lik (1847) 240 58.9 7.1e-08 NP_115797 (OMIM: 603432) CAP-Gly domain-containing (1011) 229 56.5 2e-07 NP_003379 (OMIM: 603432) CAP-Gly domain-containing (1046) 229 56.5 2.1e-07 NP_056341 (OMIM: 607382) CAP-Gly domain-containing ( 547) 202 50.9 5.5e-06 NP_001186499 (OMIM: 607382) CAP-Gly domain-contain ( 547) 202 50.9 5.5e-06 NP_003184 (OMIM: 241410,244460,604934) tubulin-spe ( 527) 155 41.3 0.0041 NP_001072983 (OMIM: 241410,244460,604934) tubulin- ( 527) 155 41.3 0.0041 NP_001274730 (OMIM: 241410,244460,604934) tubulin- ( 578) 155 41.3 0.0045 >>NP_001272 (OMIM: 601303) tubulin-folding cofactor B is (244 aa) initn: 1610 init1: 1610 opt: 1610 Z-score: 1787.7 bits: 338.2 E(85289): 7.8e-93 Smith-Waterman score: 1610; 100.0% identity (100.0% similar) in 244 aa overlap (1-244:1-244) 10 20 30 40 50 60 pF1KE5 MEVTGVSAPTVTVFISSSLNTFRSEKRYSRSLTIAEFKCKLELLVGSPASCMELELYGVD :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: NP_001 MEVTGVSAPTVTVFISSSLNTFRSEKRYSRSLTIAEFKCKLELLVGSPASCMELELYGVD 10 20 30 40 50 60 70 80 90 100 110 120 pF1KE5 DKFYSKLDQEDALLGSYPVDDGCRIHVIDHSGARLGEYEDVSRVEKYTISQEAYDQRQDT :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: NP_001 DKFYSKLDQEDALLGSYPVDDGCRIHVIDHSGARLGEYEDVSRVEKYTISQEAYDQRQDT 70 80 90 100 110 120 130 140 150 160 170 180 pF1KE5 VRSFLKRSKLGRYNEEERAQQEAEAAQRLAEEKAQASSIPVGSRCEVRAAGQSPRRGTVM :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: NP_001 VRSFLKRSKLGRYNEEERAQQEAEAAQRLAEEKAQASSIPVGSRCEVRAAGQSPRRGTVM 130 140 150 160 170 180 190 200 210 220 230 240 pF1KE5 YVGLTDFKPGYWIGVRYDEPLGKNDGSVNGKRYFECQAKYGAFVKPAVVTVGDFPEEDYG :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: NP_001 YVGLTDFKPGYWIGVRYDEPLGKNDGSVNGKRYFECQAKYGAFVKPAVVTVGDFPEEDYG 190 200 210 220 230 240 pF1KE5 LDEI :::: NP_001 LDEI >>NP_001287900 (OMIM: 601303) tubulin-folding cofactor B (193 aa) initn: 1294 init1: 1294 opt: 1294 Z-score: 1440.5 bits: 273.6 E(85289): 1.7e-73 Smith-Waterman score: 1294; 100.0% identity (100.0% similar) in 193 aa overlap (52-244:1-193) 30 40 50 60 70 80 pF1KE5 FRSEKRYSRSLTIAEFKCKLELLVGSPASCMELELYGVDDKFYSKLDQEDALLGSYPVDD :::::::::::::::::::::::::::::: NP_001 MELELYGVDDKFYSKLDQEDALLGSYPVDD 10 20 30 90 100 110 120 130 140 pF1KE5 GCRIHVIDHSGARLGEYEDVSRVEKYTISQEAYDQRQDTVRSFLKRSKLGRYNEEERAQQ :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: NP_001 GCRIHVIDHSGARLGEYEDVSRVEKYTISQEAYDQRQDTVRSFLKRSKLGRYNEEERAQQ 40 50 60 70 80 90 150 160 170 180 190 200 pF1KE5 EAEAAQRLAEEKAQASSIPVGSRCEVRAAGQSPRRGTVMYVGLTDFKPGYWIGVRYDEPL :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: NP_001 EAEAAQRLAEEKAQASSIPVGSRCEVRAAGQSPRRGTVMYVGLTDFKPGYWIGVRYDEPL 100 110 120 130 140 150 210 220 230 240 pF1KE5 GKNDGSVNGKRYFECQAKYGAFVKPAVVTVGDFPEEDYGLDEI ::::::::::::::::::::::::::::::::::::::::::: NP_001 GKNDGSVNGKRYFECQAKYGAFVKPAVVTVGDFPEEDYGLDEI 160 170 180 190 >>NP_001128512 (OMIM: 105400,168605,601143,607641) dynac (1253 aa) initn: 235 init1: 206 opt: 241 Z-score: 266.0 bits: 59.0 E(85289): 4.4e-08 Smith-Waterman score: 241; 46.0% identity (67.8% similar) in 87 aa overlap (145-231:14-96) 120 130 140 150 160 170 pF1KE5 DQRQDTVRSFLKRSKLGRYNEEERAQQEAEAAQRLAEEKAQASSIPVGSRCEVRAAGQSP ...:.. : :.: . :::: :: . :. NP_001 MAQSKRHVYSRTPSGSRMSAE-ASARPLRVGSRVEVIGKGH-- 10 20 30 40 180 190 200 210 220 230 pF1KE5 RRGTVMYVGLTDFKPGYWIGVRYDEPLGKNDGSVNGKRYFECQAKYGAFVKPAVVTVGDF :::: ::: : : : :.:: :: :::::.:.:..:: :. .: ::. . . : NP_001 -RGTVAYVGATLFATGKWVGVILDEAKGKNDGTVQGRKYFTCDEGHGIFVRQSQIQVFED 50 60 70 80 90 240 pF1KE5 PEEDYGLDEI NP_001 GADTTSPETPDSSASKVLKREGTDTTAKTSKLPTRPASTGVAGASSSLGPSGSASAGELS 100 110 120 130 140 150 >>NP_001177766 (OMIM: 105400,168605,601143,607641) dynac (1271 aa) initn: 235 init1: 206 opt: 241 Z-score: 265.9 bits: 59.0 E(85289): 4.5e-08 Smith-Waterman score: 241; 46.0% identity (67.8% similar) in 87 aa overlap (145-231:14-96) 120 130 140 150 160 170 pF1KE5 DQRQDTVRSFLKRSKLGRYNEEERAQQEAEAAQRLAEEKAQASSIPVGSRCEVRAAGQSP ...:.. : :.: . :::: :: . :. NP_001 MAQSKRHVYSRTPSGSRMSAE-ASARPLRVGSRVEVIGKGH-- 10 20 30 40 180 190 200 210 220 230 pF1KE5 RRGTVMYVGLTDFKPGYWIGVRYDEPLGKNDGSVNGKRYFECQAKYGAFVKPAVVTVGDF :::: ::: : : : :.:: :: :::::.:.:..:: :. .: ::. . . : NP_001 -RGTVAYVGATLFATGKWVGVILDEAKGKNDGTVQGRKYFTCDEGHGIFVRQSQIQVFED 50 60 70 80 90 240 pF1KE5 PEEDYGLDEI NP_001 GADTTSPETPDSSASKVLKREGTDTTAKTSKLAPTARKTTTRRPKPTRPASTGVAGASSS 100 110 120 130 140 150 >>NP_004073 (OMIM: 105400,168605,601143,607641) dynactin (1278 aa) initn: 235 init1: 206 opt: 241 Z-score: 265.9 bits: 59.0 E(85289): 4.5e-08 Smith-Waterman score: 241; 46.0% identity (67.8% similar) in 87 aa overlap (145-231:14-96) 120 130 140 150 160 170 pF1KE5 DQRQDTVRSFLKRSKLGRYNEEERAQQEAEAAQRLAEEKAQASSIPVGSRCEVRAAGQSP ...:.. : :.: . :::: :: . :. NP_004 MAQSKRHVYSRTPSGSRMSAE-ASARPLRVGSRVEVIGKGH-- 10 20 30 40 180 190 200 210 220 230 pF1KE5 RRGTVMYVGLTDFKPGYWIGVRYDEPLGKNDGSVNGKRYFECQAKYGAFVKPAVVTVGDF :::: ::: : : : :.:: :: :::::.:.:..:: :. .: ::. . . : NP_004 -RGTVAYVGATLFATGKWVGVILDEAKGKNDGTVQGRKYFTCDEGHGIFVRQSQIQVFED 50 60 70 80 90 240 pF1KE5 PEEDYGLDEI NP_004 GADTTSPETPDSSASKVLKREGTDTTAKTSKLRGLKPKKAPTARKTTTRRPKPTRPASTG 100 110 120 130 140 150 >>NP_001177765 (OMIM: 105400,168605,601143,607641) dynac (1236 aa) initn: 235 init1: 206 opt: 239 Z-score: 263.9 bits: 58.6 E(85289): 5.8e-08 Smith-Waterman score: 239; 50.0% identity (67.5% similar) in 80 aa overlap (154-231:3-79) 130 140 150 160 170 180 pF1KE5 FLKRSKLGRYNEEERAQQEAEAAQRLAEEKAQASSIP--VGSRCEVRAAGQSPRRGTVMY :.::. : :::: :: . :. :::: : NP_001 MSAEASARPLRVGSRVEVIGKGH---RGTVAY 10 20 190 200 210 220 230 240 pF1KE5 VGLTDFKPGYWIGVRYDEPLGKNDGSVNGKRYFECQAKYGAFVKPAVVTVGDFPEEDYGL :: : : : :.:: :: :::::.:.:..:: :. .: ::. . . : NP_001 VGATLFATGKWVGVILDEAKGKNDGTVQGRKYFTCDEGHGIFVRQSQIQVFEDGADTTSP 30 40 50 60 70 80 pF1KE5 DEI NP_001 ETPDSSASKVLKREGTDTTAKTSKLPTRPASTGVAGASSSLGPSGSASAGELSSSEPSTP 90 100 110 120 130 140 >>XP_016875281 (OMIM: 179838) PREDICTED: CAP-Gly domain- (1355 aa) initn: 233 init1: 233 opt: 239 Z-score: 263.3 bits: 58.6 E(85289): 6.3e-08 Smith-Waterman score: 239; 46.5% identity (65.1% similar) in 86 aa overlap (154-235:200-284) 130 140 150 160 170 pF1KE5 FLKRSKLGRYNEEERAQQEAEAAQRLAEEKAQASSIPVGSRCEVRAAGQ----SPRRGTV ..:.:: : : :.. . . . . :.: XP_016 KPSQPAAKEPSATPPISNLTKTASESISNLSEAGSIKKGER-ELKIGDRVLVGGTKAGVV 170 180 190 200 210 220 180 190 200 210 220 230 pF1KE5 MYVGLTDFKPGYWIGVRYDEPLGKNDGSVNGKRYFECQAKYGAFVKPAVVTVGDFPEEDY ..: ::: : : ::. :::::::::.: : :::.:: ::: :. :: :: XP_016 RFLGETDFAKGEWCGVELDEPLGKNDGAVAGTRYFQCQPKYGLFAPVHKVTKIGFPSTTP 230 240 250 260 270 280 240 pF1KE5 GLDEI XP_016 AKAKANAVRRVMATTSASLKRSPSASSLSSMSSVASSVSSRPSRTGLLTETSSRYARKIS 290 300 310 320 330 340 >-- initn: 199 init1: 156 opt: 224 Z-score: 246.7 bits: 55.6 E(85289): 5.3e-07 Smith-Waterman score: 224; 47.1% identity (70.0% similar) in 70 aa overlap (161-230:60-125) 140 150 160 170 180 190 pF1KE5 GRYNEEERAQQEAEAAQRLAEEKAQASSIPVGSRCEVRAAGQSPRRGTVMYVGLTDFKPG :: : : . :..: : ....: :.: :: XP_016 VAPVEKTISSEKASSTPSSETQEEFVDDFRVGER--VWVNGNKP--GFIQFLGETQFAPG 30 40 50 60 70 80 200 210 220 230 240 pF1KE5 YWIGVRYDEPLGKNDGSVNGKRYFECQAKYGAFVKPAVVTVGDFPEEDYGLDEI : :. :::.::::::: : :::.:. : :..:. .: XP_016 QWAGIVLDEPIGKNDGSVAGVRYFQCEPLKGIFTRPSKLTRKVQAEDEANGLQTTPASRA 90 100 110 120 130 140 XP_016 TSPLCTSTASMVSSSPSTPSNIPQKPSQPAAKEPSATPPISNLTKTASESISNLSEAGSI 150 160 170 180 190 200 >>XP_016875280 (OMIM: 179838) PREDICTED: CAP-Gly domain- (1390 aa) initn: 233 init1: 233 opt: 239 Z-score: 263.1 bits: 58.6 E(85289): 6.4e-08 Smith-Waterman score: 239; 46.5% identity (65.1% similar) in 86 aa overlap (154-235:200-284) 130 140 150 160 170 pF1KE5 FLKRSKLGRYNEEERAQQEAEAAQRLAEEKAQASSIPVGSRCEVRAAGQ----SPRRGTV ..:.:: : : :.. . . . . :.: XP_016 KPSQPAAKEPSATPPISNLTKTASESISNLSEAGSIKKGER-ELKIGDRVLVGGTKAGVV 170 180 190 200 210 220 180 190 200 210 220 230 pF1KE5 MYVGLTDFKPGYWIGVRYDEPLGKNDGSVNGKRYFECQAKYGAFVKPAVVTVGDFPEEDY ..: ::: : : ::. :::::::::.: : :::.:: ::: :. :: :: XP_016 RFLGETDFAKGEWCGVELDEPLGKNDGAVAGTRYFQCQPKYGLFAPVHKVTKIGFPSTTP 230 240 250 260 270 280 240 pF1KE5 GLDEI XP_016 AKAKANAVRRVMATTSASLKRSPSASSLSSMSSVASSVSSRPSRTGLLTETSSRYARKIS 290 300 310 320 330 340 >-- initn: 199 init1: 156 opt: 224 Z-score: 246.6 bits: 55.6 E(85289): 5.4e-07 Smith-Waterman score: 224; 47.1% identity (70.0% similar) in 70 aa overlap (161-230:60-125) 140 150 160 170 180 190 pF1KE5 GRYNEEERAQQEAEAAQRLAEEKAQASSIPVGSRCEVRAAGQSPRRGTVMYVGLTDFKPG :: : : . :..: : ....: :.: :: XP_016 VAPVEKTISSEKASSTPSSETQEEFVDDFRVGER--VWVNGNKP--GFIQFLGETQFAPG 30 40 50 60 70 80 200 210 220 230 240 pF1KE5 YWIGVRYDEPLGKNDGSVNGKRYFECQAKYGAFVKPAVVTVGDFPEEDYGLDEI : :. :::.::::::: : :::.:. : :..:. .: XP_016 QWAGIVLDEPIGKNDGSVAGVRYFQCEPLKGIFTRPSKLTRKVQAEDEANGLQTTPASRA 90 100 110 120 130 140 XP_016 TSPLCTSTASMVSSSPSTPSNIPQKPSQPAAKEPSATPPISNLTKTASESISNLSEAGSI 150 160 170 180 190 200 >>NP_937883 (OMIM: 179838) CAP-Gly domain-containing lin (1392 aa) initn: 233 init1: 233 opt: 239 Z-score: 263.1 bits: 58.6 E(85289): 6.5e-08 Smith-Waterman score: 239; 46.5% identity (65.1% similar) in 86 aa overlap (154-235:200-284) 130 140 150 160 170 pF1KE5 FLKRSKLGRYNEEERAQQEAEAAQRLAEEKAQASSIPVGSRCEVRAAGQ----SPRRGTV ..:.:: : : :.. . . . . :.: NP_937 KPSQPAAKEPSATPPISNLTKTASESISNLSEAGSIKKGER-ELKIGDRVLVGGTKAGVV 170 180 190 200 210 220 180 190 200 210 220 230 pF1KE5 MYVGLTDFKPGYWIGVRYDEPLGKNDGSVNGKRYFECQAKYGAFVKPAVVTVGDFPEEDY ..: ::: : : ::. :::::::::.: : :::.:: ::: :. :: :: NP_937 RFLGETDFAKGEWCGVELDEPLGKNDGAVAGTRYFQCQPKYGLFAPVHKVTKIGFPSTTP 230 240 250 260 270 280 240 pF1KE5 GLDEI NP_937 AKAKANAVRRVMATTSASLKRSPSASSLSSMSSVASSVSSRPSRTGLLTETSSRYARKIS 290 300 310 320 330 340 >-- initn: 199 init1: 156 opt: 224 Z-score: 246.6 bits: 55.6 E(85289): 5.4e-07 Smith-Waterman score: 224; 47.1% identity (70.0% similar) in 70 aa overlap (161-230:60-125) 140 150 160 170 180 190 pF1KE5 GRYNEEERAQQEAEAAQRLAEEKAQASSIPVGSRCEVRAAGQSPRRGTVMYVGLTDFKPG :: : : . :..: : ....: :.: :: NP_937 VAPVEKTISSEKASSTPSSETQEEFVDDFRVGER--VWVNGNKP--GFIQFLGETQFAPG 30 40 50 60 70 80 200 210 220 230 240 pF1KE5 YWIGVRYDEPLGKNDGSVNGKRYFECQAKYGAFVKPAVVTVGDFPEEDYGLDEI : :. :::.::::::: : :::.:. : :..:. .: NP_937 QWAGIVLDEPIGKNDGSVAGVRYFQCEPLKGIFTRPSKLTRKVQAEDEANGLQTTPASRA 90 100 110 120 130 140 NP_937 TSPLCTSTASMVSSSPSTPSNIPQKPSQPAAKEPSATPPISNLTKTASESISNLSEAGSI 150 160 170 180 190 200 >>XP_016875279 (OMIM: 179838) PREDICTED: CAP-Gly domain- (1392 aa) initn: 233 init1: 233 opt: 239 Z-score: 263.1 bits: 58.6 E(85289): 6.5e-08 Smith-Waterman score: 239; 46.5% identity (65.1% similar) in 86 aa overlap (154-235:200-284) 130 140 150 160 170 pF1KE5 FLKRSKLGRYNEEERAQQEAEAAQRLAEEKAQASSIPVGSRCEVRAAGQ----SPRRGTV ..:.:: : : :.. . . . . :.: XP_016 KPSQPAAKEPSATPPISNLTKTASESISNLSEAGSIKKGER-ELKIGDRVLVGGTKAGVV 170 180 190 200 210 220 180 190 200 210 220 230 pF1KE5 MYVGLTDFKPGYWIGVRYDEPLGKNDGSVNGKRYFECQAKYGAFVKPAVVTVGDFPEEDY ..: ::: : : ::. :::::::::.: : :::.:: ::: :. :: :: XP_016 RFLGETDFAKGEWCGVELDEPLGKNDGAVAGTRYFQCQPKYGLFAPVHKVTKIGFPSTTP 230 240 250 260 270 280 240 pF1KE5 GLDEI XP_016 AKAKANAVRRVMATTSASLKRSPSASSLSSMSSVASSVSSRPSRTGLLTETSSRYARKIS 290 300 310 320 330 340 >-- initn: 199 init1: 156 opt: 224 Z-score: 246.6 bits: 55.6 E(85289): 5.4e-07 Smith-Waterman score: 224; 47.1% identity (70.0% similar) in 70 aa overlap (161-230:60-125) 140 150 160 170 180 190 pF1KE5 GRYNEEERAQQEAEAAQRLAEEKAQASSIPVGSRCEVRAAGQSPRRGTVMYVGLTDFKPG :: : : . :..: : ....: :.: :: XP_016 VAPVEKTISSEKASSTPSSETQEEFVDDFRVGER--VWVNGNKP--GFIQFLGETQFAPG 30 40 50 60 70 80 200 210 220 230 240 pF1KE5 YWIGVRYDEPLGKNDGSVNGKRYFECQAKYGAFVKPAVVTVGDFPEEDYGLDEI : :. :::.::::::: : :::.:. : :..:. .: XP_016 QWAGIVLDEPIGKNDGSVAGVRYFQCEPLKGIFTRPSKLTRKVQAEDEANGLQTTPASRA 90 100 110 120 130 140 XP_016 TSPLCTSTASMVSSSPSTPSNIPQKPSQPAAKEPSATPPISNLTKTASESISNLSEAGSI 150 160 170 180 190 200 244 residues in 1 query sequences 60827320 residues in 85289 library sequences Tcomplib [36.3.4 Apr, 2011] (8 proc) start: Mon Nov 7 22:31:15 2016 done: Mon Nov 7 22:31:16 2016 Total Scan time: 4.310 Total Display time: -0.020 Function used was FASTA [36.3.4 Apr, 2011]