Bed file of union of all transcripts for a set of genes
1
0
Entering edit mode
5.2 years ago
Floydian_slip ▴ 170

Hi, I have a set of 100 genes and I want to obtain a bed file containing a union of all transcripts for each of those genes with hg19. IOW, for each gene, I want a bed file that contains the union of all reported transcripts for that gene. What is the fastest/best way to do so? Thanks!

transcripts gene exons • 1.0k views
ADD COMMENT
0
Entering edit mode
ADD REPLY
0
Entering edit mode

Example input: a list of genes such as BRCA1 and ALDH1. Desired output: a bed file containing their transcripts. I understand there will be multiple lines for each for each of the transcript regions. Thanks!

ADD REPLY
0
Entering edit mode

I was referring to pasting some real data examples, like a table that shows how things should look, not just words. You already explained with words in your question, but no coded solution yet because people do not know how exactly the data look like.

ADD REPLY
0
Entering edit mode
5.2 years ago
michael.ante ★ 3.9k

I'd grep a gtf with the gene list and use bedtools merge on the result.

ADD COMMENT

Login before adding your answer.

Traffic: 2518 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

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

Powered by the version 2.3.6