Rna-Seq Tag Count
4
1
Entering edit mode
9.6 years ago
Rna-Seq ▴ 50

Hi,

I've RNA-seq illumina data and i would like to find differentially expressed genes,I mapped the data and now i should write a script for tag counting. I wonder is there a ready made function for tag counting that i can use or i have to write my own script..

Thank you in advance and Best Regards.

rna • 3.7k views
ADD COMMENT
3
Entering edit mode
9.6 years ago
Vitis ★ 2.4k

HTSeq package does that, it's based on python. It takes a gff3 as part of the input to define tags. Also, R package IRanges has a function called countOverlap, which also does the same thing, counting reads mapped to defined tags. Here, you may use a gff3 as input, convert it into a Granges object and count overlapped reads on top of that. I'd recommend the HTSeq solution if your bam/sam files are very large, as R tends to take the entire alignment at once. If you don't have a huge memory, you're out of luck.

ADD COMMENT
0
Entering edit mode

I think HTSeq-count is a superior soultion to countOverlap because it has options to deal different ambiguities scenarios one would have to code oneself.

ADD REPLY
2
Entering edit mode
9.6 years ago
Marina Manrique ★ 1.3k

Have you tried Cufflinks? It's pretty straightforward

http://cufflinks.cbcb.umd.edu/

ADD COMMENT
0
Entering edit mode
9.6 years ago
Travis ★ 2.8k

You could also look up RUM (RNASeq Unified Mapper) as a complete mapping/counting solution for Illumina data.

ADD COMMENT
0
Entering edit mode

edgeR, DESeq, NOISeq don't do tag counting, but take as input a tag counted file.

ADD REPLY
0
Entering edit mode

Comment edited, thanks.

ADD REPLY
0
Entering edit mode
9.6 years ago

coverageBed from bedtools can do it too.

ADD COMMENT

Login before adding your answer.

Traffic: 1174 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