alexdobin / STAR

RNA-seq aligner

Geek Repo:Geek Repo

Github PK Tool:Github PK Tool

Issue aligning PacBio reads with STARlong

diazdc opened this issue · comments

Good day,

I have PacBio library that I would like to run using STARlong. However, when attempting to run STARlong I get an error immediately that states the quality score does not match the sequence length.

EXITING because of FATAL ERROR in reads input: quality string length is not equal to sequence length
@m84039_240419_200920_s3/249238533/ccs/18_1224
...(Some 650bp sequence)...

Not surprisingly, the error was thrown for the very first sequence in the file, so I took a look at the number of characters in both the sequence and quality scores for the fastq record listed in the error. Both had exactly 1076 characters. Given the error only listed 650 bases, is there a max limit for read length that STARlong can handle? I searched the manual and Googled around for an answer, but couldn't find any information on this. I can also share a sample of the data.

Thanks