Question: Tophat prep_read error
0
gravatar for Yuka Takemon
3.4 years ago by
Yuka Takemon20
The Jackson Laboratory
Yuka Takemon20 wrote:

Hello,

I am running a Tophat script that has worked previously, however with the new data set I am trying to map the sequence length and my quality length in all of my fastq files are incongruent (I am guess that quality length and sequence length are usually the same).

I am getting the following error message:

[2015-12-18 22:09:55] Checking for Bowtie

                  Bowtie version:        2.2.3.0

[2015-12-18 22:09:55] Checking for Bowtie index files (genome)..

[2015-12-18 22:09:55] Checking for reference FASTA file

[2015-12-18 22:09:55] Generating SAM header for /pathto/genome_index/Mus_musculus/UCSC/mm10/Sequence/Bowtie2Index/genome

[2015-12-18 22:10:04] Preparing reads

        [FAILED]

Error running 'prep_reads'

Error: qual length (47) differs from seq length (100) for fastq record !

Each fastq file has 100 seq length but varying qual lengths. 

When running Fastqc, my reads seem normal with good quality score with no indication of different quality and seq lengths.

Am I understanding this error correctly and how do I fix this?

Thanks in advance!

 

rna-seq tophat illumina • 1.5k views
ADD COMMENTlink modified 2.3 years ago by lprivatjogja10 • written 3.4 years ago by Yuka Takemon20
2
gravatar for Brian Bushnell
3.4 years ago by
Walnut Creek, USA
Brian Bushnell16k wrote:

If the error message is correct, your fastq file is corrupt; you should redownload it.

ADD COMMENTlink written 3.4 years ago by Brian Bushnell16k
0
gravatar for Yuka Takemon
3.4 years ago by
Yuka Takemon20
The Jackson Laboratory
Yuka Takemon20 wrote:

Cheers Brian, 

I think you're right. I ran into a common bug on Chrome, where users downloading large files will often have incomplete files. I am re-downloading all the data again using a different browser. 

Yuka

ADD COMMENTlink written 3.4 years ago by Yuka Takemon20
0
gravatar for lprivatjogja1
2.3 years ago by
lprivatjogja10 wrote:

i wish i could help..

ADD COMMENTlink modified 2.3 years ago • written 2.3 years ago by lprivatjogja10
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 707 users visited in the last hour