bcftools merge single-sample vcf files
1
0
Entering edit mode
4.0 years ago
j.lunger18 ▴ 30

Hi. I'm trying to merge .vcf files from TCGA. Each of these files is from a single sample, which is indicated by the file name. However, the vcf file header contains the following for each sample:

#CHROM  POS ID  REF ALT QUAL    FILTER  INFO    FORMAT  NORMAL  TUMOR

When I try to merge, I get the error:

Error: Duplicate sample names (NORMAL), use --force-samples to proceed anyway.

I could replace NORMAL with the file name so that I can associate each annotation with the correct ID, but I'm wondering if I have to do this. Will --force-samples remove information?

bcftools vcf SNP • 895 views
ADD COMMENT
1
Entering edit mode
4.0 years ago

I would suggest doing bcftools reheader to change all sample names per file to ID-N and ID-T before merging. If not you won't be able to query which Tumor and Normal belong to the same individual.

ADD COMMENT
0
Entering edit mode

Thanks, I'll look into it.

ADD REPLY

Login before adding your answer.

Traffic: 2609 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6