annotate variant_effect_predictor/Bio/EnsEMBL/Utils/VegaCuration/Gene.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 =head1 SYNOPSIS
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
24
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
25 =head1 DESCRIPTION
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
26
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
27 =head1 METHODS
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
28
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
29 =cut
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
30
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
31 package Bio::EnsEMBL::Utils::VegaCuration::Gene;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
32
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
33 use strict;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
34 use warnings;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
35 use vars qw(@ISA);
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
36
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
37 use Bio::EnsEMBL::Utils::ConversionSupport;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
38
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
39 @ISA = qw(Bio::EnsEMBL::Utils::ConversionSupport);
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
40
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
41
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
42 =head2 find_gaps
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
43
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
44 Args : arrayref of B::E::Transcripts
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
45 Example : my $gaps = find_gaps($all_transcripts)
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
46 Description: identifies regions of a gene that are not covered by any transcript
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
47 Returntype : int
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
48 Exceptions : none
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
49 Caller : internal
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
50
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
51 =cut
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
52
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
53 sub find_gaps {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
54 my $self = shift;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
55 my ($all_transcripts) = @_;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
56 my $gaps = 0;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
57 my @sorted_transcripts = sort {$a->start <=> $b->start || $b->end <=> $a->end} @{$all_transcripts};
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
58 if ( my $first_transcript = shift @sorted_transcripts ) {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
59 my $pos = $first_transcript->end;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
60 foreach my $transcript (@sorted_transcripts) {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
61 next if ($transcript->end < $pos );
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
62 if ($transcript->start < $pos && $transcript->end > $pos ) {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
63 $pos = $transcript->end;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
64 next;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
65 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
66 elsif ($transcript->end > $pos) {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
67 $gaps++;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
68 $pos = $transcript->end;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
69 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
70 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
71 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
72 return $gaps;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
73 }