Question: how to get Homo_sapiens.GRCh38.84.bed file
1
gravatar for yuxinghai
20 months ago by
yuxinghai10
yuxinghai10 wrote:

I'm using geneBody_coverage2.py in RseQC package to calculate Read coverage over gene body, which need a bed12 file. I did alignment based on Homo_sapiens.GRCh38.84.gtf , so i need a Homo_sapiens.GRCh38.84.bed. But gtfToBed can't produce a BED12 file, how can I get it? anyone can help me?

annotiation rna-seq gene • 1.0k views
ADD COMMENTlink modified 20 months ago by Devon Ryan85k • written 20 months ago by yuxinghai10
8
gravatar for Devon Ryan
20 months ago by
Devon Ryan85k
Freiburg, Germany
Devon Ryan85k wrote:

awk '{if($3 != "gene") print $0}' Homo_sapiens.GRCh38.84.gtf | grep -v "^#" | gtfToGenePred /dev/stdin /dev/stdout | genePredToBed stdin Homo_sapiens.GRCh38.84.bed

You can download the two tools from UCSC.

If you're curious, this is from here and is essentially the code I use to generate the reference data for our institute.

ADD COMMENTlink modified 20 months ago • written 20 months ago by Devon Ryan85k
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: 1687 users visited in the last hour