SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
bowtie number of mismatches and multiple aligned reads carolW Bioinformatics 17 12-23-2014 05:03 PM
Ideas on how to count the number of specific amplicons in aligned amplicon data swNGS Bioinformatics 6 10-15-2012 11:57 PM
Does Picard remove reads aligned to multiple locations? Heisman Bioinformatics 0 04-05-2012 10:14 PM
How to count aligned RNA-seq reads after sequenced and aligned by Illumina? IceWater Illumina/Solexa 5 04-05-2012 09:18 AM
count A/T/C/G/N in a aligned site johnsequence Bioinformatics 6 11-17-2011 12:02 PM

Reply
 
Thread Tools
Old 01-27-2015, 12:44 AM   #1
a.kmg
Member
 
Location: Internet

Join Date: Aug 2014
Posts: 15
Default A great tool to count multiple aligned reads.

Dear members of seqanswers,

Do you know a tool for counting multiple aligned reads?

Currently I use the htseq-count tool, but it doesn't take into account this type of read.

I know cufflinks allows this kind of thing, but I do not like this tool, because for me it is a "black box".

There are some time ERANGE allowed to do this assigning multiple aligned reads in the proportion of single aligned reads in each respective region. I find it effective, unfortunately ERANGE is no longer maintained.

I thank you in advance for your help.

Best,
a.kmg
a.kmg is offline   Reply With Quote
Old 01-27-2015, 01:00 AM   #2
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,480
Default

The appropriate tool will depend on how you want the reads counted. featureCounts is one option, though obviously if you're counting over genes and using those integer counts to look at differential expression then this is a bad idea. Other options would be things like eXpress or any other tool meant to generate expected counts using a expectation maximization method. This will yield fractional counts, which is appropriate for multimapped reads.

Last edited by dpryan; 01-27-2015 at 01:08 AM.
dpryan is offline   Reply With Quote
Old 02-20-2015, 07:07 AM   #3
a.kmg
Member
 
Location: Internet

Join Date: Aug 2014
Posts: 15
Default

Thank you dpryan.

I currently testing Rcount, a new tool where alignments of multireads are weighted based on the number of uniquely aligned reads in the neighborhood.

And I tested featureCount, that actually allows me to treat multiple aligned reads, but the fact of counting the reads for each chromosomal locations on which it aligns leads to a biais, right?
a.kmg is offline   Reply With Quote
Old 02-22-2015, 05:05 AM   #4
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,480
Default

Quote:
Originally Posted by a.kmg View Post
And I tested featureCount, that actually allows me to treat multiple aligned reads, but the fact of counting the reads for each chromosomal locations on which it aligns leads to a biais, right?
Not really. The reads arose from the neighboring reads, so the counting would accurately represent the biological reality.
dpryan is offline   Reply With Quote
Reply

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 10:04 PM.


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