Question: Map RNA-Seq reads in 3' UTR region
0
gravatar for luca
14 months ago by
luca20
luca20 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 • 617 views
ADD COMMENTlink written 14 months ago by luca20
0
gravatar for ATpoint
14 months ago by
ATpoint38k
Germany
ATpoint38k 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 14 months ago • written 14 months ago by ATpoint38k

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

ADD REPLYlink modified 14 months ago • written 14 months ago by luca20

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

ADD REPLYlink written 14 months ago by ATpoint38k

Thanks! I will have a look at it!

ADD REPLYlink written 14 months ago by luca20
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: 1242 users visited in the last hour