2016年5月18日
extracting fasta records from a multi-fasta file based on a list using awk
摘要： ``` for i in $(cat gene_list) do awk -v RS=">" '($1==a){print ">"$0}' a=$i input.fasta >> out.fa sed -i '/^$/d' out.fa done ``` 阅读全文
posted @ 2016-05-18 18:48 liuhui_pine 阅读(121) 评论(0) 推荐(0)