Biostar Beta. Not for public use.
Question: counting the reads that mapped to a part of transcriptome
0
Entering edit mode

I have RNAseq data and trying to count the reads that mapped only to a part of 5'UTR (in fact the whole 5'UTR except the first 50 nts). do you know how to do that?

ADD COMMENTlink 2.4 years ago Sara • 30 • updated 2.4 years ago glihm • 600
0
Entering edit mode

Hi Sara,

Did you align the raw reads to a reference genome? if so what is the current format of your aligned reads?

Tarek

ADD COMMENTlink 2.4 years ago tarek.mohamed • 250
Entering edit mode
0

I can actually align to the part of transcriptome that I want (first I have to make fasta file and use it as reference file to align to). I can also make gtf file using the same file and count using that. but in both cases it returns 0 counts. so I thought maybe I should align to the genome and count on the reads that mapped to that part. but here I do not know how to do that.

ADD REPLYlink 2.4 years ago
Sara
• 30
0
Entering edit mode

Hello Sara,

once you have mapped the data against the genome of your interest, you can simply use FeatureCounts program for instance:

featureCounts -a yourGTF.gtf -t UTR5 -g transcript_id -o counts_UTR5_transcripts.txt

The -t option allows you to select a type of features to be counted. The -g option allows you to choose "how to group your data". Using transcript_id (or something equivalent in your GTF/GFF), you tell to FeatureCounts: "ouput the counts for each transcript".

ADD COMMENTlink 2.4 years ago glihm • 600
Entering edit mode
0

Also if you prefer to use R, you can use "Rsubread" package to count your reads using featureCounts () as mention by glihm.

T

ADD REPLYlink 2.4 years ago
tarek.mohamed
• 250

Login before adding your answer.

Similar Posts
Loading Similar Posts
Powered by the version 2.0