coverage of the reads over coding sequence only for one gene
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?
• 1,609 views
•
link
written
by
alireza346 •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Changing length of x axis on ggplot
written by salman_96 •
-
How to create coverage graph/plot for an entire Illumina run
written by SaltedPork
-
Plotting .tab/.bam file
written by zack.henning •
-
dispersionPlot in cummeRbund package
written by manjumoorthy95 •
-
Plot line graph when counts are not automatically given
written by kristina.mahan
-
making the gene specific bam file
written by alireza346 •
-
How to have frequency displayed on the Y axis in weblogo
written by Aishwarya Kulkarni
-
How to interpret and trim this plots?
written by eyonesi •
-
How to plot motif footprint?
written by tamadomi •
-
out put figure of a NGS tool
written by kanwarjag
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)
I read it as the % is relative to everything, including non-coding regions of the gene.
try coveragebed -hist option from bedtools.