I have ONT direct RNA sequencing data consisting of 958 fastq.gz files per sample. I need to merge all of these files into a single file. When attempting to combine them using the command zcat PAW*.fastq.gz > fld.fastq.gz, I encountered an error indicating that the files are not in gzip format. However, the files were basecalled by the ONT machine and have the .gz extension.
How can I resolve this issue and ensure that the files are in the correct gzip format? I tried to unzip them using gunzip but again error says the same
2 answers
You can check the validity of gzipped files using gzip -v -t file.gz.
Should look something like this:
$ gzip -v -t test_file.gz
test_file.gz: OK
If they fail, you can check if they are uncompressed with head file.gz, and it is human readable then it's not compressed. If it's not readable, then maybe it's not compressed with gzip but another program.
You should use cat on .gz files, part of the magic of .gz ...I guess if you use zcat they are unzipped and this is why you get the error - they have the right extension but they are not gzipped. Proof: if you enter head fld.fastq.gz does it return readable text or gibberish?
cat PAW*.fastq.gz > fld.fastq.gz
Log in to answer this question.
-
Error in converting vcf file to plink format
written by Madhawa •
-
How can I makedatabase for env_nt.tar.gz (from NCBI ftp site) in Diamond
written by nonaddldy •
-
readFilesCommand not working
written by vinishavvenugopal •
-
gzip: merged.fastq.gz not in gzip format
written by Kai_Qi
-
bedtools bamtofastq generates fastq which cannot be gzip correctly?
written by francois •
-
Size difference using gunzip-cat-gzip and zcat
written by dazhudou1122
-
STAR Mapping muliple files
written by ipalmisa •
-
Error when unzipping .gz files from GATK?
written by vctrm67 •
-
invalid compressed data--crc error
written by microbez.uaf •
-
fastq.gz goes into gz cpgz loop
written by S0phia •
