Question: DEXseq count reads tool obtaining only zero counts for all exons
gravatar for f.carlisle
5 months ago by
f.carlisle0 wrote:

I'd like to run DEXseq on my RNA-seq dataset to identify any potential differences in exon inclusion between samples. I used HISAT2 to generate BAM files of aligned reads, and then DEXseq-Count mode-prepare to generate an annotation GFF file. My reference genome for alignments was hg38, and I used the ENSEMBL hg38 gene annotation GTF to prepare my GFF file. However, whenever I try to run DEXseq-Count in count mode, I obtain zero counts for all exons. I have separately run transcript-level analyses on these BAM files, so I know they're okay and I should be getting non-zero exon counts. Any ideas what I might be doing wrong?

Thanks for your help.

ADD COMMENTlink modified 5 months ago by Jennifer Hillman Jackson25k • written 5 months ago by f.carlisle0
gravatar for Jennifer Hillman Jackson
5 months ago by
United States
Jennifer Hillman Jackson25k wrote:


Compare your analysis methods to those in our tutorial utilizing the DEXseq tool and it should reveal where the problem is:

A data mismatch could also be a factor (mismatched chromosome identifiers between the inputs) - even if the prior transcript-level analysis was successful. How to check:

For reference, should you need more help:

Thanks! Jen, Galaxy team

ADD COMMENTlink written 5 months ago by Jennifer Hillman Jackson25k

Hi Jen, Thanks for the suggestions. My analysis methods are identical, so it looks like a data mismatch may be the problem.

ADD REPLYlink written 5 months ago by f.carlisle0
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 115 users visited in the last hour