Genes near SNP within certain window frame and distance
Entering edit mode
10 days ago

Hi everyone, I have SNPs and genes in bed format in two files, I can see the genes that map to SNP within certain windoframe lets say 200Kb using the following command:

 windowBed -a SNPs.bed -b Gene.bed -w 200000 > snps.withGenesWithin.bed

This command identifies the set of genes within certain windowframe but doesn't provide a distance, whereas the following command provides the genes distance from SNP but wasn't able to find option to cannot specify the window frame.

closest-features --delim '\t' --dist SNPs.bed Genes.bed > answers.bed

Is there any command which I can use to specify the window frame as well as it estimates the genes distance from SNP.

bedops bedtools • 214 views
Entering edit mode

Not sure if I am following the question, or what are you looking after, but maybe you could just combine both commands? Something like:

windowBed -a Gene.bed -b SNPs.bed -w 200000 -u > Genes_nearby_SNPs.bed

And then:

closest-features --delim '\t' --dist SNPs.bed Genes_nearby_SNPs.bed > answers.bed


Login before adding your answer.

Traffic: 2271 users visited in the last hour
Help About
Access RSS

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

Powered by the version 2.3.6