Merge two vcf's, keep only intersection of REF/ALT alleles

Question: Merge two vcf's, keep only intersection of REF/ALT alleles

6 weeks ago by

Hi all,

I would like to merge two vcf files chr1g.vcf.gz and chr1hk.vcf.gz. I would like the resulting file to have only the intersection of the two files upon merge.

chr1g.vcf.gz excerpt:

CHROM .POS ....ID... REF... ALT ... QUAL... FILTER....... INFO

chr1 ..... 10031 .. NA .. T ....... C .. NA .. AS_VQSR ..AC=0,etc

chr1 ..... 10055 .. NA .. T ....... C .. NA .. AS_VQSR ..AC=0,etc

chr1 ..... 10061... NA... T........ C .. NA .. AS_VQSR ..AC=0,etc

chr1 ..... 10061... NA ... T ...... TAACC.. .. NA .. AS_VQSR ...AC=0,etc

chr1 .....10109 ... NA ... A....... T .. NA .. AS_VQSR ..AC=0,etc

chr1 .... 10109... NA... AACCCT A .. NA .. AS_VQSR .. AC=0,etc

chr1hk.vcf.gz excerpt:

CHROM .POS ....ID... REF... ALT ... QUAL... FILTER....... INFO

chr1 ..... 10055 .. NA .. T ....... C .. NA .. AS_VQSR ..AC=0,etc

chr1 ..... 10061... NA... T........ C .. NA .. AS_VQSR ..AC=0,etc

chr1 .....10109 ... NA ... A....... T .. NA .. AS_VQSR ..AC=0,etc

chr1 .... 10109... NA... AACCCT A .. NA .. AS_VQSR .. AC=0,etc

Merging goal ex:

CHROM .POS ....ID... REF... ALT ... QUAL... FILTER....... INFO

chr1 ..... 10055 .. NA .. T ....... C .. NA .. AS_VQSR ..AC=0,etc

chr1 ..... 10061... NA... T........ C .. NA .. AS_VQSR ..AC=0,etc

chr1 .....10109 ... NA ... A....... T .. NA .. AS_VQSR ..AC=0,etc

chr1 .... 10109... NA... AACCCT A .. NA .. AS_VQSR .. AC=0,etc

.. .. The code I have been working with is as follows: "bcftools merge --merge none chr1g.vcf.gz chr1hk.vcf.gz > chr1merge.vcf" This code works to merge based off of REF/ALT allele matches, but is the union of the two original files. How can I tweak it to keep only the intersection?

Thank you!

vcftools bcftools merge intersection vcf • 167 views

ADD COMMENT • link •

modified 16 days ago by Elucidata • 120 • written 6 weeks ago by just learning • 0

6 weeks ago by

Pierre Lindenbaum ♦ 134k

France/Nantes/Institut du Thorax - INSERM UMR1087

Pierre Lindenbaum ♦ 134k wrote:

How can I tweak it to keep only the intersection?

process the output of bcftools merge with an invocation of bcftools isec

ADD COMMENT • link written 6 weeks ago by Pierre Lindenbaum ♦ 134k

Thank you for the reply! I was attempting to avoid using bcftools isec after the merge because it outputs four extremely large data sets- however if this is the only option I will work with it.

ADD REPLY • link written 6 weeks ago by just learning • 0

16 days ago by

Elucidata • 120

Elucidata • 120 wrote:

One of the common tools used to merge and intersect vcf files based on the REF/ALT alleles is bedtools intersect (also can be used as intersectbed). One can use this tool to find overlapping entries between files or exclusive entries between files by mentioning the corresponding flags of the tool.

You can look up the variety of uses the tool offers here.

In your case, to get the overlapping (intersecting) entries of the two files with the output files containing the entries of the chr1g.vcf.gz file (including the REF/ALT entries) you can use the following command:

bedtools intersect -header -wa -a chr1g.vcf.gz -b chr1hk.vcf.gz > intersect.vcf OR intersectBed -header -wa -a chr1g.vcf.gz -b chr1hk.vcf.gz > intersect.vcf

You can change the -wa flag to -wb flag to check the entries of file chr1hk.vcf.gz

ADD COMMENT • link written 16 days ago by Elucidata • 120

Please log in to add an answer.

Similar posts • Search »