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
_______________________________________________
EMBOSS mailing list
[email protected]
http://lists.open-bio.org/mailman/listinfo/emboss

Reply via email to