linux 中 awk命令实现将fasta文件中每个scaffold中的所有碱基转换为一行

 

001、

复制代码
[b20223040323@admin1 test]$ ls
test.fa
[b20223040323@admin1 test]$ cat test.fa     ## 测试数据
>chr1
aabbccdd
eeff
>chr2
xxyyzzqq
gg
>chr3
ddjjiill
ssff                                        ## 转换程序
[b20223040323@admin1 test]$ awk '{if(NR == 1) {print $0} else if($0 ~ />/) {printf("\n"); print $0} else {printf("%s", $0)}} END {printf("\n")}' test.fa
>chr1
aabbccddeeff
>chr2
xxyyzzqqgg
>chr3
ddjjiillssff
复制代码

 

 

posted @   小鲨鱼2018  阅读(57)  评论(0编辑  收藏  举报
相关博文:
阅读排行:
· 震惊!C++程序真的从main开始吗?99%的程序员都答错了
· 【硬核科普】Trae如何「偷看」你的代码?零基础破解AI编程运行原理
· 单元测试从入门到精通
· 上周热点回顾(3.3-3.9)
· winform 绘制太阳,地球,月球 运作规律
历史上的今天:
2021-02-22 python中实现列表的倒序排列
点击右上角即可分享
微信分享提示