SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
HTSeq - paired reads count nofar Bioinformatics 0 10-02-2014 04:02 AM
HTseq-count : how to get unique reads? ThePresident Bioinformatics 11 08-19-2013 04:11 PM
[HTSeq-Count] How are reads counted? syintel87 Bioinformatics 0 06-17-2013 01:18 PM
Reads Count - HTSeq alternative? NGSAwesome RNA Sequencing 1 10-15-2012 09:13 AM
htseq-count gets more reads? deepsea Bioinformatics 3 03-29-2012 11:27 AM

Reply
 
Thread Tools
Old 11-12-2014, 05:05 AM   #1
Fernas
Member
 
Location: Middle East

Join Date: Apr 2013
Posts: 74
Default Ambiguous Reads Treatment by HTSeq-Count

Currently, if a person uses HTSeq-Count to count reads on transcript level and the tool finds a read that overlaps two overlapped transcripts (this is very common) then the tool will identify the read as ambiguous and will not include it in the counting process. As the author assumed, this is useful for identification of differentially expressed genes. However, this approach will underestimate the expression of overlapping transcripts and overestimate unique ones.

Can we modify somehow HTSeq-Count to include such reads? any other tools that overcome such situation?
Fernas is offline   Reply With Quote
Old 11-12-2014, 05:13 AM   #2
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,480
Default

htseq-count isn't intended for that use, so while you could modify it, why both? You can get featureCounts to do this, but you really need to be careful about what you then do with these counts. If you're going to use these in something like DESeq2 or edgeR then you're simply doing it wrong (you'll just simply get wrong results).

In general, transcript abundance is better estimated with an expectation maximization algorithm (e.g. FluxCapacitor, or eXpress).
dpryan is offline   Reply With Quote
Old 11-12-2014, 05:18 AM   #3
Fernas
Member
 
Location: Middle East

Join Date: Apr 2013
Posts: 74
Default

Thank you very much @dpryan for the very informative reply. I will have a look into these tools.
Fernas is offline   Reply With Quote
Reply

Tags
htseq-count

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 02:28 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2019, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO