染色体のACGTをカウント

デバッグのために、各染色体の塩基組成を調べます。


tail +n 2  chrX.fa | sed -e "s/[^aA]//g" | tr -d '/n' | wc


awk使った方がもっと早くいけたかなぁ。
最初小文字大文字あるのを忘れてました。。

コメント

このブログの人気の投稿

再現性なんてないさ(?)

特に収穫がない日

jupytherにチャレンジ