Question: How to get the no-overlapping exome region (.bed file)?
0
gravatar for weiyuchung
3.2 years ago by
weiyuchung0
United States
weiyuchung0 wrote:

I want to get the hg38 whole exome regions (.bed) for the downstream analysis. However, some regions in the .bed file browsed from the UCSC table are overlapped. How to get the no-overlapping exome region .bed file?

My settings for the UCSC table browser as follows:

  • clade: mammal
  • genome: human
  • assembly: Dec.2013 (GRCh38/hg38)
  • group: Genes and Gens Predictions
  • track: GENCODE v24
  • table: knownGene
  • position: chr22:1-50818468
  • output format: BED

Click "Get Output", select “Exons plus 0 bases at each end”, and get the .bed file.

In this .bed file, the regions from row 29:15690025-15690709 and row 40:15690077-15690314 are overlapped. There are some other overlapped regions in this file. How to get the .bed file with no-overlapping exome regions? BEDtools?

Thank you!!!

ucsc hg38 exon bed exome • 1.8k views
ADD COMMENTlink modified 3.2 years ago by Alex Reynolds29k • written 3.2 years ago by weiyuchung0
0
gravatar for Alex Reynolds
3.2 years ago by
Alex Reynolds29k
Seattle, WA USA
Alex Reynolds29k wrote:

With BEDOPS bedops -n:

$ bedops -n 1 regions.bed regions.bed > non-overlapping-regions.bed

What comes out of UCSC is sorted correctly, usually, so I think you can skip presorting.

ADD COMMENTlink modified 3.2 years ago • written 3.2 years ago by Alex Reynolds29k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1945 users visited in the last hour