Question: Vcftools Get-Info Output Name
0
gravatar for mark.dunning
6.1 years ago by
mark.dunning200
The University of Sheffield
mark.dunning200 wrote:

Hi all,

I am using the vcftools to extract data from a VCF file into tabular form, so that I can read into R and do some computations and visualisation

e.g.

vcftools --vcf sample1.vcf --get-INFO MQ --get-INFO MQ0

The file it creates is always called out.INFO. Is there really no way of controlling the name of the output file to something like sample1.INFO? If I want to extract info from several vcfs in the same directory, they will overwrite each other.

Regards,

Mark

vcf vcftools • 2.2k views
ADD COMMENTlink modified 6.1 years ago by Matt Shirley8.7k • written 6.1 years ago by mark.dunning200
2
gravatar for Matt Shirley
6.1 years ago by
Matt Shirley8.7k
Cambridge, MA
Matt Shirley8.7k wrote:

Based on the documentation, I think you are able to specify something like this:

vcftools --vcf sample1.vcf --get-INFO MQ --get-INFO MQ0 --out sample1

You can write a shell script that will set the --out prefix based on the input filename. If you need help with this, just ask.

ADD COMMENTlink written 6.1 years ago by Matt Shirley8.7k

Brilliant, just what I was looking for...Should have read the documentation in a bit more detail first though :/

ADD REPLYlink written 6.1 years ago by mark.dunning200
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: 1640 users visited in the last hour