根据你的程序,是不是就是要一个这样的功能:
> 开头的行略去。
不是的行一肯只含有 agtc 这四个字母,数一数就行了?
这个东西能慢到哪去呀。你上传个数据实例我测测呗。
[ 本帖最后由 pangding 于 2011-7-1 16:40 编辑 ]
> 开头的行略去。
不是的行一肯只含有 agtc 这四个字母,数一数就行了?
这个东西能慢到哪去呀。你上传个数据实例我测测呗。
[ 本帖最后由 pangding 于 2011-7-1 16:40 编辑 ]
[home@cpp]cat input_file.txt >header_1 ACGTGAGAGATAGAGC >header_2 AGATGAGATGAGAGATAGA . .
[home@cpp]awk '$0!~">"{total += length($1)}END{print total}' input_file.txt > output_file.txt [home@cpp]cat output_file.txt 35
[home@cpp]cat input_file.txt >header_1 ACGTGAGAGATAGAGC >header_2 AGATGAGATGAGAGATAGA . .
[home@cpp]awk '$0!~">"{total += length($1)}END{print total}' input_file.txt > output_file.txt [home@cpp]cat output_file.txt 35