如何编写循环遍历chicken_topmotifs.fasta文件的所有偶数行？

我有一个结构如下的FASTA文件，FASTA是一种文本文件，以>开头的行是头部信息，下方对应序列文本。需要编写循环遍历名为chicken_topmotifs.fasta文件的所有偶数行。

文件内容示例： #

>gene8
ATGAATTATTATACACCTCAAATACTCTCCTCAATCTCTCCAACATTCCCCACCACAATTCTCGGTGACTTTACTACACTTCTACAATCATACACTTCT

>gene12
ATGGTAGATCTCTATTACGATTATCTTTCTTAGATCACATAATTATCACCCCCCCTTATAAATCTACACTTCTACAACCAATTACACTTCTACAAAACA
>gene18
ATGCTTTTACACTTCTACAACTACTTTTAACTCGATACTTCTACAATCTACACATATCACAATAACAAAAACAAAAAGCTACTAATATATATATATACA
>gene21
ATGTCTCAATTTCACCAATCTATAATTTACTACGCCGTACTCTTTATAACCTTACTTTCTTAAATAACATTACACTTCTACATTACATATTTTACATCA

原代码问题 #

你当前写的循环逻辑有误：

for sequence in chicken_topmotifs.fasta;
do
    echo $sequence
done

这个循环只是把文件名chicken_topmotifs.fasta作为唯一的循环变量值，只会执行一次并输出文件名，完全没有读取文件内容。

正确实现方式 #

方法一：用awk提取偶数行（高效简洁）

如果只是需要输出偶数行，直接用awk一行搞定：

awk 'NR%2==0' chicken_topmotifs.fasta

如果需要循环处理每一行内容，可结合while循环：

while read -r line; do
    echo "$line"
    # 在这里添加你对每行序列的处理逻辑
done < <(awk 'NR%2==0' chicken_topmotifs.fasta)

方法二：用for循环遍历偶数行号（适合必须用for循环的场景）

先获取文件总行数，再遍历所有偶数行号，用sed提取对应行：

total_lines=$(wc -l < chicken_topmotifs.fasta)
for line_num in $(seq 2 2 "$total_lines"); do
    current_line=$(sed -n "${line_num}p" chicken_topmotifs.fasta)
    echo "$current_line"
    # 添加你的处理逻辑
done

内容的提问来源于stack exchange，提问作者Arjun Harer