Hi all,
New to this, and I'm trying to work through our pipeline. I am having difficulty extracting my regions of interest from my merged VCF file. The code is correct, however, I get an error that I can only extract one chr at a time. Wondering if there is something I can do to this to make it extract multiple at once, or if I can use something else.
Note: do not have an available bed file, I have a .txt file with these regions in it. Also, I ran this script with just one chr, and it worked. (Which is how I know it's correct).
Thank you!
Code (imagine this with multiple chr regions to extract):
for FILE in XXX
do
vcftools --gzvcf $FILE --chr chr8 --from-bp 89933325 --to-bp 90003238 --out output --recode --recode-INFO-all
done
1 answer
You can use bedtools intesect. Keep regions of interest in a separate bed file. Intersect vcf file with bed file with regions of interest. Your text file should have 3 columns (chromosome name, start and end positions). Example code and example bed for regions of interest:
$ bedtools intersect -header -wa -a test.vcf -b test.txt
$ cat test.txt
chr12 133433173 133433373
chr12 133781605 133781705
-header will print header from VCF
-wa will print matching lines from VCF (file a)
From your code, it seems you have multiple VCF files and you are trying to extract single region of interest from each vcf file in a loop. Let us say you want to extract multiple regions of interest from each vcf file. You can use parallel:
parallel bedtools intersect -header -wa -a test.vcf -b test.txt >{.}.out ::: *.vcf
Log in to answer this question.
-
Using tabix to subset a region from a VCF file
written by michael.flower.14 •
-
VCF file
written by Nai •
-
--max-missing filter in vcftools removes most of my SNP data
written by truebeliever24 •
-
extracting SNPs from .vcf file along with CHR and POS
written by anamaria
-
Separating a region of the VCF file
written by mostafarafiepour
-
vcftools counts output not written to file
written by spider_monkey •
-
vcftools: unwanted filtering
written by Molly_K
-
Quickest way to filter vcf using bed files
written by spiral01
-
how to cut the gz. file by chromosome
written by chengyi31000 •
-
Vcftools: Filtering By Multiple Regions (--Positions Flag?)
written by Matt W •
Switch to bcftools. It's faster and has a LOT of extended options.
Can I extract more than one region at once?
Can I extract more than one region at once?
Hi someler ,
This reply should be a reply to my comment. Could you make the appropriate change please? That would involve the following steps:
Select All->Copythere).moderateback in your answer here: A: Extracting multiple regions of interest with vcftoolsDelete PostSubmitbutton.Thank you!
PS: In the future, do not add answers unless you're answering the top level question.