FastQ: quality string length is not equal to sequence length
0
0
Entering edit mode
8 weeks ago
yvonnehhe • 0

I'm using starlong to map fastq reads according to the index (SB1003index2) created with star.

My code:

STARlong --quantMode GeneCounts --genomeDir SB1003index2 \ --readFilesIn all_DE.fastq \ --outFileNamePrefix DE_ \ --outSAMtype BAM SortedByCoordinate

My error: enter image description here

I checked and the sequence and the quality scores have the same length (same number of characters).

I googled the error message, and someone suggested to remove all the run information after the first part of the identifier. I tried that but I'm still receiving the same error.

Any advice would be appreciated!

star fastq starlong • 305 views
ADD COMMENT
0
Entering edit mode

When I've had this error with STAR running seqkit sana resolved it, although this was for short read.

ADD REPLY
0
Entering edit mode

Thanks for the reply, I used seqkit sana but all the lines passed and it discarded 0 lines.

ADD REPLY
0
Entering edit mode

If it can rescue a read it will do so without discarding it, so try realigning now to see if it worked.

ADD REPLY
0
Entering edit mode

Well, what does that entry look like in the fastq?

ADD REPLY
0
Entering edit mode

Can you try aligning the data using minimap2? This may be a bug in STARlong.

ADD REPLY

Login before adding your answer.

Traffic: 1297 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