VCFtools --exclude-bed Does it include from the BED file or exclude only?
1
0
Entering edit mode
3.8 years ago
joomi • 0

Hello,

I have a question about VCFtools --exclude-bed. I have a bed file that contains <chrom> <start> and <stop> coordinates. However I want to filter my VCF files to only contain those coordinates in the bed file I have made. I know I have to use --exclude-bed however online it is mentioned constantly that --exclude-bed is used to both include and exclude. For example:

--exclude-bed <filename> Include or exclude a set of sites on the basis of a BED file. Only the first three columns (chrom, chromStart and chromEnd) are required. The BED file is expected to have a header line. A site will be kept or excluded if any part of any allele (REF or ALT) at a site is within the range of one of the BED entries.

I only want to make sure that I'm keeping the coordinates in the bed file in my new VCF file. Do I use --exclude-bed? When I try --include-bed as a flag it does not work.

vcftools --vcf $Files/PT_sanger2yTCA454RNvsB73_with_read_groups_filtered_refina2_snps.ann.gatk.vcf --out $Files_made/VCF_intervals_with_transcripts.vcf --exclude-bed $Files_made/VCF_intervals.bed
RNA-Seq • 2.4k views
ADD COMMENT
0
Entering edit mode
2.2 years ago

Hi joomi, maybe too late, but try with --bed option instead of --bed-exclude.

Greetings, Francisco Abad.

ADD COMMENT

Login before adding your answer.

Traffic: 1708 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6