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
do you know how to do that?