更改fasta标头以包含计数信息

时间:2016-03-11 07:57:52

标签: unix header fasta

我有一个包含多行的文件:

  

fg-000001 ATGGATGACCGATGCTAGC 34 23 33 45 34 23 34 35 34 35 56 47

我想将它转换为一个fasta文件,其中包含计数信息的总和以及名称基本上是一个fasta文件

  

FG-000001_x433   ATGGATGACCGATGCTAGC

这可能吗? 提前谢谢!

1 个答案:

答案 0 :(得分:0)

awk脚本应该按照您的要求执行:

awk '{sum=0; for (i=3; i<=NF; i++) {sum += $i}; print $1"_x"sum, $2}' file