Hello,
I have a VCF file with SNPs and genes in GFF file. I would like to extract only SNPs that reside on the genes or flanking +/- 5kb. I could artificially add/substract 5kb to gff file and then follow the answer in this link: https://www.biostars.org/p/333734/), but is there an option to do so? Also, is it possible to append the gene name and positions to the vcf file?
•
written
by
User000
1 answer
I'd probably just add the 5kb to the gff via bedtools slop, then intersect with bedtools intersect -wa -wb, and you'd basically have everything you want.
Bedtools doesn't preserve headers by default. You can try using the --header option (but it has to be in their specific format), or you can pipe the header to your output file first, then append your results to bedtools to that file. Or just open and tack it on afterwards, which is often easiest if you want to keep the headers from both files after an intersect.
Log in to answer this question.
-
How to download a file from zenodo link without leaving terminal?
written by SlowSD •
-
bcftools query add header
written by User000
-
How does one harmonize WGS data with Microarray data?
written by Moni •
-
Bedtools intersect/bcftools annotate dublicate the SNPs
written by User000
-
Mapping SNP position to gff file
written by shachigahoimbi •
-
How to extract genomic positions of SNPs using snp-sites?
written by pltbiotech_tkarthi
-
Extract SNPs flanking sequences based on VCF and genome Fasta files
written by Denis
-
Extract SNPs from VCFfile located in genes based on GFF file information
written by Denis
-
SNP annotation using gff3
written by reza
-
SNP file preparation for genotyping
written by hpapoli
