然后在Linux中解压缩

Question

所以我有一个解开文件的脚本：

#!/bin/bash -e 
# will unzip the data without removing the zipped version
gzip -dc $1 > RawData/unzipped/$(basename $1 .gz)

我想在该未拉紧文件上执行代码，我有

# will run fast qc on the argument passed
fastqc RawData/unzipped/$(basename $1 .gz) --outdir=fastReports/

但是第二个脚本似乎从未执行。（请注意，这些是在同一脚本中，所以我假设它将在第二个脚本之前执行初始脚本）

拉链：

14624_1#10_1.fastq.gz  14624_1#12_2.fastq.gz  14624_1#4_1.fastq.gz  14624_1#7_1.fastq.gz
14624_1#10_2.fastq.gz  14624_1#1_2.fastq.gz   14624_1#4_2.fastq.gz  14624_1#7_2.fastq.gz
14624_1#11_1.fastq.gz  14624_1#2_1.fastq.gz   14624_1#5_1.fastq.gz  14624_1#8_1.fastq.gz
14624_1#11_2.fastq.gz  14624_1#2_2.fastq.gz   14624_1#5_2.fastq.gz  14624_1#8_2.fastq.gz
14624_1#1_1.fastq.gz   14624_1#3_1.fastq.gz   14624_1#6_1.fastq.gz  14624_1#9_1.fastq.gz
14624_1#12_1.fastq.gz  14624_1#3_2.fastq.gz   14624_1#6_2.fastq.gz  14624_1#9_2.fastq.gz

提取：

14624_1#10_1.fastq  14624_1#12_1.fastq  14624_1#3_1.fastq  14624_1#5_2.fastq  14624_1#8_1.fastq
14624_1#10_2.fastq  14624_1#12_2.fastq  14624_1#3_2.fastq  14624_1#6_1.fastq  14624_1#8_2.fastq
14624_1#11_1.fastq  14624_1#1_2.fastq   14624_1#4_1.fastq  14624_1#6_2.fastq  14624_1#9_1.fastq
14624_1#11_2.fastq  14624_1#2_1.fastq   14624_1#4_2.fastq  14624_1#7_1.fastq  14624_1#9_2.fastq
14624_1#1_1.fastq   14624_1#2_2.fastq   14624_1#5_1.fastq  14624_1#7_2.fastq

Answer 1

您可能只使用zcat并即时处理文件：

fastqc <(zcat path/to/file.gz)

btw， <()语法是一个过程替代。

如果您同时需要解压缩文件和过程结果，则可以使用tee：

fastqc <(zcat path/to/file.gz | tee file)