Jump to content

  • Log in with Facebook Log in with Twitter Log in with Windows Live Log In with Google      Sign In   
  • Create Account

Submit your paper to J Biol Methods today!
- - - - -

RNA seq reads - assessing distance from start codon

  • Please log in to reply
1 reply to this topic

#1 Ahrenhase



  • Awaiting Authorisation
  • PipPipPipPipPipPipPipPipPipPip
  • 196 posts

Posted 02 December 2013 - 08:24 PM

I have RNAseq data and I'm trying to relate the reads to distance from the start codon.  Anyone have experience with this. 

#2 Mikhail Shugay

Mikhail Shugay


  • Active Members
  • Pip
  • 6 posts

Posted 09 December 2013 - 07:20 AM

Do you need to extract reads that are e.g. less than 1kb from start codon? You should first map your reads (with Tophat for example) and get a sam/bam file. Then should try bed tools (see doc here http://bedtools.read...ools-suite.html) to intersect your sam/bam with bed file that contains regions near the start codon. To get this bed file you need to do some manual work: get RefSeq gene track from UCSC Gb, find the starts of CDS (dont forget about gene strand), extend them to 1kb regions (upstream and downstream) and convert to bed format.

On the other hand you could run through all mapped reads in sam/bam using Picard Java API and process them as desired.

Home - About - Terms of Service - Privacy - Contact Us

©1999-2013 Protocol Online, All rights reserved.