gnomAD

Read vcf.bgz file from gnomAD with HTSlib

gnomAD is very convinient. However, it can only provide the variants in a gene or a region. I am going to provide a software tool that search through the whole genome and output the variants after filtering with the criteria defined by users. The under developing code can be found at https://github.com/peng-gang/ReadGnomad. It is interesting that I can use function ‘bcf_read1’ to read records in vcf.bgz file. The output of this function is a data struction with compressedd information.