A tool to compress FASTQ files with ultra-high compression ratio and high speed. repaq supports compressing the FASTQ to .rfq or .rfq.xz formats. Compressing to .rfq is ultra fast, while compressing to .rfq.xz provides very high compression ratio.
For NovaSeq data, as an example:
- the
.rfqfile can be much smaller than.fq.gz, and the compressing time is usually less than 1/5 of gzip compression. - The
.rfq.xzfile can be as small as 5% of the original FASTQ file, or smaller than 30% of the.fq.gzfile.
For paired-end FASTQ files, repaq compresses them into one single file to provide higher compression ratio.
This tool also supports non-Illumina format FASTQ (i.e. the BGI-SEQ format), but the compression ratio is not as good Illumina format FASTQ.
Citation: Chen S, Chen Y, Wang Z, Qin W, Zhang J, Nand H, Zhang J, Li J, Zhang X, Liang X and Xu M (2023) Efficient sequencing data compression and FPGA acceleration based on a two-step framework. Front. Genet. 14:1260531. doi: 10.3389/fgene.2023.1260531
Here we demonstrate the compression ratio of two paired-end NovaSeq data. You can download these files and test locally.
nova.R1.fq: 1704 MB, the original read1 file, http://opengene.org/repaq/testdata/nova.R1.fqnova.R2.fq: 1704 MB, the original read2 file, http://opengene.org/repaq/testdata/nova.R2.fqnova.R1.fq.gz: 308 MB (CR 18.08%), the gzipped read1, http://opengene.org/repaq/testdata/nova.R1.fq.gznova.R2.fq.gz: 325 MB (CR 19.07%), the gzipped read2, http://opengene.org/repaq/testdata/nova.R2.fq.gznova.rfq: 333 MB (CR 9.77%), the repacked file of read1+read2, http://opengene.org/repaq/testdata/nova.rfqnova.rfq.xz: 134 MB (CR 3.93%), the xz compressednova.rfq, http://opengene.org/repaq/testdata/nova.rfq.xz
See? The size of final nova.rfq.xz is only 3.39% of the original FASTQ files! You can decompress it and check the md5 to see whether they are identical!
Typically with one single CPU core, it takes less than 1 minute to convert nova.R1.fq + nova.R2.fq to nova.rfq, and takes less than 5 minutes to compress the nova.rfq to nova.rfq.xz by xz.
conda install -c bioconda repaqThis binary is only for Linux systems: http://opengene.org/repaq/repaq
# this binary was compiled on CentOS, and tested on CentOS/Ubuntu
wget http://opengene.org/repaq/repaq
chmod a+x ./repaq# get source (you can also use browser to download from master or releases)
git clone https://github.com/OpenGene/repaq.git
# build
cd repaq
make
# Install
sudo make installFor single-end mode:
# compress to .rfq.xz
repaq -c -i in.fq -o out.rfq.xz
# decompress from .rfq.xz
repaq -d -i in.rfq.xz -o out.fqFor paired-end mode:
# compress to .rfq.xz
repaq -c -i in.R1.fq -I in.R2.fq -o out.rfq.xz
# decompress from .rfq.xz
repaq -d -i in.rfq.xz -o out.R1.fq -O out.R2.fqTips:
-iand-Ialways denote the first and second input files, while-oand-Oalways denote the first and second output files.- the FASTQ input/output files can be gzipped if their names are ended with
.gz. - for paired-end data. the .rfq file created in paired-end mode is usually much smaller than the sum of the .rfq files created in single-end mode for R1 and R2 respectively. To obtain high A6B8 compression rate, please always use PE mode for PE data.
- if you want higher speed and are not concern with compression ratio, replace
xxx.rfq.xzwithxxx.rfq, then repaq will compress or decompress.rfqformat.
- Memory: 16G RAM
- CPU: 4 cores
repaq offers a compare mode to check the consistency of the original FASTQ file(s) and the compressed .rfq or .rfq.xz file.
- set
--compareto enable thecomparemode - specify the .rfq or .rfq.xz file by
-roption - specify the FASTQ files by
-iand-Ioptions.
Examples:
# for single-end data
repaq --compare -i original.R1.fq -r compressed.rfq.xz
# for paired-end data
repaq --compare -i original.R1.fq.gz -I original.R2.fq.gz -r compressed.rfq.xzWithout any expection, you will get an output of a JSON like:
{
"result":"passed",
"msg":"",
"fastq_reads":50000,
"rfq_reads":50000,
"fastq_bases":7419082,
"rfq_bases":7419082
}The result will be "failed" if the compressed file is not consistent with the original FASTQ files.
repaq can read the input from STDIN, and write the output to STDOUT.
- specify
--stdinif you want to read the STDIN for compression or decompression. - specify
--stdoutif you want to output to the STDOUT for compression or decompression - in decompression mode, if
--stdoutis specified, the output will be interleaved PE stream. - if the STDIN is an interleaved paired-end stream, specify
--interleaved_into indicate that. - be noted that STDIN cannot be read when the input is a .xz file, and STDOUT cannot be written when the output is a .xz file
Here gives you an example of compressing the interleaved PE output from fastp by directly using pipes:
fastp -i R1.fq -I R2.fq --stdout | repaq -c --interleaved_in --stdin -o out.rfq.xzrepaq was initially designed for compressing Illumina data, but it also works with data from other platforms, like BGI-Seq. To work with repaq, the FASTQ format should meet following condidtions:
- only has bases A/T/C/G/N.
- each FASTQ record has, and only has four lines (name, sequence, strand, quality).
- the name and strand line cannot be longer than 255 bytes.
- the number of different quality characters cannot be more than 127.
repaq works best for Illumina data directly output by bcl2fastq.
options:
-i, --in1 input file name
-o, --out1 output file name
-I, --in2 read2 input file name when encoding paired-end FASTQ files
-O, --out2 read2 output file name when decoding to paired-end FASTQ files
-c, --compress compress input to output
-d, --decompress decompress input to output
-k, --chunk the chunk size (kilo bases) for encoding, default 1000=1000kb.
--stdin input from STDIN. If the STDIN is interleaved paired-end FASTQ, please also add --interleaved_in.
--stdout write to STDOUT. When decompressing PE data, this option will result in interleaved FASTQ output for paired-end input. Disabled by defaut.
--interleaved_in indicate that <in1> is an interleaved paired-end FASTQ which contains both read1 and read2. Disabled by defaut.
# following options are used to check the consistency of the compressed data
-p, --compare compare the files read by read to check the compression consistency. <rfq_to_compare> should be specified in this mode.
-r, --rfq_to_compare the RFQ file to be compared with the input. This option is only used in compare mode.
-j, --json_compare_result the file to store the comparison result. This is optional since the result is also printed on STDOUT.
# options for .xz output
-t, --thread thread number for xz compression (default 1). When compression level (-z) is >= 4, no threading will be used.
-z, --compression compression level. Higher level means higher compression ratio, and more RAM usage (1~9), default 3.
# verify the output when compressing
-v, --verify verify the output stream to ensure compression is correct.
-f, --fast_verify only verify part (10%) of the output stream to save time.
-?, --help print this messagerepaq makes a system call in order to run the xz compression tool available on GNU/Linux systems. If xz isn't installed, repaq will fail with the message:
failed to call xz, please confirm that xz is installed in your system
Shifu Chen, Yaru Chen, Zhouyang Wang, Wenjian Qin, Jing Zhang, Heera Nand, Jishuai Zhang, Jun Li, Xiaoni Zhang, Xiaoming Liang, Mingyan Xu. Efficient sequencing data compression and FPGA acceleration based on a two-step framework. Frontiers in Genetics, 2023, https://doi.org/10.3389/fgene.2023.1260531