将多行序列变为一行|删除ID后的注释信息并将ID与seq分行>输出新文件
命令行:
awk '/^>/&&NR>1{print "";}{ printf "%s",/^>/ ? $0" ":$0 }' cds.fa|sed 's/trans=[a-zA-Z0-9]*\-[A-Z]*//'|sed 's/length=[0-9]* /\n/'>cDNA.fa
原始序列文件:
结果序列文件:
78 篇文章
如果觉得我的文章对您有用,请随意打赏。你的支持将鼓励我继续创作!