R/Python script to collect allele frequency for batch of SNPs in gnomAD database

Question

[Shicheng Guo]

Subject: Seeking R/Python Script for Allele Frequency among multiple population Data Collection from gnomAD

Message: Hello everyone,

I'm in search of R or Python script capable of fetching allele frequency data for a batch of SNPs from the gnomAD database, encompassing all populations. Ideally, this script would accept SNP IDs as input and output a data frame containing the allele frequency across different populations, along with the count of homozygous carriers.

Appreciate any pointers or resources you could share. Thank you!

Potential but not final solution: https://github.com/KalinNonchev/gnomAD_DB/blob/master/scripts/GettingStartedwithGnomAD_DB.ipynb

[Pierre Lindenbaum]

bcftools query -i xxxx -f yyyy in.vcf.gz