forked from Ensembl/VEP_plugins
-
Notifications
You must be signed in to change notification settings - Fork 0
/
CCDSFilter.pm
55 lines (32 loc) · 1.65 KB
/
CCDSFilter.pm
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
=head1 LICENSE
Copyright (c) 1999-2011 The European Bioinformatics Institute and
Genome Research Limited. All rights reserved.
This software is distributed under a modified Apache license.
For license details, please see
http://www.ensembl.org/info/about/code_licence.html
=head1 CONTACT
Graham Ritchie <grsr@ebi.ac.uk>
=cut
=head1 NAME
CCDSFilter
=head1 SYNOPSIS
mv CCDSFilter.pm ~/.vep/Plugins
perl variant_effect_predictor.pl -i variations.vcf --plugin CCDSFilter
=head1 DESCRIPTION
A simple VEP filter plugin that limits output to variants that
fall in transcripts which have CCDS coding sequences.
=cut
package CCDSFilter;
use strict;
use warnings;
use base qw(Bio::EnsEMBL::Variation::Utils::BaseVepFilterPlugin);
sub feature_types {
return ['Transcript'];
}
sub include_line {
my ($self, $tva) = @_;
my $t = $tva->transcript;
my @entries = grep {$_->database eq 'CCDS'} @{$t->get_all_DBEntries};
return scalar @entries;
}
1;