Question: Vcftools --Geno And --Hwe Filter Options -No Output File Created
1
gravatar for ttom
6.6 years ago by
ttom210
ttom210 wrote:

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 • 6.8k views
ADD COMMENTlink modified 6.0 years ago by Biostar ♦♦ 20 • written 6.6 years ago by ttom210
5
gravatar for ttom
6.6 years ago by
ttom210
ttom210 wrote:
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 COMMENTlink modified 6.6 years ago • written 6.6 years ago by ttom210

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

ADD REPLYlink written 6.6 years ago by Giovanni M Dall'Olio26k

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 REPLYlink modified 6.6 years ago • written 6.6 years ago by ttom210
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: 1407 users visited in the last hour