Question: Filtering a vcf file with vcftools, where is my output?
gravatar for devenvyas
6.2 years ago by
Stony Brook
devenvyas630 wrote:

Hello, I am a trying to filter some vcf files based off a set of rs#s ($FILE corresponds to chromosome #, so for now that is just 21, 22, and Y)

vcftools --gzvcf "sorted_AltaiNea.hg19_1000g."$FILE".mod.vcf.gz" --snps 330k.txt --out "filtered_AltaiNea.hg19_1000g."$FILE"_"

I run that code on the cluster, but then all I get out are


Which contain text as such:

VCFtools - v0.1.11
(C) Adam Auton 2009

Parameters as interpreted:
    --gzvcf sorted_AltaiNea.hg19_1000g.21.mod.vcf.gz
    --out filtered_AltaiNea.hg19_1000g.21_
    --snps 330k.txt

Using zlib version: 1.2.3
Versions of zlib >= 1.2.4 will be *much* faster when reading zipped VCF files.
Index file is older than variant file. Will regenerate.
Building new index file.
    Scanning Chromosome: 21
    Warning - file contains entries with the same position. These entries will be processed separately.

Writing Index file.
File contains 35104060 entries and 1 individuals.
Applying Required Filters.
Keeping sites by user-supplied list
After filtering, kept 1 out of 1 Individuals
After filtering, kept 5364 out of a possible 35104060 Sites
Run Time = 176.00 seconds

I do not have any vcf output, just those log files. Is there something I am doing wrong? Thanks!


snp vcftools vcf • 15k views
ADD COMMENTlink modified 6.2 years ago by chefer300 • written 6.2 years ago by devenvyas630
gravatar for chefer
6.2 years ago by
Pretoria, ZA
chefer300 wrote:

You should add the --recode flag to produce a new vcf file.

See example 2 on

ADD COMMENTlink modified 9 months ago by RamRS27k • written 6.2 years ago by chefer300

You must be kidding me!

ADD REPLYlink written 2.5 years ago by rmf940

so what happened if we don't add --recode? if just delete the SNPs from the raw input?

ADD REPLYlink modified 9 months ago by RamRS27k • written 23 months ago by Shicheng Guo8.3k
Please log in to add an answer.


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