Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
DEXSeq error in estimateDispersions: match.arg(start.method, c("log(y)", "mean")) fpadilla Bioinformatics 14 07-03-2013 03:11 PM
Relatively large proportion of "LOWDATA", "FAIL" of FPKM_status running cufflink ruben6um Bioinformatics 3 10-12-2011 01:39 AM
The position file formats ".clocs" and "_pos.txt"? Ist there any difference? elgor Illumina/Solexa 0 06-27-2011 08:55 AM
Bfast output and "Empty Sequence Dictionary" in .sam output aiden Bioinformatics 1 05-28-2010 07:50 PM
"Systems biology and administration" & "Genome generation: no engineering allowed" seb567 Bioinformatics 0 05-25-2010 01:19 PM

Thread Tools
Old 07-09-2014, 06:45 AM   #1
Junior Member
Location: Switzerland

Join Date: Mar 2013
Posts: 9
Default Vcftools output in "012" format is empty

Hi everyone,
I've been trying to get vcftools to alter my vcf files (97 samples, around 1,000,000 variants) into matrix format, however the ".012" file is always empty .
The other files (FILENAME.012.indv, FILENAME.012.pos) are fine (i.e. the IDs of the 97 are listed in FILENAME.012.indv and the variant positions are listed in FILENAME.012.pos). There is no error message, and the log file looks good.

Here's my command line:
vcftools --gzvcf TEST.vcf.gz --012 --out TEST

I did notice that there are 97 emtpy lines in the FILENAME.012 file, coinciding with the number of samples I've got in my original .vcf file.

When I try this with the --plink option, I don't get any output either, but the --plink-tped seems to work fine.

Any ideas whats going on here and how I can get my .vcf file into 0 1 2 format?

Help very much appreciated,
CowGirl is offline   Reply With Quote
Old 07-09-2014, 01:24 PM   #2
Location: Mountain View, CA

Join Date: Jun 2013
Posts: 15

If a slightly different format is ok (main body is 012, but marker/sample IDs are saved differently), try

plink --vcf TEST.vcf.gz --recode A --out TEST

(this requires plink 1.9)
chrchang is offline   Reply With Quote
Old 07-10-2014, 01:24 AM   #3
Junior Member
Location: Switzerland

Join Date: Mar 2013
Posts: 9

Hi Chrchang,

super, thanks for the tip!

Works like a charm, but I had to include "--cow" to adjust the number of chromosomes PLINK should read.

It would still be good to know whats happening with the vcftools results - anyone have an idea whats going on there? And why there's no error message?


Last edited by CowGirl; 07-10-2014 at 01:27 AM.
CowGirl is offline   Reply With Quote

empty 012 output, vcftools

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 05:07 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2019, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO