DEXSeq python script error message

[shilez]

When I try to run dexseq_prepare_annotation.py, I got a python error at line 29:
exons[f.iv] += ( f.attr['gene_id'], f.attr['transcript_id'] )
TypeError: unsupported operand type(s) for +=: 'StepVector' and 'tuple'
Anyone has any idea how to fix it?

[Jane M]

I have no idea, but I managed to run this script, so there is no error in it. Nevertheless, if I remeber well, a certain version of python is needed...

[Simon Anders]

Could it be that you have an outdated version of HTSeq installed?

[shilez]

Thank you, guys.
Once I update the HTSeq, it works fine.
However, I do got the error message "Same name found on two chromosomes".
Does it mean we need to remove the duplicated genes from the annotation gtf file first before running dexseq_prepare_annotation.py?

*and also an error from HTSeq
ValueError: start is larger than end

[areyes]

The first error can be solved by adding an additional ID tag to these genes. (e.g. if the gene is called A, you could write A_dup1 and A_dup2 foreach one so that the gene names become unique). They are usually not to many in the annotation files. (I have done it manually for the ensembl mouse gtf annotation).

The second I have also seen it, for some genes, the start value is higher to the end value. I guess this is a mistake in the gtf files. Can be corrected just by shifting the columns in this cases. Also, there are usually not to many.