View Single Post
Old 08-05-2010, 11:52 PM   #1
Junior Member
Location: Vienna

Join Date: Jul 2010
Posts: 2
Default DESeq: Read counts vs. BP counts

Hi Simon,

is it possible to use raw bp counts (how often a bp is sequenced within a region) of an exon/transcript for differential expression analysis with DESeq?
First I think, that if we consider the expression of an individual exon the bp count would avoid double counting of split reads (especially if it belongs to a major part to one exon). It would only assign a portion of the read to one exon. Second, it would be nice to use with cufflinks where cov*length = bp count. Also for cufflinks the bp count divided by the read length is not necessarily equal the read count, especially for alternative transcripts or if you provide a gtf of known genes (even though its an approximation).

Thanks a lot and best regards,
burkard is offline   Reply With Quote