Hi all :) !

I need help: I want to intersect 2 vcf file (fileA and fileB) in order to have three file:

1) common variants fileA and fileB

2) variants in fileA that are not present in fileB

3) variants in fileB that are not present in fileA

I tried with this command:

bcftools isec fileA fileB -p directoryOutput

In my directoryOuput, I found 4 file: 0000.vcf, 0001.vcf, 0002.vcf, 0003.vcf , README.txt and sites.txt.

README.txt reported:

Using the following file names:

0000.vcf for records private to fileA.vcf.gz

0001.vcf for records private to fileB.vcf.gz

0002.vcf for records from fileA.vcf.gz shared by both fileA.vcf.gz fileB.vcf.gz

0003.vcf for records from fileB.vcf.gz shared by both fileA.vcf.gz fileB.vcf.gz

I was exspecting to have three file:

1) for records private to fileA,

2) for records private to fileB

3) for records shared by both fileA and fileB

what did I do wrong?

Thank you in advance

Best