提取以P1结尾的全部序列
awk '/^>.*P1$/{flag=1;print;next} />/{flag=0} flag' xxx.fa > P1_xxx.fa
这个命令使用awk来查找包含"以>开头,以P1结尾"的行,并将标志位"flag"设置为1以指示接下来的行需要被输出。然后,当下一个包含">"的行被找到时,标志位将被重置为0。最后,输出结果被重定向到一个名为"P1_sequences.txt"的新文件中
提取以P1结尾的全部序列
awk '/^>.*P1$/{flag=1;print;next} />/{flag=0} flag' xxx.fa > P1_xxx.fa
这个命令使用awk来查找包含"以>开头,以P1结尾"的行,并将标志位"flag"设置为1以指示接下来的行需要被输出。然后,当下一个包含">"的行被找到时,标志位将被重置为0。最后,输出结果被重定向到一个名为"P1_sequences.txt"的新文件中
本文标题:2023-07-20 提取序列中包含某个字符的全部序列
本文链接:https://www.haomeiwen.com/subject/lxvoudtx.html
网友评论