Error in HTSeq
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?
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
HTseq-count error
written by sanatbhadsavle •
-
htseq gene atrribute error
written by Balatheskulter •
-
Error with HTseq RNAseq read count
written by Bioinfonext
-
HTseq code error repeatedly
written by bnayer26 •
-
Aligning miRNA reads to miRbase
written by glady
-
HTSeq error when reading BAM file
written by Mike
-
HTSeq analysis for Illumina tissue samples?
written by pbio
-
Htseq-count error "start is larger than end"
written by ejmcmaho •
-
genes count in sam file using htseq-count
written by bolbolman2000 •
-
Why Does The Htseq Program Indicate A Malformed Sam File With Error Mrnm==*
written by joaslucas •
Not a solution, rather recommendation - better use featurecounts, faster and more flexible.
Somehow read on line
2305is 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.