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
bcftools query -i xxxx -f yyyy in.vcf.gz