2023-07-20 提取序列中包含某个字符的全部序列

提取以P1结尾的全部序列

awk '/^>.*P1$/{flag=1;print;next} />/{flag=0} flag'  xxx.fa > P1_xxx.fa

这个命令使用awk来查找包含"以>开头，以P1结尾"的行，并将标志位"flag"设置为1以指示接下来的行需要被输出。然后，当下一个包含">"的行被找到时，标志位将被重置为0。最后，输出结果被重定向到一个名为"P1_sequences.txt"的新文件中