如果您喜欢我们的文章,请设“关注”哦~
七、通配符:(glob):*、?、[]
* 匹配任意数量的字符:ls *.fastq.gz ?匹配单个字符:ls sample?.bam ls sample_[12].fastq.gz for f in *.fastq.gz; do echo "Processing $f" done 花括号可将一组命令作为整体进行重定向或逻辑控制。注意 {后需留空格:
{ echo "Header"; echo "Body";} > out.txt
{ echo "Sample QC:"; fastqc sample.fastq.gz;} > qc_report.txt
echo "Price is \$20" echo '$HOME' " "会保留变量替换:echo "$HOME"
bcftools view -i 'QUAL > 30' bwa mem \ -t 8 \ ref.fa R1.fq R2.fq \ > aln.sam 这种写法在运行 GATK 等长命令时非常常见。
zcat sample.fastq.gz \ | fastp -i - -o clean.fastq \ 2> fastp.log \ && echo "QC done" \ & 该命令通过管道将解压的 FASTQ 直接传给 fastp,同时将错误输出写入日志。在 fastp 成功后输出提示信息,并将整个任务放入后台执行。
@H 开头的 header、统计剩余行数,并分别保存标准输出与错误输出。参考命令如下:zcat reads.fastq.gz \ | grep -v '^@H' \ | wc -l > count.txt \ 2> err.txt
for f in *.fastq.gz; do fastqc "$f" >> stdout.log 2>> stderr.log &done
(bcftools view -i 'DP>10 && QUAL>30' input.vcf || echo "filter failed") \ > filtered.vcf 2> log.txt 该命令在子 Shell 中运行过滤操作;若过滤失败,则输出提示信息。随后将标准输出写入filtered.vcf,并将标准错误写入log.txt。
⭐ SDbioinfo-往期回顾 ⭐
【栏目1:播客分享会】
【栏目2:前沿信息早知道】
【栏目3:知识解读】
【栏目4:最新招聘】

【?我是谁?】
SDbioinfo,专注「计算生物学」领域的学习智库
我们汇聚了来自哈佛剑桥港大、北清复交等高校科研工作者。我们的愿景:链接领域最强大脑,用数据和分析改变生命。
不管你是:
⭐ 找实习、找工作、出国留学
⭐ 找校友、找搭子、找资源
⭐ 找科研方法、前沿信息、大咖分享
统统可以帮你,加入SDbioinfo,和全球生信人一起成长!
(添加微信:SDbioinfo_2023)
