Hi all,
I have a vcf file that I perform the BED file by vcftools to subset the mutations in my interest genes. I have downloaded the BED file from the UCSC (selecting the BED file for the whole gene of my gene list). However, after performing the vcftools, I have found some mutations in genes that weren't in my gene list. How to subset the mutations that occur just in my interest genes?
Thanks for any help.
1 answer
Extract Sub-Set Of Regions From Vcf File
Use tabix, http://www.htslib.org/doc/tabix.html, and see its -R option to only retrieve VCF entries that overlap the BED file you provide for that argument. That would need to be a BED file with the coordinates you are interested in, e.g. the entire gene interval, or its exons. You could get the coordinates directly from a reference GTF file.
Log in to answer this question.
-
How to extract mutations in a particular gene from the VCF file of whole genome
written by saamar.rajput
-
Does PCAWG have any advantage compared to TCGA for studying protein-coding regions?
written by Zahra •
-
Is the promoter part of whole of gene in the UCSC database?
written by Zahra •
-
How to find the position of genes that their variants called by hs37d5?
written by Zahra •
-
A question about BED file and finding mutation in vcf file
written by Zahra •
-
How is better perform the analyze the somatic mutations? (the mutations of my interest gene)
written by Sonia •
-
Can non-driver mutations disrupt the function of proteins and further affect a cellular pathway?
written by Zahra •
-
Study variation of some genes by VCF file
written by Zahra •
-
BAM file download for just a few genes
written by Zahra •
-
intersecting bed and vcf files for polyploid samples
written by hpapoli
