Permalink
Browse files

README

  • Loading branch information...
1 parent 3f43d93 commit 75321bf41c74b7028f9991053bf106b44fd614e5 @pjotrp committed Oct 27, 2013
Showing with 7 additions and 1 deletion.
  1. +7 −1 README.md
View
@@ -98,8 +98,14 @@ Rewrite the FASTA descriptors
fasta_filter.rb --rewrite 'rec.descr =~ /gene=(\S+)/; rec.descr = $1' test.fa
```
+Filters and rewrites can be combined. The rest is up to your
+imagination! One final example to remove low quality sequences from an
+amino acid file (one amino acid dominates):
+
+```sh
+fasta_filter.rb --filter "count = {} ; rec.seq.each_char { |c| count[c] ||= 1 ; count[c] += 1 }; count.values.max.to_f/rec.seq.size < 0.30" < aa.fa > aa1.fa
+```
-Filters and rewrites can be combined. The rest is up to your imagination!
# API Examples

0 comments on commit 75321bf

Please sign in to comment.