在处理fasta序列时，我们常常需要知道每条序列的长度。今天小编将与大家分享四种获取fasta序列长度的方法。
一、使用awk
以下是使用awk获取fasta序列长度的代码：
该命令的输出如下：
二、使用bioawk
安装bioawk并使用以下命令获取fasta序列长度：
输出结果如下：
三、使用samtools
首先生成.fai文件，然后提取前两列：
Docker Desktop（linux）
当前 Docker 最新稳定版本之一，主要针对稳定性和兼容性进行了修复优化，适合生产环境与日常开发使用。该版本继续强化 AI 开发支持、容器日志管理以及 Docker Engine 的安全能力，对 Windows/macOS/Linux 平台兼容性进行了进一步优化。
下载
生成的.fai文件如下，前两列显示fasta序列的名字和长度。
.fai文件的每一列具体含义：
四、使用seqkit
安装seqkit并使用以下命令获取fasta序列长度：
awk '/^>/{if (l!="") print l; print; l=0; next}{l+=length($0)}END{print l}' test.fasta>Chr115704606>Chr108327059>Chr118696443>Chr127853492>Chr134340075>Chr146288713>Chr29509017>Chr311222275>Chr47647452>Chr57499939>Chr64872821>Chr78973615>Chr88275968>Chr98318069conda install bioawkbioawk -c fastx '{ print $name, length($seq) }' Chr1 15704606Chr10 8327059Chr11 8696443Chr12 7853492Chr13 4340075Chr14 6288713Chr2 9509017Chr3 11222275Chr4 7647452Chr5 7499939Chr6 4872821Chr7 8973615Chr8 8275968Chr9 8318069#生成.fai文件samtools faidx test.fasta#提取前两列cut -f1-2 test.fasta.fai第一列 NAME : 序列的名称，只保留“>”后，第一个空白之前的内容；第二列 LENGTH: 序列的长度， 单位为bp；第三列 OFFSET : 第一个碱基的偏移量， 从0开始计数，换行符也统计进行；第四列 LINEBASES : 除了最后一行外， 其他代表序列的行的碱基数， 单位为bp；第五列 LINEWIDTH : 行宽， 除了最后一行外， 其他代表序列的行的长度， 包括换行符，conda install seqkitseqkit fx2tab --length --name --header-line test.fasta