我一直在尝试在 fasta 文件中查找每个物种的 1 秒数量,如下所示:
>111
1100101010
>102
1110000001
Run Code Online (Sandbox Code Playgroud)
期望的输出是:
>111
5
>102
4
Run Code Online (Sandbox Code Playgroud)
我知道如何获取文件中 1 的数量:
grep -c 1 file
Run Code Online (Sandbox Code Playgroud)
我的问题是我无法找到跟踪每个物种的 1 数量(而不是文件中的总数)的方法。