This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Haplotype Frequency Of A Range From 1000 Genomes Data Using Tabix And Other Utils

Aim: Download public data in a range, calculate the frequency of haplotypes in that region for overall and each ethnic population.

I want to download the region around BRCA1 from the 1000genomes data. 

tabix -fh http://ftp.1000genomes.ebi.ac.uk/vol1/ftp/release/20110521/ALL.chr17.phase1_release_v3.20101123.snps_indels_svs.genotypes.vcf.gz 17:41,196,312-41,277,340 >BRCA1_1000g_20101123.vcf

So I have my BRCA1 genotype data and I want to check the frequency just as a QC measure.

vcftools --gzvcf BRCA1_1000g_20101123.vcf.gz \
    --freq \
    --out BRCA1Copy_1000g_20101123.vcf.freq

Now, I want to now find the common and "all" haplotype blocks and the frequency of haplotypes in this region.
What filters should be applied on allele frequencies . Any help is very much appreciated.

haplotype tabix plink haploview
• 3,404 views • link
written by jvijai • jvijai

Can you make an example of the output that you would expect to see? Do you want the frequency of all the possible haplotypes, or only of the ancestral alleles haplotypes?

• 0 views • link
written by Giovanni M Dall'Olio Giovanni M Dall'Olio

0 answers

No answers yet.

Log in to answer this question.

More posts like this
  • perl, vcf-subset: Broken VCF header, no column names?
    written by severalorks severalorks
    EDIT: I found I can also use Data Slicer from 1000Genomes browser tools to do what I need I ran: tabix -h /data/ALL.chrX.phase3_shapeit2_mvncall_integrated_v1b.20130502.genotypes.vcf.gz | perl …
  • extract specific population (vcf) from 1000 genomes - keep option
    written by muralinmars • muralinmars
    hi all, how can I extract 1000 genome data in vcf format of a specific population using vcftools. Currently I am using tabix to download …
  • [Error] getting fasta file from 1000genomes
    written by Mari • Mari
    Hi, I'd like to get a fasta file of haplotypes from 1000 genomes. I ran the commands below: http://www.1000genomes.org/faq/are-there-any-fasta-files-containing-1000-genomes-variants-or-haplotypes ``` tabix -h ftp://ftp.1000genomes.ebi.ac.uk/vol1/ftp/release/20110521/ALL.chr17.phase1_release_v3.20101123.snps_indels_svs.genotypes.vcf.gz 17:1471000-1472000 | …
  • Haplotype Frequency Calculator For Vcf Files
    written by diviya.smith • diviya.smith
    <p>Aim: Download public data in a range, calculate the haplotype frequency for SNPs in the region for each ethnic population.</p> <p>I want to compute haplotype …
  • Variant Effect And The Impact Of Samples
    written by User 1933 • User 1933
    <p>I have extracted a region from 1KG webserver like </p> <pre><code>mysql --user=genome --host=genome-mysql.cse.ucsc.edu -A -D hg19 -N -e 'select concat("ftp://ftp-trace.ncbi.nih.gov/1000genomes/ftp/release/20110521/ALL.",K.chrom,".phase1_release_v3.20101123.snps_indels_svs.genotypes.vcf.gz"),concat(RIGHT(K.chrom,LENGTH(chrom)-3),":",MIN(K.txStart)+1,"-",MAX(K.txEnd)) from knownGene as K,kgXref as …
  • Variant Effect Of Each Individual In The 1Kg Project
    written by User 1933 • User 1933
    <p><strong>UPDATE : I added another relevant question to this post instead of generating a new post.</strong> So, I should start my post by saying that …
  • What Type Of Information Does The Data Downloaded From The 1000 Genomes Contain?
    written by Pappu Pappu
    <p>With this command, I only get a vcf file containing numbers:</p> <pre><code>tabix -fh ftp://ftp-trace.ncbi.nih.gov/1000genomes/ftp/release/20110521/ALL.chr11.phase1_release_v3.20101123.snps_indels_svs.genotypes.vcf.gz 11:62379194-62382592 &gt;1000.vcf </code></pre> <p>Could you tell me what is going on …
  • 1000 Genomes - Sample Subset Of A Vcf File
    written by pablo.riesgo • pablo.riesgo
    <p>Hi there,</p> <p>I'm trying to retrieve genotype data for a given sample from the 1000 genomes FTP repository, following their guidelines: <a href="http://www.1000genomes.org/faq/how-do-i-get-sub-section-vcf-file">http://www.1000genomes.org/faq/how-do-i-get-sub-section-vcf-file</a></p> <p>So I …
  • How To Interpret Ciend And Cipos Fields In Vcf Files From 1000 Genomes?
    written by agnieszka • agnieszka
    <p>Here is an example variant for genome NA20317. It comes from 1000 genomes data (file: <a href='ftp://ftp.ncbi.nih.gov///1000genomes/ftp/release/20110521/ALL.chr8.phase1_release_v3.20101123.snps_indels_svs.genotypes.vcf.gz'>ftp://ftp.ncbi.nih.gov///1000genomes/ftp/release/20110521/ALL.chr8.phase1_release_v3.20101123.snps_indels_svs.genotypes.vcf.gz</a>) and was obtained with tabix and vcf-subset script:</p> …
  • How To Filter Vcf Files From 1000 Genomes Release V3.2010-11 (Alternative Source)?
    written by user56 user56
    <p>I want to use VCF files from WGS to arrive at pharmacogenomics clinical recommendations (relevant to a single patient, not a population).</p> <p>I decided to …