Hi, I downloaded some data from gnomad - https://gnomad.broadinstitute.org/downloads. it comes in the form of VCF.bgz file and I would like to read it as a vcf file. found some code to do it in python :

#!/usr/bin/env python3
import gzip
ifile = gzip.GzipFile("gnomad.genomes.r2.1.1.sites.2.vcf.bgz")
ofile = open("truncated.vcf", "wb")
LINES_TO_EXTRACT = 100000

for line in range(LINES_TO_EXTRACT):
    ofile.write(ifile.readline())

ifile.close()
ofile.close()

i tried it on my data :

import gzip
ifile = gzip.GzipFile("gnomad.exomes.r2.1.1.sites.4.vcf.bgz")
ofile = open("truncated.vcf", "wb")
LINES_TO_EXTRACT = 1000

for line in range(LINES_TO_EXTRACT):
    ofile.write(ifile.readline())

ifile.close()
ofile.close()

but it returned an error :

FileNotFoundError: [Errno 2] No such file or directory: 'gnomad.exomes.r2.1.1.sites.4.vcf.bgz'

although I downloaded it to my computer.

is there any way to open vcf.bgz file as vcf file in python?