Mercurial > repos > mahtabm > ensembl
comparison variant_effect_predictor/Bio/UpdateableSeqI.pm @ 0:1f6dce3d34e0
Uploaded
author | mahtabm |
---|---|
date | Thu, 11 Apr 2013 02:01:53 -0400 |
parents | |
children |
comparison
equal
deleted
inserted
replaced
-1:000000000000 | 0:1f6dce3d34e0 |
---|---|
1 # $Id: UpdateableSeqI.pm,v 1.4 2002/10/22 07:45:09 lapp Exp $ | |
2 # | |
3 # BioPerl module for Bio::UpdateableSeqI | |
4 # | |
5 # Cared for by David Block <dblock@gene.pbi.nrc.ca> | |
6 # | |
7 # Copyright David Block | |
8 # | |
9 # You may distribute this module under the same terms as perl itself | |
10 | |
11 # POD documentation - main docs before the code | |
12 | |
13 =head1 NAME | |
14 | |
15 Bio::UpdateableSeqI - Descendant of Bio::SeqI that allows updates | |
16 | |
17 =head1 SYNOPSIS | |
18 | |
19 See Bio::SeqI for most of the documentation. | |
20 See the documentation of the methods for further details. | |
21 | |
22 =head1 DESCRIPTION | |
23 | |
24 Bio::UpdateableSeqI is an interface for Sequence objects which are | |
25 expected to allow users to perform basic editing functions (update/delete) | |
26 on their component SeqFeatures. | |
27 | |
28 =head1 FEEDBACK | |
29 | |
30 =head2 Mailing Lists | |
31 | |
32 User feedback is an integral part of the evolution of this and other | |
33 Bioperl modules. Send your comments and suggestions preferably to | |
34 the Bioperl mailing list. Your participation is much appreciated. | |
35 | |
36 bioperl-l@bioperl.org - General discussion | |
37 http://bioperl.org/MailList.shtml - About the mailing lists | |
38 | |
39 =head2 Reporting Bugs | |
40 | |
41 Report bugs to the Bioperl bug tracking system to help us keep track | |
42 of the bugs and their resolution. Bug reports can be submitted via | |
43 email or the web: | |
44 | |
45 bioperl-bugs@bioperl.org | |
46 http://bugzilla.bioperl.org/ | |
47 | |
48 =head1 AUTHOR - David Block | |
49 | |
50 Email dblock@gene.pbi.nrc.ca | |
51 | |
52 =head1 CONTRIBUTORS | |
53 | |
54 Ewan Birney forced me to this... | |
55 | |
56 =head1 APPENDIX | |
57 | |
58 The rest of the documentation details each of the object methods. | |
59 Internal methods are usually preceded with a _ | |
60 | |
61 =cut | |
62 | |
63 | |
64 # Let the code begin... | |
65 | |
66 | |
67 package Bio::UpdateableSeqI; | |
68 use vars qw(@ISA); | |
69 use strict; | |
70 use Carp; | |
71 | |
72 # Object preamble - inherits from Bio::Root::Root | |
73 | |
74 use Bio::SeqI; | |
75 | |
76 | |
77 @ISA = qw(Bio::SeqI); | |
78 | |
79 | |
80 =head2 delete_feature | |
81 | |
82 Title : delete_feature | |
83 Usage : my $orphanlist=$self->delete_feature($feature,$transcript,$gene); | |
84 Function: deletes the specified $feature from the given transcript, if $transcript is sent and exists and $feature is a feature of $transcript, | |
85 or from $gene if the $feature is a feature of $gene, or from $self if $transcript and $gene are not sent. Keeps track of the features | |
86 of the $gene object that may be left as orphans and returns them as a listref. | |
87 Example : I want to delete transcript 'abc' of gene 'def', with three exons, leaving only transcript 'ghi' with two exons. | |
88 This will leave exons 1 and 3 part of 'ghi', but exon 2 will become an orphan. | |
89 my $orphanlist=$seq->delete_feature($transcript{'abc'},undef,$gene{'def'}); | |
90 $orphanlist is a reference to a list containing $exon{'2'}; | |
91 Returns : a listref of orphaned features after the deletion of $feature (optional) | |
92 Args : $feature - the feature to be deleted | |
93 $transcript - the transcript containing the $feature, so that a $feature can be removed from only one transcript when there are multiple | |
94 transcripts in a gene. | |
95 $gene - the gene containing the $transcript and/or the $feature | |
96 | |
97 | |
98 =cut | |
99 | |
100 sub delete_feature{ | |
101 my ($self,$feature,$transcript,$gene) = @_; | |
102 | |
103 $self->throw_not_implemented(); | |
104 } | |
105 | |
106 | |
107 1; |