Error in HTSeq
0
0
Entering edit mode
6.1 years ago
glady ▴ 320

Hello, I have performed the mapping of miRNA reads to the mirbase precursor file by SHRiMP. But I'm facing a probem in HTSeq while taking the read counts from the output.sam file.

The error which I'm trying to solve is;

Error occured when processing SAM input (line 2305 of file HBI_10.sam): ("'seq' and 'qualstr' do not have the same length.", 'line 2305 of file HBI_10.sam') [Exception type: ValueError, raised in _HTSeq.pyx:808]

Has anyone of you come across this error? What might be the solution?

RNA-Seq • 1.4k views
ADD COMMENT
0
Entering edit mode

Not a solution, rather recommendation - better use featurecounts, faster and more flexible.

ADD REPLY
0
Entering edit mode

Somehow read on line 2305 is mangled (if the error message is correct) such that it has sequence and quality strings that are of different length. You could backtrack and see where/what step of the processing that happened.

ADD REPLY
0
Entering edit mode

How can we correct this? Because the SAM file was a output of SHRiMP algorithm.

ADD REPLY

Login before adding your answer.

Traffic: 2541 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6