bcftools sec for multiple files
Hello,
I am comparing five vcf files each with multiple samples using bcftools isec:
bcftools isec -p dir -n=5 a.vcf.gz b.vcf.gz c.vcf.gz d.vcf.gz e.vcf.gz
It results in these files:
0000.vcf 0001.vcf 0002.vcf 0003.vcf 0004.vcf README.txt sites.txt
However, the output does not show intersection sites among all five vcf files. I am not sure if something is wrong in the command line. I will appreciate any suggestions. Thank you.
1 answer
Please see this How bcftools isec works ? . why do you have -n=5 and not 7 if you have 7 files?
• 0 views
•
link
written
by
onestop_data
Log in to answer this question.
More posts like this
-
PLINK | Updating sex information issue
written by S •
-
Help me understand "for stripped" in bcftools isec output
written by Axzd •
-
How to get intersected variants between multiple vcf files?
written by Dana •
-
Interpreting bcftools isec output
written by pabe •
-
bcftools isec not producing sites.txt
written by steve
-
Find SNPs common to mutants, absent in WT (vcf files)
written by charlesh •
-
bcftools isec: tandem mutations
written by whb •
-
Using BCFtools to intersect multiple VCF files
written by Hadi M •
-
Intersection of 3 multi-sampled vcf files
written by michelle.zl.kee •
-
bcftools isec: weird number of output file
written by pingEde •
