alexdobin / STAR

RNA-seq aligner
MIT License
1.77k stars 495 forks source link

Issue aligning PacBio reads with STARlong #2131

Open diazdc opened 2 months ago

diazdc commented 2 months ago

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