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

2022-12-15  本文已影响0人  田小田_8afd

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

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

注解如下:

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

上一篇下一篇

猜你喜欢

热点阅读