如何使用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
答案 0 :(得分:2)
$ awk '/>/{sub(/\|/,"-"++i"|")}1' file
>chromosome1-1|2199-2200
----------------------
>chromosome1-2|3546-3548
----------------------
>chromosome1-3|6489-6548
----------------------