This is a test version of Biostars. For the public version, visit https://www.biostars.org.
CDS read counts per gene

Hi Guys,

I have aligned RNA-seq data and have bam files. I am goint to coundt the number of reads per gene but I need only the number of reads per gene only for coding sequence part of genes. do you guys know how to get those counts? tHANKS

rna-seq

Dear Behzad, Hi

Does this post can offer any helps to you ?

And I guess if you have GTF file there would be some transcripts coordinates such as CDS in it.

~Best

If I understand correctly, he needs counts, not reads. This is exactly what featureCounts or htseq-counts does. "Counting reads per annotation element".

Dear @Wouter, Hi.

Do you mean this and this ?

Yes, as such he can "count the number of reads per gene" -> as he asked for ;-)

0 answers

No answers yet.

Log in to answer this question.