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

Hello,

Compare your analysis methods to those in our tutorial utilizing the DEXseq tool and it should reveal where the problem is: https://galaxyproject.github.io/training-material/topics/transcriptomics/tutorials/ref-based/tutorial.html#inference-of-the-differential-exon-usage

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: https://galaxyproject.org/support/chrom-identifiers/

For reference, should you need more help:

Thanks! Jen, Galaxy team

ADD COMMENTlink written 11 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 11 months ago by f.carlisle0
Please log in to add an answer.

Help
Access

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