DESeq2 error using HTseq count

1

Hello!

I am trying to run DESeq2 using HTseq count as input but I am running in the following error every time, I already checked and input data is correct.

Error in Ops.factor(a$V1, l[[1]]$V1) :
level sets of factors are different
Calls: DESeqDataSetFromHTSeqCount -> sapply -> sapply -> lapply -> FUN -> Ops.factor
Warning message:
In is.na(e1) | is.na(e2) :
longer object length is not a multiple of shorter object length

I have 4 replicates per treatment (two treatments) and they have different number of lines (outpu htseq) but I think this is normal since htseq counts the number of exons and they could be different among replicates.

I read that the problem could be the gff3 format but the annotation of the genome is not in gtf format. Is there a tool in Galaxy to do the conversion?.

Thank you very much in advance,
Maria

2

Hi,
I have the same problem with HTseq outputs and DESeq2: my error message said ‘’ a problem with inputs’’. Please In addition of HTseq outputs, factor and factor levels, do I need to provide other information to run DESeq2 in galaxy?
Just in case, you can use gffread Filters and/or converts GFF3/GTF2 records from Cufflinks for converting GFF3 into GTF.

Thank you,
Ferdi