"select distinct(vf.variation_name) from $SNPDB.transcript_variation as tv, $SNPDB.variation_feature as vf where vf.variation_feature_id = tv.variation_feature_id and tv.transcript_stable_id in(?)"
);
};
# my $snp_sth = eval {
# $dbh->prepare(
#"select distinct(vf.variation_name) from $SNPDB.transcript_variation as tv, $SNPDB.variation_feature as vf where vf.variation_feature_id = tv.variation_feature_id and tv.transcript_stable_id in(?)"
# );
# };
my$haplotypes=$dbh->selectall_hashref(
"select gene_id from gene g, assembly_exception ae where
...
...
@@ -517,12 +517,13 @@ g.seq_region_id=ae.seq_region_id and ae.exc_type='HAP'", [qw(gene_id)]
if($old{'gene_id'}!=$gene_id){
if($old{'gene_id'}){
if($SNPDB&&$DBeq'core'){
my@transcript_stable_ids=
keys%{$old{transcript_stable_ids}};
$snp_sth->execute("@transcript_stable_ids");
$old{snps}=$snp_sth->fetchall_arrayref;
}
# if ( $SNPDB && $DB eq 'core' ) {
# my @transcript_stable_ids =
# keys %{ $old{transcript_stable_ids} };
# $snp_sth->execute("@transcript_stable_ids");
# $old{snps} = $snp_sth->fetchall_arrayref;
# }
if($want_species_orthologs){
$old{orthologs}=
...
...
@@ -643,11 +644,11 @@ g.seq_region_id=ae.seq_region_id and ae.exc_type='HAP'", [qw(gene_id)]