使用序号

时间:2016-06-16 12:21:47

标签: bash awk sed gsub fasta

如何使用awk / sed / grep更改文件中更多fasta序列的标题,如:

>chromosome1|2199-2200
----------------------

>chromosome1|3546-3548
----------------------

>chromosome1|6489-6548
----------------------

到此:

>chromosome1-1|2199-2200
------------------------

>chromosome1-2|3546-3548
------------------------

>chromosome1-3|6489-6548
------------------------

目前,我可以在名称的开头或结尾添加一个序号,如:

awk '/^>/{gsub(/^>/,">Seq"i++" ");}1'

>Seq1chromosome1|2199-2200

1 个答案:

答案 0 :(得分:2)

$ awk '/>/{sub(/\|/,"-"++i"|")}1' file
>chromosome1-1|2199-2200
----------------------

>chromosome1-2|3546-3548
----------------------

>chromosome1-3|6489-6548
----------------------