Entering edit mode
10.8 years ago
blueskypie
▴
70
I have the following files:
A.vcfwhich contains one sample s1B.vcfwhich contains one sample s1region.bed
I'd like to merge A.vcf and B.vcf for sites in region.bed, and use A.vcf for sites present in both files. Here is my command but it does not work as I expected:
bcftools isec A.vcf.gz B.vcf.gz -R region.bed -n+1 -o tmp.vcf.gz -O z -c all -w1
Could anyone help?
Thanks a lot!
0
Entering edit mode
- Do you want to get sites that are common for A, B and region.bed?
- Or you want to merge A and B files and find common sites with region.bed?
- And find A sites that are common for B and region.bed?
ADD REPLY
• link
updated 3.6 years ago by
Ram
45k
•
written 10.4 years ago by
Evgeniia Golovina
★
1.3k