Citrus Sinensis ID: 034451

---

Local Sequence Feature Prediction

Residue Number Marker
--------10--------20--------30--------40--------50--------60--------70--------80--------90----
Protein Sequence
MELPLKSIALTIVIVTVLTWAWRVLNWVWLRPKKLEKFLRQQGLKGNSYRLLFGDLKENSIELKEAKARPLSLDDNIAIRVNPFLHKLVNDYGM
Secondary Structure (PSIPRED)
ccccHHHHHHHHHHHHHHHHHHHHHHHHHHcHHHHHHHHHHccccccccccccccHHHHHHHHHHHccccccccccHHHHHccHHHHHHHHccc
Secondary Structure Prediction (SSPRO)
ccHHHHHHHHHHHHHHHHHHHHHHHHHHHHcHHHHHHHHHHccccccccEEccccHHHHHHHHHHHHcccccccccHHHHHHHHHHHHHHHHcc
Coil and Loop (DISEMBL)
MELPLKSIALTIVIVTVLTWAWRVLNWVWLRPKKLEKFLRQqglkgnsyrLLFGDLKENSIELKeakarplslddniairVNPFLHKLVNDYGM
Flexible Loop (DISEMBL)
melplksiALTIVIVTVLTWAWRVLNWVWLRPKKLEKFLrqqglkgnsyrLLFGDLKENSIELKeakarplslddniairvNPFLHKLVNDYGM
Low Complexity Region (SEG)
MELPLKSIALTIVIVTVLTWAWRVLNWVWLRPKKLEKFLRQQGLKGNSYRLLFGDLKENSIELKEAKARPLSLDDNIAIRVNPFLHKLVNDYGM
Disordered region (IsUnstruct)
*****KSIALTIVIVTVLTWAWRVLNWVWLRPKKLEKFLRQQGLKGNSYRLLFGDLKENSIELKEAKARPLSLDDNIAIRVNPFLHKLVN****
Disordered Region (DISOPRED)
**LPLKSIALTIVIVTVLTWAWRVLNWVWLRPKKLEKFLRQQGLKGNSYRLLFGDLK*****************DNIAIRVNPFLHKLVNDYGM
Disordered Region (DISEMBL)
MELPLKSIALTIVIVTVLTWAWRVLNWVWLRPKKLEKFLRQQGLKGNSYRLLFGDLKENSIELKEAKARPLSLDDNIAIRVNPFLHKLVNDYGM
Disordered Region (DISPRO)
*ELPLKSIALTIVIVTVLTWAWRVLNWVWLRPKKLEKFLRQQGLKGNSYRLLFGDLKENSIELKEAKARPLSLDDNIAIRVNPFLHKLVNDYG*
Transmembrane Helix (TMHMM)
iiiiiiHHHHHHHHHHHHHHHHHHHHHHHooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo
Transmembrane Helix (HMMTOP)
oooooooHHHHHHHHHHHHHHHHHHHHHHHiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiii
Transmembrane Helix (MEMSAT)
ooooooooooHHHHHHHHHHHHHHHHHHHHHHHHHiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiii
TM Helix, Signal Peptide (MEMSAT_SVM)
oooooooooooooooHHHHHHHHHHHHHHHHiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiii
TM Helix, Signal Peptide (Phobius)
oooooHHHHHHHHHHHHHHHHHHHHHHHHHiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiii
Signal Peptide (SignalP HMM Mode)
SSSSSSSSSSSSSSSSSSSSSSoooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo
Signal Peptide (SignalP NN Mode)
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
Coiled Coils (COILS)
MELPLKSIALTIVIVTVLTWAWRVLNWVWLRPKKLEKFLRQQGLKGNSYRLLFGDLKENSIELKEAKARPLSLDDNIAIRVNPFLHKLVNDYGM
Positional Conservation
no confident homologs detected

---

Close Homologs for Annotation Transfer

Close Homologs in SWISS-PROT Database Detected by BLAST

Original result of BLAST against SWISS-PROT Database

ID	Length	Definition	RBH(Q2H)	RBH(H2Q)	Q cover	H cover	Identity	E-value
Query	94	2.2.26 [Sep-21-2011]
Q05047	524	Secologanin synthase OS=C	N/A	no	0.978	0.175	0.451	4e-19
Q9SHG5	519	Cytochrome P450 72C1 OS=A	no	no	0.946	0.171	0.449	9e-19
O48786	520	Cytochrome P450 734A1 OS=	no	no	0.851	0.153	0.310	3e-06
Q6Z6D6	557	Cytochrome P450 734A2 OS=	no	no	0.776	0.131	0.346	2e-05
Q69XM6	538	Cytochrome P450 734A4 OS=	no	no	0.755	0.131	0.369	3e-05
B9X287	542	Cytochrome P450 734A6 OS=	no	no	0.755	0.130	0.324	0.0003

>sp|Q05047|C72A1_CATRO Secologanin synthase OS=Catharanthus roseus GN=CYP72A1 PE=2 SV=1

 Score = 93.2 bits (230), Expect = 4e-19,   Method: Compositional matrix adjust.
 Identities = 42/93 (45%), Positives = 63/93 (67%), Gaps = 1/93 (1%)

Query: 1  MELPLKSIALTIVIVTVLTWAWRVLNWVWLRPKKLEKFLRQQGLKGNSYRLLFGDLKENS 60
          M+   K+IA TI  + V+ WAWRVL+W W  PK++EK LRQQG +GN YR L GD+KE+ 
Sbjct: 5  MDTIRKAIAATIFAL-VMAWAWRVLDWAWFTPKRIEKRLRQQGFRGNPYRFLVGDVKESG 63

Query: 61 IELKEAKARPLSLDDNIAIRVNPFLHKLVNDYG 93
             +EA ++P+  +++I  R+ P ++  +N YG
Sbjct: 64 KMHQEALSKPMEFNNDIVPRLMPHINHTINTYG 96

Converts loganin into secologanin.
Catharanthus roseus (taxid: 4058)
EC: 1EC: .EC: 3EC: .EC: 3EC: .EC: 9

>sp|Q9SHG5|C72C1_ARATH Cytochrome P450 72C1 OS=Arabidopsis thaliana GN=CYP72C1 PE=2 SV=2

 Score = 92.0 bits (227), Expect = 9e-19,   Method: Compositional matrix adjust.
 Identities = 40/89 (44%), Positives = 61/89 (68%)

Query: 5  LKSIALTIVIVTVLTWAWRVLNWVWLRPKKLEKFLRQQGLKGNSYRLLFGDLKENSIELK 64
          ++ + L   ++ +L W WR +NWVWLRPK+LEK+L++QG  GNSYR+L GD++E++   +
Sbjct: 7  VRKVFLIGFLILILNWVWRAVNWVWLRPKRLEKYLKKQGFSGNSYRILMGDMRESNQMDQ 66

Query: 65 EAKARPLSLDDNIAIRVNPFLHKLVNDYG 93
           A + PL LD +   R+ PFLH  V  +G
Sbjct: 67 VAHSLPLPLDADFLPRMMPFLHHTVLKHG 95

Atypical cytochrome P450 involved in brassinosteroids (BRs) inactivation and regulation of BRs homeostasis. Does not possess carbon 26 hydroxylase activity and may inactivate BRs by hydroxylation of carbons other than C-26. Acts in association with CYP734A1 to inactivate BRs and modulate photomorphogenesis.

Arabidopsis thaliana (taxid: 3702)

EC: 1

EC: .

EC: 1

EC: 4

EC: .

EC: -

EC: .

EC: -

>sp|O48786|C734A_ARATH Cytochrome P450 734A1 OS=Arabidopsis thaliana GN=CYP734A1 PE=2 SV=1

 Score = 50.4 bits (119), Expect = 3e-06,   Method: Compositional matrix adjust.
 Identities = 27/87 (31%), Positives = 46/87 (52%), Gaps = 7/87 (8%)

Query: 7  SIALTIVIVTVLTWAWRVLNWVWLRPKKLEKFLRQQGLKGNSYRLLFGDLKENSIELKEA 66
          S+ L++VIV       + ++ +W RP+K+E+   +QG++G  Y    G++KE    + +A
Sbjct: 17 SVILSLVIV-------KGMSLLWWRPRKIEEHFSKQGIRGPPYHFFIGNVKELVGMMLKA 69

Query: 67 KARPLSLDDNIAIRVNPFLHKLVNDYG 93
           + P+    NI  RV  F H     YG
Sbjct: 70 SSHPMPFSHNILPRVLSFYHHWRKIYG 96

Cytochrome P450 involved in brassinosteroids (BRs) inactivation and regulation of BRs homeostasis. Inactivates the BRs castasterone and brassinolide through carbon 26 hydroxylation. Acts in association with CYP72C1 to inactivate BRs and modulate photomorphogenesis.

Arabidopsis thaliana (taxid: 3702)

EC: 1

EC: .

EC: 1

EC: 4

EC: .

EC: -

EC: .

EC: -

>sp|Q6Z6D6|C7342_ORYSJ Cytochrome P450 734A2 OS=Oryza sativa subsp. japonica GN=CYP734A2 PE=2 SV=1

 Score = 47.8 bits (112), Expect = 2e-05,   Method: Composition-based stats.
 Identities = 26/75 (34%), Positives = 39/75 (52%), Gaps = 2/75 (2%)

Query: 21  AWRVLNWVWLRPKKLEKFLRQQGLKGNSYRLLFGDLKENSIELKEAKARPLS--LDDNIA 78
           A R+ + +W RP++LE     QG++G  YR L G ++E    + EA ++P+S     N  
Sbjct: 34  AARMADALWWRPRRLEAHFAAQGVRGPPYRFLLGSVREMVALMAEASSKPMSPPTSHNAL 93

Query: 79  IRVNPFLHKLVNDYG 93
            RV  F H     YG
Sbjct: 94  PRVLAFYHYWRKIYG 108

Cytochrome P450 involved in brassinosteroids (BRs) inactivation and regulation of BRs homeostasis. Is a multifunctional and multisubstrate enzyme that controls the endogenous bioactive BR content both by direct inactivation of castasterone and by decreasing the levels of BR precursors. Catalyzes the oxidation of carbon 22 hydroxylated BR intermediates to produce C26 oxidized metabolites.

Oryza sativa subsp. japonica (taxid: 39947)

EC: 1

EC: .

EC: 1

EC: 4

EC: .

EC: -

EC: .

EC: -

>sp|Q69XM6|C7344_ORYSJ Cytochrome P450 734A4 OS=Oryza sativa subsp. japonica GN=CYP734A4 PE=2 SV=1

 Score = 47.0 bits (110), Expect = 3e-05,   Method: Composition-based stats.
 Identities = 27/73 (36%), Positives = 38/73 (52%), Gaps = 2/73 (2%)

Query: 23 RVLNWVWLRPKKLEKFLRQQGLKGNSYRLLFGDLKENSIELKEAKARPL--SLDDNIAIR 80
          RV + VW RP++LE     QG++G  YR L G ++E    + EA A+P+  +   N   R
Sbjct: 20 RVADAVWWRPRRLEAHFAGQGVRGPPYRFLVGCVREMVALMAEATAKPMPPAAPHNALPR 79

Query: 81 VNPFLHKLVNDYG 93
          V  F H     YG
Sbjct: 80 VLAFYHYWRKIYG 92

Cytochrome P450 involved in brassinosteroids (BRs) inactivation and regulation of BRs homeostasis. Is a multifunctional and multisubstrate enzyme that controls the endogenous bioactive BR content both by direct inactivation of castasterone and by decreasing the levels of BR precursors. Catalyzes the oxidation of carbon 22 hydroxylated BR intermediates to produce C26 oxidized metabolites.

Oryza sativa subsp. japonica (taxid: 39947)

EC: 1

EC: .

EC: 1

EC: 4

EC: .

EC: -

EC: .

EC: -

>sp|B9X287|C7346_ORYSJ Cytochrome P450 734A6 OS=Oryza sativa subsp. japonica GN=CYP734A6 PE=2 SV=1

 Score = 43.5 bits (101), Expect = 3e-04,   Method: Composition-based stats.
 Identities = 24/74 (32%), Positives = 38/74 (51%), Gaps = 3/74 (4%)

Query: 23 RVLNWVWLRPKKLEKFLRQQGLKGNSYRLLFGDLKENSIELKEAKARPLS---LDDNIAI 79
          +VL  +W RP+++E+   +QG+ G  YR L G ++E    +  A A+P+       N+  
Sbjct: 20 KVLEVLWWRPRRVEEHFARQGITGPRYRFLVGCVREMVALMVAASAKPMPPPYRSHNVLP 79

Query: 80 RVNPFLHKLVNDYG 93
          RV  F H     YG
Sbjct: 80 RVLAFYHHWKKIYG 93

Cytochrome P450 involved in brassinosteroids (BRs) inactivation and regulation of BRs homeostasis. Is a multifunctional and multisubstrate enzyme that controls the endogenous bioactive BR content both by direct inactivation of castasterone and by decreasing the levels of BR precursors. Catalyzes the oxidation of carbon 22 hydroxylated BR intermediates to produce C26 oxidized metabolites.

Oryza sativa subsp. japonica (taxid: 39947)

EC: 1

EC: .

EC: 1

EC: 4

EC: .

EC: -

EC: .

EC: -

Close Homologs in the Non-Redundant Database Detected by BLAST

Original result of BLAST against Nonredundant Database

GI	Length	Definition	Q cover	H cover	Identity	E-value
Query	94
224105979	518	cytochrome P450 [Populus trichocarpa] gi	0.989	0.179	0.623	2e-26
224113635	104	cytochrome P450 [Populus trichocarpa] gi	1.0	0.903	0.606	6e-26
224119126	518	cytochrome P450 [Populus trichocarpa] gi	0.989	0.179	0.612	7e-26
224119134	518	cytochrome P450 [Populus trichocarpa] gi	0.989	0.179	0.612	7e-26
224119146	96	predicted protein [Populus trichocarpa]	0.989	0.968	0.602	1e-25
359494197	526	PREDICTED: secologanin synthase-like [Vi	0.989	0.176	0.559	7e-24
359494591	391	PREDICTED: secologanin synthase-like [Vi	0.989	0.237	0.559	8e-24
147795635	518	hypothetical protein VITISV_015554 [Viti	0.914	0.166	0.569	1e-23
147858656	193	hypothetical protein VITISV_020336 [Viti	0.968	0.471	0.560	3e-23
359494208	515	PREDICTED: secologanin synthase-like [Vi	0.968	0.176	0.560	4e-23

>gi|224105979|ref|XP_002333744.1| cytochrome P450 [Populus trichocarpa] gi|222838385|gb|EEE76750.1| cytochrome P450 [Populus trichocarpa]

 Score =  122 bits (307), Expect = 2e-26,   Method: Compositional matrix adjust.
 Identities = 58/93 (62%), Positives = 69/93 (74%)

Query: 1  MELPLKSIALTIVIVTVLTWAWRVLNWVWLRPKKLEKFLRQQGLKGNSYRLLFGDLKENS 60
          M + + SI ++IV V VL WAWRVLNWVW RPKK+E+ LRQQG  G  YRLLFGD KENS
Sbjct: 1  MTVTVTSILVSIVYVAVLRWAWRVLNWVWFRPKKVERCLRQQGFAGKPYRLLFGDWKENS 60

Query: 61 IELKEAKARPLSLDDNIAIRVNPFLHKLVNDYG 93
            LKEA+ +P+ L D +  RV PFLH+LV DYG
Sbjct: 61 DMLKEARTKPIGLSDALLPRVMPFLHQLVKDYG 93

Source: Populus trichocarpa

Species: Populus trichocarpa

Genus: Populus

Family: Salicaceae

Order: Malpighiales

Phylum: Streptophyta

Superkingdom: Eukaryota

>gi|224113635|ref|XP_002332522.1| cytochrome P450 [Populus trichocarpa] gi|222832634|gb|EEE71111.1| cytochrome P450 [Populus trichocarpa]

 Score =  121 bits (303), Expect = 6e-26,   Method: Compositional matrix adjust.
 Identities = 57/94 (60%), Positives = 69/94 (73%)

Query: 1  MELPLKSIALTIVIVTVLTWAWRVLNWVWLRPKKLEKFLRQQGLKGNSYRLLFGDLKENS 60
          M + + SI ++IV V VL WAWRVLNWVW RPKK+E+ LRQQG  G  YRLLFGD KE S
Sbjct: 1  MSITVTSILVSIVYVAVLRWAWRVLNWVWFRPKKVERCLRQQGFAGKPYRLLFGDWKEKS 60

Query: 61 IELKEAKARPLSLDDNIAIRVNPFLHKLVNDYGM 94
            LKEA+ +P+ L D +  RV PFLH+LV DYG+
Sbjct: 61 DMLKEARTKPIGLSDALLPRVMPFLHQLVKDYGI 94

Source: Populus trichocarpa

Species: Populus trichocarpa

Genus: Populus

Family: Salicaceae

Order: Malpighiales

Phylum: Streptophyta

Superkingdom: Eukaryota

>gi|224119126|ref|XP_002331331.1| cytochrome P450 [Populus trichocarpa] gi|222873914|gb|EEF11045.1| cytochrome P450 [Populus trichocarpa]

 Score =  121 bits (303), Expect = 7e-26,   Method: Compositional matrix adjust.
 Identities = 57/93 (61%), Positives = 69/93 (74%)

Query: 1  MELPLKSIALTIVIVTVLTWAWRVLNWVWLRPKKLEKFLRQQGLKGNSYRLLFGDLKENS 60
          M + + SI ++IV V VL WAWRVLNWVW RPKK+E+ LRQQG  G  YRLLFGD KE+S
Sbjct: 1  MTVTVSSILVSIVYVAVLRWAWRVLNWVWFRPKKVERCLRQQGFAGKPYRLLFGDWKESS 60

Query: 61 IELKEAKARPLSLDDNIAIRVNPFLHKLVNDYG 93
            LKEA+ +P+ L D +  RV PFLH+LV DYG
Sbjct: 61 DMLKEARTKPIGLSDALLPRVMPFLHQLVKDYG 93

Source: Populus trichocarpa

Species: Populus trichocarpa

Genus: Populus

Family: Salicaceae

Order: Malpighiales

Phylum: Streptophyta

Superkingdom: Eukaryota

>gi|224119134|ref|XP_002331333.1| cytochrome P450 [Populus trichocarpa] gi|222873916|gb|EEF11047.1| cytochrome P450 [Populus trichocarpa]

 Score =  120 bits (302), Expect = 7e-26,   Method: Compositional matrix adjust.
 Identities = 57/93 (61%), Positives = 69/93 (74%)

Query: 1  MELPLKSIALTIVIVTVLTWAWRVLNWVWLRPKKLEKFLRQQGLKGNSYRLLFGDLKENS 60
          M + + SI ++IV V VL WAWRVLNWVW RPKK+E+ LRQQG  G  YRLLFGD KE+S
Sbjct: 1  MTVTVTSILVSIVYVAVLRWAWRVLNWVWFRPKKVERCLRQQGFAGKPYRLLFGDWKESS 60

Query: 61 IELKEAKARPLSLDDNIAIRVNPFLHKLVNDYG 93
            LKEA+ +P+ L D +  RV PFLH+LV DYG
Sbjct: 61 DMLKEARTKPIGLSDALLPRVMPFLHQLVKDYG 93

Source: Populus trichocarpa

Species: Populus trichocarpa

Genus: Populus

Family: Salicaceae

Order: Malpighiales

Phylum: Streptophyta

Superkingdom: Eukaryota

>gi|224119146|ref|XP_002331336.1| predicted protein [Populus trichocarpa] gi|222873919|gb|EEF11050.1| predicted protein [Populus trichocarpa]

 Score =  120 bits (300), Expect = 1e-25,   Method: Compositional matrix adjust.
 Identities = 56/93 (60%), Positives = 69/93 (74%)

Query: 1  MELPLKSIALTIVIVTVLTWAWRVLNWVWLRPKKLEKFLRQQGLKGNSYRLLFGDLKENS 60
          M + + SI ++IV V V+ WAWRVLNWVW RPKK+E+ LRQQG  G  YRLLFGD KE+S
Sbjct: 1  MTVTVSSILVSIVYVAVVRWAWRVLNWVWFRPKKVERCLRQQGFAGKPYRLLFGDWKESS 60

Query: 61 IELKEAKARPLSLDDNIAIRVNPFLHKLVNDYG 93
            LKEA+ +P+ L D +  RV PFLH+LV DYG
Sbjct: 61 DMLKEARTKPIGLSDALLPRVMPFLHQLVKDYG 93

Source: Populus trichocarpa

Species: Populus trichocarpa

Genus: Populus

Family: Salicaceae

Order: Malpighiales

Phylum: Streptophyta

Superkingdom: Eukaryota

>gi|359494197|ref|XP_002265784.2| PREDICTED: secologanin synthase-like [Vitis vinifera]

 Score =  114 bits (285), Expect = 7e-24,   Method: Compositional matrix adjust.
 Identities = 52/93 (55%), Positives = 68/93 (73%)

Query: 1   MELPLKSIALTIVIVTVLTWAWRVLNWVWLRPKKLEKFLRQQGLKGNSYRLLFGDLKENS 60
           +E+ L S+A++ V +T+L +AWR+LNWVWLRPKKLE+ LRQQGL GNSYRLL GD +E S
Sbjct: 10  VEMKLSSVAISFVFITLLIYAWRLLNWVWLRPKKLERCLRQQGLTGNSYRLLHGDFREMS 69

Query: 61  IELKEAKARPLSLDDNIAIRVNPFLHKLVNDYG 93
             + EA +RP+S  D+I  RV PF    +  YG
Sbjct: 70  RMINEANSRPISFSDDIVQRVLPFHDHSIQKYG 102

Source: Vitis vinifera

Species: Vitis vinifera

Genus: Vitis

Family: Vitaceae

Order: Vitales

Phylum: Streptophyta

Superkingdom: Eukaryota

>gi|359494591|ref|XP_003634812.1| PREDICTED: secologanin synthase-like [Vitis vinifera]

 Score =  114 bits (285), Expect = 8e-24,   Method: Compositional matrix adjust.
 Identities = 52/93 (55%), Positives = 68/93 (73%)

Query: 1   MELPLKSIALTIVIVTVLTWAWRVLNWVWLRPKKLEKFLRQQGLKGNSYRLLFGDLKENS 60
           +E+ L S+A++ V +T+L +AWR+LNWVWLRPKKLE+ LRQQGL GNSYRLL GD +E S
Sbjct: 10  VEMKLSSVAISFVFITLLIYAWRLLNWVWLRPKKLERCLRQQGLTGNSYRLLHGDFREMS 69

Query: 61  IELKEAKARPLSLDDNIAIRVNPFLHKLVNDYG 93
             + EA +RP+S  D+I  RV PF    +  YG
Sbjct: 70  RMINEANSRPISFSDDIVQRVLPFHDHSIQKYG 102

Source: Vitis vinifera

Species: Vitis vinifera

Genus: Vitis

Family: Vitaceae

Order: Vitales

Phylum: Streptophyta

Superkingdom: Eukaryota

>gi|147795635|emb|CAN67740.1| hypothetical protein VITISV_015554 [Vitis vinifera]

 Score =  114 bits (284), Expect = 1e-23,   Method: Compositional matrix adjust.
 Identities = 49/86 (56%), Positives = 66/86 (76%)

Query: 8  IALTIVIVTVLTWAWRVLNWVWLRPKKLEKFLRQQGLKGNSYRLLFGDLKENSIELKEAK 67
          +AL+   +T+L +AWRVLNW+WLRPK+LE+ L+QQGL GNSYRLL+GD KE S+ +KEA 
Sbjct: 9  VALSFAFITILIYAWRVLNWMWLRPKRLERCLKQQGLAGNSYRLLYGDFKEMSMMIKEAT 68

Query: 68 ARPLSLDDNIAIRVNPFLHKLVNDYG 93
          +RP+S  D+I  RV PF +  +  YG
Sbjct: 69 SRPISFSDDILQRVAPFHYHSIKKYG 94

Source: Vitis vinifera

Species: Vitis vinifera

Genus: Vitis

Family: Vitaceae

Order: Vitales

Phylum: Streptophyta

Superkingdom: Eukaryota

>gi|147858656|emb|CAN80407.1| hypothetical protein VITISV_020336 [Vitis vinifera]

 Score =  112 bits (280), Expect = 3e-23,   Method: Compositional matrix adjust.
 Identities = 51/91 (56%), Positives = 66/91 (72%)

Query: 3  LPLKSIALTIVIVTVLTWAWRVLNWVWLRPKKLEKFLRQQGLKGNSYRLLFGDLKENSIE 62
          + L S+A++ V +T+L +AWR+LNWVWLRPKKLE+ LRQQGL GNSYRLL GD +E S  
Sbjct: 1  MKLSSVAISFVFITLLIYAWRLLNWVWLRPKKLERCLRQQGLTGNSYRLLHGDFREMSXM 60

Query: 63 LKEAKARPLSLDDNIAIRVNPFLHKLVNDYG 93
          + EA +RP+S  D+I  RV PF    +  YG
Sbjct: 61 INEANSRPISFSDDIVQRVLPFHDHSIQKYG 91

Source: Vitis vinifera

Species: Vitis vinifera

Genus: Vitis

Family: Vitaceae

Order: Vitales

Phylum: Streptophyta

Superkingdom: Eukaryota

>gi|359494208|ref|XP_002270501.2| PREDICTED: secologanin synthase-like [Vitis vinifera]

 Score =  112 bits (279), Expect = 4e-23,   Method: Compositional matrix adjust.
 Identities = 51/91 (56%), Positives = 66/91 (72%)

Query: 3  LPLKSIALTIVIVTVLTWAWRVLNWVWLRPKKLEKFLRQQGLKGNSYRLLFGDLKENSIE 62
          + L S+A++ V +T+L +AWR+LNWVWLRPKKLE+ LRQQGL GNSYRLL GD +E S  
Sbjct: 1  MKLSSVAISFVFITLLIYAWRLLNWVWLRPKKLERCLRQQGLTGNSYRLLHGDFREMSRM 60

Query: 63 LKEAKARPLSLDDNIAIRVNPFLHKLVNDYG 93
          + EA +RP+S  D+I  RV PF    +  YG
Sbjct: 61 INEANSRPISFSDDIVQRVLPFHDHSIQKYG 91

Source: Vitis vinifera

Species: Vitis vinifera

Genus: Vitis

Family: Vitaceae

Order: Vitales

Phylum: Streptophyta

Superkingdom: Eukaryota

---

Prediction of Gene Ontology (GO) Terms

Close Homologs with Gene Ontology terms Detected by BLAST

Original result of BLAST against Gene Ontology (AMIGO)

ID	Length	Definition	Q cover	H cover	Identity	E-value
Query	94
TAIR|locus:2089521	512	CYP72A15 ""cytochrome P450, fa	0.989	0.181	0.451	8.8e-19
TAIR|locus:2089586	512	CYP72A11 ""cytochrome P450, fa	0.989	0.181	0.408	1.8e-17
TAIR|locus:2089596	512	CYP72A13 ""cytochrome P450, fa	0.989	0.181	0.419	6.3e-17
TAIR|locus:2089561	514	CYP72A10 ""cytochrome P450, fa	0.989	0.180	0.419	6.3e-17
TAIR|locus:2089526	512	CYP72A7 ""cytochrome P450, fam	0.904	0.166	0.4	4.6e-16
TAIR|locus:2089531	515	CYP72A8 ""cytochrome P450, fam	0.925	0.168	0.409	1.6e-15
TAIR|locus:2089546	508	CYP72A9 ""cytochrome P450, fam	0.936	0.173	0.414	2e-15
TAIR|locus:2089621	512	CYP72A14 ""cytochrome P450, fa	0.989	0.181	0.397	9.2e-15
TAIR|locus:2041399	572	CYP709B2 ""cytochrome P450, fa	0.914	0.150	0.366	5e-07
UNIPROTKB|Q69XM6	538	CYP734A4 "Cytochrome P450 734A	0.946	0.165	0.329	9.7e-07

TAIR|locus:2089521 CYP72A15 ""cytochrome P450, family 72, subfamily A, polypeptide 15"" [Arabidopsis thaliana (taxid:3702)]

 Score = 233 (87.1 bits), Expect = 8.8e-19, P = 8.8e-19
 Identities = 42/93 (45%), Positives = 61/93 (65%)

Query:     1 MELPLKSIALTIVIVTVLTWAWRVLNWVWLRPKKLEKFLRQQGLKGNSYRLLFGDLKENS 60
             ME+ + S+ +++V+  V  W WR L WVW +PK LE +LR+QGL G  Y  L GDLK+N 
Sbjct:     1 MEISVASVTISVVLAVVSWWIWRTLQWVWFKPKMLEHYLRRQGLAGTPYTPLVGDLKKNF 60

Query:    61 IELKEAKARPLSLDDNIAIRVNPFLHKLVNDYG 93
               L EA+++PL L D+I+ RV P+  ++   YG
Sbjct:    61 TMLSEARSKPLKLTDDISPRVVPYPLQMFKTYG 93

GO:0005506 "iron ion binding" evidence=IEA

GO:0009055 "electron carrier activity" evidence=IEA

GO:0016705 "oxidoreductase activity, acting on paired donors, with incorporation or reduction of molecular oxygen" evidence=IEA

GO:0019825 "oxygen binding" evidence=ISS

GO:0020037 "heme binding" evidence=IEA

GO:0055114 "oxidation-reduction process" evidence=IEA

TAIR|locus:2089586 CYP72A11 ""cytochrome P450, family 72, subfamily A, polypeptide 11"" [Arabidopsis thaliana (taxid:3702)]

 Score = 221 (82.9 bits), Expect = 1.8e-17, P = 1.8e-17
 Identities = 38/93 (40%), Positives = 61/93 (65%)

Query:     1 MELPLKSIALTIVIVTVLTWAWRVLNWVWLRPKKLEKFLRQQGLKGNSYRLLFGDLKENS 60
             ME+ + S+ +++ +V V  W WR L WVW +PK LE +LR+QGL G  Y  L GDLK+N 
Sbjct:     1 MEISVASVTVSVAVVVVSWWVWRTLQWVWFKPKMLESYLRRQGLAGTPYTPLVGDLKKNF 60

Query:    61 IELKEAKARPLSLDDNIAIRVNPFLHKLVNDYG 93
                 EA+++P++L D+I  R+ P+  +++  +G
Sbjct:    61 SMRAEARSKPINLTDDITPRIVPYPLQMLKTHG 93

GO:0005506 "iron ion binding" evidence=IEA

GO:0009055 "electron carrier activity" evidence=IEA

GO:0016705 "oxidoreductase activity, acting on paired donors, with incorporation or reduction of molecular oxygen" evidence=IEA

GO:0019825 "oxygen binding" evidence=ISS

GO:0020037 "heme binding" evidence=IEA

GO:0055114 "oxidation-reduction process" evidence=IEA

TAIR|locus:2089596 CYP72A13 ""cytochrome P450, family 72, subfamily A, polypeptide 13"" [Arabidopsis thaliana (taxid:3702)]

 Score = 216 (81.1 bits), Expect = 6.3e-17, P = 6.3e-17
 Identities = 39/93 (41%), Positives = 61/93 (65%)

Query:     1 MELPLKSIALTIVIVTVLTWAWRVLNWVWLRPKKLEKFLRQQGLKGNSYRLLFGDLKENS 60
             ME+ + S+ +++ +V V  W WR L  VWL+PK LE +LR+QGL G  Y  L GDLK N 
Sbjct:     1 MEISVASVTVSVAVVVVSWWVWRTLQRVWLKPKMLESYLRRQGLAGTPYTPLVGDLKRNF 60

Query:    61 IELKEAKARPLSLDDNIAIRVNPFLHKLVNDYG 93
               L EA+++P++L D+I  R+ P+  +++  +G
Sbjct:    61 SMLAEARSKPINLTDDITPRIVPYPLQMLKTHG 93

GO:0005506 "iron ion binding" evidence=IEA

GO:0009055 "electron carrier activity" evidence=IEA

GO:0016705 "oxidoreductase activity, acting on paired donors, with incorporation or reduction of molecular oxygen" evidence=IEA

GO:0019825 "oxygen binding" evidence=ISS

GO:0020037 "heme binding" evidence=IEA

GO:0055114 "oxidation-reduction process" evidence=IEA

TAIR|locus:2089561 CYP72A10 ""cytochrome P450, family 72, subfamily A, polypeptide 10"" [Arabidopsis thaliana (taxid:3702)]

 Score = 216 (81.1 bits), Expect = 6.3e-17, P = 6.3e-17
 Identities = 39/93 (41%), Positives = 58/93 (62%)

Query:     1 MELPLKSIALTIVIVTVLTWAWRVLNWVWLRPKKLEKFLRQQGLKGNSYRLLFGDLKENS 60
             ME+ +  + +++V+V V  W WR L WVW +PK LE +LR+QGL G  Y  L GDLK N 
Sbjct:     1 MEISVACVTVSVVVVVVSWWVWRTLQWVWFKPKMLESYLRRQGLAGTPYTPLIGDLKRNV 60

Query:    61 IELKEAKARPLSLDDNIAIRVNPFLHKLVNDYG 93
               L EA ++P+ L ++I  RV P   +++  +G
Sbjct:    61 NMLTEATSKPIKLTEDITPRVLPHPFQMLKTHG 93

GO:0005506 "iron ion binding" evidence=IEA

GO:0009055 "electron carrier activity" evidence=IEA

GO:0016705 "oxidoreductase activity, acting on paired donors, with incorporation or reduction of molecular oxygen" evidence=IEA

GO:0019825 "oxygen binding" evidence=ISS

GO:0020037 "heme binding" evidence=IEA

GO:0055114 "oxidation-reduction process" evidence=IEA

TAIR|locus:2089526 CYP72A7 ""cytochrome P450, family 72, subfamily A, polypeptide 7"" [Arabidopsis thaliana (taxid:3702)]

 Score = 208 (78.3 bits), Expect = 4.6e-16, P = 4.6e-16
 Identities = 34/85 (40%), Positives = 57/85 (67%)

Query:     9 ALTIVIVTVLTWAWRVLNWVWLRPKKLEKFLRQQGLKGNSYRLLFGDLKENSIELKEAKA 68
             AL +++  V+ W WR++ WVW++PK LE  L++QGL G  Y  L GD+K N   + EA++
Sbjct:     8 ALPVLVAVVVLWTWRIVKWVWIKPKMLESSLKRQGLTGTPYTPLVGDIKRNVDMMMEARS 67

Query:    69 RPLSLDDNIAIRVNPFLHKLVNDYG 93
             +P+++ D+I  R+ P   K++N +G
Sbjct:    68 KPINVTDDITPRLLPLALKMLNSHG 92

GO:0005506 "iron ion binding" evidence=IEA

GO:0009055 "electron carrier activity" evidence=IEA

GO:0016705 "oxidoreductase activity, acting on paired donors, with incorporation or reduction of molecular oxygen" evidence=IEA

GO:0019825 "oxygen binding" evidence=ISS

GO:0020037 "heme binding" evidence=IEA

GO:0055114 "oxidation-reduction process" evidence=IEA

GO:0005783 "endoplasmic reticulum" evidence=IDA

TAIR|locus:2089531 CYP72A8 ""cytochrome P450, family 72, subfamily A, polypeptide 8"" [Arabidopsis thaliana (taxid:3702)]

 Score = 203 (76.5 bits), Expect = 1.6e-15, P = 1.6e-15
 Identities = 36/88 (40%), Positives = 58/88 (65%)

Query:     7 SIALTIVIVTVLT-WAWRVLNWVWLRPKKLEKFLRQQGLKGNSYRLLFGDLKENSIELKE 65
             ++A  +V+VT +T W W+ LN  WLRPKK E +L++QGL G  +  L GD+K  +  +++
Sbjct:    10 AVAAAVVVVTTVTVWIWKGLNVAWLRPKKNEAYLKRQGLSGTPFTFLVGDIKREASMVEQ 69

Query:    66 AKARPLSLDDNIAIRVNPFLHKLVNDYG 93
              K+RP++L D+   RV P + + V D+G
Sbjct:    70 EKSRPINLTDDYTHRVMPLIQQTVKDHG 97

GO:0005506 "iron ion binding" evidence=IEA

GO:0009055 "electron carrier activity" evidence=IEA

GO:0016705 "oxidoreductase activity, acting on paired donors, with incorporation or reduction of molecular oxygen" evidence=IEA

GO:0019825 "oxygen binding" evidence=ISS

GO:0020037 "heme binding" evidence=IEA

GO:0055114 "oxidation-reduction process" evidence=IEA

GO:0007623 "circadian rhythm" evidence=IEP

TAIR|locus:2089546 CYP72A9 ""cytochrome P450, family 72, subfamily A, polypeptide 9"" [Arabidopsis thaliana (taxid:3702)]

 Score = 202 (76.2 bits), Expect = 2.0e-15, P = 2.0e-15
 Identities = 39/94 (41%), Positives = 60/94 (63%)

Query:     1 MELPLKSIALTIVIVTVLTWA-WRVLNWVWLRPKKLEKFLRQQGLKGNSYRLLFGDLKEN 59
             ME+ + S+AL +V+     W  WR+L WVWL+PK LE +LR+QGL G  Y  L GD++ +
Sbjct:     1 MEIVIASLALVVVL-----WCIWRILEWVWLKPKMLESYLRRQGLVGTRYTPLVGDVRRS 55

Query:    60 SIELKEAKARPLSLDDNIAIRVNPFLHKLVNDYG 93
                LKEA+++P+   D++   V P+   ++N YG
Sbjct:    56 FSMLKEARSKPMKPTDDLISLVMPYSFHMLNTYG 89

GO:0005506 "iron ion binding" evidence=IEA

GO:0009055 "electron carrier activity" evidence=IEA

GO:0016705 "oxidoreductase activity, acting on paired donors, with incorporation or reduction of molecular oxygen" evidence=IEA

GO:0019825 "oxygen binding" evidence=ISS

GO:0020037 "heme binding" evidence=IEA

GO:0055114 "oxidation-reduction process" evidence=IEA

TAIR|locus:2089621 CYP72A14 ""cytochrome P450, family 72, subfamily A, polypeptide 14"" [Arabidopsis thaliana (taxid:3702)]

 Score = 196 (74.1 bits), Expect = 9.2e-15, P = 9.2e-15
 Identities = 37/93 (39%), Positives = 54/93 (58%)

Query:     1 MELPLKSIALTIVIVTVLTWAWRVLNWVWLRPKKLEKFLRQQGLKGNSYRLLFGDLKENS 60
             ME+ + S+  ++ +V V  W WR L WVW  PK LE+ LR+QGL G SY  L GD K+  
Sbjct:     1 MEISVSSVTFSLAVVVVSWWVWRTLKWVWFTPKMLERSLRRQGLSGTSYTPLIGDFKKMI 60

Query:    61 IELKEAKARPLSLDDNIAIRVNPFLHKLVNDYG 93
                 EA ++P+   D+I  RV P   +++  +G
Sbjct:    61 SMFIEATSKPIKPTDDITPRVMPHPLQMLKTHG 93

GO:0005506 "iron ion binding" evidence=IEA

GO:0009055 "electron carrier activity" evidence=IEA

GO:0016705 "oxidoreductase activity, acting on paired donors, with incorporation or reduction of molecular oxygen" evidence=IEA

GO:0019825 "oxygen binding" evidence=ISS

GO:0020037 "heme binding" evidence=IEA

GO:0055114 "oxidation-reduction process" evidence=IEA

GO:0000041 "transition metal ion transport" evidence=RCA

TAIR|locus:2041399 CYP709B2 ""cytochrome P450, family 709, subfamily B, polypeptide 2"" [Arabidopsis thaliana (taxid:3702)]

 Score = 125 (49.1 bits), Expect = 5.0e-07, P = 5.0e-07
 Identities = 33/90 (36%), Positives = 50/90 (55%)

Query:     5 LKSIALTIVIVTVLTWAWRVLNWVWLRPKKLEKFLRQQGLKGNSYRLLFGDLKENSIELK 64
             L +IAL +++V  +  A R+L W   RP  L +  ++QG+ G  YR+L+G+L+E      
Sbjct:    64 LLAIALVLLVVPKIYGACRILVW---RPWMLSRRFKKQGISGPKYRILYGNLREIRKMKN 120

Query:    65 EAKARPLSLDDN-IAIRVNPFLHKLVNDYG 93
             EAK   L  + N I  RV P L +  + YG
Sbjct:   121 EAKLMVLDPNSNDIVPRVLPHLQQWKSQYG 150

GO:0005506 "iron ion binding" evidence=IEA

GO:0009055 "electron carrier activity" evidence=IEA

GO:0009507 "chloroplast" evidence=ISM

GO:0016705 "oxidoreductase activity, acting on paired donors, with incorporation or reduction of molecular oxygen" evidence=IEA

GO:0019825 "oxygen binding" evidence=ISS

GO:0020037 "heme binding" evidence=IEA

GO:0055114 "oxidation-reduction process" evidence=IEA

UNIPROTKB|Q69XM6 CYP734A4 "Cytochrome P450 734A4" [Oryza sativa Japonica Group (taxid:39947)]

 Score = 122 (48.0 bits), Expect = 9.7e-07, P = 9.7e-07
 Identities = 30/91 (32%), Positives = 49/91 (53%)

Query:     5 LKSIALTIVIVTVLTWAWRVLNWVWLRPKKLEKFLRQQGLKGNSYRLLFGDLKENSIELK 64
             ++++A+   ++ +L  A RV + VW RP++LE     QG++G  YR L G ++E    + 
Sbjct:     2 MEAVAVAAAVLLLLHVAARVADAVWWRPRRLEAHFAGQGVRGPPYRFLVGCVREMVALMA 61

Query:    65 EAKARPL--SLDDNIAIRVNPFLHKLVNDYG 93
             EA A+P+  +   N   RV  F H     YG
Sbjct:    62 EATAKPMPPAAPHNALPRVLAFYHYWRKIYG 92

GO:0004497 "monooxygenase activity" evidence=IDA

GO:0010268 "brassinosteroid homeostasis" evidence=IDA

GO:0016131 "brassinosteroid metabolic process" evidence=IDA

---

Prediction of Enzyme Commission (EC) Number

EC Number Prediction by Annotation Transfer from SWISS-PROT Entries

Original result of BLAST against SWISS-PROT

No confident hit for EC number transfering in SWISSPROT detected by BLAST

EC Number Prediction by Ezypred Server

Original result from Ezypred Server

1st Layer

is NOT an Enzyme protein

EC Number Prediction by EFICAz Software

No EC number assignment, probably not an enzyme!

---

Prediction of Functionally Associated Proteins

Functionally Associated Proteins Detected by STRING

Original result from the STRING server

Your Input: CYP72A47 cytochrome P450 (518 aa) (Populus trichocarpa)

Predicted Functional Partners:

Sorry, there are no predicted associations at the current settings.

---

Conserved Domains and Related Protein Families

Conserved Domains Detected by RPS-BLAST

Original result of RPS-BLAST against CDD database part I

ID	Length	Definition	E-value
Query	94
PLN02290	516	PLN02290, PLN02290, cytokinin trans-hydroxylase	3e-05

>gnl|CDD|215164 PLN02290, PLN02290, cytokinin trans-hydroxylase

 Score = 40.6 bits (95), Expect = 3e-05
 Identities = 18/75 (24%), Positives = 44/75 (58%), Gaps = 1/75 (1%)

Query: 10 LTIVIVTVLTWAWRVLNWVWLRPKKLEKFLRQQGLKGNSYRLLFGDLKENSIELKEAKAR 69
          L I +  +L  A+  ++  +L P++++K + +QG++G   R L G++ + S  + ++ ++
Sbjct: 10 LVIFLTLLLRVAYDTISCYFLTPRRIKKIMERQGVRGPKPRPLTGNILDVSALVSQSTSK 69

Query: 70 PL-SLDDNIAIRVNP 83
           + S+  +I  R+ P
Sbjct: 70 DMDSIHHDIVGRLLP 84

Length = 516

Conserved Domains Detected by HHsearch

Original result of HHsearch against CDD database

ID	Length	Definition	Probability
Query	94
PLN02290	516	cytokinin trans-hydroxylase	99.35
PTZ00404	482	cytochrome P450; Provisional	97.37
KOG0157	497	consensus Cytochrome P450 CYP4/CYP19/CYP26 subfami	96.88
KOG0158	499	consensus Cytochrome P450 CYP3/CYP5/CYP6/CYP9 subf	96.53
PLN02302	490	ent-kaurenoic acid oxidase	96.08
PLN02183	516	ferulate 5-hydroxylase	95.69
KOG0156	489	consensus Cytochrome P450 CYP2 subfamily [Secondar	95.4
PLN00168	519	Cytochrome P450; Provisional	95.2
PLN03234	499	cytochrome P450 83B1; Provisional	94.86
COG3898	531	Uncharacterized membrane-bound protein [Function u	94.82
PLN02687	517	flavonoid 3'-monooxygenase	94.4
PLN02738	633	carotene beta-ring hydroxylase	94.3
PLN03112	514	cytochrome P450 family protein; Provisional	93.65
PLN02394	503	trans-cinnamate 4-monooxygenase	93.37
PLN02966	502	cytochrome P450 83A1	93.09
PLN02774	463	brassinosteroid-6-oxidase	92.82
PLN02196	463	abscisic acid 8'-hydroxylase	92.77
PLN02169	500	fatty acid (omega-1)-hydroxylase/midchain alkane h	90.65
PLN02971	543	tryptophan N-hydroxylase	90.49
PF07219	108	HemY_N: HemY protein N-terminus; InterPro: IPR0108	90.24
PLN03195	516	fatty acid omega-hydroxylase; Provisional	88.28
PF00067	463	p450: Cytochrome P450 p450 superfamily signature b	86.84
PLN02500	490	cytochrome P450 90B1	86.26
PLN02655	466	ent-kaurene oxidase	82.08

>PLN02290 cytokinin trans-hydroxylase

Probab=99.35  E-value=1.5e-11  Score=96.70  Aligned_cols=78  Identities=22%  Similarity=0.478  Sum_probs=65.8

Q ss_pred             HHHHHHHHHHHHHhchHHHHHHHHHcCCCCCCCcCCCCCHHHHHHHHHHhhcCCCC-CCCCcccccchHHHHHHHhhCC
Q 034451           17 VLTWAWRVLNWVWLRPKKLEKFLRQQGLKGNSYRLLFGDLKENSIELKEAKARPLS-LDDNIAIRVNPFLHKLVNDYGM   94 (94)
Q Consensus        17 ~~~~~~~~~~~~w~~P~rl~~~lr~QGI~GP~~~~l~Gn~~e~~~~~~~a~~~~~~-~shdi~prv~P~~~~w~k~YGk   94 (94)
                      +.-++++.++.+.+.|+|++++++++|++||++.+++||+.++.....++..++++ ..+|+..+..+++.+|.++||+
T Consensus        17 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~PGP~~~P~iGnl~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~yG~   95 (516)
T PLN02290         17 LLRVAYDTISCYFLTPRRIKKIMERQGVRGPKPRPLTGNILDVSALVSQSTSKDMDSIHHDIVGRLLPHYVAWSKQYGK   95 (516)
T ss_pred             HHHHHHHHHHHHhccHHHHHHHHHHcCCCCCCCCcCCCCHHHHHHHHHHhhcCCCCCCCcccccccchHHHHHHHHhCC
Confidence            34466788899999999999999999999999999999999997765565566665 4677777788899999999995

>PTZ00404 cytochrome P450; Provisional

Probab=97.37  E-value=0.00028  Score=55.14  Aligned_cols=47  Identities=23%  Similarity=0.216  Sum_probs=34.6

Q ss_pred             hchHHHHHHHHHcCCCCCCCcCCCCCHHHHHHHHHHhhcCCCCCCCCcccccchHHHHHHHhhCC
Q 034451           30 LRPKKLEKFLRQQGLKGNSYRLLFGDLKENSIELKEAKARPLSLDDNIAIRVNPFLHKLVNDYGM   94 (94)
Q Consensus        30 ~~P~rl~~~lr~QGI~GP~~~~l~Gn~~e~~~~~~~a~~~~~~~shdi~prv~P~~~~w~k~YGk   94 (94)
                      +.+.+..+++++++|+||++.+++||..++.+                  ...-++++|.++||+
T Consensus        17 ~~~~~~~~~~~~~~~pgp~~~p~~G~~~~~~~------------------~~~~~~~~~~~~yG~   63 (482)
T PTZ00404         17 HNAYKKYKKIHKNELKGPIPIPILGNLHQLGN------------------LPHRDLTKMSKKYGG   63 (482)
T ss_pred             HHHHHHhhhccCCCCCCCCCCCeeccHhhhcc------------------cHHHHHHHHHHHhCC
Confidence            33555556799999999999999999866632                  112367888999984

>KOG0157 consensus Cytochrome P450 CYP4/CYP19/CYP26 subfamilies [Secondary metabolites biosynthesis, transport and catabolism; Lipid transport and metabolism]

Probab=96.88  E-value=0.0037  Score=49.89  Aligned_cols=39  Identities=26%  Similarity=0.341  Sum_probs=30.7

Q ss_pred             HHcCCCCCCCcCCCCCHHHHHHHHHHhhcCCCCCCCCcccccchHHHHHHHhhC
Q 034451           40 RQQGLKGNSYRLLFGDLKENSIELKEAKARPLSLDDNIAIRVNPFLHKLVNDYG   93 (94)
Q Consensus        40 r~QGI~GP~~~~l~Gn~~e~~~~~~~a~~~~~~~shdi~prv~P~~~~w~k~YG   93 (94)
                      +.++.+||++.+++||..|+.++.               .....++.++..+||
T Consensus        33 ~~~~~~gp~~~P~iG~~~~~~~~~---------------~~~~~~~~~~~~~~~   71 (497)
T KOG0157|consen   33 KKKLPPGPPGWPLIGNLLEFLKPL---------------EEILDFVTELLSRYG   71 (497)
T ss_pred             HhccCCCCCCCCcccchHHhhcch---------------hHHHHHHHHHHHHcC
Confidence            788999999999999999987631               334457777777776

>KOG0158 consensus Cytochrome P450 CYP3/CYP5/CYP6/CYP9 subfamilies [Secondary metabolites biosynthesis, transport and catabolism]

Probab=96.53  E-value=0.0054  Score=49.84  Aligned_cols=29  Identities=17%  Similarity=0.244  Sum_probs=25.1

Q ss_pred             HHHHHHHHcCCCCCCCcCCCCCHHHHHHH
Q 034451           34 KLEKFLRQQGLKGNSYRLLFGDLKENSIE   62 (94)
Q Consensus        34 rl~~~lr~QGI~GP~~~~l~Gn~~e~~~~   62 (94)
                      +...+|++.||+||+|.+++||+..+.+.
T Consensus        23 ~~~~yw~rrGi~~~~p~p~~Gn~~~~~~~   51 (499)
T KOG0158|consen   23 WTYSYWRRRGIPGPKPLPFLGNLPGMLKR   51 (499)
T ss_pred             hhhhhhccCCCCCCCCCCcEecHHHHHhc
Confidence            45568889999999999999999998773

>PLN02302 ent-kaurenoic acid oxidase

Probab=96.08  E-value=0.024  Score=44.07  Aligned_cols=37  Identities=16%  Similarity=0.310  Sum_probs=24.4

Q ss_pred             CCCCCcCCCCCHHHHHHHHHHhhcCCCCCCCCcccccchHHHHHHHhhCC
Q 034451           45 KGNSYRLLFGDLKENSIELKEAKARPLSLDDNIAIRVNPFLHKLVNDYGM   94 (94)
Q Consensus        45 ~GP~~~~l~Gn~~e~~~~~~~a~~~~~~~shdi~prv~P~~~~w~k~YGk   94 (94)
                      +||++-+++||+.++.....          +   .....++.+|.++||+
T Consensus        45 Pgp~~~PilG~l~~~~~~~~----------~---~~~~~~~~~~~~kyG~   81 (490)
T PLN02302         45 PGDLGWPVIGNMWSFLRAFK----------S---SNPDSFIASFISRYGR   81 (490)
T ss_pred             CCCCCCCccccHHHHHHhcc----------c---CCcHHHHHHHHHHhCC
Confidence            67888899999877754100          0   1112467899999984

>PLN02183 ferulate 5-hydroxylase

Probab=95.69  E-value=0.019  Score=45.54  Aligned_cols=46  Identities=15%  Similarity=0.043  Sum_probs=31.6

Q ss_pred             HHHHhchHHHHHHHHHcCCCCCCCcCCCCCHHHHHHHHHHhhcCCCCCCCCcccccchHHHHHHHhhCC
Q 034451           26 NWVWLRPKKLEKFLRQQGLKGNSYRLLFGDLKENSIELKEAKARPLSLDDNIAIRVNPFLHKLVNDYGM   94 (94)
Q Consensus        26 ~~~w~~P~rl~~~lr~QGI~GP~~~~l~Gn~~e~~~~~~~a~~~~~~~shdi~prv~P~~~~w~k~YGk   94 (94)
                      -.+|.++.|-     .+-++||++.+++||+.++.+              +    ...++.+|.++||+
T Consensus        25 ~~~~~~~~~~-----~~~ppgp~~~Pl~G~l~~~~~--------------~----~~~~~~~~~~~yG~   70 (516)
T PLN02183         25 LGLISRLRRR-----LPYPPGPKGLPIIGNMLMMDQ--------------L----THRGLANLAKQYGG   70 (516)
T ss_pred             HHHHhhccCC-----CCCCcCCCCCCeeccHHhcCC--------------c----chHHHHHHHHHhCC
Confidence            3446666552     356889999999999865521              1    12468999999995

>KOG0156 consensus Cytochrome P450 CYP2 subfamily [Secondary metabolites biosynthesis, transport and catabolism]

Probab=95.40  E-value=0.023  Score=45.88  Aligned_cols=34  Identities=21%  Similarity=0.238  Sum_probs=26.1

Q ss_pred             CCCCCCcCCCCCHHHHHHHHHHhhcCCCCCCCCcccccchHHHHHHHhhCC
Q 034451           44 LKGNSYRLLFGDLKENSIELKEAKARPLSLDDNIAIRVNPFLHKLVNDYGM   94 (94)
Q Consensus        44 I~GP~~~~l~Gn~~e~~~~~~~a~~~~~~~shdi~prv~P~~~~w~k~YGk   94 (94)
                      -+||++-+++||+.++..                .| ..-.|++|+++||+
T Consensus        28 PPGP~~lPiIGnl~~l~~----------------~~-~h~~~~~ls~~yGp   61 (489)
T KOG0156|consen   28 PPGPPPLPIIGNLHQLGS----------------LP-PHRSFRKLSKKYGP   61 (489)
T ss_pred             CcCCCCCCccccHHHcCC----------------Cc-hhHHHHHHHHHhCC
Confidence            469999999999977754                11 12378999999995

>PLN00168 Cytochrome P450; Provisional

Probab=95.20  E-value=0.027  Score=44.67  Aligned_cols=53  Identities=19%  Similarity=0.176  Sum_probs=35.3

Q ss_pred             HHHhchHHHHHHHHHcCCCCCCCcCCCCCHHHHHHHHHHhhcCCCCCCCCcccccchHHHHHHHhhCC
Q 034451           27 WVWLRPKKLEKFLRQQGLKGNSYRLLFGDLKENSIELKEAKARPLSLDDNIAIRVNPFLHKLVNDYGM   94 (94)
Q Consensus        27 ~~w~~P~rl~~~lr~QGI~GP~~~~l~Gn~~e~~~~~~~a~~~~~~~shdi~prv~P~~~~w~k~YGk   94 (94)
                      .+|+.|.+..-..+.+=.+||++.+++||..++..              + .+....++.+|.++||+
T Consensus        20 ~~~~~~~~~~~~~~~~lpPgp~~~pl~G~l~~~~~--------------~-~~~~~~~~~~~~~~yG~   72 (519)
T PLN00168         20 LLLGKHGGRGGKKGRRLPPGPPAVPLLGSLVWLTN--------------S-SADVEPLLRRLIARYGP   72 (519)
T ss_pred             HhhhhhhccCCCCCCCCCcCCCCCcccccHHhhcc--------------c-cccHHHHHHHHHHHhCC
Confidence            44566666655555555679999999999865521              0 11233578999999984

>PLN03234 cytochrome P450 83B1; Provisional

Probab=94.86  E-value=0.042  Score=43.11  Aligned_cols=35  Identities=26%  Similarity=0.262  Sum_probs=25.5

Q ss_pred             CCCCCCCcCCCCCHHHHHHHHHHhhcCCCCCCCCcccccchHHHHHHHhhCC
Q 034451           43 GLKGNSYRLLFGDLKENSIELKEAKARPLSLDDNIAIRVNPFLHKLVNDYGM   94 (94)
Q Consensus        43 GI~GP~~~~l~Gn~~e~~~~~~~a~~~~~~~shdi~prv~P~~~~w~k~YGk   94 (94)
                      .++||++.+++||+.++..             +|    ...++.+|.++||+
T Consensus        29 ~pPgp~~~P~iG~~~~~~~-------------~~----~~~~~~~~~~~yG~   63 (499)
T PLN03234         29 LPPGPKGLPIIGNLHQMEK-------------FN----PQHFLFRLSKLYGP   63 (499)
T ss_pred             CCcCCCCCCeeccHHhcCC-------------CC----ccHHHHHHHHHcCC
Confidence            4899999999999866521             11    12477889999984

>COG3898 Uncharacterized membrane-bound protein [Function unknown]

Probab=94.82  E-value=0.1  Score=42.22  Aligned_cols=50  Identities=20%  Similarity=0.409  Sum_probs=39.7

Q ss_pred             HHHHHHHHHHHHhchHHHHHHHH----HcCCCCCCCcCC---CCCHHHHHHHHHHhh
Q 034451           18 LTWAWRVLNWVWLRPKKLEKFLR----QQGLKGNSYRLL---FGDLKENSIELKEAK   67 (94)
Q Consensus        18 ~~~~~~~~~~~w~~P~rl~~~lr----~QGI~GP~~~~l---~Gn~~e~~~~~~~a~   67 (94)
                      +..+|.+++++|-.|++++|++|    .||++--+--++   -||-.+-.+|.+++.
T Consensus        55 v~llwwlv~~iw~sP~t~~Ryfr~rKRdrgyqALStGliAagAGda~lARkmt~~~~  111 (531)
T COG3898          55 VLLLWWLVRSIWESPYTARRYFRERKRDRGYQALSTGLIAAGAGDASLARKMTARAS  111 (531)
T ss_pred             HHHHHHHHHHHHhCcHHHHHHHHHHHhhhHHHHHhhhhhhhccCchHHHHHHHHHHH
Confidence            34678999999999999999986    468876666666   489988888876654

>PLN02687 flavonoid 3'-monooxygenase

Probab=94.40  E-value=0.044  Score=43.46  Aligned_cols=34  Identities=18%  Similarity=0.180  Sum_probs=24.2

Q ss_pred             CCCCCCCcCCCCCHHHHHHHHHHhhcCCCCCCCCcccccchHHHHHHHhhCC
Q 034451           43 GLKGNSYRLLFGDLKENSIELKEAKARPLSLDDNIAIRVNPFLHKLVNDYGM   94 (94)
Q Consensus        43 GI~GP~~~~l~Gn~~e~~~~~~~a~~~~~~~shdi~prv~P~~~~w~k~YGk   94 (94)
                      ..+||++-+++||..++.+              +    ...++.+|.++||+
T Consensus        35 ~pPgp~~~P~iG~~~~~~~--------------~----~~~~~~~~~~~yG~   68 (517)
T PLN02687         35 LPPGPRGWPVLGNLPQLGP--------------K----PHHTMAALAKTYGP   68 (517)
T ss_pred             CCccCCCCCccccHHhcCC--------------c----hhHHHHHHHHHhCC
Confidence            4568888899999865521              1    22477899999984

>PLN02738 carotene beta-ring hydroxylase

Probab=94.30  E-value=0.044  Score=45.35  Aligned_cols=21  Identities=24%  Similarity=0.286  Sum_probs=18.2

Q ss_pred             chHHHHHHHHHcCCCCCCCcC
Q 034451           31 RPKKLEKFLRQQGLKGNSYRL   51 (94)
Q Consensus        31 ~P~rl~~~lr~QGI~GP~~~~   51 (94)
                      .|..++..|+++||+||+...
T Consensus       104 ~~~~~~~~~~~~~~pgp~laa  124 (633)
T PLN02738        104 FPATLRNGLAKLGPPGELLAF  124 (633)
T ss_pred             chHHHHhhhhhCCCCCchHHH
Confidence            589999999999999996543

>PLN03112 cytochrome P450 family protein; Provisional

Probab=93.65  E-value=0.073  Score=42.00  Aligned_cols=33  Identities=21%  Similarity=0.188  Sum_probs=23.1

Q ss_pred             CCCCCCcCCCCCHHHHHHHHHHhhcCCCCCCCCcccccchHHHHHHHhhCC
Q 034451           44 LKGNSYRLLFGDLKENSIELKEAKARPLSLDDNIAIRVNPFLHKLVNDYGM   94 (94)
Q Consensus        44 I~GP~~~~l~Gn~~e~~~~~~~a~~~~~~~shdi~prv~P~~~~w~k~YGk   94 (94)
                      .+||++.+++||+.++.+              +  +  .-++.+|.++||+
T Consensus        34 ppgp~~~pl~G~~~~~~~--------------~--~--~~~~~~~~~kyG~   66 (514)
T PLN03112         34 PPGPPRWPIVGNLLQLGP--------------L--P--HRDLASLCKKYGP   66 (514)
T ss_pred             ccCCCCCCeeeeHHhcCC--------------c--h--HHHHHHHHHHhCC
Confidence            579999999999865521              0  1  1356888888884

>PLN02394 trans-cinnamate 4-monooxygenase

Probab=93.37  E-value=0.14  Score=40.25  Aligned_cols=35  Identities=17%  Similarity=0.191  Sum_probs=25.2

Q ss_pred             CCCCCCCcCCCCCHHHHHHHHHHhhcCCCCCCCCcccccchHHHHHHHhhCC
Q 034451           43 GLKGNSYRLLFGDLKENSIELKEAKARPLSLDDNIAIRVNPFLHKLVNDYGM   94 (94)
Q Consensus        43 GI~GP~~~~l~Gn~~e~~~~~~~a~~~~~~~shdi~prv~P~~~~w~k~YGk   94 (94)
                      -++||++..++||..++.+              |   .....+++|.++||+
T Consensus        31 ~pPgp~~~p~~g~l~~~~~--------------~---~~~~~~~~~~~~yG~   65 (503)
T PLN02394         31 LPPGPAAVPIFGNWLQVGD--------------D---LNHRNLAEMAKKYGD   65 (503)
T ss_pred             CCcCCCCCCeeeeHHhcCC--------------C---chhHHHHHHHHHhCC
Confidence            3679999999999866522              1   012477899999984

>PLN02966 cytochrome P450 83A1

Probab=93.09  E-value=0.063  Score=42.39  Aligned_cols=35  Identities=20%  Similarity=0.232  Sum_probs=25.4

Q ss_pred             CCCCCCCcCCCCCHHHHHHHHHHhhcCCCCCCCCcccccchHHHHHHHhhCC
Q 034451           43 GLKGNSYRLLFGDLKENSIELKEAKARPLSLDDNIAIRVNPFLHKLVNDYGM   94 (94)
Q Consensus        43 GI~GP~~~~l~Gn~~e~~~~~~~a~~~~~~~shdi~prv~P~~~~w~k~YGk   94 (94)
                      =++||++.+++||+.++..                 .....++.+|.++||+
T Consensus        30 ~ppgp~~~p~~G~l~~l~~-----------------~~~~~~~~~~~~~yG~   64 (502)
T PLN02966         30 LPPGPSPLPVIGNLLQLQK-----------------LNPQRFFAGWAKKYGP   64 (502)
T ss_pred             CCcCCCCCCeeccHHhcCC-----------------CChhHHHHHHHHHhCC
Confidence            3789999999999865521                 1123577999999995

>PLN02774 brassinosteroid-6-oxidase

Probab=92.82  E-value=0.13  Score=40.20  Aligned_cols=41  Identities=27%  Similarity=0.299  Sum_probs=28.3

Q ss_pred             HHHHHHcCC-CCCCCcCCCCCHHHHHHHHHHhhcCCCCCCCCcccccchHHHHHHHhhCC
Q 034451           36 EKFLRQQGL-KGNSYRLLFGDLKENSIELKEAKARPLSLDDNIAIRVNPFLHKLVNDYGM   94 (94)
Q Consensus        36 ~~~lr~QGI-~GP~~~~l~Gn~~e~~~~~~~a~~~~~~~shdi~prv~P~~~~w~k~YGk   94 (94)
                      ++++++.|. +||++-+++||..++.+              |  +  .-++++|.++||+
T Consensus        24 ~~~~~r~~~ppgp~~~P~~G~~~~~~~--------------~--~--~~~~~~~~~~yG~   65 (463)
T PLN02774         24 EVRYSKKGLPPGTMGWPLFGETTEFLK--------------Q--G--PDFMKNQRLRYGS   65 (463)
T ss_pred             hhccCCCCCCCCCCCCCchhhHHHHHH--------------h--h--HHHHHHHHHHhcc
Confidence            344566687 48888899999877643              1  1  1267888899985

>PLN02196 abscisic acid 8'-hydroxylase

Probab=92.77  E-value=0.2  Score=39.19  Aligned_cols=36  Identities=17%  Similarity=0.102  Sum_probs=23.2

Q ss_pred             cCCCCCCCcCCCCCHHHHHHHHHHhhcCCCCCCCCcccccchHHHHHHHhhCC
Q 034451           42 QGLKGNSYRLLFGDLKENSIELKEAKARPLSLDDNIAIRVNPFLHKLVNDYGM   94 (94)
Q Consensus        42 QGI~GP~~~~l~Gn~~e~~~~~~~a~~~~~~~shdi~prv~P~~~~w~k~YGk   94 (94)
                      .+.+||++.+++||..++.+             .|  |  .-++.+|.++||+
T Consensus        35 ~~Ppgp~~~P~iG~~~~~~~-------------~~--~--~~~~~~~~~~yG~   70 (463)
T PLN02196         35 PLPPGTMGWPYVGETFQLYS-------------QD--P--NVFFASKQKRYGS   70 (463)
T ss_pred             CCCCCCCCCCccchHHHHHh-------------cC--H--HHHHHHHHHHhhh
Confidence            35556666788999876532             11  1  1357889999984

>PLN02169 fatty acid (omega-1)-hydroxylase/midchain alkane hydroxylase

Probab=90.65  E-value=1.2  Score=35.51  Aligned_cols=23  Identities=13%  Similarity=-0.010  Sum_probs=20.0

Q ss_pred             HHHHcCCCCCCCcCCCCCHHHHH
Q 034451           38 FLRQQGLKGNSYRLLFGDLKENS   60 (94)
Q Consensus        38 ~lr~QGI~GP~~~~l~Gn~~e~~   60 (94)
                      .||+.|++||++.+++||+.++.
T Consensus        27 ~~~~~~~p~p~~~pl~G~~~~~~   49 (500)
T PLN02169         27 HKKPHGQPILKNWPFLGMLPGML   49 (500)
T ss_pred             HhccCCCCCCCCCCcccchHHHH
Confidence            56778999999999999997764

>PLN02971 tryptophan N-hydroxylase

Probab=90.49  E-value=0.48  Score=37.99  Aligned_cols=34  Identities=21%  Similarity=0.201  Sum_probs=24.5

Q ss_pred             CCCCCCcCCCCCHHHHHHHHHHhhcCCCCCCCCcccccchHHHHHHHhhC
Q 034451           44 LKGNSYRLLFGDLKENSIELKEAKARPLSLDDNIAIRVNPFLHKLVNDYG   93 (94)
Q Consensus        44 I~GP~~~~l~Gn~~e~~~~~~~a~~~~~~~shdi~prv~P~~~~w~k~YG   93 (94)
                      .+||++.+++||+.++.+              | .+ ...++++|.++||
T Consensus        59 PPGP~~lPiiGnl~~l~~--------------~-~~-~~~~l~~~~~~yg   92 (543)
T PLN02971         59 PPGPTGFPIVGMIPAMLK--------------N-RP-VFRWLHSLMKELN   92 (543)
T ss_pred             CcCCCCCCcccchHHhcc--------------C-Cc-HhHHHHHHHHHhC
Confidence            569999999999876632              1 11 2246789999998

>PF07219 HemY_N: HemY protein N-terminus; InterPro: IPR010817 This entry represents the N terminus (approximately 150 residues) of bacterial HemY porphyrin biosynthesis proteins

Probab=90.24  E-value=1.3  Score=28.51  Aligned_cols=52  Identities=21%  Similarity=0.269  Sum_probs=33.2

Q ss_pred             HHHHHHHHHHHHHHHhchHHHHHHHHHcCCCCC-------CCcCCCCCHHHHHHHHHHh
Q 034451           15 VTVLTWAWRVLNWVWLRPKKLEKFLRQQGLKGN-------SYRLLFGDLKENSIELKEA   66 (94)
Q Consensus        15 ~~~~~~~~~~~~~~w~~P~rl~~~lr~QGI~GP-------~~~~l~Gn~~e~~~~~~~a   66 (94)
                      +++++.+++++..+|--|.++.+.++++.-+--       =-.+.-||..+-.+...++
T Consensus        27 ~~~l~ll~~ll~~~~~~p~~~~~~~~~rr~~ka~~al~~Gl~al~~G~~~~A~k~~~~a   85 (108)
T PF07219_consen   27 FVVLYLLLRLLRRLLSLPSRVRRWRRRRRRRKAQRALSRGLIALAEGDWQRAEKLLAKA   85 (108)
T ss_pred             HHHHHHHHHHHHHHHhChHHHHHHHHHHHHHHHHHHHHHHHHHHHCCCHHHHHHHHHHH
Confidence            344556789999999999988765544432110       0124568988877765554

These are membrane protein involved in a late step of protoheme IX synthesis [].

>PLN03195 fatty acid omega-hydroxylase; Provisional

Probab=88.28  E-value=1.5  Score=34.62  Aligned_cols=21  Identities=24%  Similarity=0.253  Sum_probs=16.6

Q ss_pred             HHHcCCCCCCCcCCCCCHHHH
Q 034451           39 LRQQGLKGNSYRLLFGDLKEN   59 (94)
Q Consensus        39 lr~QGI~GP~~~~l~Gn~~e~   59 (94)
                      |++++++||++-+++||..++
T Consensus        27 ~~~~~~pgp~~~p~~G~~~~~   47 (516)
T PLN03195         27 WSQRNRKGPKSWPIIGAALEQ   47 (516)
T ss_pred             hhccccCCCCCCCeecchHHH
Confidence            445679999999999997544

>PF00067 p450: Cytochrome P450 p450 superfamily signature b-class p450 signature mitochondrial p450 signature E-class p450 group I signature E-class p450 group II signature E-class p450 group IV signature; InterPro: IPR001128 Cytochrome P450 enzymes are a superfamily of haem-containing mono-oxygenases that are found in all kingdoms of life, and which show extraordinary diversity in their reaction chemistry

Probab=86.84  E-value=0.59  Score=34.51  Aligned_cols=34  Identities=21%  Similarity=0.294  Sum_probs=25.5

Q ss_pred             CCCCCcCCCCCHHHHHHHHHHhhcCCCCCCCCcccccchHHHHHHHhhCC
Q 034451           45 KGNSYRLLFGDLKENSIELKEAKARPLSLDDNIAIRVNPFLHKLVNDYGM   94 (94)
Q Consensus        45 ~GP~~~~l~Gn~~e~~~~~~~a~~~~~~~shdi~prv~P~~~~w~k~YGk   94 (94)
                      +||++-+++||..++.+                .....-++++|.++||+
T Consensus         2 pgp~~~p~~G~~~~~~~----------------~~~~~~~~~~~~~kyG~   35 (463)
T PF00067_consen    2 PGPPPLPILGNLLQFRR----------------KGNPHEFFRELHKKYGP   35 (463)
T ss_dssp             SCSSSBTTTBTHHHHHT----------------THHHHHHHHHHHHHHTS
T ss_pred             cCCCCcCceeEHHHhcC----------------CCcHHHHHHHHHHHhCC
Confidence            58999999999988863                01122367999999995

In mammals, these proteins are found primarily in microsomes of hepatocytes and other cell types, where they oxidise steroids, fatty acids and xenobiotics, and are important for the detoxification and clearance of various compounds, as well as for hormone synthesis and breakdown, cholesterol synthesis and vitamin D metabolism. In plants, these proteins are important for the biosynthesis of several compounds such as hormones, defensive compounds and fatty acids. In bacteria, they are important for several metabolic processes, such as the biosynthesis of antibiotic erythromycin in Saccharopolyspora erythraea (Streptomyces erythraeus). Cytochrome P450 enzymes use haem to oxidise their substrates, using protons derived from NADH or NADPH to split the oxygen so a single atom can be added to a substrate. They also require electrons, which they receive from a variety of redox partners. In certain cases, cytochrome P450 can be fused to its redox partner to produce a bi-functional protein, such as with P450BM-3 from Bacillus megaterium [], which has haem and flavin domains. Organisms produce many different cytochrome P450 enzymes (at least 58 in humans), which together with alternative splicing can provide a wide array of enzymes with different substrate and tissue specificities. Individual cytochrome P450 proteins follow the nomenclature: CYP, followed by a number (family), then a letter (subfamily), and another number (protein); e.g. CYP3A4 is the fourth protein in family 3, subfamily A. In general, family members should share >40% identity, while subfamily members should share >55% identity. Cytochrome P450 proteins can also be grouped by two different schemes. One scheme was based on a taxonomic split: class I (prokaryotic/mitochondrial) and class II (eukaryotic microsomes). The other scheme was based on the number of components in the system: class B (3-components) and class E (2-components). These classes merge to a certain degree. Most prokaryotes and mitochondria (and fungal CYP55) have 3-component systems (class I/class B) - a FAD-containing flavoprotein (NAD(P)H-dependent reductase), an iron-sulphur protein and P450. Most eukaryotic microsomes have 2-component systems (class II/class E) - NADPH:P450 reductase (FAD and FMN-containing flavoprotein) and P450. There are exceptions to this scheme, such as 1-component systems that resemble class E enzymes [, , ]. The class E enzymes can be further subdivided into five sequence clusters, groups I-V, each of which may contain more than one cytochrome P450 family (eg, CYP1 and CYP2 are both found in group I). The divergence of the cytochrome P450 superfamily into B- and E-classes, and further divergence into stable clusters within the E-class, appears to be very ancient, occurring before the appearance of eukaryotes. More information about these proteins can be found at Protein of the Month: Cytochrome P450 [].; GO: 0005506 iron ion binding, 0009055 electron carrier activity, 0016705 oxidoreductase activity, acting on paired donors, with incorporation or reduction of molecular oxygen, 0020037 heme binding, 0055114 oxidation-reduction process; PDB: 2RFC_B 2RFB_A 3EJB_H 3EJE_H 3EJD_H 1N6B_A 1NR6_A 1DT6_A 3EL3_A 3DBG_B ....

>PLN02500 cytochrome P450 90B1

Probab=86.26  E-value=2.4  Score=33.28  Aligned_cols=38  Identities=16%  Similarity=0.275  Sum_probs=23.2

Q ss_pred             CCCCCCcCCCCCHHHHHHHHHHhhcCCCCCCCCcccccchHHHHHHHhhCC
Q 034451           44 LKGNSYRLLFGDLKENSIELKEAKARPLSLDDNIAIRVNPFLHKLVNDYGM   94 (94)
Q Consensus        44 I~GP~~~~l~Gn~~e~~~~~~~a~~~~~~~shdi~prv~P~~~~w~k~YGk   94 (94)
                      .+||++.+++||..++.+.          ..|+   ...-++.+|.++||+
T Consensus        40 PPgp~~~PiiGn~~~~~~~----------~~~~---~~~~~~~~~~~~yG~   77 (490)
T PLN02500         40 PPGNMGWPFLGETIGYLKP----------YSAT---SIGEFMEQHISRYGK   77 (490)
T ss_pred             CCCCcCCCchhhHHHHHhh----------cccC---ChHHHHHHHHHHhcc
Confidence            3499999999998554320          0011   112356888888984

>PLN02655 ent-kaurene oxidase

Probab=82.08  E-value=1.1  Score=35.05  Aligned_cols=33  Identities=15%  Similarity=0.116  Sum_probs=24.4

Q ss_pred             CCCCCcCCCCCHHHHHHHHHHhhcCCCCCCCCcccccchHHHHHHHhhCC
Q 034451           45 KGNSYRLLFGDLKENSIELKEAKARPLSLDDNIAIRVNPFLHKLVNDYGM   94 (94)
Q Consensus        45 ~GP~~~~l~Gn~~e~~~~~~~a~~~~~~~shdi~prv~P~~~~w~k~YGk   94 (94)
                      +||++.+++||+.++..              +   ....++.+|.++||+
T Consensus         2 pgp~~lP~iG~l~~~~~--------------~---~~~~~~~~~~~~yG~   34 (466)
T PLN02655          2 PAVPGLPVIGNLLQLKE--------------K---KPHRTFTKWSEIYGP   34 (466)
T ss_pred             cCCCCCCccccHHHcCC--------------C---chhHHHHHHHHHhCC
Confidence            69999999999876632              0   112478999999994

---

Homologous Structure Templates

Structure Templates Detected by BLAST

Original result of BLAST against Protein Data Bank

No homologous structure with e-value below 0.005

Structure Templates Detected by RPS-BLAST

Original result of RPS-BLAST against PDB70 database

No hit with e-value below 0.005

Structure Templates Detected by HHsearch

Original result of HHsearch against PDB70 database

ID	Length	Definition	Probability
Query	94
3nxu_A	485	Cytochrome P450 3A4; alpha beta protein, cytochrom	97.43
3s79_A	503	Cytochrome P450 19A1; oxidoreductase; HET: HEM ASD	96.15
3k9v_A	482	1,25-dihydroxyvitamin D(3) 24-hydroxylase, mitocho	93.78
1r9o_A	477	Cytochrome P450 2C9; monooxygenase, drug metaboliz	93.19
1po5_A	476	Cytochrome P450 2B4; oxidoreductase, membrane prot	92.64
3tbg_A	479	Cytochrome P450 2D6; monooxygenase, thioridazine,	92.12
3e6i_A	476	CYPIIE1, P450-J, cytochrome P450 2E1; CYP2E1, mono	91.32
2fdv_A	476	Cytochrome P450 2A6; CYP2A6, monooxygenase, drug m	91.0
2cd8_A	436	Cytochrome P450 monooxygenase; oxidoreductase, PIK	90.92
2ij2_A	470	Cytochrome P450 BM3; monoxygenase, heme binding pr	89.99
2hi4_A	495	Cytochrome P450 1A2; CYP1A2, monooxygenase, drug m	89.82
3swz_A	494	Steroid 17-alpha-hydroxylase/17,20 lyase; cytochro	89.08
3ld6_A	461	Lanosterol 14-alpha demethylase; cytochrome P450,	88.39
3pm0_A	507	Cypib1, cytochrome P450 1B1; CYP1B1, monooxygenase	87.57
3czh_A	481	Cytochrome P450 2R1; vitamin D, vitamin S 25-hydro	87.46
3dan_A	473	Cytochrome P450 74A2; AOS heme cytochrome P450 str	86.28
2ve3_A	444	Putative cytochrome P450 120; oxidoreductase, mono	84.94
3i3k_A	461	Lanosterol 14-alpha demethylase; cytochrome P450,	84.24
3b98_A	475	Prostaglandin I2 synthase; prostacyclin synthase,	83.6
3gw9_A	450	Sterol 14alpha-demethylase; CYP51, cytochrome P450	82.45
3b6h_A	498	Prostacyclin synthase; enzyme-inhibitor complex, C	81.71
3dsk_A	495	Cytochrome P450 74A, chloroplast; P450 fold, fatty	80.95
2cib_A	455	Cytochrome P450 51; heme, heme lipid synthesis, me	80.12

>3nxu_A Cytochrome P450 3A4; alpha beta protein, cytochrome P450 fold, hemoprotein, monoo cytochrome P450 reductase, endoplasmic reticulum; HET: HEM RIT; 2.00A {Homo sapiens} SCOP: a.104.1.1 PDB: 1w0e_A* 1w0g_A* 2j0d_A* 2v0m_A* 1w0f_A* 1tqn_A* 3ua1_A* 3tjs_A*

Probab=97.43  E-value=7.6e-05  Score=55.79  Aligned_cols=41  Identities=15%  Similarity=0.269  Sum_probs=32.1

Q ss_pred             HHHHHHcCCCCCCCcCCCCCHHHHHHHHHHhhcCCCCCCCCcccccchHHHHHHHhhCC
Q 034451           36 EKFLRQQGLKGNSYRLLFGDLKENSIELKEAKARPLSLDDNIAIRVNPFLHKLVNDYGM   94 (94)
Q Consensus        36 ~~~lr~QGI~GP~~~~l~Gn~~e~~~~~~~a~~~~~~~shdi~prv~P~~~~w~k~YGk   94 (94)
                      ..+++++|++||++.+++||+.++.                  .....++.+|.++||+
T Consensus         8 ~~~~k~~~~PGP~~~PliGn~~~~~------------------~~~~~~~~~~~~~yG~   48 (485)
T 3nxu_A            8 HGLFKKLGIPGPTPLPFLGNILSYH------------------KGFCMFDMECHKKYGK   48 (485)
T ss_dssp             TTHHHHHTCCCCCCBTTTBTGGGGG------------------GCHHHHHHHHHHHHCS
T ss_pred             chHHhhCCCCCCCCcCeecCcHHhh------------------cChHHHHHHHHHHcCC
Confidence            3689999999999999999985442                  2223467899999995

>3s79_A Cytochrome P450 19A1; oxidoreductase; HET: HEM ASD; 2.75A {Homo sapiens} PDB: 3eqm_A* 3s7s_A* 4gl5_A* 4gl7_A*

Probab=96.15  E-value=0.0018  Score=48.43  Aligned_cols=42  Identities=17%  Similarity=0.141  Sum_probs=28.4

Q ss_pred             HHHcCCCCCCCcCCCCCHHHHHHHHHHhhcCCCCCCCCcccccchHHHHHHHhhCC
Q 034451           39 LRQQGLKGNSYRLLFGDLKENSIELKEAKARPLSLDDNIAIRVNPFLHKLVNDYGM   94 (94)
Q Consensus        39 lr~QGI~GP~~~~l~Gn~~e~~~~~~~a~~~~~~~shdi~prv~P~~~~w~k~YGk   94 (94)
                      .++++++||++-+++||+.++.+...           +   .....+.+|.++||+
T Consensus        42 ~~~~~pPGP~~~P~iGnl~~~~~~~~-----------~---~~~~~~~~~~~~yG~   83 (503)
T 3s79_A           42 EGTSSIPGPGYCMGIGPLISHGRFLW-----------M---GIGSACNYYNRVYGE   83 (503)
T ss_dssp             ---CCCCSCCCCSSSHHHHHHHHHHH-----------H---CHHHHHHHHHHHSCS
T ss_pred             hccCCCCCCCCCceeeehhccccccc-----------c---chhHHHHHHHHHhCC
Confidence            34579999999999999988754211           0   111256889999995

>3k9v_A 1,25-dihydroxyvitamin D(3) 24-hydroxylase, mitochondrial; mitochondrial cytochrome P450, monotopic membrane protein, monooxygenase; HET: HEM CPS; 2.50A {Rattus norvegicus} PDB: 3k9y_A*

Probab=93.78  E-value=0.063  Score=39.59  Aligned_cols=48  Identities=21%  Similarity=0.106  Sum_probs=31.0

Q ss_pred             HHHHhchHHHHHHHHHcCCCCCCCcCCCCCHHHHHHHHHHhhcCCCCCCCCcccccchHHHHHHHhhCC
Q 034451           26 NWVWLRPKKLEKFLRQQGLKGNSYRLLFGDLKENSIELKEAKARPLSLDDNIAIRVNPFLHKLVNDYGM   94 (94)
Q Consensus        26 ~~~w~~P~rl~~~lr~QGI~GP~~~~l~Gn~~e~~~~~~~a~~~~~~~shdi~prv~P~~~~w~k~YGk   94 (94)
                      .....+|.+        .++||++.+++||+.++..             ..-.....-++.+|.++||+
T Consensus        16 ~~~~~~pl~--------~~PGP~~~p~iG~~~~~~~-------------~~~~~~~~~~~~~l~~~YG~   63 (482)
T 3k9v_A           16 TDGETRNVT--------DLPGPTNWPLLGSLLEIFW-------------KGGLKKQHDTLAEYHKKYGQ   63 (482)
T ss_dssp             ---CCEEGG--------GCCCSCCCTTTBTHHHHHH-------------TTCGGGHHHHHHHHHHHHCS
T ss_pred             cccccCCCC--------CCCCCCCCCccccHHHHhc-------------cCCcccHHHHHHHHHHHcCC
Confidence            444556643        4999999999999988754             00011222467899999995

>1r9o_A Cytochrome P450 2C9; monooxygenase, drug metabolizing enzyme, oxidoreductas; HET: HEM FLP; 2.00A {Homo sapiens} SCOP: a.104.1.1 PDB: 1og5_A* 1og2_A* 2nnj_A* 1pq2_A* 2nni_A* 2nnh_A* 2vn0_A* 1nr6_A* 1dt6_A* 1n6b_A*

Probab=93.19  E-value=0.041  Score=41.01  Aligned_cols=42  Identities=17%  Similarity=0.153  Sum_probs=20.5

Q ss_pred             HHHHHHc--CCCCCCCcCCCCCHHHHHHHHHHhhcCCCCCCCCcccccchHHHHHHHhhCC
Q 034451           36 EKFLRQQ--GLKGNSYRLLFGDLKENSIELKEAKARPLSLDDNIAIRVNPFLHKLVNDYGM   94 (94)
Q Consensus        36 ~~~lr~Q--GI~GP~~~~l~Gn~~e~~~~~~~a~~~~~~~shdi~prv~P~~~~w~k~YGk   94 (94)
                      .+++++.  .++||++.+++||+.++.                 .....+++.+|.++||+
T Consensus         3 ~~~~~~~~~~pPgP~~~p~iG~~~~~~-----------------~~~~~~~~~~l~~~yG~   46 (477)
T 1r9o_A            3 KKTSSGRGKLPPGPTPLPVIGNILQIG-----------------IKDISKSLTNLSKVYGP   46 (477)
T ss_dssp             ------CCBCCCCSSSCC-----CCBC-----------------HHHHHHHHHHHHHHHCS
T ss_pred             cccccCCCCCCCCCCCCceeccHhhcC-----------------CCChHHHHHHHHHHhCC
Confidence            3556666  578999999999984431                 11123567888888984

>1po5_A Cytochrome P450 2B4; oxidoreductase, membrane protein, CYP 2B4, CYP LM2, cytochro monooxygenase; HET: HEM; 1.60A {Oryctolagus cuniculus} SCOP: a.104.1.1 PDB: 3mvr_A* 2bdm_A* 3g5n_A* 3g93_A* 3kw4_A* 3me6_A* 1suo_A* 3r1a_A* 3r1b_A* 2q6n_A* 3tk3_A* 3ibd_A* 3qoa_A* 3qu8_A*

Probab=92.64  E-value=0.063  Score=40.02  Aligned_cols=35  Identities=20%  Similarity=0.159  Sum_probs=25.1

Q ss_pred             CCCCCCCcCCCCCHHHHHHHHHHhhcCCCCCCCCcccccchHHHHHHHhhCC
Q 034451           43 GLKGNSYRLLFGDLKENSIELKEAKARPLSLDDNIAIRVNPFLHKLVNDYGM   94 (94)
Q Consensus        43 GI~GP~~~~l~Gn~~e~~~~~~~a~~~~~~~shdi~prv~P~~~~w~k~YGk   94 (94)
                      .++||++.+++||+.++.                 ......++.+|.++||+
T Consensus        11 ~pPgP~~~p~iG~~~~~~-----------------~~~~~~~~~~l~~~yG~   45 (476)
T 1po5_A           11 LPPGPSPLPVLGNLLQMD-----------------RKGLLRSFLRLREKYGD   45 (476)
T ss_dssp             CCCCSCCBTTTBTGGGCC-----------------TTCHHHHHHHHHHHHCS
T ss_pred             CCcCCCCCCccccHHhcc-----------------CCcHHHHHHHHHHHhCC
Confidence            578999999999985441                 11123577899999985

>3tbg_A Cytochrome P450 2D6; monooxygenase, thioridazine, oxidoreductase; HET: RTZ HEM; 2.10A {Homo sapiens} PDB: 3qm4_A* 2f9q_A*

Probab=92.12  E-value=0.089  Score=38.74  Aligned_cols=41  Identities=22%  Similarity=0.234  Sum_probs=23.7

Q ss_pred             HHHHHcC-CC-CCCCcCCCCCHHHHHHHHHHhhcCCCCCCCCcccccchHHHHHHHhhCC
Q 034451           37 KFLRQQG-LK-GNSYRLLFGDLKENSIELKEAKARPLSLDDNIAIRVNPFLHKLVNDYGM   94 (94)
Q Consensus        37 ~~lr~QG-I~-GP~~~~l~Gn~~e~~~~~~~a~~~~~~~shdi~prv~P~~~~w~k~YGk   94 (94)
                      |.-+++| ++ ||++-+++||+.++..             .+    ..-.+.+|.++||+
T Consensus         3 ~~~ss~~kLPPGP~~lP~iGn~~~~~~-------------~~----~~~~~~~~~~kYG~   45 (479)
T 3tbg_A            3 KKTSSKGKLPPGPLPLPGLGNLLHVDF-------------QN----TPYCFDQLRRRFGD   45 (479)
T ss_dssp             ------CCCCCCSCCBTTTBTGGGCCT-------------TS----HHHHHHHHHHHHCS
T ss_pred             CCCCCCCCCCCCCCCcCcccchHhhcC-------------CC----HHHHHHHHHHHhCC
Confidence            3445556 54 9999999999844311             11    11246899999995

>3e6i_A CYPIIE1, P450-J, cytochrome P450 2E1; CYP2E1, monooxygenase, acetaminophen, oxidoreductase, heme, endoplasmic reticulum, iron, membrane; HET: HEM; 2.20A {Homo sapiens} PDB: 3e4e_A* 3gph_A* 3koh_A* 3lc4_A* 3t3z_A*

Probab=91.32  E-value=0.1  Score=38.42  Aligned_cols=34  Identities=15%  Similarity=0.184  Sum_probs=24.7

Q ss_pred             CCCCCCcCCCCCHHHHHHHHHHhhcCCCCCCCCcccccchHHHHHHHhhCC
Q 034451           44 LKGNSYRLLFGDLKENSIELKEAKARPLSLDDNIAIRVNPFLHKLVNDYGM   94 (94)
Q Consensus        44 I~GP~~~~l~Gn~~e~~~~~~~a~~~~~~~shdi~prv~P~~~~w~k~YGk   94 (94)
                      .+||++.+++||+.++.                 ......++.+|.++||+
T Consensus        12 pPgP~~~PliG~~~~~~-----------------~~~~~~~~~~~~~~yG~   45 (476)
T 3e6i_A           12 PPGPFPLPIIGNLFQLE-----------------LKNIPKSFTRLAQRFGP   45 (476)
T ss_dssp             CCCCCCBTTTBTGGGCC-----------------TTCHHHHHHHHHHHHCS
T ss_pred             CcCCCCcccccChhhhc-----------------cccHhHHHHHHHHHhCC
Confidence            68999999999985431                 11222467999999995

>2fdv_A Cytochrome P450 2A6; CYP2A6, monooxygenase, drug metabolizing enzyme, coumarin 7-hydroxylase, nicotine oxidase, oxidoreductase; HET: HEM D2G; 1.65A {Homo sapiens} PDB: 1z11_A* 1z10_A* 2fdu_A* 2fdw_A* 2fdy_A* 3t3r_A* 2pg5_A* 2pg7_A* 2pg6_A* 3t3q_A* 3ebs_A* 2p85_A* 3t3s_A*

Probab=91.00  E-value=0.11  Score=38.61  Aligned_cols=41  Identities=17%  Similarity=0.202  Sum_probs=25.2

Q ss_pred             HHHHHc-C-CCCCCCcCCCCCHHHHHHHHHHhhcCCCCCCCCcccccchHHHHHHHhhCC
Q 034451           37 KFLRQQ-G-LKGNSYRLLFGDLKENSIELKEAKARPLSLDDNIAIRVNPFLHKLVNDYGM   94 (94)
Q Consensus        37 ~~lr~Q-G-I~GP~~~~l~Gn~~e~~~~~~~a~~~~~~~shdi~prv~P~~~~w~k~YGk   94 (94)
                      ++|++. + ++||++.+++||+.++..             .|    ...++.+|.++||+
T Consensus         3 ~~~~~~~~~pPgP~~~p~iG~~~~~~~-------------~~----~~~~~~~l~~~yG~   45 (476)
T 2fdv_A            3 KKTSSKGKLPPGPTPLPFIGNYLQLNT-------------EQ----MYNSLMKISERYGP   45 (476)
T ss_dssp             -----CCBCCCCCCCBTTTBTGGGCCT-------------TB----HHHHHHHHHHHHCS
T ss_pred             ccccccCCCCCCCCCCcccccHhhcCC-------------cc----hHHHHHHHHHHhCC
Confidence            455555 4 589999999999854310             11    23577899999985

>2cd8_A Cytochrome P450 monooxygenase; oxidoreductase, PIKC, macrolide monooxygenase, antibiotic biosynthesis, heme, iron, metal-binding; HET: HEM PXI; 1.7A {Streptomyces venezuelae} PDB: 2c6h_A* 2bvj_A* 2ca0_A* 2c7x_A* 2vzm_A* 2vz7_A* 2vsj_A* 2wi9_A* 2whw_A*

Probab=90.92  E-value=0.07  Score=39.94  Aligned_cols=35  Identities=26%  Similarity=0.291  Sum_probs=17.6

Q ss_pred             HcCCCCCCCcCCCCCHH-HHHHHHHHhhcCCCCCCCCcccccchHHHHHHHhhCC
Q 034451           41 QQGLKGNSYRLLFGDLK-ENSIELKEAKARPLSLDDNIAIRVNPFLHKLVNDYGM   94 (94)
Q Consensus        41 ~QGI~GP~~~~l~Gn~~-e~~~~~~~a~~~~~~~shdi~prv~P~~~~w~k~YGk   94 (94)
                      ++|++||++.+++||.. ++.+              |    -.++++++ ++||+
T Consensus        25 ~~~~PGP~~~p~lG~~~~~~~~--------------~----p~~~~~~l-~~yGp   60 (436)
T 2cd8_A           25 QQGTTASPPVLDLGALGQDFAA--------------D----PYPTYARL-RAEGP   60 (436)
T ss_dssp             ---------CCBHHHHHHHHHH--------------C----CHHHHHHH-HTTCS
T ss_pred             ccCCCCCCccccCCCCCccccc--------------C----hHHHHHHH-HHhCC
Confidence            57899999999999986 3322              2    13577888 88985

>2ij2_A Cytochrome P450 BM3; monoxygenase, heme binding protein, atomic resolution, oxidoreductase; HET: HEM; 1.20A {Bacillus megaterium} SCOP: a.104.1.1 PDB: 2hpd_A* 1fag_A* 1jpz_A* 1zo9_A* 1zo4_A* 1zoa_A* 3m4v_A* 3ekb_A* 3ben_A* 1fah_A* 2nnb_A* 3kx3_A* 3ekd_A* 3ekf_A* 1smi_A* 1smj_A* 3kx4_A* 2ij3_A* 2ij4_A* 3hf2_A* ...

Probab=89.99  E-value=0.16  Score=37.98  Aligned_cols=36  Identities=14%  Similarity=0.086  Sum_probs=25.5

Q ss_pred             cCCCCCCCcCCCCCHHHHHHHHHHhhcCCCCCCCCcccccchHHHHHHHhhCC
Q 034451           42 QGLKGNSYRLLFGDLKENSIELKEAKARPLSLDDNIAIRVNPFLHKLVNDYGM   94 (94)
Q Consensus        42 QGI~GP~~~~l~Gn~~e~~~~~~~a~~~~~~~shdi~prv~P~~~~w~k~YGk   94 (94)
                      .+++||++.+++||+.++..             .    .....+.+|.++||+
T Consensus         3 ~~~PGP~~~p~iG~l~~~~~-------------~----~~~~~~~~~~~~yG~   38 (470)
T 2ij2_A            3 KEMPQPKTFGELKNLPLLNT-------------D----KPVQALMKIADELGE   38 (470)
T ss_dssp             CCCCCCCCCGGGTTGGGGCS-------------S----CHHHHHHHHHHHHCS
T ss_pred             CCCCCCCCCCccccHHHHhc-------------c----cchHHHHHHHHHhCC
Confidence            47899999999999855411             1    112367889999995

>2hi4_A Cytochrome P450 1A2; CYP1A2, monooxygenase, drug metabolizing enzyme, alpha-naphthoflavone, benzo(H)flavone, 7,8- benzoflavone, oxidoreductase; HET: HEM BHF; 1.95A {Homo sapiens}

Probab=89.82  E-value=0.14  Score=38.32  Aligned_cols=34  Identities=18%  Similarity=0.112  Sum_probs=25.3

Q ss_pred             CCCCCCCcCCCCCHHHHHHHHHHhhcCCCCCCCCcccccchHHHHHHHhhCC
Q 034451           43 GLKGNSYRLLFGDLKENSIELKEAKARPLSLDDNIAIRVNPFLHKLVNDYGM   94 (94)
Q Consensus        43 GI~GP~~~~l~Gn~~e~~~~~~~a~~~~~~~shdi~prv~P~~~~w~k~YGk   94 (94)
                      +++||++.+++||+.++.+              |    ...++.+|.++||+
T Consensus        17 ~~PgP~~~p~~G~~~~~~~--------------~----~~~~~~~l~~~yG~   50 (495)
T 2hi4_A           17 SPPEPWGWPLLGHVLTLGK--------------N----PHLALSRMSQRYGD   50 (495)
T ss_dssp             CCCCCCCBTTTBTHHHHTT--------------C----HHHHHHHHHHHHCS
T ss_pred             CCCCCCCCcceeeHHhcCc--------------c----HHHHHHHHHHHhCC
Confidence            4599999999999976622              1    12467899999995

>3swz_A Steroid 17-alpha-hydroxylase/17,20 lyase; cytochrome P450, CYP17A1, P450C17, P450 17A1, monooxyg 17A-hydroxylase, heme protein; HET: HEM TOK; 2.40A {Homo sapiens} PDB: 3ruk_A*

Probab=89.08  E-value=0.074  Score=39.67  Aligned_cols=34  Identities=24%  Similarity=0.139  Sum_probs=21.1

Q ss_pred             CCCCCcCCCCCHHHHHHHHHHhhcCCCCCCCCcccccchHHHHHHHhhCC
Q 034451           45 KGNSYRLLFGDLKENSIELKEAKARPLSLDDNIAIRVNPFLHKLVNDYGM   94 (94)
Q Consensus        45 ~GP~~~~l~Gn~~e~~~~~~~a~~~~~~~shdi~prv~P~~~~w~k~YGk   94 (94)
                      +||++-+++||+.++.+             +   ....-.+.+|.++||+
T Consensus        11 PgP~~lPliGnl~~l~~-------------~---~~~~~~~~~~~~kYG~   44 (494)
T 3swz_A           11 KSLLSLPLVGSLPFLPR-------------H---GHMHNNFFKLQKKYGP   44 (494)
T ss_dssp             -CCBCCCEEEEESSCTT-------------S---SCHHHHHHHTHHHHCS
T ss_pred             CCCCCCCeEcchHHhCC-------------C---CchhHHHHHHHHHcCC
Confidence            56666899999833211             0   1122367899999995

>3ld6_A Lanosterol 14-alpha demethylase; cytochrome P450, ketoconazole, S genomics, structural genomics consortium, SGC; HET: HEM KKK BCD; 2.80A {Homo sapiens} PDB: 3juv_A* 3jus_A*

Probab=88.39  E-value=0.37  Score=35.63  Aligned_cols=34  Identities=24%  Similarity=0.352  Sum_probs=23.3

Q ss_pred             CCCCCCC----cCCCCCHHHHHHHHHHhhcCCCCCCCCcccccchHHHHHHHhhCC
Q 034451           43 GLKGNSY----RLLFGDLKENSIELKEAKARPLSLDDNIAIRVNPFLHKLVNDYGM   94 (94)
Q Consensus        43 GI~GP~~----~~l~Gn~~e~~~~~~~a~~~~~~~shdi~prv~P~~~~w~k~YGk   94 (94)
                      |.++||.    -+++||+.++.+              |  |  ..++.+|.++||+
T Consensus         9 g~~~PP~~~~~lP~iG~~~~~~~--------------~--~--~~~~~~~~~kYG~   46 (461)
T 3ld6_A            9 GVKSPPYIFSPIPFLGHAIAFGK--------------S--P--IEFLENAYEKYGP   46 (461)
T ss_dssp             -CCCCCBCCCSSTTTBTHHHHHH--------------C--H--HHHHHHHHHHHCS
T ss_pred             CCCCCCCCCCCcCeeeeHHHhhh--------------C--H--HHHHHHHHHHhCC
Confidence            6666653    357899988765              2  1  2478999999995

>3pm0_A Cypib1, cytochrome P450 1B1; CYP1B1, monooxygenase, alpha-naphthoflavone, 17BETA-estradiol, oxidoreductase; HET: HEM BHF; 2.70A {Homo sapiens}

Probab=87.57  E-value=0.18  Score=37.43  Aligned_cols=34  Identities=18%  Similarity=0.165  Sum_probs=10.7

Q ss_pred             CCCCCCCcCCCCCHHHHHHHHHHhhcCCCCCCCCcccccchHHHHHHHhhCC
Q 034451           43 GLKGNSYRLLFGDLKENSIELKEAKARPLSLDDNIAIRVNPFLHKLVNDYGM   94 (94)
Q Consensus        43 GI~GP~~~~l~Gn~~e~~~~~~~a~~~~~~~shdi~prv~P~~~~w~k~YGk   94 (94)
                      -.+||++.+++||+.++.+              |    ...++.+|.++||+
T Consensus        10 ~pPgP~~~P~iG~~~~~~~--------------~----~~~~~~~~~~~yG~   43 (507)
T 3pm0_A           10 KPPGPFAWPLIGNAAAVGQ--------------A----AHLSFARLARRYGD   43 (507)
T ss_dssp             ------------------C--------------C----HHHHHHHHHHHHCS
T ss_pred             CCcCCCCCCeeCchhhcCc--------------c----HHHHHHHHHHHhCC
Confidence            4789999999999866521              1    12367999999995

>3czh_A Cytochrome P450 2R1; vitamin D, vitamin S 25-hydroxylase, drug metabolism, structural genomics, structural genomics consortium, SGC; HET: BCD HEM D2V; 2.30A {Homo sapiens} SCOP: a.104.1.1 PDB: 2ojd_A* 3c6g_A* 3dl9_A*

Probab=87.46  E-value=0.38  Score=35.95  Aligned_cols=36  Identities=14%  Similarity=0.170  Sum_probs=26.2

Q ss_pred             CCCCCCcCCCCCHHHHHHHHHHhhcCCCCCCCCcccccchHHHHHHHhhCC
Q 034451           44 LKGNSYRLLFGDLKENSIELKEAKARPLSLDDNIAIRVNPFLHKLVNDYGM   94 (94)
Q Consensus        44 I~GP~~~~l~Gn~~e~~~~~~~a~~~~~~~shdi~prv~P~~~~w~k~YGk   94 (94)
                      ++||++.+++||+.++..             .+  ....+++.+|.++||+
T Consensus        14 pPgP~~~p~~G~~~~~~~-------------~~--~~~~~~~~~l~~~yG~   49 (481)
T 3czh_A           14 PPGPPGLPFIGNIYSLAA-------------SS--ELPHVYMRKQSQVYGE   49 (481)
T ss_dssp             CCCCCCBTTTBHHHHHHH-------------CS--SCHHHHHHHHHHHHCS
T ss_pred             CCCCCCCcccccHhhcCc-------------cc--CcHHHHHHHHHHHhCC
Confidence            589999999999987744             11  0123578899999995

>3dan_A Cytochrome P450 74A2; AOS heme cytochrome P450 structure, fatty acid biosynthesis, heme, iron, lipid synthesis, lyase, metal-binding; HET: HEM; 1.80A {Parthenium argentatum} PDB: 3dam_A* 3dbm_A*

Probab=86.28  E-value=0.22  Score=36.67  Aligned_cols=39  Identities=10%  Similarity=0.136  Sum_probs=28.1

Q ss_pred             CCCCCCCcCCCCCHHHHHHHHHHhhcCCCCCCCCccc-ccchHHHHHHHhhCC
Q 034451           43 GLKGNSYRLLFGDLKENSIELKEAKARPLSLDDNIAI-RVNPFLHKLVNDYGM   94 (94)
Q Consensus        43 GI~GP~~~~l~Gn~~e~~~~~~~a~~~~~~~shdi~p-rv~P~~~~w~k~YGk   94 (94)
                      .|+||.+-+++||+.++.+-             .... ...+++.+|.++||+
T Consensus        10 ~iPGp~g~P~iG~~~~~~~~-------------~~~~g~~~~~~~~~~~kyG~   49 (473)
T 3dan_A           10 EIPGSYGIPFFQPIKDRLEY-------------FYGTGGRDEYFRSRMQKYQS   49 (473)
T ss_dssp             CCCCCCCSTTHHHHHHHHHH-------------HHSTTHHHHHHHHHHHHHTC
T ss_pred             CCCCCCCCcchhhHHHHHHH-------------HHhhcCchHHHHhHHHHhCC
Confidence            68999988999999877541             1111 223688999999995

>2ve3_A Putative cytochrome P450 120; oxidoreductase, monooxygenase, metal-binding, heme, iron; HET: HEM REA; 2.10A {Synechocystis SP} PDB: 2ve4_A*

Probab=84.94  E-value=0.92  Score=33.44  Aligned_cols=33  Identities=15%  Similarity=0.015  Sum_probs=25.7

Q ss_pred             CCCCCCCcCCCCCHHHHHHHHHHhhcCCCCCCCCcccccchHHHHHHHhhCC
Q 034451           43 GLKGNSYRLLFGDLKENSIELKEAKARPLSLDDNIAIRVNPFLHKLVNDYGM   94 (94)
Q Consensus        43 GI~GP~~~~l~Gn~~e~~~~~~~a~~~~~~~shdi~prv~P~~~~w~k~YGk   94 (94)
                      .++||++.+++||..++.+              |  |.   ++.+|.++||+
T Consensus        13 ~pPgp~~~P~iG~~~~~~~--------------~--~~---~~~~~~~~yg~   45 (444)
T 2ve3_A           13 IPPGDFGLPWLGETLNFLN--------------D--GD---FGKKRQQQFGP   45 (444)
T ss_dssp             CCCCCCCBTTTBTHHHHHH--------------C--TT---HHHHHHHHHCS
T ss_pred             CCCCCCCCCccccHHHHhc--------------C--cH---HHHHHHHHcCC
Confidence            5799999999999987754              2  22   56789999985

>3b98_A Prostaglandin I2 synthase; prostacyclin synthase, cytochrome P450 8A1, CYP8A1, isomerase; HET: HEM; 2.08A {Danio rerio} PDB: 3b99_A*

Probab=83.60  E-value=0.3  Score=36.02  Aligned_cols=34  Identities=15%  Similarity=0.072  Sum_probs=25.8

Q ss_pred             CCCCCCC-cCCCCCHHHHHHHHHHhhcCCCCCCCCcccccchHHHHHHHhhCC
Q 034451           43 GLKGNSY-RLLFGDLKENSIELKEAKARPLSLDDNIAIRVNPFLHKLVNDYGM   94 (94)
Q Consensus        43 GI~GP~~-~~l~Gn~~e~~~~~~~a~~~~~~~shdi~prv~P~~~~w~k~YGk   94 (94)
                      -.+||++ -+++||+.++.+              |  |  .+++.+|.++||+
T Consensus        19 ~pPgp~~~~P~iG~~~~~~~--------------~--~--~~~~~~l~~~yG~   53 (475)
T 3b98_A           19 EPPLDKGMIPWLGHALEFGK--------------D--A--AKFLTRMKEKHGD   53 (475)
T ss_dssp             CCCEECCSSTTTBTHHHHHH--------------C--H--HHHHHHHHHHHCS
T ss_pred             CCCCCCCCcchHHhHHHHhh--------------C--H--HHHHHHHHHHhCC
Confidence            4689988 689999987743              2  1  3578999999995

>3gw9_A Sterol 14alpha-demethylase; CYP51, cytochrome P450, heme, oxidoreductase, monooxygenase, sterol biosynthesis, lipids, endoplasmic reticulum; HET: HEM VNI; 1.87A {Trypanosoma brucei} PDB: 3tik_A* 3g1q_A* 3p99_A* 2wv2_A* 2x2n_A* 3khm_A* 3k1o_A* 3ksw_A* 2wx2_A* 2wuz_A* 3l4d_A*

Probab=82.45  E-value=0.93  Score=32.89  Aligned_cols=32  Identities=6%  Similarity=-0.002  Sum_probs=23.2

Q ss_pred             CCCCCC-cCCCCCHHHHHHHHHHhhcCCCCCCCCcccccchHHHHHHHhhC
Q 034451           44 LKGNSY-RLLFGDLKENSIELKEAKARPLSLDDNIAIRVNPFLHKLVNDYG   93 (94)
Q Consensus        44 I~GP~~-~~l~Gn~~e~~~~~~~a~~~~~~~shdi~prv~P~~~~w~k~YG   93 (94)
                      -+|||+ -+++||+.++.+              |  +  ..++.+|.++||
T Consensus         5 PPg~p~~~P~iG~~~~~~~--------------~--~--~~~~~~~~~~yG   37 (450)
T 3gw9_A            5 PPVYPVTVPILGHIIQFGK--------------S--P--LGFMQECKRQLK   37 (450)
T ss_dssp             CCBCCCCSTTTBTHHHHHH--------------C--H--HHHHHHHHHHHT
T ss_pred             CCCCCCCcchhccHHHHcc--------------C--H--HHHHHHHHHHhC
Confidence            468776 678999988754              1  1  136799999999

>3b6h_A Prostacyclin synthase; enzyme-inhibitor complex, CYP8A1, cytochrome P450, endoplasmic reticulum, fatty acid biosynthesis, heme, iron, isomerase; HET: BOG MXD HEM; 1.62A {Homo sapiens} PDB: 2iag_A*

Probab=81.71  E-value=0.37  Score=35.89  Aligned_cols=34  Identities=12%  Similarity=0.020  Sum_probs=25.3

Q ss_pred             CCCCCCC-cCCCCCHHHHHHHHHHhhcCCCCCCCCcccccchHHHHHHHhhCC
Q 034451           43 GLKGNSY-RLLFGDLKENSIELKEAKARPLSLDDNIAIRVNPFLHKLVNDYGM   94 (94)
Q Consensus        43 GI~GP~~-~~l~Gn~~e~~~~~~~a~~~~~~~shdi~prv~P~~~~w~k~YGk   94 (94)
                      -.+||++ -+++||+.++.+              |  +  ..++.+|.++||+
T Consensus        19 ~PPgp~~~~P~iG~~~~~~~--------------~--~--~~~~~~l~~kYG~   53 (498)
T 3b6h_A           19 EPPLDLGSIPWLGYALDFGK--------------D--A--ASFLTRMKEKHGD   53 (498)
T ss_dssp             CCCEECCSSTTTBTHHHHHH--------------C--H--HHHHHHHHHHHCS
T ss_pred             CCCCCCCCCcchhhHHHhcc--------------C--H--HHHHHHHHHHcCC
Confidence            3689988 588999977743              1  1  2478999999995

>3dsk_A Cytochrome P450 74A, chloroplast; P450 fold, fatty acid biosynthesis, heme, iron, synthesis, lyase, metal-binding, oxylipin biosynthesis; HET: HEM T25; 1.55A {Arabidopsis thaliana} PDB: 2rcm_A* 3dsj_A* 3dsi_A* 2rcl_A* 2rch_A* 3cli_A*

Probab=80.95  E-value=0.57  Score=34.79  Aligned_cols=39  Identities=15%  Similarity=0.224  Sum_probs=27.9

Q ss_pred             CCCCCCCcCCCCCHHHHHHHHHHhhcCCCCCCCCcccccchHHHHHHHhhCC
Q 034451           43 GLKGNSYRLLFGDLKENSIELKEAKARPLSLDDNIAIRVNPFLHKLVNDYGM   94 (94)
Q Consensus        43 GI~GP~~~~l~Gn~~e~~~~~~~a~~~~~~~shdi~prv~P~~~~w~k~YGk   94 (94)
                      -++||++-+++||+.++.+-..             ...-..++.++.++||+
T Consensus        29 ~pPGp~g~P~iG~~~~~~~~~~-------------~~~~~~f~~~~~~kyG~   67 (495)
T 3dsk_A           29 NIPGNYGLPIVGPIKDRWDYFY-------------DQGAEEFFKSRIRKYNS   67 (495)
T ss_dssp             CCCCCCCSTTHHHHHHHHHHHT-------------TSCHHHHHHHHHHHHTC
T ss_pred             CCCCCCCCCccchHHHHHHHHH-------------hcCcHHHHHHHHHHhCC
Confidence            5889988899999988754110             11223588999999995

>2cib_A Cytochrome P450 51; heme, heme lipid synthesis, metal-binding, monooxygenase, NADP, oxidoreductase, protein-inhibitor complex; HET: HEM CM6; 1.50A {Mycobacterium tuberculosis} SCOP: a.104.1.1 PDB: 2bz9_A* 1x8v_A* 2ci0_A* 2vku_A* 2w09_A* 2w0b_A* 2w0a_A* 1h5z_A* 1ea1_A* 1e9x_A* 1u13_A*

Probab=80.12  E-value=0.75  Score=34.05  Aligned_cols=35  Identities=14%  Similarity=0.159  Sum_probs=24.3

Q ss_pred             cCCCC-CCCcCCCCCHHHHHHHHHHhhcCCCCCCCCcccccchHHHHHHHhhCC
Q 034451           42 QGLKG-NSYRLLFGDLKENSIELKEAKARPLSLDDNIAIRVNPFLHKLVNDYGM   94 (94)
Q Consensus        42 QGI~G-P~~~~l~Gn~~e~~~~~~~a~~~~~~~shdi~prv~P~~~~w~k~YGk   94 (94)
                      ...+| |++.+++||..++.+              |  |  ..++.+|.++||+
T Consensus         4 ~~PPg~p~~~P~iG~~~~~~~--------------~--~--~~~~~~l~~~yG~   39 (455)
T 2cib_A            4 VALPRVSGGHDEHGHLEEFRT--------------D--P--IGLMQRVRDELGD   39 (455)
T ss_dssp             -CCCBCSCCCBTTBTHHHHTT--------------C--H--HHHHHHHHHHHCS
T ss_pred             CCCCCCCCCCCCccCHHHHhh--------------C--h--HHHHHHHHHHcCC
Confidence            34577 667889999977632              2  1  2577899999984

---

Homologous Structure Domains

Structure Domains Detected by RPS-BLAST

Original result of RPS-BLAST against SCOP70(version1.75) database

No hit with e-value below 0.005

Homologous Domains Detected by HHsearch

Original result of HHsearch against SCOP70(version1.75) database

ID	Length	Definition	Probability
Query	94
d1tqna_	472	Mammalian cytochrome P450 3a4 {Human (Homo sapiens	97.16
d3czha1	463	Vitamin D 25-hydroxylase Cyp2R1 {Human (Homo sapie	92.09
d2ij2a1	453	Cytochrome P450 bm-3 {Bacillus megaterium [TaxId:	90.03
d1r9oa_	467	Mammalian cytochrome p450 2c9 {Human (Homo sapiens	85.69

>d1tqna_ a.104.1.1 (A:) Mammalian cytochrome P450 3a4 {Human (Homo sapiens) [TaxId: 9606]}

class: All alpha proteins
fold: Cytochrome P450
superfamily: Cytochrome P450
family: Cytochrome P450
domain: Mammalian cytochrome P450 3a4
species: Human (Homo sapiens) [TaxId: 9606]

Probab=97.16  E-value=0.00012  Score=51.98  Aligned_cols=40  Identities=15%  Similarity=0.272  Sum_probs=31.5

Q ss_pred             HHHHHcCCCCCCCcCCCCCHHHHHHHHHHhhcCCCCCCCCcccccchHHHHHHHhhCC
Q 034451           37 KFLRQQGLKGNSYRLLFGDLKENSIELKEAKARPLSLDDNIAIRVNPFLHKLVNDYGM   94 (94)
Q Consensus        37 ~~lr~QGI~GP~~~~l~Gn~~e~~~~~~~a~~~~~~~shdi~prv~P~~~~w~k~YGk   94 (94)
                      ..|+|.||+||++-+++||+-++.+                  ...-++.+|.++||+
T Consensus         4 ~~~~~~~iPGP~~~P~iG~~~~~~~------------------~~~~~~~~~~~kyG~   43 (472)
T d1tqna_           4 GLFKKLGIPGPTPLPFLGNILSYHK------------------GFCMFDMECHKKYGK   43 (472)
T ss_dssp             THHHHTTCCCCCCBTTTBTGGGGGG------------------CHHHHHHHHHHHHCS
T ss_pred             cchhhcCCCCCCCcCceeEHHHhhC------------------CHHHHHHHHHHHhCC
Confidence            4689999999999999999855532                  223477899999995

>d3czha1 a.104.1.1 (A:40-502) Vitamin D 25-hydroxylase Cyp2R1 {Human (Homo sapiens) [TaxId: 9606]}

class: All alpha proteins
fold: Cytochrome P450
superfamily: Cytochrome P450
family: Cytochrome P450
domain: Vitamin D 25-hydroxylase Cyp2R1
species: Human (Homo sapiens) [TaxId: 9606]

Probab=92.09  E-value=0.045  Score=37.59  Aligned_cols=35  Identities=14%  Similarity=0.192  Sum_probs=25.8

Q ss_pred             CCCCCcCCCCCHHHHHHHHHHhhcCCCCCCCCcccccchHHHHHHHhhCC
Q 034451           45 KGNSYRLLFGDLKENSIELKEAKARPLSLDDNIAIRVNPFLHKLVNDYGM   94 (94)
Q Consensus        45 ~GP~~~~l~Gn~~e~~~~~~~a~~~~~~~shdi~prv~P~~~~w~k~YGk   94 (94)
                      +||+.-+++||+.++..             +...+.  -+++++.++||+
T Consensus         2 PGP~~~P~iG~~~~~~~-------------~~~~~~--~~~~~~~~kyG~   36 (463)
T d3czha1           2 PGPPGLPFIGNIYSLAA-------------SSELPH--VYMRKQSQVYGE   36 (463)
T ss_dssp             CCCCCBTTTBHHHHHHH-------------CSSCHH--HHHHHHHHHHCS
T ss_pred             cCCCCcCccccHHHhhC-------------CCCcHH--HHHHHHHHHhCC
Confidence            69998899999988865             111122  278999999995

>d2ij2a1 a.104.1.1 (A:3-455) Cytochrome P450 bm-3 {Bacillus megaterium [TaxId: 1404]}

class: All alpha proteins
fold: Cytochrome P450
superfamily: Cytochrome P450
family: Cytochrome P450
domain: Cytochrome P450 bm-3
species: Bacillus megaterium [TaxId: 1404]

Probab=90.03  E-value=0.056  Score=37.93  Aligned_cols=36  Identities=14%  Similarity=0.068  Sum_probs=26.2

Q ss_pred             cCCCCCCCcCCCCCHHHHHHHHHHhhcCCCCCCCCcccccchHHHHHHHhhCC
Q 034451           42 QGLKGNSYRLLFGDLKENSIELKEAKARPLSLDDNIAIRVNPFLHKLVNDYGM   94 (94)
Q Consensus        42 QGI~GP~~~~l~Gn~~e~~~~~~~a~~~~~~~shdi~prv~P~~~~w~k~YGk   94 (94)
                      +.|+||++-+++||+.++.+             .+    ..-++.+|.++||+
T Consensus         1 r~iPGP~~~p~lG~l~~l~~-------------~~----~~~~~~~~~~kyG~   36 (453)
T d2ij2a1           1 KEMPQPKTFGELKNLPLLNT-------------DK----PVQALMKIADELGE   36 (453)
T ss_dssp             CCCCCCCCCGGGTTGGGGCS-------------SC----HHHHHHHHHHHHCS
T ss_pred             CCCccCCCcchhhCHHHhCC-------------CC----HHHHHHHHHHHhCC
Confidence            46899999999999855421             12    22467999999995

>d1r9oa_ a.104.1.1 (A:) Mammalian cytochrome p450 2c9 {Human (Homo sapiens) [TaxId: 9606]}

class: All alpha proteins
fold: Cytochrome P450
superfamily: Cytochrome P450
family: Cytochrome P450
domain: Mammalian cytochrome p450 2c9
species: Human (Homo sapiens) [TaxId: 9606]

Probab=85.69  E-value=0.17  Score=34.86  Aligned_cols=33  Identities=21%  Similarity=0.358  Sum_probs=18.0

Q ss_pred             CCCCCcCCCCCHHHHHHHHHHhhcCCCCCCCCcccccchHHHHHHHhhCC
Q 034451           45 KGNSYRLLFGDLKENSIELKEAKARPLSLDDNIAIRVNPFLHKLVNDYGM   94 (94)
Q Consensus        45 ~GP~~~~l~Gn~~e~~~~~~~a~~~~~~~shdi~prv~P~~~~w~k~YGk   94 (94)
                      +||++-+++||+.++..             +|    ...++.+|.++||+
T Consensus         6 PGP~~~P~lG~~~~l~~-------------~~----~~~~~~~~~~kyG~   38 (467)
T d1r9oa_           6 PGPTPLPVIGNILQIGI-------------KD----ISKSLTNLSKVYGP   38 (467)
T ss_dssp             CCSSSCC-----CCBCH-------------HH----HHHHHHHHHHHHCS
T ss_pred             cCCCCCCccccHHHhCC-------------cC----HHHHHHHHHHHhCC
Confidence            49999999999854421             12    22467889999985