def range_filter(reference_genome, variants, reference_name, start, end):
"""Filters variants by range depending on the reference_genome"""
if 'chr' in reference_name:
reference_name = reference_name.replace('chr', '')
variants = variants.filter(Chr=reference_name)
if reference_genome == 'hg36':
variants = variants.order_by('Hg36_Start')
variants = variants.filter(Hg36_Start__lt=end, Hg36_End__gt=start)
elif reference_genome == 'hg37':
variants = variants.order_by('Hg37_Start')
variants = variants.filter(Hg37_Start__lt=end, Hg37_End__gt=start)
elif reference_genome == 'hg38':
variants = variants.order_by('Hg38_Start')
variants = variants.filter(Hg38_Start__lt=end, Hg38_End__gt=start)
return variants
def range_filter(reference_genome, variants, reference_name, start, end):
"""Filters variants by range depending on the reference_genome"""
if 'chr' in reference_name:
reference_name = reference_name.replace('chr', '')
variants = variants.filter(Chr=reference_name)
if reference_genome == 'hg36':
variants = variants.order_by('Hg36_Start')
variants = variants.filter(Hg36_Start__lt=end, Hg36_End__gt=start)
elif reference_genome == 'hg37':
variants = variants.order_by('Hg37_Start')
variants = variants.filter(Hg37_Start__lt=end, Hg37_End__gt=start)
elif reference_genome == 'hg38':
variants = variants.order_by('Hg38_Start')
variants = variants.filter(Hg38_Start__lt=end, Hg38_End__gt=start)
return variants