Question: Mapping aligned reads to UTRs
gravatar for jpcchoy
11 months ago by
jpcchoy0 wrote:

Hi everyone,

I'm trying to work out which of my aligned reads in a sorted BAM file map to UTRs.
I have a text file with transcript and CDS coordinates as follows:

name chrom strand txStart txEnd cdsStart cdsEnd
ENST00000371007.6 chr1 - 67092164 67231852 67093004 67127240
ENST00000371006.5 chr1 - 67092175 67127261 67093004 67127240
ENST00000475209.6 chr1 - 67092175 67127261 67093579 67127240

My reads have been aligned to hg38, and the features above are also from hg38.
For each feature (transcript) I think that the region txStart - cdsStart is the 5' UTR, and the region cdsEnd - txEnd is the 3' UTR, but am not sure how to do this mapping.

Does anyone know any packages that might serve this purpose?
Please let me know! Thanks so much.

rna-seq utr alignment cds • 261 views
ADD COMMENTlink modified 11 months ago by Ace70 • written 11 months ago by jpcchoy0
gravatar for Pierre Lindenbaum
11 months ago by
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum128k wrote:

convert the transcript file to a BED file using awk and then use samtools view with option -L

ADD COMMENTlink written 11 months ago by Pierre Lindenbaum128k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1662 users visited in the last hour