rsem-prepare-reference for human
リファンレンスゲノムとGTFの両方Ensemblから調達。実行時点で最新のEnsembl89を使って検索用のリファレンス作成。 [shell] time rsem-prepare-reference --gtf Homo_sapiens.GRCh38.89.chr.gtf --bowtie2 --bowtie-path /usr/local/bin Homo_sapiens.GRCh38.dna.toplevel.fa rsem-bowtie2/human [/shell] 実行時間はわりとかかる。
2406.83s user 40.73s system 99% cpu 40:51.46 total
これと並行してやっていた実行するためのサンプル配列を取得して、FASTQをdumpして、トリミングするほうがずっと時間結構かかるという罠。