How to solve STAR tools error : quality string length is not equal to sequence length ?
1
0
Entering edit mode
16 months ago

I have used STAR tool to align short read (150-160bp) transcript against hg38. Although, few suggested to use STARlong which resulted in increasing mapping percentage and mapped few transcripts as well. As STARlong used for long RNA reads. Therefore, It seems unfavourable to use for short reads. Kindly, suggest the possible solution or approach.

command line:

STAR --runThreadN 16 --runMode alignReads --genomeDir grch38_STARindex --readFilesCommand zcat --readFilesIn C_S30_1P.fastq.gz C_S30_2P.fastq.gz --outSAMtype BAM SortedByCoordinate --quantMode GeneCounts --outFileNamePrefix C_S30_

Output :

EXITING because of FATAL ERROR in reads input: quality string length is not equal to sequence length
    @A00573:363:HGYWYDSX7:2:2677:23728:19132
    +
    7:2:2654:20591:23594 24831472 N 0
    SOLUTION: fix your fastq file

    @A00573:363:HGYWYDSX7:2:2677:23728:19132
    +
    7:2:2654:20591:23594 24831472 N 0
    SOLUTION: fix your fastq file
STAR illumina sequencing RNAseq Alignment • 2.1k views
ADD COMMENT
0
Entering edit mode
16 months ago
GenoMax 147k

What program are you running that gives you this error?

In general it looks like your fastq file was corrupted in pre-processing. Tell us about what you did starting with original fastq data.

ADD COMMENT
0
Entering edit mode

Hey GenoMax, I have updated the details.

ADD REPLY
0
Entering edit mode

You will need to check C_S30_1P.fastq.gz C_S30_2P.fastq.gz using a fastq validation program such as fastqvalidator (LINK). If you have done nothing to these datafiles then perhaps they were corrupted during download from provider.

ADD REPLY
0
Entering edit mode

I executed fastqvalidator in this file including others as well. As such no error was detected in the fastq file. So I guess the files are not corrupted at least.

ADD REPLY
0
Entering edit mode

Hey GenoMax, Yes the files were corrupted. Thanks for the help.

ADD REPLY

Login before adding your answer.

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