awk批量提取序列

awk批量提取序列

在提取前需保证序列文件仅有一列!

awk '{print$1}' input.fa > ouput.fa#就可将ID后面的其余注释信息去掉,仅保留ID
awk -F '>' 'NR=FNR{a[$0]};NR!=FNR;NF==2{b=($2 in a)} b'  id.txt all.pep.fa > oupt_result.fa
smyang2018

跳至工具栏