convert and merge somatic vcf file to plink format
I want to convert vcf files to plink format, but the vcf are missing IDs (annotated as "." for all variants) instead of the usual rsXXX that we usually see for germline mutations. After I convert I try to merge all the plink files using this command:
plink --bfile sample1.vcf --merge-list files.bed.txt --make-bed --out Merged.Plink
But I get a lot of warnings:
Warning: Multiple chromosomes seen for variant '.'.
Is it possible to work with somatic vcf files using plink?
• 910 views
•
link
updated
by
LChart
•
written
by
Palgrave •
•
written
by
Palgrave •
1 answer
You can use --set-missing-var-ids "@:#" to re-set the (missing) variant IDs to chromosome:position nomenclature. See https://www.cog-genomics.org/plink/1.9/data#set_missing_var_ids.
Log in to answer this question.
More posts like this
-
Plink merge errors.
written by Dinmukhamed •
-
Problem merging plink files for somatic variants.
written by Palgrave •
-
Processing WES VCF for case control GWAS analysis
written by sonsunjirachote •
-
Merging multiple PLINK files
written by Nejla •
-
Plink genome flag error: Pedigree graph is cyclic
written by rwilliams7653 •
-
PLINK:Error: Variant '.' is not biallelic. To obtain a full list of merge failures
written by williamsbrian5064
-
TreeMix: ERROR: cannot open file treemix.freq.gz
written by k.prata •
-
multi sample vcf to data matrix
written by Omics data mining
-
Using Plink to calculate LD between loci on different chromosomes
written by jr.azaria •
-
Best way to convert VCF to PLINK file format and merge chromosomes?
written by nchuang
