multiBamCov or htseq-count to count read per feature ?

NicoBxl · 06-13-2012, 05:10 AM

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.

dawe · 07-03-2012, 03:05 AM

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?

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

06-13-2012, 05:10 AM	#1
NicoBxl not just another member Location: Belgium Join Date: Aug 2010 Posts: 264	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.

07-03-2012, 03:05 AM	#2
dawe Senior Member Location: 45°30'25.22"N / 9°15'53.00"E Join Date: Apr 2009 Posts: 258	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?