Question: coverage of the reads over coding sequence only for one gene
0
gravatar for alireza346
10 months ago by
alireza3460
alireza3460 wrote:

I have RNAseq data and trying to get the coverage of the reads over coding sequence but only for one gene for example GAPDH. and plot the coverage in this way: in fact the y-axis would be coverage or normalized read count and x-axis would be the coding sequence. but the coding sequence which would be the x-axis will be shown between 0 to 100 (%). I have made the bam file for gapdh. do you know how to do that?

rna-seq • 296 views
ADD COMMENTlink written 10 months ago by alireza3460

You can use bedtools coverage to compute the coverage at every base of your gene, then use R to output a graph of coverage over the length of your gene. (I don't really get the "0 to 100%" part)

ADD REPLYlink written 10 months ago by guillaume.rbt560

I read it as the % is relative to everything, including non-coding regions of the gene.

ADD REPLYlink written 10 months ago by Eric Lim1.4k

try coveragebed -hist option from bedtools.

ADD REPLYlink modified 10 months ago • written 10 months ago by cpad011211k
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: 1171 users visited in the last hour