标签: grep
我需要从CARD数据库中分离出特定的基因(基本上是fasta文件)。我创建了一个感兴趣的基因名称列表
grep -wA1 -f list.file database.file
(A1也得到各自的序列)
问题是,即使命令正确识别感兴趣的基因,输出仍然是完整的数据库文件,只是标记了搜索模式。那当然没有太多帮助......
有人知道我做错了什么吗?我知道,这个问题可以通过循环解决,但我认为-f选项正是我需要的,而且它曾经很好地工作。