Counting reads in CDS - but skipping few first/last codons
0
0
Entering edit mode
6.3 years ago
PocketChange ▴ 10

I've been struggling with this for some time now, and finally decided to ask you guys for advice: I'm trying to analyze some ribosome profiling data, and would like to count reads (e.g. using HT-Seq) in cds, but at the same time I'd like to exclude several first and last codons (I think in general people skip first 15 and last 5 codons). I suppose I could potentially achieve this by modifying my gtf? Are there any tools you could recommend I use to get this done? Or do you have any other ideas/suggestions on how I start with this? Thank you for your help.

rna-seq htseq • 1.5k views
ADD COMMENT
0
Entering edit mode

I think bedtools can be used to change the gtf.

ADD REPLY
0
Entering edit mode

More specifically: bedtools slop is able to increase the size of features. However, using negative integers for -l and -r shrinks features at the desired end.

ADD REPLY
0
Entering edit mode

Thanks! So in this case proper way would be to extract CDS exons containing start codons, resize, and then merge back? Thanks again!

ADD REPLY

Login before adding your answer.

Traffic: 2011 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6