Vcftools --Geno And --Hwe Filter Options -No Output File Created
1
1
Entering edit mode
8.3 years ago
ttom ▴ 210

I used the following vcftools filter commands --geno and --hwe to filter VCF file

vcftools --vcf sample.vcf --hwe 0.001

vcftools --vcf sample.vcf --geno 0.9

Both created log files named out.log. But could not find the output file. Like for example when we use vcftools --relatedness --vcf generates default output file named out.relatedness

However the log file (out.log) for --geno option has the following lines

File contains 259375 entries and 435 individuals. Applying Required Filters. Filtering sites by allele frequency and call rate After filtering, kept 435 out of 435 Individuals After filtering, kept 197076 out of a possible 259375 Sites Run Time = 41.00 seconds

Any help much appreciated

Thanks, Tinu

vcftools • 8.1k views
ADD COMMENT
5
Entering edit mode
8.3 years ago
ttom ▴ 210
vcftools –vcf sample.vcf –hwe 0.001 --recode

Adding recode option will produce an output file named out.recode.vcf which will have filtered data.

Adam Auton helped by replying to the same question posted in VCFtools help discussions

http://sourceforge.net/mailarchive/forum.php?thread_name=4F2CFD25-89F8-46A3-A97D-0C000870E78C%40gmail.com&forum_name=vcftools-help

ADD COMMENT
0
Entering edit mode

To specify the suffix for the output file (instead of "out"), use the --out option

ADD REPLY
0
Entering edit mode

Just using --recode would remove contents from the INFO field

--recode-INFO-all would retains INFO field

Ideally it should be

vcftools --vcf sample.vcf --geno 0.9 --recode --recode-INFO-all --out genotypefiltered_sample

ADD REPLY

Login before adding your answer.

Traffic: 1139 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6