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?
Not a solution, rather recommendation - better use featurecounts, faster and more flexible.
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.How can we correct this? Because the SAM file was a output of SHRiMP algorithm.