From ee79b8675e99fe1622d061bfd088ba352331c21b Mon Sep 17 00:00:00 2001 From: kojix2 <2xijok@gmail.com> Date: Sun, 13 Dec 2020 09:39:00 +0900 Subject: [PATCH] Update README.md --- README.md | 8 +++++++- 1 file changed, 7 insertions(+), 1 deletion(-) diff --git a/README.md b/README.md index 0c409cb..f4dfb9c 100644 --- a/README.md +++ b/README.md @@ -89,8 +89,14 @@ cat gencode.v35.annotation.gff3 \ ![count](https://user-images.githubusercontent.com/5798442/101999832-30b1ae80-3d24-11eb-96fe-e5000bed1f5c.png) -Note: `count` is slower than other Unix commands because it runs in a Ruby script. +Note: `count` command is not very fast because it runs in a Ruby script. +This is fine if the data is small, that is, in most cases. However, if you want to visualize huge data, it is faster to use a combination of common Unix commands as shown below. +```sh +cat gencode.v35.annotation.gff3 | grep -v '#' | grep 'gene' | cut -f1 \ +|sort | uniq -c | sort -nrk2 | awk '{print $2,$1}' \ +| uplot bar -d ' ' -t "The number of human gene annotations per chromosome" -c blue +``` ## Usage