Skip to content

Commit

Permalink
Merge pull request #307 from kshefchek/flybase-fix-nt
Browse files Browse the repository at this point in the history
Fix bnode generation in Genotype models
  • Loading branch information
kshefchek committed May 12, 2016
2 parents 0f62b6b + 6e9b281 commit b29821e
Show file tree
Hide file tree
Showing 2 changed files with 12 additions and 5 deletions.
1 change: 1 addition & 0 deletions dipper/models/Genotype.py
Original file line number Diff line number Diff line change
Expand Up @@ -476,6 +476,7 @@ def addReagentTargetedGene(
# akin to a variant locus
if targeted_gene_id is None:
targeted_gene_id = '_' + gene_id + '-' + reagent_id
targeted_gene_id = targeted_gene_id.replace(":", "")
self.gu.addIndividualToGraph(
self.graph, targeted_gene_id, targeted_gene_label,
self.genoparts['reagent_targeted_gene'], description)
Expand Down
16 changes: 11 additions & 5 deletions dipper/sources/FlyBase.py
Original file line number Diff line number Diff line change
Expand Up @@ -1589,7 +1589,8 @@ def _process_feature_relationship(self, limit):
# subject to type_id = 219,33 object type_id 219 #??? TEC
# subject = variation
# object = gene
allele_id = gene_id = None
allele_id = None
gene_id = None
if subject_id in self.idhash['allele']:
allele_id = self.idhash['allele'][subject_id]
if object_id in self.idhash['gene']:
Expand Down Expand Up @@ -1634,7 +1635,10 @@ def _process_feature_relationship(self, limit):
"this thing %s is not a gene", feature_id)
elif int(type_id) == 59983: # associated_with

allele_id = gene_id = reagent_id = ti_id = None
allele_id = None
gene_id = None
reagent_id = None
ti_id = None

if object_id in self.idhash['allele']:
allele_id = self.idhash['allele'][object_id]
Expand Down Expand Up @@ -1683,7 +1687,7 @@ def _process_feature_relationship(self, limit):
# note that this relationship is only specified between
# an allele and a tp. therefore we know the FBal should be
# a transgenic_insertion
allele_id = tp_id = None
allele_id = None

if subject_id in self.idhash['allele']:
allele_id = self.idhash['allele'][subject_id]
Expand All @@ -1705,7 +1709,8 @@ def _process_feature_relationship(self, limit):
# a reagent-targeted-gene (FBal) and
# the reagent that targetes it (FBsf)

allele_id = reagent_id = None
allele_id = None
reagent_id = None

if subject_id in self.idhash['allele']:
allele_id = self.idhash['allele'][subject_id]
Expand All @@ -1722,7 +1727,8 @@ def _process_feature_relationship(self, limit):
elif int(type_id) == 27:
# i'm looking just for the relationships between
# ti and tp features... so doing a bit of a hack
ti_id = tp_id = None
ti_id = None
tp_id = None
if subject_id in self.idhash['feature']:
ti_id = self.idhash['feature'][subject_id]
if not re.search(r'FBti', ti_id):
Expand Down

0 comments on commit b29821e

Please sign in to comment.