This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Confirm if a genome position falls into a coding exon

Good morning,

I do have few human positions of interest and I want to check if they fall inside a coding exon. I know genome browser (https://genome.ucsc.edu/cgi-bin/hgGateway) is very useful for visualization, but I am more interesting in identify different overlaps. \Would it be possible to download it from somewhere? Thank you!

exon ucsc coding
• 1,241 views • link
updated by Pierre Lindenbaum Pierre Lindenbaum • written by Lila M Lila M

1 answer

get a gtf, convert to bed, sort and merge with bcftools merge. e.g: awk -F '\t' '($3=="exon") {printf("%s\t%d\t%d\n",$1,int($4)-1,$5);}' in.gtf

put your position(s) in a bed and sort

use bedtools intersect.

• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

Thanks both. Will be $3=="exon" enough to select only the 'coding' exons (and to exclude any non-coding exon)?

• 0 views • link
written by Lila M Lila M

Log in to answer this question.

More posts like this