Mapped Reads Into A Gene Level Count
3
0
Entering edit mode
11.6 years ago
joaslucas ▴ 90

Hi everyone, This is a couple of simple questions

After mapping reads (RNA-seq) against a reference bacterial genome I need to know how many reads I got for each gene. How do I do that, which program should I use? And finally, big genes will have much more reads than smaller genes, and that cannot be accounted for DE, Is this problem solved in the normalization step?

Thanks, I am a biologist tying to understand RNA-seq. I ve read some papers but I still far from understanding the whole process.

Thanks for the help.

Thanks.

rna-seq counts • 3.4k views
ADD COMMENT
1
Entering edit mode
11.6 years ago

A simple strategy:

Make a .bed of genes coordinates

Use BEDTools to find out how many reads hit each gene

Correct for gene length

Correct for differing # of reads between samples

ADD COMMENT
0
Entering edit mode

Sorry, from where can I get the gene coordinates? I have a bam file then how I can have a gene coordinates?

ADD REPLY
0
Entering edit mode
ADD COMMENT
0
Entering edit mode
11.6 years ago
Asaf 10k

Another great place to start is in this Galaxy exercise: https://main.g2.bx.psu.edu/u/jeremy/p/galaxy-rna-seq-analysis-exercise Since you're dealing with bacteria you should probably use Bowtie and not Tophat but I think it won't have any effect.

ADD COMMENT

Login before adding your answer.

Traffic: 1040 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6