Skip to content

Commit

Permalink
removing AF
Browse files Browse the repository at this point in the history
  • Loading branch information
costero-e committed Feb 12, 2024
1 parent bcbbc0c commit 5be7476
Showing 1 changed file with 9 additions and 16 deletions.
25 changes: 9 additions & 16 deletions genomicVariations_vcf.py
Original file line number Diff line number Diff line change
Expand Up @@ -27,16 +27,11 @@ def generate(dict_properties):

pbar = tqdm(total = num_rows)
for v in vcf:
'''
try:
if v.INFO["AF"] > 0.1: continue
except Exception:
pass
try:
if v.INFO.get('VT') == 'SV': continue
except Exception:
pass
'''

#print(v)
dict_to_xls={}
ref=v.REF
Expand All @@ -52,16 +47,14 @@ def generate(dict_properties):
except Exception:
pass
#print(v.INFO.get('ANN'))
try:
if v.INFO.get('ANN') is not None:
annot = v.INFO.get('ANN')
annotations = annot.split("|")
dict_to_xls['molecularAttributes|molecularEffects|label'] = annotations[1]
dict_to_xls['molecularAttributes|molecularEffects|id'] = "ENSGLOSSARY:0000174"
dict_to_xls['molecularAttributes|aminoacidChanges'] = "."
dict_to_xls['molecularAttributes|geneIds'] = annotations[3]
except Exception:
pass
if v.INFO.get('ANN') is not None:
annot = v.INFO.get('ANN')
annotations = annot.split("|")
dict_to_xls['molecularAttributes|molecularEffects|label'] = annotations[1]
dict_to_xls['molecularAttributes|molecularEffects|id'] = "ENSGLOSSARY:0000174"
dict_to_xls['molecularAttributes|aminoacidChanges'] = "."
dict_to_xls['molecularAttributes|geneIds'] = annotations[3]



zigosity={}
Expand Down

0 comments on commit 5be7476

Please sign in to comment.