This is a test version of Biostars. For the public version, visit https://www.biostars.org.
coding sequnce read coun

Hi Guys, I asked a question yesterday but did not get my answer and I really need it. I have RNA-seq data and its bam file. I want to get the read count per gene. but the problem is that I need only the read counts in the coding sequence part (I don't need the read counts at 5' and 3' UTRs). if I use Htseq And Feature Count I would get the read count for the transcript. do you guys know how I can get the read count per gene only at the CODING SEQUENCE? thanks

rna-seq
• 1,775 views • link
updated by WouterDeCoster WouterDeCoster • written by ashkan ashkan

2 answers

You could leave only the entries for stuff you want counts on and take the rest of the lines out of your GTF (is that what you are using) file?

• 0 views • link
written by GenoMax GenoMax

You could just remove the UTRs from your gtf file... grep -v

• 0 views • link
written by WouterDeCoster WouterDeCoster

Log in to answer this question.

More posts like this