通过注释查找相关基因并导出基因序列

查找TPS相关基因

grep -wi --color terpene.*synthase combined_function_annotation.xls >> TPS1.xls

查看搜索到的基因id并导入到txt文件

awk '{print $1}' TPS1.xls >> TPS1.txt | less -S TPS1.txt

根据基因ID循环语句搜索基因序列并导出到TPScds.txt

awk '{print $1}' TPS1.txt | while read id; do grep -A1 ${id} final.cds >> TPScds.txt; done