Hey guys,
Is there any way to use bedtools intersect on multiple VCF files to get only variants that are common for all those files? From the manual I counclude that in standard use output would be variants that are present in query and at least one database file, while I need those that would be present in all databse files.
All the best, Agata
•
written
by
ciemanek
1 answer
I think you need to use bcftools isec
Creates intersections, unions and complements of VCF files. Depending on the options, the program can output records from one (or more) files which have (or do not have) corresponding records with the same position in the other files.
assuming that you have 3 files
bcftools isec -p dir -n=3 A.vcf.gz B.vcf.gz C.vcf.gz
Note: vcf-isec do the same work
Log in to answer this question.
-
Bedtools loop over subfolders
written by cwang130 •
-
How to find variants which are not present in second file
written by shubhamkumbhar420 •
-
Tutorial: Real World BEDtools. Practical Examples for Variant Filtering and Gene Annotations
written by sdsubaru33 •
-
Genomics regions that common to atleast 10 cell types
written by rohitsatyam102
-
Categorize patients based on continuous variables
written by ciemanek
-
Combine VCF files with multiple alternative alleles in GATK
written by ciemanek
-
VCFtools intersect vs. GATK intersect
written by ciemanek
-
Intersection of a BED genes with multiple VCFs: report only if at least N intersection.
written by 乙
-
Intersect multiple BED files
written by int11ap1
-
Integrative Bioinformatics: How To Integrate My Queries?
written by Chong Tang
