如何将样品名添加到contig名称的开头位置

需要将样品名称添加到contig 名称的前面,用"__"分隔,并删除序列的长度及其他信息,可采用如下方法:

sed '/^>/s/k/SRR1035_k/' | awk -F '[ ]' '{print $1} > SRRSRR1035.contig.fasta'


注解如下:

但是目前仅实现了单个文件的替换,下一步摸索采用循环替换。自学的路是道阻且长,加油吧!

©著作权归作者所有,转载或内容合作请联系作者
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

推荐阅读更多精彩内容