gravatar for Peter Chung

2 hours ago by

Hong Kong

I have a paired end fastq file and when I ran trim galore, the error said lengths of bases and qualities do not match,
so I searched the solution and recommended to using bbtools to discard reads that have mismatching lengths of bases and qualities: in=pair_1.fq.gz in2=pair_2.fq.gz out=fixed_1.fq.gz out2=fixed_2.fq.gz tossbrokenreads=t

The error
Set INTERLEAVED to false
Input is being processed as paired
pigz: abort: read error on pair_1.fq.gz (Input/output error)
There appear to be different numbers of reads in the paired input files.
The pairing may have been corrupted by an upstream process. It may be fixable by running
at stream.ConcurrentGenericReadInputStream.pair(
at stream.ConcurrentGenericReadInputStream.readLists(
at stream.ConcurrentGenericReadInputStream.run0(

so I tried to repair it by using bbtools in1=pair_1.fq.gz in2=pair_2.fq.gz out1=fixed_1.fq.gz out2=fixed_2.fq.gz outs=singletons.fq repair

Set INTERLEAVED to false
Started output stream.
pigz: abort: read error on pair_1.fq.gz (Input/output error)
Mismatch between length of bases and qualities for read 107893745 (id=ST-E00126:1085:HF3YVCCX2:1:2106:16620:58339 1:N:0:TAAGCTCC+AGATCTCG).


This can be bypassed with the flag 'tossbrokenreads' or 'nullifybrokenquality'
at shared.KillSwitch.kill(
at stream.Read.validateQualityLength(
at stream.Read.validate(
at stream.Read.<init>(
at stream.Read.<init>(
at stream.FASTQ.quadToRead_slow(
at stream.FASTQ.toReadList(
at stream.FastqReadInputStream.fillBuffer(
at stream.FastqReadInputStream.nextList(
at stream.ConcurrentGenericReadInputStream$ReadThread.readLists(
at stream.ConcurrentGenericReadInputStream$

but the error said tossbrokenreads again, it kinda fall into a loop.

Anyone has experience this, please advice. Thanks.


modified 1 hour ago

2 hours ago

Peter Chung120

Source link