# /hgtech/tools/solaris8/bin/fasta34_t -T 8 -b50 -d10 -E0.01 -H -Off03515.fasta.nr -Q ff03515.ptfa /cdna2/lib/nr/nr 2 FASTA searches a protein or DNA sequence data bank version 34.26.5 April 26, 2007 Please cite: W.R. Pearson & D.J. Lipman PNAS (1988) 85:2444-2448 ff03515, 192 aa vs /cdna2/lib/nr/nr library 2355815254 residues in 6825066 sequences statistics sampled from 60000 to 6824095 sequences Expectation_n fit: rho(ln(x))= 4.9250+/-0.00018; mu= 8.5743+/- 0.010 mean_var=62.0650+/-12.065, 0's: 44 Z-trim: 52 B-trim: 12 in 1/65 Lambda= 0.162799 FASTA (3.5 Sept 2006) function [optimized, BL50 matrix (15:-5)] ktup: 2 join: 36, opt: 24, open/ext: -10/-2, width: 16 The best scores are: opt bits E(6825066) gi|62088560|dbj|BAD92727.1| H2A histone family, me ( 192) 1234 297.7 5.9e-79 gi|71042448|pdb|1ZR5|A Chain A, Crystal Structure ( 214) 1124 271.9 3.9e-71 gi|32492946|gb|AAN08620.1| medulloblastoma antigen ( 366) 1124 272.1 6e-71 gi|3341992|gb|AAC39908.1| histone macroH2A1.2 [Hom ( 371) 1124 272.1 6e-71 gi|90110023|sp|O75367|H2AY_HUMAN Core histone macr ( 372) 1124 272.1 6e-71 gi|114601789|ref|XP_001168507.1| PREDICTED: simila ( 380) 1124 272.1 6.1e-71 gi|73971330|ref|XP_538635.2| PREDICTED: similar to ( 412) 1124 272.1 6.6e-71 gi|87578157|gb|AAI13288.1| H2A histone family, mem ( 371) 1123 271.9 7.1e-71 gi|81869559|sp|Q9QZQ8|H2AY_MOUSE Core histone macr ( 372) 1116 270.2 2.2e-70 gi|90110024|sp|Q02874|H2AY_RAT Core histone macro- ( 371) 1115 270.0 2.6e-70 gi|149726345|ref|XP_001504412.1| PREDICTED: simila ( 372) 1115 270.0 2.6e-70 gi|63100748|gb|AAH95406.1| H2A histone family, mem ( 372) 1115 270.0 2.6e-70 gi|10435335|dbj|BAB14565.1| unnamed protein produc ( 372) 1115 270.0 2.6e-70 gi|3493529|gb|AAC33433.1| histone macroH2A1.2 [Hom ( 372) 1115 270.0 2.6e-70 gi|126290042|ref|XP_001365244.1| PREDICTED: simila ( 372) 1108 268.3 8.2e-70 gi|126290045|ref|XP_001365313.1| PREDICTED: simila ( 373) 1108 268.3 8.2e-70 gi|126290048|ref|XP_001365386.1| PREDICTED: simila ( 381) 1108 268.3 8.3e-70 gi|114601791|ref|XP_001168478.1| PREDICTED: simila ( 371) 1087 263.4 2.5e-68 gi|82190792|sp|O93327|H2AY_CHICK Core histone macr ( 372) 1086 263.2 2.9e-68 gi|119582622|gb|EAW62218.1| H2A histone family, me ( 173) 972 236.2 1.8e-60 gi|90108649|pdb|1ZR3|A Chain A, Crystal Structure ( 211) 972 236.2 2.1e-60 gi|119582624|gb|EAW62220.1| H2A histone family, me ( 368) 972 236.4 3.3e-60 gi|90076556|dbj|BAE87958.1| unnamed protein produc ( 368) 972 236.4 3.3e-60 gi|3493531|gb|AAC33434.1| histone macroH2A1.1 [Hom ( 369) 972 236.4 3.3e-60 gi|114601787|ref|XP_001168453.1| PREDICTED: simila ( 370) 972 236.4 3.4e-60 gi|86373534|gb|AAA41561.2| histone H2A.1 [Rattus n ( 368) 963 234.3 1.4e-59 gi|78100987|pdb|1YD9|A Chain A, 1.6a Crystal Struc ( 193) 959 233.1 1.7e-59 gi|13905316|gb|AAH06955.1| H2afy protein [Mus musc ( 369) 961 233.8 2e-59 gi|126290039|ref|XP_001365186.1| PREDICTED: simila ( 369) 955 232.4 5.3e-59 gi|3395780|gb|AAC28846.1| histone macroH2A1.1 [Gal ( 369) 930 226.5 3.1e-57 gi|120538711|gb|AAI29543.1| LOC100036858 protein [ ( 369) 916 223.2 3e-56 gi|58403354|gb|AAH89240.1| H2afy protein [Xenopus ( 369) 911 222.1 6.9e-56 gi|92098326|gb|AAI15315.1| Zgc:136891 [Danio rerio ( 357) 861 210.3 2.3e-52 gi|49900191|gb|AAH76893.1| H2A histone [Xenopus tr ( 372) 816 199.7 3.6e-49 gi|33416694|gb|AAH56065.1| H2afy2-prov protein [Xe ( 371) 813 199.0 5.9e-49 gi|66911347|gb|AAH97115.1| Zgc:114051 [Danio rerio ( 367) 811 198.6 8.1e-49 gi|49115274|gb|AAH73272.1| MGC80637 protein [Xenop ( 371) 809 198.1 1.1e-48 gi|74229913|gb|ABA00475.1| H2A histone family memb ( 318) 792 194.1 1.6e-47 gi|73952817|ref|XP_546142.2| PREDICTED: similar to ( 306) 776 190.3 2.1e-46 gi|73952823|ref|XP_861172.1| PREDICTED: similar to ( 363) 776 190.3 2.4e-46 gi|73952825|ref|XP_861204.1| PREDICTED: similar to ( 372) 776 190.4 2.4e-46 gi|73952821|ref|XP_850298.1| PREDICTED: similar to ( 372) 776 190.4 2.4e-46 gi|114630951|ref|XP_001170464.1| PREDICTED: simila ( 306) 775 190.1 2.4e-46 gi|73952819|ref|XP_861114.1| PREDICTED: similar to ( 388) 776 190.4 2.5e-46 gi|114630945|ref|XP_001170543.1| PREDICTED: simila ( 363) 775 190.1 2.8e-46 gi|194042746|ref|XP_001927061.1| PREDICTED: simila ( 372) 775 190.1 2.9e-46 gi|12585260|sp|Q9P0M6.3|H2AW_HUMAN Core histone ma ( 372) 775 190.1 2.9e-46 gi|114630947|ref|XP_001170515.1| PREDICTED: simila ( 372) 775 190.1 2.9e-46 gi|94534822|gb|AAI16137.1| H2A histone family, mem ( 372) 775 190.1 2.9e-46 gi|118599884|gb|AAH22725.1| H2afy2 protein [Mus mu ( 192) 768 188.3 5.3e-46 >>gi|62088560|dbj|BAD92727.1| H2A histone family, member (192 aa) initn: 1234 init1: 1234 opt: 1234 Z-score: 1572.6 bits: 297.7 E(): 5.9e-79 Smith-Waterman score: 1234; 100.000% identity (100.000% similar) in 192 aa overlap (1-192:1-192) 10 20 30 40 50 60 ff0351 DFPALLLFPTSCFSLQLNLIHSEISNLAGFEVEAIINPTNADIDLKDDLGNTLEKKGGKE :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: gi|620 DFPALLLFPTSCFSLQLNLIHSEISNLAGFEVEAIINPTNADIDLKDDLGNTLEKKGGKE 10 20 30 40 50 60 70 80 90 100 110 120 ff0351 FVEAVLELRKKNGPLEVAGAAVSAGHGLPAKFVIHCNSPVWGADKCEELLEKTVKNCLAL :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: gi|620 FVEAVLELRKKNGPLEVAGAAVSAGHGLPAKFVIHCNSPVWGADKCEELLEKTVKNCLAL 70 80 90 100 110 120 130 140 150 160 170 180 ff0351 ADDKKLKSIAFPSIGSGRNGFPKQTAAQLILKAISSYFVSTMSSSIKTVYFVLFDSESIG :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: gi|620 ADDKKLKSIAFPSIGSGRNGFPKQTAAQLILKAISSYFVSTMSSSIKTVYFVLFDSESIG 130 140 150 160 170 180 190 ff0351 IYVQEMAKLDAN :::::::::::: gi|620 IYVQEMAKLDAN 190 >>gi|71042448|pdb|1ZR5|A Chain A, Crystal Structure Of T (214 aa) initn: 1124 init1: 1124 opt: 1124 Z-score: 1432.3 bits: 271.9 E(): 3.9e-71 Smith-Waterman score: 1124; 99.435% identity (100.000% similar) in 177 aa overlap (16-192:38-214) 10 20 30 40 ff0351 DFPALLLFPTSCFSLQLNLIHSEISNLAGFEVEAIINPTNADIDL .::::::::::::::::::::::::::::: gi|710 SKAASADSTTEGTPADGFTVLSTKSLFLGQKLNLIHSEISNLAGFEVEAIINPTNADIDL 10 20 30 40 50 60 50 60 70 80 90 100 ff0351 KDDLGNTLEKKGGKEFVEAVLELRKKNGPLEVAGAAVSAGHGLPAKFVIHCNSPVWGADK :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: gi|710 KDDLGNTLEKKGGKEFVEAVLELRKKNGPLEVAGAAVSAGHGLPAKFVIHCNSPVWGADK 70 80 90 100 110 120 110 120 130 140 150 160 ff0351 CEELLEKTVKNCLALADDKKLKSIAFPSIGSGRNGFPKQTAAQLILKAISSYFVSTMSSS :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: gi|710 CEELLEKTVKNCLALADDKKLKSIAFPSIGSGRNGFPKQTAAQLILKAISSYFVSTMSSS 130 140 150 160 170 180 170 180 190 ff0351 IKTVYFVLFDSESIGIYVQEMAKLDAN ::::::::::::::::::::::::::: gi|710 IKTVYFVLFDSESIGIYVQEMAKLDAN 190 200 210 >>gi|32492946|gb|AAN08620.1| medulloblastoma antigen MU- (366 aa) initn: 1124 init1: 1124 opt: 1124 Z-score: 1429.0 bits: 272.1 E(): 6e-71 Smith-Waterman score: 1124; 99.435% identity (100.000% similar) in 177 aa overlap (16-192:190-366) 10 20 30 40 ff0351 DFPALLLFPTSCFSLQLNLIHSEISNLAGFEVEAIINPTNADIDL .::::::::::::::::::::::::::::: gi|324 SKAASADSTTEGTPADGFTVLSTKSLFLGQKLNLIHSEISNLAGFEVEAIINPTNADIDL 160 170 180 190 200 210 50 60 70 80 90 100 ff0351 KDDLGNTLEKKGGKEFVEAVLELRKKNGPLEVAGAAVSAGHGLPAKFVIHCNSPVWGADK :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: gi|324 KDDLGNTLEKKGGKEFVEAVLELRKKNGPLEVAGAAVSAGHGLPAKFVIHCNSPVWGADK 220 230 240 250 260 270 110 120 130 140 150 160 ff0351 CEELLEKTVKNCLALADDKKLKSIAFPSIGSGRNGFPKQTAAQLILKAISSYFVSTMSSS :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: gi|324 CEELLEKTVKNCLALADDKKLKSIAFPSIGSGRNGFPKQTAAQLILKAISSYFVSTMSSS 280 290 300 310 320 330 170 180 190 ff0351 IKTVYFVLFDSESIGIYVQEMAKLDAN ::::::::::::::::::::::::::: gi|324 IKTVYFVLFDSESIGIYVQEMAKLDAN 340 350 360 >>gi|3341992|gb|AAC39908.1| histone macroH2A1.2 [Homo sa (371 aa) initn: 1124 init1: 1124 opt: 1124 Z-score: 1428.9 bits: 272.1 E(): 6e-71 Smith-Waterman score: 1124; 99.435% identity (100.000% similar) in 177 aa overlap (16-192:195-371) 10 20 30 40 ff0351 DFPALLLFPTSCFSLQLNLIHSEISNLAGFEVEAIINPTNADIDL .::::::::::::::::::::::::::::: gi|334 SKAASADSTTEGTPADGFTVLSTKSLFLGQKLNLIHSEISNLAGFEVEAIINPTNADIDL 170 180 190 200 210 220 50 60 70 80 90 100 ff0351 KDDLGNTLEKKGGKEFVEAVLELRKKNGPLEVAGAAVSAGHGLPAKFVIHCNSPVWGADK :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: gi|334 KDDLGNTLEKKGGKEFVEAVLELRKKNGPLEVAGAAVSAGHGLPAKFVIHCNSPVWGADK 230 240 250 260 270 280 110 120 130 140 150 160 ff0351 CEELLEKTVKNCLALADDKKLKSIAFPSIGSGRNGFPKQTAAQLILKAISSYFVSTMSSS :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: gi|334 CEELLEKTVKNCLALADDKKLKSIAFPSIGSGRNGFPKQTAAQLILKAISSYFVSTMSSS 290 300 310 320 330 340 170 180 190 ff0351 IKTVYFVLFDSESIGIYVQEMAKLDAN ::::::::::::::::::::::::::: gi|334 IKTVYFVLFDSESIGIYVQEMAKLDAN 350 360 370 >>gi|90110023|sp|O75367|H2AY_HUMAN Core histone macro-H2 (372 aa) initn: 1124 init1: 1124 opt: 1124 Z-score: 1428.8 bits: 272.1 E(): 6e-71 Smith-Waterman score: 1124; 99.435% identity (100.000% similar) in 177 aa overlap (16-192:196-372) 10 20 30 40 ff0351 DFPALLLFPTSCFSLQLNLIHSEISNLAGFEVEAIINPTNADIDL .::::::::::::::::::::::::::::: gi|901 SKAASADSTTEGTPADGFTVLSTKSLFLGQKLNLIHSEISNLAGFEVEAIINPTNADIDL 170 180 190 200 210 220 50 60 70 80 90 100 ff0351 KDDLGNTLEKKGGKEFVEAVLELRKKNGPLEVAGAAVSAGHGLPAKFVIHCNSPVWGADK :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: gi|901 KDDLGNTLEKKGGKEFVEAVLELRKKNGPLEVAGAAVSAGHGLPAKFVIHCNSPVWGADK 230 240 250 260 270 280 110 120 130 140 150 160 ff0351 CEELLEKTVKNCLALADDKKLKSIAFPSIGSGRNGFPKQTAAQLILKAISSYFVSTMSSS :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: gi|901 CEELLEKTVKNCLALADDKKLKSIAFPSIGSGRNGFPKQTAAQLILKAISSYFVSTMSSS 290 300 310 320 330 340 170 180 190 ff0351 IKTVYFVLFDSESIGIYVQEMAKLDAN ::::::::::::::::::::::::::: gi|901 IKTVYFVLFDSESIGIYVQEMAKLDAN 350 360 370 >>gi|114601789|ref|XP_001168507.1| PREDICTED: similar to (380 aa) initn: 1124 init1: 1124 opt: 1124 Z-score: 1428.7 bits: 272.1 E(): 6.1e-71 Smith-Waterman score: 1124; 99.435% identity (100.000% similar) in 177 aa overlap (16-192:204-380) 10 20 30 40 ff0351 DFPALLLFPTSCFSLQLNLIHSEISNLAGFEVEAIINPTNADIDL .::::::::::::::::::::::::::::: gi|114 SKAASADSTTEGTPADGFTVLSTKSLFLGQKLNLIHSEISNLAGFEVEAIINPTNADIDL 180 190 200 210 220 230 50 60 70 80 90 100 ff0351 KDDLGNTLEKKGGKEFVEAVLELRKKNGPLEVAGAAVSAGHGLPAKFVIHCNSPVWGADK :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: gi|114 KDDLGNTLEKKGGKEFVEAVLELRKKNGPLEVAGAAVSAGHGLPAKFVIHCNSPVWGADK 240 250 260 270 280 290 110 120 130 140 150 160 ff0351 CEELLEKTVKNCLALADDKKLKSIAFPSIGSGRNGFPKQTAAQLILKAISSYFVSTMSSS :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: gi|114 CEELLEKTVKNCLALADDKKLKSIAFPSIGSGRNGFPKQTAAQLILKAISSYFVSTMSSS 300 310 320 330 340 350 170 180 190 ff0351 IKTVYFVLFDSESIGIYVQEMAKLDAN ::::::::::::::::::::::::::: gi|114 IKTVYFVLFDSESIGIYVQEMAKLDAN 360 370 380 >>gi|73971330|ref|XP_538635.2| PREDICTED: similar to H2A (412 aa) initn: 1124 init1: 1124 opt: 1124 Z-score: 1428.2 bits: 272.1 E(): 6.6e-71 Smith-Waterman score: 1124; 99.435% identity (100.000% similar) in 177 aa overlap (16-192:236-412) 10 20 30 40 ff0351 DFPALLLFPTSCFSLQLNLIHSEISNLAGFEVEAIINPTNADIDL .::::::::::::::::::::::::::::: gi|739 SKAASADSTTEGTPADGFTVLSTKSLFLGQKLNLIHSEISNLAGFEVEAIINPTNADIDL 210 220 230 240 250 260 50 60 70 80 90 100 ff0351 KDDLGNTLEKKGGKEFVEAVLELRKKNGPLEVAGAAVSAGHGLPAKFVIHCNSPVWGADK :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: gi|739 KDDLGNTLEKKGGKEFVEAVLELRKKNGPLEVAGAAVSAGHGLPAKFVIHCNSPVWGADK 270 280 290 300 310 320 110 120 130 140 150 160 ff0351 CEELLEKTVKNCLALADDKKLKSIAFPSIGSGRNGFPKQTAAQLILKAISSYFVSTMSSS :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: gi|739 CEELLEKTVKNCLALADDKKLKSIAFPSIGSGRNGFPKQTAAQLILKAISSYFVSTMSSS 330 340 350 360 370 380 170 180 190 ff0351 IKTVYFVLFDSESIGIYVQEMAKLDAN ::::::::::::::::::::::::::: gi|739 IKTVYFVLFDSESIGIYVQEMAKLDAN 390 400 410 >>gi|87578157|gb|AAI13288.1| H2A histone family, member (371 aa) initn: 1123 init1: 1123 opt: 1123 Z-score: 1427.6 bits: 271.9 E(): 7.1e-71 Smith-Waterman score: 1123; 98.870% identity (100.000% similar) in 177 aa overlap (16-192:195-371) 10 20 30 40 ff0351 DFPALLLFPTSCFSLQLNLIHSEISNLAGFEVEAIINPTNADIDL .::::::::::::::::::::::::::::: gi|875 SKAASADSTTEGTPADGFTVLSTKSLFLGQKLNLIHSEISNLAGFEVEAIINPTNADIDL 170 180 190 200 210 220 50 60 70 80 90 100 ff0351 KDDLGNTLEKKGGKEFVEAVLELRKKNGPLEVAGAAVSAGHGLPAKFVIHCNSPVWGADK ::::::::::::::::::::::::::::::::::::::::::::::::::::::.::::: gi|875 KDDLGNTLEKKGGKEFVEAVLELRKKNGPLEVAGAAVSAGHGLPAKFVIHCNSPIWGADK 230 240 250 260 270 280 110 120 130 140 150 160 ff0351 CEELLEKTVKNCLALADDKKLKSIAFPSIGSGRNGFPKQTAAQLILKAISSYFVSTMSSS :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: gi|875 CEELLEKTVKNCLALADDKKLKSIAFPSIGSGRNGFPKQTAAQLILKAISSYFVSTMSSS 290 300 310 320 330 340 170 180 190 ff0351 IKTVYFVLFDSESIGIYVQEMAKLDAN ::::::::::::::::::::::::::: gi|875 IKTVYFVLFDSESIGIYVQEMAKLDAN 350 360 370 >>gi|81869559|sp|Q9QZQ8|H2AY_MOUSE Core histone macro-H2 (372 aa) initn: 1116 init1: 1116 opt: 1116 Z-score: 1418.7 bits: 270.2 E(): 2.2e-70 Smith-Waterman score: 1116; 97.740% identity (100.000% similar) in 177 aa overlap (16-192:196-372) 10 20 30 40 ff0351 DFPALLLFPTSCFSLQLNLIHSEISNLAGFEVEAIINPTNADIDL .::::::::::::::::::::::::::::: gi|818 SKAASADSTTEGTPTDGFTVLSTKSLFLGQKLNLIHSEISNLAGFEVEAIINPTNADIDL 170 180 190 200 210 220 50 60 70 80 90 100 ff0351 KDDLGNTLEKKGGKEFVEAVLELRKKNGPLEVAGAAVSAGHGLPAKFVIHCNSPVWGADK ::::::::::::::::::::::::::::::::::::.::::::::::::::::::::::: gi|818 KDDLGNTLEKKGGKEFVEAVLELRKKNGPLEVAGAAISAGHGLPAKFVIHCNSPVWGADK 230 240 250 260 270 280 110 120 130 140 150 160 ff0351 CEELLEKTVKNCLALADDKKLKSIAFPSIGSGRNGFPKQTAAQLILKAISSYFVSTMSSS ::::::::::::::::::.::::::::::::::::::::::::::::::::::::::::: gi|818 CEELLEKTVKNCLALADDRKLKSIAFPSIGSGRNGFPKQTAAQLILKAISSYFVSTMSSS 290 300 310 320 330 340 170 180 190 ff0351 IKTVYFVLFDSESIGIYVQEMAKLDAN ::::::.:::::::::::::::::::: gi|818 IKTVYFMLFDSESIGIYVQEMAKLDAN 350 360 370 >>gi|90110024|sp|Q02874|H2AY_RAT Core histone macro-H2A. (371 aa) initn: 1115 init1: 1115 opt: 1115 Z-score: 1417.4 bits: 270.0 E(): 2.6e-70 Smith-Waterman score: 1115; 98.305% identity (100.000% similar) in 177 aa overlap (16-192:195-371) 10 20 30 40 ff0351 DFPALLLFPTSCFSLQLNLIHSEISNLAGFEVEAIINPTNADIDL .::::::::::::::::::::::::::::: gi|901 SKAASADSTTEGAPTDGFTVLSTKSLFLGQKLNLIHSEISNLAGFEVEAIINPTNADIDL 170 180 190 200 210 220 50 60 70 80 90 100 ff0351 KDDLGNTLEKKGGKEFVEAVLELRKKNGPLEVAGAAVSAGHGLPAKFVIHCNSPVWGADK :::::.:::::::::::::::::::::::::::::::::::::::::::::::::::::: gi|901 KDDLGSTLEKKGGKEFVEAVLELRKKNGPLEVAGAAVSAGHGLPAKFVIHCNSPVWGADK 230 240 250 260 270 280 110 120 130 140 150 160 ff0351 CEELLEKTVKNCLALADDKKLKSIAFPSIGSGRNGFPKQTAAQLILKAISSYFVSTMSSS ::::::::::::::::::.::::::::::::::::::::::::::::::::::::::::: gi|901 CEELLEKTVKNCLALADDRKLKSIAFPSIGSGRNGFPKQTAAQLILKAISSYFVSTMSSS 290 300 310 320 330 340 170 180 190 ff0351 IKTVYFVLFDSESIGIYVQEMAKLDAN ::::::::::::::::::::::::::: gi|901 IKTVYFVLFDSESIGIYVQEMAKLDAN 350 360 370 192 residues in 1 query sequences 2355815254 residues in 6825066 library sequences Tcomplib [34.26] (8 proc) start: Sat Aug 9 03:06:44 2008 done: Sat Aug 9 03:09:37 2008 Total Scan time: 635.790 Total Display time: 0.030 Function used was FASTA [version 34.26.5 April 26, 2007]