forked from Ensembl/VEP_plugins
-
Notifications
You must be signed in to change notification settings - Fork 0
/
NonSynonymousFilter.pm
58 lines (35 loc) · 1.76 KB
/
NonSynonymousFilter.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
55
56
57
=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 <[email protected]>
=cut
=head1 NAME
NonSynonymousFilter
=head1 SYNOPSIS
mv NonSynonymousFilter.pm ~/.vep/Plugins
perl variant_effect_predictor.pl -i variations.vcf --plugin NonSynonymousFilter
=head1 DESCRIPTION
A simple example VEP filter plugin that limits output to non-synonymous variants
=cut
package NonSynonymousFilter;
use strict;
use warnings;
use base qw(Bio::EnsEMBL::Variation::Utils::BaseVepFilterPlugin);
sub feature_types {
return ['Transcript'];
}
sub include_line {
my ($self, $tva) = @_;
# just check if there are alternative amino acids in the
# pep_allele_string, this means we'll catch stop gained
# or lost as well
if (my $pep_alleles = $tva->pep_allele_string) {
return $pep_alleles =~ /\//;
}
return 0;
}
1;