linux 中 shell 统计每条scaffold的长度

 

001、

复制代码
[root@pc1 test]# ls
a.fa
[root@pc1 test]# cat a.fa     ## 测试fasta文件
>chr1
aatt
cc
>chr2
ttgg
ccgg
>chr3
aa
[root@pc1 test]# awk '{if($0 ~ /^>/) {print chr, sum; chr = $0; sum = 0} else {sum += length($0)}} END {print chr, sum}' a.fa

>chr1 6
>chr2 8
>chr3 2
复制代码

 

posted @   小鲨鱼2018  阅读(46)  评论(0编辑  收藏  举报
相关博文:
阅读排行:
· 震惊!C++程序真的从main开始吗?99%的程序员都答错了
· 【硬核科普】Trae如何「偷看」你的代码?零基础破解AI编程运行原理
· 单元测试从入门到精通
· 上周热点回顾(3.3-3.9)
· winform 绘制太阳,地球,月球 运作规律
历史上的今天:
2021-10-27 R语言中%*%运算符
2021-10-27 windows中如何查看端口占用情况、端口是否开启
2021-10-27 R语言中setdiff、intersect、union函数
2021-10-27 R语言中merge函数
2021-10-27 linux系统中paste命令实现列转行
2020-10-27 linux vmware虚拟机添加硬盘
点击右上角即可分享
微信分享提示