One Liners

Count number of sequences in a multi-sequence fasta file

grep -c "^>" fasta_file

and if the ”-” key is broken on your keyboard

grep "^>" fasta_file | wc

Count number of records in a multi-records GenBank file:

grep -c "^//$" genbank_file

Normalize headers in a multi-sequence fasta file

seqret -sequence myfasta.fa -outseq myfasta2.fa

via What you're doing is rather desperate, more on seqret.

 
one-liners.txt · Last modified: 2007/08/20 13:27 (external edit)
 
Recent changes RSS feed Creative Commons License Donate Powered by PHP Valid XHTML 1.0 Valid CSS Driven by DokuWiki