This is a test version of Biostars. For the public version, visit https://www.biostars.org.
VCF file format for VCFTOOLs

I used PLINK to generate a vcf and want to use vcftools to analyze it. Here is the first several line of the vcf files:

##fileformat=VCFv4.2
##fileDate=20220520
##source=PLINKv1.90
##contig=<ID=0,length=2147483645
##INFO=<ID=PR,Number=0,Type=Flag,Description="Provisional reference allele, may not be based on real reference genome"
##FORMAT=<ID=GT,Number=1,Type=String,Description="Genotype"
#CHROM  POS ID  REF ALT QUAL    FILTER  INFO    FORMAT  14LM0008_14LM0008   14...**

When I run vcftools, I get the error code:

Error: Unknown option: ##fileDate=20220520

After I delete the first several lines, the new error code occurs:

Error: Unknown option: ##POS

Here is my code:

vcftools --vcf $file --keep $list --max-missing 0.8 --hwe 0.1 --relatedness2 --out ${name}maxmissing.vcf

My vcftools version is 0.1.15

I wonder why this will happen.

plink vcftools vcf
• 1,387 views • link
updated by Pierre Lindenbaum Pierre Lindenbaum • written by syxsu • syxsu

Have you tried FixVcfHeader yo try to sanity check the VCFs?

• 0 views • link
written by iraun iraun

what is the output of

echo vcftools --vcf "$file" --keep "$list" --max-missing 0.8 --hwe 0.1 --relatedness2 --out "${name}maxmissing.vcf"
file $file
file $list
• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

0 answers

No answers yet.

Log in to answer this question.

More posts like this