SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
htseq-count paolo.kunder Bioinformatics 10 10-22-2014 04:45 AM
htseq-count on UTR emilyjia2000 Bioinformatics 9 04-06-2012 12:13 PM
htseq-count gets more reads? deepsea Bioinformatics 3 03-29-2012 11:27 AM
htseq-count output Palgrave Bioinformatics 7 03-05-2012 07:04 AM
htseq-count with warning for every read to represent all of zero counts in output hibachings2013 RNA Sequencing 10 07-15-2011 10:19 AM

Reply
 
Thread Tools
Old 06-13-2012, 04:10 AM   #1
NicoBxl
not just another member
 
Location: Belgium

Join Date: Aug 2010
Posts: 264
Default multiBamCov or htseq-count to count read per feature ?

Hi,

I'm wondering what is the best method to extract the number of reads for each feature in a gtf (or gff, bed,...) file. I tried htseq-count and multiBamCov but they gave me different results.it seems that multiBam count all the reads (complete and partial aligned) associated with each exon. It means there are many reads are count twice or more time.

After doing DE analysis (DESeq) on both read count matrix (one from htseq, one from multiBamCov), the results are quite surprising.

pval adjusted < 0.05 multiBamCov : 123 gene differentially expressed htseq : 880 gene

Intersection : 118 gene

So which one to use ? is it possible to change multiBamCov to be more strict ? maybe is it possible to use other tools from bedtools ?

Thanks,

N.
NicoBxl is offline   Reply With Quote
Old 07-03-2012, 02:05 AM   #2
dawe
Senior Member
 
Location: 4530'25.22"N / 915'53.00"E

Join Date: Apr 2009
Posts: 258
Default

AFAIK htseq-count supports several ways to perform the intersection, hence you may want to post the actual command lines with all parameters used. multiBamCov supports flag filtering.
What do you mean by strictness here?
dawe 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 02:50 PM.


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