bedtools interset doesn't return a VCF file?
I am filtering a VCF file with a bed file using Bedtools. I have carried out this successfully with
bedtools intersect -wb -a myVCF.vcf -b myBEDfile.bed > output.txt
However, what I want is to get a VCF file with the metadata and the header because after this, the application I am developing will require the format of a VCF file. I thought that is a very common filtering operation but reading bedtools interest documentation I can't find anything.
1 answer
Log in to answer this question.
More posts like this
-
Two identical reads in my bam files but with a different in CIGAR are mapped in a different locatio…
written by ManuelDB
-
bedtools intersect multiple files using for loop
written by clemaster •
-
Beedtools intersect unable to open a vcf file
written by ManuelDB
-
Determine the significance of peak overlaps (bedtools)
written by kstangline •
-
bedtools intersect correct usage
written by storm1907 •
-
Does bedtools intersect output overlaps with isoforms?
written by Kash
-
Bedtools intersect behavior under different coordinate systems?
written by edsouza •
-
extract unique variant in sample
written by Medhat
-
bedtools segmentation fault
written by denise.mauldin •
-
bedtools intersect mistakes
written by schelarina •
