kuroの覚え書き

96の個人的覚え書き

2017-05-10から1日間の記事一覧

RNAseqパイプライン3

cufflinksの結果をcuffmergeで統合し、cuffdiffでサンプル間のFPKMを比較するところまで自動化 cuffmerge_diffスクリプト #!/bin/sh #tophat_resultsディレクトリを引数として実行 # cd "$@" ls */*/transcripts.gtf > transcripts.gtf.txt #アノテーション…

RNAseqパイプライン2.1

フルパス付きとか拡張子とか気にせずにすべてドラッグ&ドロップでできるようにスクリプトをブラッシュアップ tophat_cufflinks2 #!/bin/sh #アノテーション情報を含むgtfファイル read -p "annotation infomation (full path) = " gtf #リファレンスファイ…