Error: no samples in vcf file

34 views
Skip to first unread message

HA_63

unread,
Dec 3, 2024, 12:27:01 AMDec 3
to plink2-users
Hi, 
I have a .SNV.VCF file that I am trying to convert to plink format, however when I try to convert it I receive "Error: No samples in .vcf file." I am unsure if this is due to the different file extension or if there is a formatting issue (however the file appears to be in the correct format).

I have atached the command and top of the vcf file below. Thanks in advance for any help! 

##fileformat=VCFv4.1
##source=array-analysis-cli 2.1.0
##genomeBuild=38
##reference=file://Homo_sapiens_assembly38.fa
##FORMAT=<ID=GT,Number=1,Type=String,Description="Genotype">
##FORMAT=<ID=GS,Number=1,Type=Float,Description="GenCall score. For merged multi-assay or multi-allelic record, min GenCall score is reported.">
##FORMAT=<ID=BAF,Number=1,Type=Float,Description="B Allele Frequency">
##FORMAT=<ID=LRR,Number=1,Type=Float,Description="Log R Ratio">
##contig=<ID=1,length=248956422>
##contig=<ID=2,length=242193529>
##contig=<ID=3,length=198295559>
##contig=<ID=4,length=190214555>
##contig=<ID=5,length=181538259>
##contig=<ID=6,length=170805979>
##contig=<ID=7,length=159345973>
##contig=<ID=8,length=145138636>
##contig=<ID=9,length=138394717>
##contig=<ID=10,length=133797422>
##contig=<ID=11,length=135086622>
##contig=<ID=12,length=133275309>
##contig=<ID=13,length=114364328>
##contig=<ID=14,length=107043718>
##contig=<ID=15,length=101991189>
##contig=<ID=16,length=90338345>
##contig=<ID=17,length=83257441>
##contig=<ID=18,length=80373285>
##contig=<ID=19,length=58617616>
##contig=<ID=20,length=64444167>
##contig=<ID=21,length=46709983>
##contig=<ID=22,length=50818468>
##contig=<ID=X,length=156040895>
##contig=<ID=Y,length=57227415>
##contig=<ID=MT,length=16569>
#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT
1 632287 rs9651229 C T . PASS . GT:GS:BAF:LRR 0/0:0.3583:0.007867098:-0.049537543

HA_63

unread,
Dec 3, 2024, 12:32:36 AMDec 3
to plink2-users
Apologies! Command used is: plink --vcf pt_3/207656660134_R12C02.snv.vcf --make-bed --out pt_3/pt_3_convert

Chris Chang

unread,
Dec 3, 2024, 7:32:05 AMDec 3
to HA_63, plink2-users
The #CHROM header line is missing a sample-ID entry at the end, and this is why plink is complaining.

--
You received this message because you are subscribed to the Google Groups "plink2-users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to plink2-users...@googlegroups.com.
To view this discussion visit https://groups.google.com/d/msgid/plink2-users/99e1cddf-e7e8-4157-a256-8b266ba7798an%40googlegroups.com.
Reply all
Reply to author
Forward
0 new messages