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= Dr_sW_001_B07
(595 letters)
Database: Non-redundant SwissProt sequences
184,735 sequences; 68,354,980 total letters
Score E
Sequences producing significant alignments: (bits) Value
sp|P16237|HMDH_SCHMA 3-hydroxy-3-methylglutaryl-coenzyme A ... 122 9e-28
sp|O76819|HMDH_AGRIP 3-hydroxy-3-methylglutaryl-coenzyme A ... 105 7e-23
sp|P14773|HMDH_DROME 3-hydroxy-3-methylglutaryl-coenzyme A ... 104 2e-22
sp|P34135|HMDH1_DICDI 3-hydroxy-3-methylglutaryl-coenzyme A... 102 6e-22
sp|Q29512|HMDH_RABIT 3-hydroxy-3-methylglutaryl-coenzyme A ... 100 4e-21
sp|P04035|HMDH_HUMAN 3-hydroxy-3-methylglutaryl-coenzyme A ... 100 4e-21
sp|Q01237|HMDH_MOUSE 3-hydroxy-3-methylglutaryl-coenzyme A ... 99 8e-21
sp|P00347|HMDH_CRIGR 3-hydroxy-3-methylglutaryl-coenzyme A ... 99 8e-21
sp|P09610|HMDH_MESAU 3-hydroxy-3-methylglutaryl-coenzyme A ... 99 8e-21
sp|P51639|HMDH_RAT 3-hydroxy-3-methylglutaryl-coenzyme A re... 99 8e-21
>sp|P16237|HMDH_SCHMA 3-hydroxy-3-methylglutaryl-coenzyme A reductase (HMG-CoA reductase)
Length = 948
Score = 122 bits (305), Expect = 9e-28
Identities = 59/110 (53%), Positives = 84/110 (76%)
Frame = +1
Query: 181 TGQDLAQVVESSSCITEMEALPQKNNTNKFDLLITCTMPCLEVGVIGGGTKLKAQKSCLS 360
TGQDLAQVV+SSSC+T++E + + L+ + TMPCLEVG +GGGT+L Q++CL
Sbjct: 774 TGQDLAQVVDSSSCLTQLEV-----DLSDDSLVASVTMPCLEVGTVGGGTRLSGQRACLD 828
Query: 361 LLNLNEEKPCENFARIICATVLCSELNLMAALQTNDLVKSHLRLNRAQNN 510
LL+L+ ++P E+ +RII TVL +EL+LMAAL T+DLVK+H+ NRA+ +
Sbjct: 829 LLDLSVDRPTEHLSRIIAGTVLAAELSLMAALDTDDLVKAHMHFNRAKQS 878
>sp|O76819|HMDH_AGRIP 3-hydroxy-3-methylglutaryl-coenzyme A reductase (HMG-CoA reductase)
Length = 833
Score = 105 bits (263), Expect = 7e-23
Identities = 62/115 (53%), Positives = 75/115 (65%), Gaps = 5/115 (4%)
Frame = +1
Query: 181 TGQDLAQVVESSSCITEMEALPQKNNTNKFDLLITCTMPCLEVGVIGGGTKLKAQKSCLS 360
TGQD AQ V SS+C T ME + N DL +TCTMP LEVG +GGGT L Q +CL
Sbjct: 708 TGQDPAQNVTSSNCSTNMEVCGE----NGEDLYVTCTMPSLEVGTVGGGTILTGQGACLD 763
Query: 361 LLNLNEE--KPCENFARI---ICATVLCSELNLMAALQTNDLVKSHLRLNRAQNN 510
+L + +P EN AR+ ICATVL EL+LMAAL +DLVKSH+R NR+ N
Sbjct: 764 ILGVKGAGARPAENSARLASLICATVLAGELSLMAALVNSDLVKSHMRHNRSTVN 818
>sp|P14773|HMDH_DROME 3-hydroxy-3-methylglutaryl-coenzyme A reductase (HMG-CoA reductase)
Length = 920
Score = 104 bits (259), Expect = 2e-22
Identities = 67/170 (39%), Positives = 91/170 (53%), Gaps = 6/170 (3%)
Frame = +1
Query: 10 TRCQLTHDIVENILKTDAVSLEKLCIAKXXXXXXXXXXXXXXXXXXXXXXXXXXXXXTGQ 189
T C ++ + ++LKTDA +L + C TGQ
Sbjct: 734 TECTISAATLRSVLKTDAKTLVE-CNKLKNMGGSAMAGSIGGNNAHAANMVTAVFLATGQ 792
Query: 190 DLAQVVESSSCITEMEALPQKNNTNKFDLLITCTMPCLEVGVIGGGTKLKAQKSCLSLLN 369
D AQ V SS+C T ME + N DL +TCTMP LEVG +GGGT L Q +CL +L
Sbjct: 793 DPAQNVTSSNCSTAMECWAE----NSEDLYMTCTMPSLEVGTVGGGTGLPGQSACLEMLG 848
Query: 370 L---NEEKPCEN---FARIICATVLCSELNLMAALQTNDLVKSHLRLNRA 501
+ + +P +N A+I+CATV+ EL+LMAAL +DLVKSH+R NR+
Sbjct: 849 VRGAHATRPGDNAKKLAQIVCATVMAGELSLMAALVNSDLVKSHMRHNRS 898
>sp|P34135|HMDH1_DICDI 3-hydroxy-3-methylglutaryl-coenzyme A reductase 1 (HMG-CoA
reductase 1)
Length = 552
Score = 102 bits (255), Expect = 6e-22
Identities = 66/169 (39%), Positives = 90/169 (53%), Gaps = 6/169 (3%)
Frame = +1
Query: 22 LTHDIVENILKTDAVSLEKLCIAKXXXXXXXXXXXXXXXXXXXXXXXXXXXXXTGQDLAQ 201
+T D+V+ +LKT+ +L L IAK TGQD AQ
Sbjct: 390 ITGDVVQRVLKTNVQALVDLNIAKNLIGSAMAGSIGGFNAHASNIVTAIFLA-TGQDCAQ 448
Query: 202 VVESSSCITEMEALPQKNNTNKFDLLITCTMPCLEVGVIGGGTKLKAQKSCLSLLNL--- 372
VESS+CIT+MEA DL IT TMP +EVG +GGGT L AQ +CL ++ +
Sbjct: 449 NVESSNCITQMEACNDGQ-----DLYITVTMPSIEVGTVGGGTSLPAQSACLDIIGVKGS 503
Query: 373 NEEKP---CENFARIICATVLCSELNLMAALQTNDLVKSHLRLNRAQNN 510
+ KP + A+ I + V+ EL+LM+AL L+KSHL+ NRA+ N
Sbjct: 504 SSSKPGANADQLAKTIASAVMAGELSLMSALSAGHLMKSHLQYNRAKTN 552
>sp|Q29512|HMDH_RABIT 3-hydroxy-3-methylglutaryl-coenzyme A reductase (HMG-CoA reductase)
Length = 888
Score = 100 bits (248), Expect = 4e-21
Identities = 57/115 (49%), Positives = 74/115 (64%), Gaps = 6/115 (5%)
Frame = +1
Query: 184 GQDLAQVVESSSCITEMEALPQKNNTNKFDLLITCTMPCLEVGVIGGGTKLKAQKSCLSL 363
GQD AQ V SS+CIT MEA N DL I+CTMP +E+G +GGGT L Q++CL +
Sbjct: 765 GQDAAQNVGSSNCITLMEASGPPNE----DLYISCTMPSIEIGTVGGGTNLLPQQACLQM 820
Query: 364 LNLN---EEKPCEN---FARIICATVLCSELNLMAALQTNDLVKSHLRLNRAQNN 510
L + ++ P EN ARI+C TV+ EL+LMAAL LVKSH+ NR++ N
Sbjct: 821 LGVQGACKDSPGENARQLARIVCGTVMAGELSLMAALAAGHLVKSHMIHNRSKIN 875
>sp|P04035|HMDH_HUMAN 3-hydroxy-3-methylglutaryl-coenzyme A reductase (HMG-CoA reductase)
Length = 888
Score = 100 bits (248), Expect = 4e-21
Identities = 57/115 (49%), Positives = 74/115 (64%), Gaps = 6/115 (5%)
Frame = +1
Query: 184 GQDLAQVVESSSCITEMEALPQKNNTNKFDLLITCTMPCLEVGVIGGGTKLKAQKSCLSL 363
GQD AQ V SS+CIT MEA N DL I+CTMP +E+G +GGGT L Q++CL +
Sbjct: 765 GQDAAQNVGSSNCITLMEASGPTNE----DLYISCTMPSIEIGTVGGGTNLLPQQACLQM 820
Query: 364 LNLN---EEKPCEN---FARIICATVLCSELNLMAALQTNDLVKSHLRLNRAQNN 510
L + ++ P EN ARI+C TV+ EL+LMAAL LVKSH+ NR++ N
Sbjct: 821 LGVQGACKDNPGENARQLARIVCGTVMAGELSLMAALAAGHLVKSHMIHNRSKIN 875
>sp|Q01237|HMDH_MOUSE 3-hydroxy-3-methylglutaryl-coenzyme A reductase (HMG-CoA reductase)
Length = 887
Score = 99.0 bits (245), Expect = 8e-21
Identities = 56/115 (48%), Positives = 74/115 (64%), Gaps = 6/115 (5%)
Frame = +1
Query: 184 GQDLAQVVESSSCITEMEALPQKNNTNKFDLLITCTMPCLEVGVIGGGTKLKAQKSCLSL 363
GQD AQ V SS+CIT MEA N DL I+CTMP +E+G +GGGT L Q++CL +
Sbjct: 764 GQDAAQNVGSSNCITLMEASGPTNE----DLYISCTMPSIEIGTVGGGTNLLPQQACLQM 819
Query: 364 LNLN---EEKPCEN---FARIICATVLCSELNLMAALQTNDLVKSHLRLNRAQNN 510
L + ++ P EN ARI+C TV+ EL+LMAAL LV+SH+ NR++ N
Sbjct: 820 LGVQGACKDNPGENARQLARIVCGTVMAGELSLMAALAAGHLVRSHMVHNRSKIN 874
>sp|P00347|HMDH_CRIGR 3-hydroxy-3-methylglutaryl-coenzyme A reductase (HMG-CoA reductase)
Length = 887
Score = 99.0 bits (245), Expect = 8e-21
Identities = 56/115 (48%), Positives = 74/115 (64%), Gaps = 6/115 (5%)
Frame = +1
Query: 184 GQDLAQVVESSSCITEMEALPQKNNTNKFDLLITCTMPCLEVGVIGGGTKLKAQKSCLSL 363
GQD AQ V SS+CIT MEA N DL I+CTMP +E+G +GGGT L Q++CL +
Sbjct: 764 GQDAAQNVGSSNCITLMEASGPTNE----DLYISCTMPSIEIGTVGGGTNLLPQQACLQM 819
Query: 364 LNLN---EEKPCEN---FARIICATVLCSELNLMAALQTNDLVKSHLRLNRAQNN 510
L + ++ P EN ARI+C TV+ EL+LMAAL LV+SH+ NR++ N
Sbjct: 820 LGVQGACKDNPGENARQLARIVCGTVMAGELSLMAALAAGHLVRSHMVHNRSKIN 874
>sp|P09610|HMDH_MESAU 3-hydroxy-3-methylglutaryl-coenzyme A reductase (HMG-CoA reductase)
Length = 887
Score = 99.0 bits (245), Expect = 8e-21
Identities = 56/115 (48%), Positives = 74/115 (64%), Gaps = 6/115 (5%)
Frame = +1
Query: 184 GQDLAQVVESSSCITEMEALPQKNNTNKFDLLITCTMPCLEVGVIGGGTKLKAQKSCLSL 363
GQD AQ V SS+CIT MEA N DL I+CTMP +E+G +GGGT L Q++CL +
Sbjct: 764 GQDAAQNVGSSNCITLMEASGPTNE----DLYISCTMPSIEIGTVGGGTNLLPQQACLQM 819
Query: 364 LNLN---EEKPCEN---FARIICATVLCSELNLMAALQTNDLVKSHLRLNRAQNN 510
L + ++ P EN ARI+C TV+ EL+LMAAL LV+SH+ NR++ N
Sbjct: 820 LGVQGACKDNPGENARQLARIVCGTVMAGELSLMAALAAGHLVRSHMVHNRSKIN 874
>sp|P51639|HMDH_RAT 3-hydroxy-3-methylglutaryl-coenzyme A reductase (HMG-CoA reductase)
Length = 887
Score = 99.0 bits (245), Expect = 8e-21
Identities = 56/115 (48%), Positives = 74/115 (64%), Gaps = 6/115 (5%)
Frame = +1
Query: 184 GQDLAQVVESSSCITEMEALPQKNNTNKFDLLITCTMPCLEVGVIGGGTKLKAQKSCLSL 363
GQD AQ V SS+CIT MEA N DL I+CTMP +E+G +GGGT L Q++CL +
Sbjct: 764 GQDAAQNVGSSNCITLMEASGPTNE----DLYISCTMPSIEIGTVGGGTNLLPQQACLQM 819
Query: 364 LNLN---EEKPCEN---FARIICATVLCSELNLMAALQTNDLVKSHLRLNRAQNN 510
L + ++ P EN ARI+C TV+ EL+LMAAL LV+SH+ NR++ N
Sbjct: 820 LGVQGACKDNPGENARQLARIVCGTVMAGELSLMAALAAGHLVRSHMVHNRSKIN 874
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: 58,502,663
Number of Sequences: 369166
Number of extensions: 1020916
Number of successful extensions: 1996
Number of sequences better than 10.0: 10
Number of HSP's better than 10.0 without gapping: 1879
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): 1909
length of database: 68,354,980
effective HSP length: 105
effective length of database: 48,957,805
effective search space used: 4504118060
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)