Reference:
Altschul, Stephen F., Thomas L. Madden, Alejandro A. Schäffer,
Jinghui Zhang, Zheng Zhang, Webb Miller, and David J. Lipman (1997),
"Gapped BLAST and PSI-BLAST: a new generation of protein database search
programs", Nucleic Acids Res. 25:3389-3402.
Query= DrC_02165
(769 letters)
Database: Non-redundant SwissProt sequences
184,735 sequences; 68,354,980 total letters
Score E
Sequences producing significant alignments: (bits) Value
sp|O35855|BCAT2_MOUSE Branched-chain-amino-acid aminotransf... 321 2e-87
sp|O35854|BCAT2_RAT Branched-chain-amino-acid aminotransfer... 318 8e-87
sp|P54687|BCAT1_HUMAN Branched-chain-amino-acid aminotransf... 318 1e-86
sp|Q9GKM4|BCAT1_SHEEP Branched-chain-amino-acid aminotransf... 315 7e-86
sp|P24288|BCAT1_MOUSE Branched-chain-amino-acid aminotransf... 314 1e-85
sp|P54690|BCAT1_RAT Branched-chain-amino-acid aminotransfer... 313 2e-85
sp|O15382|BCAT2_HUMAN Branched-chain-amino-acid aminotransf... 313 4e-85
sp|P54688|BCAT_CAEEL Branched-chain-amino-acid aminotransfe... 288 9e-78
sp|P47176|BCA2_YEAST Branched-chain-amino-acid aminotransfe... 278 2e-74
sp|P38891|BCA1_YEAST Branched-chain-amino-acid aminotransfe... 267 3e-71
>sp|O35855|BCAT2_MOUSE Branched-chain-amino-acid aminotransferase, mitochondrial precursor
(BCAT(m))
Length = 393
Score = 321 bits (822), Expect = 2e-87
Identities = 148/257 (57%), Positives = 186/257 (72%), Gaps = 1/257 (0%)
Frame = +2
Query: 2 KQKPTDSDKISFGSVFTDHMLRIDWNIKTGWSKPLIKPVENFSLHPGCIVFHYATEIFEG 181
++KP S + FG FTDHML ++WN K GW P I+P +N +LHP C HY+ ++FEG
Sbjct: 46 QKKPAPSQALLFGKTFTDHMLMVEWNNKAGWGPPRIQPFQNLTLHPACSGLHYSLQLFEG 105
Query: 182 LKAFYGVDDTVRLFRADKNIDRLLNSAKRAMLPEFDKTELLSCIKELIHIDKEWVPKKRG 361
LKA+ G D VRLFR N+DR+L SA+R LP+FDK ELL CI++LI +DK+WVP G
Sbjct: 106 LKAYKGGDQQVRLFRPWLNMDRMLRSARRLCLPDFDKQELLECIRQLIEVDKDWVPDGNG 165
Query: 362 HSLYVRPTMIAVDPKLGVSKSDSAILYVLLCPVGAYYP-TGFKPVSLLADTKYVRAWPGG 538
SLYVRP +I +P LGV A+LYV+LCPVG+Y+P PVSLLAD +VRAW GG
Sbjct: 166 TSLYVRPVLIGNEPSLGVGMVTQALLYVILCPVGSYFPGDSMTPVSLLADPSFVRAWIGG 225
Query: 539 CGDSKMGSNYAPTIYVQNEAIKQGCQQVLWLFGKDHRLSEVGTMNIFIYMVSKSGKREVV 718
GD K+G NY PT+ VQ EA K+GC+QVLWL+G DH+L+EVGTMNIF+Y + G E+V
Sbjct: 226 VGDCKLGGNYGPTVAVQREAQKRGCEQVLWLYGPDHQLTEVGTMNIFVYWTHEDGVLELV 285
Query: 719 TPPLNGTILPGITRMSL 769
TPPLNG ILPG+ R SL
Sbjct: 286 TPPLNGVILPGVVRQSL 302
>sp|O35854|BCAT2_RAT Branched-chain-amino-acid aminotransferase, mitochondrial precursor
(BCAT(m))
Length = 393
Score = 318 bits (816), Expect = 8e-87
Identities = 146/257 (56%), Positives = 187/257 (72%), Gaps = 1/257 (0%)
Frame = +2
Query: 2 KQKPTDSDKISFGSVFTDHMLRIDWNIKTGWSKPLIKPVENFSLHPGCIVFHYATEIFEG 181
++KP S + FG FTDHML ++WN KTGW P I+P +N +LHP C HY+ ++FEG
Sbjct: 46 QKKPAPSQPLLFGKTFTDHMLMVEWNSKTGWGPPRIQPFQNLTLHPACSGLHYSLQLFEG 105
Query: 182 LKAFYGVDDTVRLFRADKNIDRLLNSAKRAMLPEFDKTELLSCIKELIHIDKEWVPKKRG 361
LKA+ G D VRLFR N+DR+L SA+R LP+FDK ELL CI++LI +DK+WVP G
Sbjct: 106 LKAYKGRDKQVRLFRPWLNMDRMLRSARRLCLPDFDKQELLECIRQLIEVDKDWVPDGNG 165
Query: 362 HSLYVRPTMIAVDPKLGVSKSDSAILYVLLCPVGAYYP-TGFKPVSLLADTKYVRAWPGG 538
SLYVRP +I +P LGV A+L+V+LCPVG+Y+P PVSLLAD +VRAW GG
Sbjct: 166 TSLYVRPVLIGNEPSLGVGMVTQALLFVILCPVGSYFPGDSMTPVSLLADPSFVRAWIGG 225
Query: 539 CGDSKMGSNYAPTIYVQNEAIKQGCQQVLWLFGKDHRLSEVGTMNIFIYMVSKSGKREVV 718
GD K+G NY PT+ VQ EA K+GC+QVLWL+G DH+L+EVGTMNIF+Y + G+ E+
Sbjct: 226 VGDCKLGGNYGPTVAVQQEAQKKGCEQVLWLYGPDHQLTEVGTMNIFVYWTHEDGELELA 285
Query: 719 TPPLNGTILPGITRMSL 769
TPPL+G ILPG+ R SL
Sbjct: 286 TPPLDGIILPGVVRQSL 302
>sp|P54687|BCAT1_HUMAN Branched-chain-amino-acid aminotransferase, cytosolic (BCAT(c))
(ECA39 protein)
Length = 386
Score = 318 bits (815), Expect = 1e-86
Identities = 145/254 (57%), Positives = 192/254 (75%), Gaps = 1/254 (0%)
Frame = +2
Query: 2 KQKPTDSDKISFGSVFTDHMLRIDWNIKTGWSKPLIKPVENFSLHPGCIVFHYATEIFEG 181
K+KP D + + FG+VFTDHML ++W+ + GW KP IKP++N SLHPG HYA E+FEG
Sbjct: 39 KEKP-DPNNLVFGTVFTDHMLTVEWSSEFGWEKPHIKPLQNLSLHPGSSALHYAVELFEG 97
Query: 182 LKAFYGVDDTVRLFRADKNIDRLLNSAKRAMLPEFDKTELLSCIKELIHIDKEWVPKKRG 361
LKAF GVD+ +RLF+ + N+DR+ SA RA LP FDK ELL CI++L+ +D+EWVP
Sbjct: 98 LKAFRGVDNKIRLFQPNLNMDRMYRSAVRATLPVFDKEELLECIQQLVKLDQEWVPYSTS 157
Query: 362 HSLYVRPTMIAVDPKLGVSKSDSAILYVLLCPVGAYYPTG-FKPVSLLADTKYVRAWPGG 538
SLY+RPT I +P LGV K A+L+VLL PVG Y+ +G F PVSL A+ KYVRAW GG
Sbjct: 158 ASLYIRPTFIGTEPSLGVKKPTKALLFVLLSPVGPYFSSGTFNPVSLWANPKYVRAWKGG 217
Query: 539 CGDSKMGSNYAPTIYVQNEAIKQGCQQVLWLFGKDHRLSEVGTMNIFIYMVSKSGKREVV 718
GD KMG NY +++ Q EA+ GCQQVLWL+G+DH+++EVGTMN+F+Y +++ G+ E+
Sbjct: 218 TGDCKMGGNYGSSLFAQCEAVDNGCQQVLWLYGEDHQITEVGTMNLFLYWINEDGEEELA 277
Query: 719 TPPLNGTILPGITR 760
TPPL+G ILPG+TR
Sbjct: 278 TPPLDGIILPGVTR 291
>sp|Q9GKM4|BCAT1_SHEEP Branched-chain-amino-acid aminotransferase, cytosolic (BCAT(c))
Length = 385
Score = 315 bits (808), Expect = 7e-86
Identities = 142/257 (55%), Positives = 193/257 (75%), Gaps = 1/257 (0%)
Frame = +2
Query: 2 KQKPTDSDKISFGSVFTDHMLRIDWNIKTGWSKPLIKPVENFSLHPGCIVFHYATEIFEG 181
K+KP D + FG+VFTDHML ++W+++ GW KP IKP++N SLHPG HYA E+FEG
Sbjct: 38 KEKP-DPSTLVFGTVFTDHMLTVEWSLELGWEKPRIKPLQNLSLHPGSSALHYAVELFEG 96
Query: 182 LKAFYGVDDTVRLFRADKNIDRLLNSAKRAMLPEFDKTELLSCIKELIHIDKEWVPKKRG 361
LKAF GVD+ +RLFR + N+DR+ SA RA LP FDK ELL CI++L+ +D+EWVP
Sbjct: 97 LKAFRGVDNKIRLFRPNLNMDRMYRSAMRATLPAFDKKELLECIQQLVKLDEEWVPYSTS 156
Query: 362 HSLYVRPTMIAVDPKLGVSKSDSAILYVLLCPVGAYYPTG-FKPVSLLADTKYVRAWPGG 538
SLY+RPT I +P LGV K A+L+V+L PVG Y+ +G F PVSL A+ KYVRAW GG
Sbjct: 157 ASLYIRPTFIGTEPSLGVKKPTKALLFVILSPVGPYFSSGSFNPVSLWANPKYVRAWKGG 216
Query: 539 CGDSKMGSNYAPTIYVQNEAIKQGCQQVLWLFGKDHRLSEVGTMNIFIYMVSKSGKREVV 718
GD KMG NY +++ Q EA++ CQQVLWL+G++++++EVGTMN+F+Y +++ G+ E+
Sbjct: 217 TGDCKMGGNYGSSLFAQCEAVENACQQVLWLYGEENQITEVGTMNLFLYWINEDGEEELA 276
Query: 719 TPPLNGTILPGITRMSL 769
TPPL+G ILPG+ R S+
Sbjct: 277 TPPLDGIILPGVMRQSI 293
>sp|P24288|BCAT1_MOUSE Branched-chain-amino-acid aminotransferase, cytosolic (BCAT(c))
(ECA39 protein)
Length = 386
Score = 314 bits (805), Expect = 1e-85
Identities = 144/257 (56%), Positives = 188/257 (73%), Gaps = 1/257 (0%)
Frame = +2
Query: 2 KQKPTDSDKISFGSVFTDHMLRIDWNIKTGWSKPLIKPVENFSLHPGCIVFHYATEIFEG 181
K+KP D D + FG+ FTDHML ++W+ +GW KP IKP N +HP V HYA E+FEG
Sbjct: 39 KEKP-DPDSLVFGATFTDHMLTVEWSSASGWEKPHIKPFGNLPIHPAASVLHYAVELFEG 97
Query: 182 LKAFYGVDDTVRLFRADKNIDRLLNSAKRAMLPEFDKTELLSCIKELIHIDKEWVPKKRG 361
LKAF GVD+ +RLFR D N+DR+ SA R LP FDK ELL CI +L+ ID+EWVP
Sbjct: 98 LKAFRGVDNKIRLFRPDLNMDRMCRSAVRTTLPMFDKEELLKCILQLLQIDQEWVPYSTS 157
Query: 362 HSLYVRPTMIAVDPKLGVSKSDSAILYVLLCPVGAYYPTG-FKPVSLLADTKYVRAWPGG 538
SLY+RPT I +P LGV K A+L+V+L PVG Y+ +G F PVSL A+ KY+RAW GG
Sbjct: 158 ASLYIRPTFIGTEPSLGVKKPSKALLFVILSPVGPYFSSGSFTPVSLWANPKYIRAWKGG 217
Query: 539 CGDSKMGSNYAPTIYVQNEAIKQGCQQVLWLFGKDHRLSEVGTMNIFIYMVSKSGKREVV 718
GD KMG NY ++ Q EA++ GCQQVLWL+GKD++++EVGTMN+F+Y +++ G+ E+
Sbjct: 218 TGDCKMGGNYGASLLAQCEAVENGCQQVLWLYGKDNQITEVGTMNLFLYWINEDGEEELA 277
Query: 719 TPPLNGTILPGITRMSL 769
TPPL+G ILPG+TR S+
Sbjct: 278 TPPLDGIILPGVTRQSI 294
>sp|P54690|BCAT1_RAT Branched-chain-amino-acid aminotransferase, cytosolic (BCAT(c))
Length = 411
Score = 313 bits (803), Expect = 2e-85
Identities = 144/257 (56%), Positives = 187/257 (72%), Gaps = 1/257 (0%)
Frame = +2
Query: 2 KQKPTDSDKISFGSVFTDHMLRIDWNIKTGWSKPLIKPVENFSLHPGCIVFHYATEIFEG 181
K+KP D + FG+ FTDHML ++W K GW KP IKP EN S+HP V HYA E+FEG
Sbjct: 64 KKKP-DPSSLVFGASFTDHMLMVEWTSKYGWDKPHIKPFENLSIHPAASVLHYAVELFEG 122
Query: 182 LKAFYGVDDTVRLFRADKNIDRLLNSAKRAMLPEFDKTELLSCIKELIHIDKEWVPKKRG 361
LKAF GVD+ +RLFR D N+ R+ SA R LPEFDK ELL C+ +LI +D+EWVP
Sbjct: 123 LKAFRGVDNKIRLFRPDLNMKRMCRSAVRTTLPEFDKEELLQCVLQLIQLDREWVPYSTS 182
Query: 362 HSLYVRPTMIAVDPKLGVSKSDSAILYVLLCPVGAYYPTG-FKPVSLLADTKYVRAWPGG 538
SLY+RPT I ++P LGV K A+L+V+L PVG+Y+ G F PVSL A+ K+VR+W GG
Sbjct: 183 ASLYIRPTFIGIEPSLGVKKPSKALLFVILSPVGSYFSNGTFSPVSLWANPKFVRSWKGG 242
Query: 539 CGDSKMGSNYAPTIYVQNEAIKQGCQQVLWLFGKDHRLSEVGTMNIFIYMVSKSGKREVV 718
GD KMG NY ++ Q EA + GC QVLWL+GK++R++EVGTMN+F+Y ++K G+ E+
Sbjct: 243 TGDFKMGCNYGSSLLAQCEAAENGCHQVLWLYGKENRITEVGTMNLFLYWINKDGEEELA 302
Query: 719 TPPLNGTILPGITRMSL 769
TPPL+G ILPG+TR S+
Sbjct: 303 TPPLDGVILPGVTRQSI 319
>sp|O15382|BCAT2_HUMAN Branched-chain-amino-acid aminotransferase, mitochondrial precursor
(BCAT(m)) (Placental protein 18) (PP18)
Length = 392
Score = 313 bits (801), Expect = 4e-85
Identities = 146/256 (57%), Positives = 186/256 (72%), Gaps = 1/256 (0%)
Frame = +2
Query: 5 QKPTDSDKISFGSVFTDHMLRIDWNIKTGWSKPLIKPVENFSLHPGCIVFHYATEIFEGL 184
+KP + + FG FTDHML ++WN K GW +P I+P +N +LHP HY+ ++FEG+
Sbjct: 47 KKPGPGEPLVFGKTFTDHMLMVEWNDK-GWGQPRIQPFQNLTLHPASSSLHYSLQLFEGM 105
Query: 185 KAFYGVDDTVRLFRADKNIDRLLNSAKRAMLPEFDKTELLSCIKELIHIDKEWVPKKRGH 364
KAF G D VRLFR N+DR+L SA R LP FDK ELL CI+ LI +DK+WVP G
Sbjct: 106 KAFKGKDQQVRLFRPWLNMDRMLRSAMRLCLPSFDKLELLECIRRLIEVDKDWVPDAAGT 165
Query: 365 SLYVRPTMIAVDPKLGVSKSDSAILYVLLCPVGAYYPTG-FKPVSLLADTKYVRAWPGGC 541
SLYVRP +I +P LGVS+ A+L+V+LCPVGAY+P G PVSLLAD ++RAW GG
Sbjct: 166 SLYVRPVLIGNEPSLGVSQPTRALLFVILCPVGAYFPGGSVTPVSLLADPAFIRAWVGGV 225
Query: 542 GDSKMGSNYAPTIYVQNEAIKQGCQQVLWLFGKDHRLSEVGTMNIFIYMVSKSGKREVVT 721
G+ K+G NY PT+ VQ EA+K+GC+QVLWL+G DH+L+EVGTMNIF+Y + G E+VT
Sbjct: 226 GNYKLGGNYGPTVLVQQEALKRGCEQVLWLYGPDHQLTEVGTMNIFVYWTHEDGVLELVT 285
Query: 722 PPLNGTILPGITRMSL 769
PPLNG ILPG+ R SL
Sbjct: 286 PPLNGVILPGVVRQSL 301
>sp|P54688|BCAT_CAEEL Branched-chain-amino-acid aminotransferase, cytosolic (BCAT) (ECA39
protein)
Length = 415
Score = 288 bits (738), Expect = 9e-78
Identities = 131/257 (50%), Positives = 186/257 (72%), Gaps = 1/257 (0%)
Frame = +2
Query: 2 KQKPTDSDKISFGSVFTDHMLRIDWNIKTGWSKPLIKPVENFSLHPGCIVFHYATEIFEG 181
K KP D K+ FG + D+M+ DW+ + GW P I+P+ +HPG V HYA+E+FEG
Sbjct: 61 KTKPLDPTKLKFGHTYADYMMTCDWDAERGWHHPKIEPIGELKIHPGAKVLHYASELFEG 120
Query: 182 LKAFYGVDDTVRLFRADKNIDRLLNSAKRAMLPEFDKTELLSCIKELIHIDKEWVPKKRG 361
+KA+ G+D+ +R+FR + N+ R+ +A RA LP+FD E+++ + EL+ +D+EWVP
Sbjct: 121 MKAYRGIDNKIRMFRPEMNMARMKRTALRAALPDFDSEEMINVLTELLRLDQEWVPNSDV 180
Query: 362 HSLYVRPTMIAVDPKLGVSKSDSAILYVLLCPVGAYYPTGFKPVSLLADTKYVRAWPGGC 541
SLY+RPT+I DP LGV + A ++V+ PVGAYY TGF+PVSLLAD++++RA+PGG
Sbjct: 181 CSLYLRPTLIGTDPTLGVGCATEAKMFVITGPVGAYYSTGFQPVSLLADSRFIRAFPGGV 240
Query: 542 GDSKMGSNYAPTIYVQNEAIKQGCQQVLWLFGKDHRLSEVGTMNIFIYMVSKSGKREVVT 721
G KMG NYAPTI+V EA + CQQVLWL+G++ L+EVGTMNIF++ ++ G E++T
Sbjct: 241 GAYKMGCNYAPTIWVGKEAASKNCQQVLWLYGENEDLTEVGTMNIFLFWKNEEGDMELIT 300
Query: 722 PPLN-GTILPGITRMSL 769
PPL+ G ILPG+TR SL
Sbjct: 301 PPLHRGLILPGVTRDSL 317
>sp|P47176|BCA2_YEAST Branched-chain-amino-acid aminotransferase, cytosolic (BCAT) (TWT2
protein)
Length = 376
Score = 278 bits (710), Expect = 2e-74
Identities = 134/255 (52%), Positives = 177/255 (69%), Gaps = 1/255 (0%)
Frame = +2
Query: 8 KPTDSDKISFGSVFTDHMLRIDWNIKTGWSKPLIKPVENFSLHPGCIVFHYATEIFEGLK 187
KP + ++ FG FTDHML +W + GW P IKP +N SL P +VFHYA E+FEG+K
Sbjct: 21 KPKPNSELVFGKSFTDHMLTAEWTAEKGWGTPEIKPYQNLSLDPSAVVFHYAFELFEGMK 80
Query: 188 AFYGVDDTVRLFRADKNIDRLLNSAKRAMLPEFDKTELLSCIKELIHIDKEWVPKKRGHS 367
A+ VD+ + +FR D N+ R+ SA+R LP FD EL++ I +LI DK VP+ +G+S
Sbjct: 81 AYRTVDNKITMFRPDMNMKRMNKSAQRICLPTFDPEELITLIGKLIQQDKCLVPEGKGYS 140
Query: 368 LYVRPTMIAVDPKLGVSKSDSAILYVLLCPVGAYYPTGFKPVSLLADTKYVRAWPGGCGD 547
LY+RPT+I LGVS D A+LYV+ CPVG YY TGFK V L A RAWPGGCGD
Sbjct: 141 LYIRPTLIGTTAGLGVSTPDRALLYVICCPVGPYYKTGFKAVRLEATDYATRAWPGGCGD 200
Query: 548 SKMGSNYAPTIYVQNEAIKQGCQQVLWLFGKDHRLSEVGTMN-IFIYMVSKSGKREVVTP 724
K+G+NYAP + Q +A +G QQ LWLFG ++ ++EVGTMN F++ SK+GK+E+VT
Sbjct: 201 KKLGANYAPCVLPQLQAASRGYQQNLWLFGPNNNITEVGTMNAFFVFKDSKTGKKELVTA 260
Query: 725 PLNGTILPGITRMSL 769
PL+GTIL G+TR S+
Sbjct: 261 PLDGTILEGVTRDSI 275
>sp|P38891|BCA1_YEAST Branched-chain-amino-acid aminotransferase, mitochondrial precursor
(BCAT) (TWT1 protein) (ECA39 protein)
Length = 393
Score = 267 bits (682), Expect = 3e-71
Identities = 133/255 (52%), Positives = 169/255 (66%), Gaps = 1/255 (0%)
Frame = +2
Query: 8 KPTDSDKISFGSVFTDHMLRIDWNIKTGWSKPLIKPVENFSLHPGCIVFHYATEIFEGLK 187
KP ++++ FG FTDHML I W+ K GW P IKP N SL P VFHYA E+FEGLK
Sbjct: 38 KPRPNEELVFGQTFTDHMLTIPWSAKEGWGTPHIKPYGNLSLDPSACVFHYAFELFEGLK 97
Query: 188 AFYGVDDTVRLFRADKNIDRLLNSAKRAMLPEFDKTELLSCIKELIHIDKEWVPKKRGHS 367
A+ +T+ +FR DKN+ R+ SA R LP F+ EL+ +LI DK VP+ G+S
Sbjct: 98 AYRTPQNTITMFRPDKNMARMNKSAARICLPTFESEELIKLTGKLIEQDKHLVPQGNGYS 157
Query: 368 LYVRPTMIAVDPKLGVSKSDSAILYVLLCPVGAYYPTGFKPVSLLADTKYVRAWPGGCGD 547
LY+RPTMI LGV A+LYV+ PVG YY TGFK V L A RAWPGG GD
Sbjct: 158 LYIRPTMIGTSKGLGVGTPSEALLYVITSPVGPYYKTGFKAVRLEATDYATRAWPGGVGD 217
Query: 548 SKMGSNYAPTIYVQNEAIKQGCQQVLWLFGKDHRLSEVGTMNIFIYMVSK-SGKREVVTP 724
K+G+NYAP I Q +A K+G QQ LWLFG + ++EVGTMN+F ++K +GK+E+VT
Sbjct: 218 KKLGANYAPCILPQLQAAKRGYQQNLWLFGPEKNITEVGTMNVFFVFLNKVTGKKELVTA 277
Query: 725 PLNGTILPGITRMSL 769
PL+GTIL G+TR S+
Sbjct: 278 PLDGTILEGVTRDSV 292
Database: Non-redundant SwissProt sequences
Posted date: Dec 6, 2005 7:40 AM
Number of letters in database: 68,354,980
Number of sequences in database: 184,735
Database: swissprot.01
Posted date: Dec 6, 2005 8:18 AM
Number of letters in database: 66,202,850
Number of sequences in database: 184,431
Lambda K H
0.318 0.134 0.401
Gapped
Lambda K H
0.267 0.0410 0.140
Matrix: BLOSUM62
Gap Penalties: Existence: 11, Extension: 1
Number of Hits to DB: 90,345,938
Number of Sequences: 369166
Number of extensions: 1900002
Number of successful extensions: 4054
Number of sequences better than 10.0: 10
Number of HSP's better than 10.0 without gapping: 3882
Number of HSP's successfully gapped in prelim test: 0
Number of HSP's that attempted gapping in prelim test: 0
Number of HSP's gapped (non-prelim): 4001
length of database: 68,354,980
effective HSP length: 108
effective length of database: 48,403,600
effective search space used: 7115329200
frameshift window, decay const: 40, 0.1
T: 12
A: 40
X1: 16 ( 7.3 bits)
X2: 38 (14.6 bits)
X3: 64 (24.7 bits)
S1: 41 (21.7 bits)