Using vcftools to find # mappable positions and # of segregating polymorphisms sites

I just installed vcftools and I'm trying to find a way to use it to tell me the # of mappable positions and # of segregating bi-allelic polymorphism sites. Does anybody know how to do this?

The vcf file I'm using is from the 1000 genomes project: ftp://ftp-trace.ncbi.nih.gov/1000genomes/ftp/release/20110521/ALL.chr21.phase1_release_v3.20101123.snps_indels_svs.genotypes.vcf.gz

vcf snp polymorphism

• 2,833 views • link

updated 4.4 years ago by Ram

• written 11.6 years ago by jolespin

0 answers

No answers yet.

More posts like this

perl, vcf-subset: Broken VCF header, no column names?
written 10.0 years ago by 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 …
How to get individual chromosome sequence in fasta format from vcf.gz and its vcf.gz.tbi file of 10…
written 10.6 years ago by David Simmons •

Hi, I am new to the field of Bioinformatics. I have downloaded files `ALL.chr1.phase1_release_v3.20101123.snps_indels_svs.genotypes.vcf.gz` and its tabix file `ALL.chr1.phase1_release_v3.20101123.snps_indels_svs.genotypes.vcf.gz.tbi` from ftp://ftp-trace.ncbi.nih.gov/1000genomes/ftp/release/20110521/. 1) The above vcf.gz …
Variant in earlier builds of 1000 genomes not found in later builds
written 11.0 years ago by jasonlouisstein •

Hi Biostars, I would like to calculate LD of genome-wide significant variants found within the most recent schizophrenia GWAS paper (see [Supp Table 2][1]) using …
extract specific population (vcf) from 1000 genomes - keep option
written 11.9 years ago by 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 12.1 years ago by 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 | …
Generate vcf.gz file and its index file vcf.gz.tbi
written 13.6 years ago by lyz10302012

Can anyone tell me how to generate `vcf.gz` file and its index file `vcf.gz.tbi` in 1000 Genomes Project? <a href='ftp://ftp-trace.ncbi.nih.gov/1000genomes/ftp/release/20110521/'>ftp://ftp-trace.ncbi.nih.gov/1000genomes/ftp/release/20110521/
Vcf File With Variants And Genotypes For A "Single Individual"
written 13.9 years ago by Nasrin •

Hi there, I am working with a software that needs a VCF file with variants and genotypes for a **single individual** and a single chromosome. …
1000 Genomes - Sample Subset Of A Vcf File
written 14.0 years ago by pablo.riesgo •

Hi there, 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> So I …
How To Interpret Ciend And Cipos Fields In Vcf Files From 1000 Genomes?
written 14.2 years ago by agnieszka •

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: …
How To Filter Vcf Files From 1000 Genomes Release V3.2010-11 (Alternative Source)?
written 14.2 years ago by user56

I want to use VCF files from WGS to arrive at pharmacogenomics clinical recommendations (relevant to a single patient, not a population). I decided to …