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

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

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


注解如下:

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

你可能感兴趣的:(如何将样品名添加到contig名称的开头位置)