Calculate the coverage of genes

Question

[Alex]

Hi,everyone , There I want to calculate the RNAseq reads coverage the numbers of genes and sequences from reference genome and I have bam files ,did anyone knows how to calculate it ? Thanks Alex

Answer 1

[Devon Ryan]

It's a bit hard to understand exactly what you want, but I assume you mean you want to find the number of alignments per gene. For that, use featureCounts from the subRead package.