Watch out for the infoseq header line you should use -noheading to avoid a +1
to your total
infoseq -sformat=genbank gbvrt1.seq -noheading -auto | wc -l
Andrew
-----Original Message-----
From: [email protected]
[mailto:[email protected]] On Behalf Of Mahmut Uludag
Sent: 20 July 2010 21:57
To: Peter; Peter Rice
Cc: [email protected]
Subject: Re: [EMBOSS] Counting the number of sequences in a file
>> $ seqret -filter -sformat=genbank gbvrt1.seq | grep -c '^>'
infoseq prints a separate line for each sequence, following command line may
also work.
> $ infoseq -filter -sformat=genbank gbvrt1.seq | wc -l
Mahmut
_______________________________________________
EMBOSS mailing list
[email protected]
http://lists.open-bio.org/mailman/listinfo/emboss
--
Disclaimer: This e-mail and any attachments are confidential and intended
solely for the use of the recipient(s) to whom they are addressed. If you have
received it in error, please destroy all copies and inform the sender. This
email and any attachments are believed to be free from viruses but BBSRC
accepts no liability in connection therewith.
_______________________________________________
EMBOSS mailing list
[email protected]
http://lists.open-bio.org/mailman/listinfo/emboss
