On 20/07/10 17:27, Peter C. wrote: > $ countseq -sformat=genbank gbvrt1.seq > 31065 Of course, you could just use: $ seqret -filter -sformat=genbank gbvrt1.seq | grep -c '^>' 31065 :-) Peter Rice