Dear all, I am analyzing paired-end stranded RNAseq data . I would be interested in using HISAT2 for the alignment, and HTseq-count to count the features of the aligned reads. When I used HTseq-count to count the reads:
htseq-count -r pos -f bam SRR3589957_sorted.bam ../reference/gencode.v26lift37.annotation.gtf > SRR3589956.count
However,I got this error :
Error occured when processing SAM input (record #69114540 in file SRR3589957_sorted.bam): my_showwarning() takes from 3 to 5 positional arguments but 6 were given [Exception type: TypeError, raised in warnings.py:99]
I converted sam file to bam and sorted the reads by pos using samtools
Does anyone has an idea of the problem with the reads aligned with HISAT2? Are they compatible with HTseq-count?
Thank you very much for any help. Jing Yu