Question: Alt allele count from 1000 genome project
0
gravatar for seta
16 months ago by
seta1.2k
Sweden
seta1.2k wrote:

Hi all,

I would like to compare the Alt allele frequency of a list of thousands SNPs derived from my project with the various populations of 1000 genome project. So, I need alt allele count, could you please help me out how to find this infromation based on GRCh38? here, I read this post about this issue, however, it used genotypes.bcf file of GRCh37 (I couldn't find it for GRCh38), and also as this file is large, it takes almost long time to download and process. I guess this information may be generated before and stored in a databases.

Your suggestions and comments would be highly appreciated.

ADD COMMENTlink modified 16 months ago by RamRS25k • written 16 months ago by seta1.2k
0
gravatar for RamRS
16 months ago by
RamRS25k
Houston, TX
RamRS25k wrote:

You'll need to combine a liftover operation on a 1000g VCF file (unless an hg38 version is already available) with a bcftools annotate operation. The annotate operation can copy over attributes from a reference VCF file to a target VCF file based on matching chrom, pos, ref and alt, assuming both VCFs are normalized.

The -c option lets you rename attributes while copying them over, so you can ensure names like AC don't clash/overwrite.

ADD COMMENTlink modified 16 months ago • written 16 months ago by RamRS25k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1583 users visited in the last hour