annotate variant_effect_predictor/Bio/EnsEMBL/DBSQL/PredictionExonAdaptor.pm @ 0:1f6dce3d34e0

Uploaded
author mahtabm
date Thu, 11 Apr 2013 02:01:53 -0400
parents
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
1 =head1 LICENSE
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
2
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
3 Copyright (c) 1999-2012 The European Bioinformatics Institute and
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
4 Genome Research Limited. All rights reserved.
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
5
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
6 This software is distributed under a modified Apache license.
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
7 For license details, please see
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
8
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
9 http://www.ensembl.org/info/about/code_licence.html
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
10
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
11 =head1 CONTACT
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
12
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
13 Please email comments or questions to the public Ensembl
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
14 developers list at <dev@ensembl.org>.
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
15
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
16 Questions may also be sent to the Ensembl help desk at
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
17 <helpdesk@ensembl.org>.
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
18
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
19 =cut
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
20
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
21 =head1 NAME
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
22
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
23 Bio::EnsEMBL::DBSQL::PredictionExonAdaptor - Performs database interaction for
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
24 PredictionExons.
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
25
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
26 =head1 SYNOPSIS
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
27
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
28 $pea = $database_adaptor->get_PredictionExonAdaptor();
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
29 $pexon = $pea->fetch_by_dbID();
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
30
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
31 my $slice =
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
32 $database_adaptor->get_SliceAdaptor->fetch_by_region( 'X', 1, 1e6 );
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
33
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
34 my @pexons = @{ $pea->fetch_all_by_Slice($slice) };
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
35
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
36 =head1 METHODS
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
37
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
38 =cut
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
39
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
40 package Bio::EnsEMBL::DBSQL::PredictionExonAdaptor;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
41
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
42 use vars qw( @ISA );
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
43 use strict;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
44
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
45
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
46 use Bio::EnsEMBL::DBSQL::BaseFeatureAdaptor;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
47 use Bio::EnsEMBL::PredictionExon;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
48 use Bio::EnsEMBL::Utils::Exception qw( warning throw deprecate );
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
49
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
50
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
51 @ISA = qw( Bio::EnsEMBL::DBSQL::BaseFeatureAdaptor );
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
52
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
53
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
54 #_tables
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
55 #
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
56 # Arg [1] : none
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
57 # Example : none
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
58 # Description: PROTECTED implementation of superclass abstract method
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
59 # returns the names, aliases of the tables to use for queries
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
60 # Returntype : list of listrefs of strings
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
61 # Exceptions : none
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
62 # Caller : internal
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
63 #
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
64
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
65 sub _tables {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
66 return ([ 'prediction_exon', 'pe' ] );
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
67 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
68
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
69
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
70
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
71 #_columns
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
72 #
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
73 # Arg [1] : none
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
74 # Example : none
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
75 # Description: PROTECTED implementation of superclass abstract method
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
76 # returns a list of columns to use for queries
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
77 # Returntype : list of strings
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
78 # Exceptions : none
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
79 # Caller : internal
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
80
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
81 sub _columns {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
82 my $self = shift;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
83
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
84 return qw( pe.prediction_exon_id
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
85 pe.seq_region_id
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
86 pe.seq_region_start
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
87 pe.seq_region_end
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
88 pe.seq_region_strand
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
89 pe.start_phase
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
90 pe.score
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
91 pe.p_value );
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
92 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
93
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
94
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
95 # _final_clause
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
96 #
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
97 # Arg [1] : none
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
98 # Example : none
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
99 # Description: PROTECTED implementation of superclass abstract method
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
100 # returns a default end for the SQL-query (ORDER BY)
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
101 # Returntype : string
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
102 # Exceptions : none
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
103 # Caller : internal
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
104
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
105 sub _final_clause {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
106 return "ORDER BY pe.prediction_transcript_id, pe.exon_rank";
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
107 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
108
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
109
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
110 =head2 fetch_all_by_PredictionTranscript
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
111
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
112 Arg [1] : Bio::EnsEMBL::PredcitionTranscript $transcript
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
113 Example : none
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
114 Description: Retrieves all Exons for the Transcript in 5-3 order
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
115 Returntype : listref Bio::EnsEMBL::Exon on Transcript slice
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
116 Exceptions : throws if transcript does not have a slice
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
117 Caller : Transcript->get_all_Exons()
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
118 Status : Stable
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
119
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
120 =cut
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
121
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
122 sub fetch_all_by_PredictionTranscript {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
123 my ( $self, $transcript ) = @_;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
124 my $constraint = "pe.prediction_transcript_id = ".$transcript->dbID();
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
125
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
126 # use 'keep_all' option to keep exons that are off end of slice
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
127
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
128 my $tslice = $transcript->slice();
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
129 my $slice;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
130
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
131 if(!$tslice) {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
132 throw("Transcript must have attached slice to retrieve exons.");
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
133 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
134
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
135 # use a small slice the same size as the prediction transcript
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
136 $slice = $self->db->get_SliceAdaptor->fetch_by_Feature($transcript);
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
137
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
138 my $exons = $self->fetch_all_by_Slice_constraint($slice, $constraint);
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
139
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
140 # remap exon coordinates if necessary
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
141 if($slice->name() ne $tslice->name()) {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
142 my @out;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
143 foreach my $ex (@$exons) {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
144 push @out, $ex->transfer($tslice);
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
145 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
146 $exons = \@out;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
147 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
148
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
149 return $exons;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
150 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
151
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
152
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
153
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
154 =head2 store
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
155
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
156 Arg [1] : Bio::EnsEMBL::PredictionExon $exon
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
157 The exon to store in this database
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
158 Arg [2] : int $prediction_transcript_id
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
159 The internal identifier of the prediction exon that that this
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
160 exon is associated with.
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
161 Arg [3] : int $rank
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
162 The rank of the exon in the transcript (starting at 1)
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
163 Example : $pexon_adaptor->store($pexon, 1211, 2);
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
164 Description: Stores a PredictionExon in the database
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
165 Returntype : none
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
166 Exceptions : thrown if exon does not have a slice attached
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
167 or if $exon->start, $exon->end, $exon->strand, or $exon->phase
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
168 are not defined or if $exon is not a Bio::EnsEMBL::PredictionExon
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
169 Caller : general
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
170 Status : Stable
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
171
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
172 =cut
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
173
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
174 sub store {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
175 my ( $self, $pexon, $pt_id, $rank ) = @_;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
176
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
177 if(!ref($pexon) || !$pexon->isa('Bio::EnsEMBL::PredictionExon') ) {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
178 throw("Expected PredictionExon argument");
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
179 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
180
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
181 throw("Expected PredictionTranscript id argument.") if(!$pt_id);
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
182 throw("Expected rank argument.") if(!$rank);
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
183
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
184 my $db = $self->db();
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
185
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
186 if($pexon->is_stored($db)) {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
187 warning('PredictionExon is already stored in this DB.');
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
188 return $pexon->dbID();
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
189 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
190
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
191 if( ! $pexon->start || ! $pexon->end ||
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
192 ! $pexon->strand || ! defined $pexon->phase ) {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
193 throw("PredictionExon does not have all attributes to store.\n" .
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
194 "start, end, strand and phase attributes must be set.");
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
195 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
196
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
197 #maintain reference to original passed-in prediction exon
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
198 my $original = $pexon;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
199 my $seq_region_id;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
200 ($pexon, $seq_region_id) = $self->_pre_store($pexon);
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
201
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
202 my $sth = $db->dbc->prepare
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
203 ("INSERT into prediction_exon (prediction_transcript_id, exon_rank, " .
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
204 "seq_region_id, seq_region_start, seq_region_end, " .
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
205 "seq_region_strand, start_phase, score, p_value) " .
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
206 "VALUES ( ?, ?, ?, ?, ?, ?, ?, ?, ? )");
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
207
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
208 $sth->bind_param(1,$pt_id,SQL_INTEGER);
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
209 $sth->bind_param(2,$rank,SQL_SMALLINT);
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
210 $sth->bind_param(3,$seq_region_id,SQL_INTEGER);
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
211 $sth->bind_param(4,$pexon->start,SQL_INTEGER);
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
212 $sth->bind_param(5,$pexon->end,SQL_INTEGER);
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
213 $sth->bind_param(6,$pexon->strand,SQL_TINYINT);
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
214 $sth->bind_param(7,$pexon->phase,SQL_TINYINT);
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
215 $sth->bind_param(8,$pexon->score,SQL_DOUBLE);
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
216 $sth->bind_param(9,$pexon->p_value,SQL_DOUBLE);
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
217
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
218 $sth->execute();
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
219
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
220 my $dbID = $sth->{'mysql_insertid'};
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
221
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
222 #set the adaptor and dbID of the object they passed in
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
223 $original->dbID($dbID);
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
224 $original->adaptor($self);
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
225
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
226 return $dbID;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
227 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
228
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
229
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
230
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
231 =head2 remove
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
232
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
233 Arg [1] : Bio::EnsEMBL::PredictionExon $exon
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
234 the exon to remove from the database
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
235 Example : $exon_adaptor->remove($exon);
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
236 Description: Removes an exon from the database
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
237 Returntype : none
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
238 Exceptions : none
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
239 Caller : general
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
240 Status : Stable
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
241
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
242 =cut
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
243
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
244 sub remove {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
245 my $self = shift;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
246 my $pexon = shift;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
247
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
248 my $db = $self->db();
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
249
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
250 if(!$pexon->is_stored($db)) {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
251 warning('PredictionExon is not in this DB - not removing');
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
252 return undef;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
253 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
254
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
255 my $sth = $self->prepare(
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
256 "DELETE FROM prediction_exon WHERE prediction_exon_id = ?");
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
257 $sth->bind_param( 1, $pexon->dbID, SQL_INTEGER );
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
258 $sth->execute();
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
259
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
260 $pexon->dbID(undef);
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
261 $pexon->adaptor(undef);
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
262 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
263
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
264
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
265
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
266 =head2 list_dbIDs
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
267
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
268 Arg [1] : none
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
269 Example : @exon_ids = @{$exon_adaptor->list_dbIDs()};
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
270 Description: Gets an array of internal ids for all exons in the current db
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
271 Arg[1] : <optional> int. not 0 for the ids to be sorted by the seq_region.
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
272 Returntype : list of ints
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
273 Exceptions : none
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
274 Caller : ?
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
275 Status : Stable
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
276
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
277 =cut
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
278
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
279 sub list_dbIDs {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
280 my ($self,$ordered) = @_;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
281
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
282 return $self->_list_dbIDs("prediction_exon",undef, $ordered);
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
283 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
284
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
285
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
286
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
287 #_objs_from_sth
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
288
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
289 # Arg [1] : Hashreference $hashref
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
290 # Example : none
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
291 # Description: PROTECTED implementation of abstract superclass method.
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
292 # responsible for the creation of Genes
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
293 # Returntype : listref of Bio::EnsEMBL::Genes in target coordinate system
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
294 # Exceptions : none
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
295 # Caller : internal
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
296 #
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
297
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
298 sub _objs_from_sth {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
299 my ($self, $sth, $mapper, $dest_slice) = @_;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
300
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
301 #
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
302 # This code is ugly because an attempt has been made to remove as many
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
303 # function calls as possible for speed purposes. Thus many caches and
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
304 # a fair bit of gymnastics is used.
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
305 #
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
306 my $sa = $self->db()->get_SliceAdaptor();
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
307
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
308 my @exons;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
309 my %slice_hash;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
310 my %sr_name_hash;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
311 my %sr_cs_hash;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
312
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
313 my($prediction_exon_id,$seq_region_id,
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
314 $seq_region_start, $seq_region_end, $seq_region_strand,
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
315 $start_phase, $score, $p_value);
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
316
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
317 $sth->bind_columns(\$prediction_exon_id,\$seq_region_id,
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
318 \$seq_region_start, \$seq_region_end, \$seq_region_strand,
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
319 \$start_phase, \$score, \$p_value);
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
320
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
321 my $asm_cs;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
322 my $cmp_cs;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
323 my $asm_cs_vers;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
324 my $asm_cs_name;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
325 my $cmp_cs_vers;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
326 my $cmp_cs_name;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
327 if($mapper) {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
328 $asm_cs = $mapper->assembled_CoordSystem();
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
329 $cmp_cs = $mapper->component_CoordSystem();
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
330 $asm_cs_name = $asm_cs->name();
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
331 $asm_cs_vers = $asm_cs->version();
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
332 $cmp_cs_name = $cmp_cs->name();
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
333 $cmp_cs_vers = $cmp_cs->version();
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
334 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
335
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
336 my $dest_slice_start;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
337 my $dest_slice_end;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
338 my $dest_slice_strand;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
339 my $dest_slice_length;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
340 my $dest_slice_cs;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
341 my $asma;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
342 my $dest_slice_sr_name;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
343 my $dest_slice_sr_id;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
344
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
345 if($dest_slice) {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
346 $dest_slice_start = $dest_slice->start();
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
347 $dest_slice_end = $dest_slice->end();
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
348 $dest_slice_strand = $dest_slice->strand();
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
349 $dest_slice_length = $dest_slice->length();
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
350 $dest_slice_cs = $dest_slice->coord_system;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
351 $dest_slice_sr_name = $dest_slice->seq_region_name();
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
352 $dest_slice_sr_id = $dest_slice->get_seq_region_id();
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
353 $asma = $self->db->get_AssemblyMapperAdaptor();
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
354 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
355
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
356 FEATURE: while($sth->fetch()) {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
357 #need to get the internal_seq_region, if present
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
358 $seq_region_id = $self->get_seq_region_id_internal($seq_region_id);
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
359 my $slice = $slice_hash{"ID:".$seq_region_id};
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
360 my $dest_mapper = $mapper;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
361
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
362
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
363 if(!$slice) {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
364 $slice = $sa->fetch_by_seq_region_id($seq_region_id);
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
365 $slice_hash{"ID:".$seq_region_id} = $slice;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
366 $sr_name_hash{$seq_region_id} = $slice->seq_region_name();
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
367 $sr_cs_hash{$seq_region_id} = $slice->coord_system();
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
368 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
369
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
370 #obtain a mapper if none was defined, but a dest_seq_region was
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
371 if(!$dest_mapper && $dest_slice &&
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
372 !$dest_slice_cs->equals($slice->coord_system)) {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
373 $dest_mapper = $asma->fetch_by_CoordSystems($dest_slice_cs,
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
374 $slice->coord_system);
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
375 $asm_cs = $dest_mapper->assembled_CoordSystem();
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
376 $cmp_cs = $dest_mapper->component_CoordSystem();
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
377 $asm_cs_name = $asm_cs->name();
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
378 $asm_cs_vers = $asm_cs->version();
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
379 $cmp_cs_name = $cmp_cs->name();
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
380 $cmp_cs_vers = $cmp_cs->version();
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
381 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
382
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
383 my $sr_name = $sr_name_hash{$seq_region_id};
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
384 my $sr_cs = $sr_cs_hash{$seq_region_id};
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
385
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
386 #
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
387 # remap the feature coordinates to another coord system
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
388 # if a mapper was provided
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
389 #
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
390 if($dest_mapper) {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
391
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
392 if (defined $dest_slice && $dest_mapper->isa('Bio::EnsEMBL::ChainedAssemblyMapper') ) {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
393 ( $seq_region_id, $seq_region_start,
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
394 $seq_region_end, $seq_region_strand )
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
395 =
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
396 $dest_mapper->map( $sr_name, $seq_region_start, $seq_region_end,
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
397 $seq_region_strand, $sr_cs, 1, $dest_slice);
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
398
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
399 } else {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
400
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
401 ( $seq_region_id, $seq_region_start,
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
402 $seq_region_end, $seq_region_strand )
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
403 = $dest_mapper->fastmap( $sr_name, $seq_region_start,
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
404 $seq_region_end, $seq_region_strand,
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
405 $sr_cs );
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
406 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
407
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
408 #skip features that map to gaps or coord system boundaries
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
409 next FEATURE if(!defined($seq_region_id));
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
410
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
411 #get a slice in the coord system we just mapped to
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
412 # if($asm_cs == $sr_cs || ($cmp_cs != $sr_cs && $asm_cs->equals($sr_cs))) {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
413 $slice = $slice_hash{"ID:".$seq_region_id} ||=
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
414 $sa->fetch_by_seq_region_id($seq_region_id);
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
415 # } else {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
416 # $slice = $slice_hash{"NAME:$sr_name:$asm_cs_name:$asm_cs_vers"} ||=
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
417 # $sa->fetch_by_region($asm_cs_name, $sr_name, undef, undef, undef,
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
418 # $asm_cs_vers);
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
419 # }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
420 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
421
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
422 #
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
423 # If a destination slice was provided convert the coords
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
424 # If the dest_slice starts at 1 and is foward strand, nothing needs doing
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
425 #
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
426 if($dest_slice) {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
427 if($dest_slice_start != 1 || $dest_slice_strand != 1) {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
428 if($dest_slice_strand == 1) {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
429 $seq_region_start = $seq_region_start - $dest_slice_start + 1;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
430 $seq_region_end = $seq_region_end - $dest_slice_start + 1;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
431 } else {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
432 my $tmp_seq_region_start = $seq_region_start;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
433 $seq_region_start = $dest_slice_end - $seq_region_end + 1;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
434 $seq_region_end = $dest_slice_end - $tmp_seq_region_start + 1;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
435 $seq_region_strand *= -1;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
436 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
437 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
438
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
439 #throw away features off the end of the requested slice
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
440 if($seq_region_end < 1 || $seq_region_start > $dest_slice_length ||
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
441 ( $dest_slice_sr_id ne $seq_region_id )) {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
442 next FEATURE;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
443 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
444
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
445 $slice = $dest_slice;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
446 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
447
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
448 # Finally, create the new PredictionExon.
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
449 push( @exons,
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
450 $self->_create_feature( 'Bio::EnsEMBL::PredictionExon', {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
451 '-start' => $seq_region_start,
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
452 '-end' => $seq_region_end,
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
453 '-strand' => $seq_region_strand,
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
454 '-adaptor' => $self,
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
455 '-slice' => $slice,
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
456 '-dbID' => $prediction_exon_id,
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
457 '-phase' => $start_phase,
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
458 '-score' => $score,
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
459 '-p_value' => $p_value
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
460 } ) );
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
461
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
462 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
463
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
464 return \@exons;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
465 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
466
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
467
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
468 1;