可以直接cat 多個fq.gz壓縮文件


案例描述:

需要將Sample_test1_R1.fastq.gz和Sample_test2_R2.fastq.gz合並為test.fastq.gz

操作方法1:

先zcat再gzip

 zcat Sample_test_1.R1.fastq.gz Sample_test_2.R2.fastq.gz | gzip - > test.fastq.gz
  • 1

操作方法2:

直接cat

cat Sample_test_1.R1.fastq.gz Sample_test_2.R2.fastq.gz > test2.fastq.gz
  • 1

結果對比

壓縮文件大小

ll test*.gz

-rw-r–r– 1 An Lau 197121 1321311 6月 13 10:14 test.fastq.gz 
-rw-r–r– 1 An Lau 197121 1321742 6月 13 10:15 test2.fastq.gz

解壓縮文件大小

gzip -cd test.fastq.gz > test.fastq 
gzip -cd test2.fastq.gz > test2.fastq 
ll test *.fastq

-rw-r–r– 1 An Lau 197121 7195186 6月 13 10:16 test.fastq 
-rw-r–r– 1 An Lau 197121 7195186 6月 13 10:16 test2.fastq

對比文件內容

zcat Sample_test_1.R1.fastq.gz Sample_test_2.R2.fastq.gz | wc

80000 100000 7195186

wc test*.fastq

80000 100000 7195186 test.fastq 
80000 100000 7195186 test2.fastq 
160000 200000 14390372 total

總結

雖然壓縮文件大小不一樣,應該是壓縮率不同導致的, 
但是不同方法出來的文件大小和內容都是一樣的。 
直接cat的速度肯定會比zcat … | gzip …快! 
所有以后直接cat就可以了。


免責聲明!

本站轉載的文章為個人學習借鑒使用,本站對版權不負任何法律責任。如果侵犯了您的隱私權益,請聯系本站郵箱yoyou2525@163.com刪除。



 
粵ICP備18138465號   © 2018-2025 CODEPRJ.COM