Question: DEXseq count reads tool obtaining only zero counts for all exons
gravatar for f.carlisle
10 weeks 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 10 weeks ago by Jennifer Hillman Jackson24k • written 10 weeks ago by f.carlisle0
gravatar for Jennifer Hillman Jackson
10 weeks ago by
United States
Jennifer Hillman Jackson24k 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 10 weeks ago by Jennifer Hillman Jackson24k

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 10 weeks 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: 81 users visited in the last hour