Skip to content
Snippets Groups Projects
Commit ab216c1d authored by Pontus Larsson's avatar Pontus Larsson
Browse files

Changed script to accommodate schema changes

parent 5969bfb9
No related branches found
No related tags found
No related merge requests found
......@@ -33,8 +33,8 @@ use Getopt::Long;
# A hard-coded hash for consequence type => attrib_type.code mapping
my %ATTRIBUTE_CODE = (
'STOP_LOST' => 'StopLost',
'STOP_GAINED' => 'StopGained'
'stop_lost' => 'StopLost',
'stop_gained' => 'StopGained'
);
# The minimum frequency of a STOP_*-causing allele in a population to be included
......@@ -244,10 +244,6 @@ if (!defined($file)) {
my $sth = $vdb->dbc->prepare($stmt) or die("Error preparing statement $stmt");
$sth->execute();
# Get a TranscriptVariationAdaptor
my $tv_adaptor = $vdb->get_TranscriptVariationAdaptor() or throw("Could not get a TranscriptVariationAdaptor from " . $vdb->dbc->dbname());
my $stop_codons = qq{TGA/TAA/TAG};
# For each variation_feature, check that the source, population and minor allele frequencies are 'dbSNP', 'CSHL-HapMap-*' and '>= $FREQUENCY_CUTOFF'
while (my ($transcript_stable_id,$variation_id,$variation_name,$allele_string,$consequence_string) = $sth->fetchrow_array()) {
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment