Question: Map RNA-Seq reads in 3' UTR region
0
gravatar for luca
3 months ago by
luca0
luca0 wrote:

Hi folks, I just have a quick question which hopefully has a simple answer. I downloaded from NCBI a set of fastq files containing mouse RNA-seq data which was generated using "classic" TruSeq RNA-library preparation kit. Hence it contains reads from the whole mRNA molecule. I want to align (and then count) reads that do map only the 3' UTR region. I am currently using STAR to map the reads and the star index was built with the latest Ensembl .gtf and .fasta files.

Any suggestion will be much appreciated Thanks Luca

rna-seq alignment • 186 views
ADD COMMENTlink written 3 months ago by luca0
0
gravatar for ATpoint
3 months ago by
ATpoint25k
Germany
ATpoint25k wrote:

The 3'UTRs should be annotated in the reference GTF/GFF3 file which you can filter for the coordinates. Once you have the coordinates use something like samtools view with the aligned bam file:

samtools view -L BEDfile_with_3UTRcoordinates.bed -o data_3UTR_reads.bam aligned.bam
ADD COMMENTlink modified 3 months ago • written 3 months ago by ATpoint25k

Thanks ATpoint for your reply. How do I extract the coordinates of UTR regions from the GTF file?

ADD REPLYlink modified 3 months ago • written 3 months ago by luca0

See for general inspiration: A: how to get intronic and intergenic sequences based on gff file?

ADD REPLYlink written 3 months ago by ATpoint25k

Thanks! I will have a look at it!

ADD REPLYlink written 3 months ago by luca0
Please log in to add an answer.

Help
Access

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