bioawk

https://github.com/lh3/bioawk

1、基本思想

使用：

usage: bioawk [-F fs] [-v var=value] [-c fmt] [-tH] [-f progfile | 'prog'] [file ...]

bioawk基本思想是把组成不同类型的文件（sam、bam、fasta、fastq、vcf）的基本元素封装成变量，直接调用即可。

上面出现的名称即可引用其变量。

2、实际例子

打印fasta序列ID、序列、长度、GC含量：

bioawk -c fastx '{print "ID: "$name"\tlength: "length($seq)"\tGC: "gc($seq)"\t"$seq}' demo.fa

引用：https://blog.csdn.net/qq_42491125/article/details/92849378