2011-05-04

Bioinformatics

• Velvet 安裝紀錄(2009-11-11)
• 測試 Velvet 1.1.03 (支援 OpenMP）
• 系統：Debian Squeeze

  ~$ sudo apt-get update; sudo apt-get upgrade; sudo apt-get dist-upgrade
  ~$ sudo apt-get install build-essential 
  ~$ wget http://www.ebi.ac.uk/~zerbino/velvet/velvet_1.1.03.tgz
  ~$ tar zxvf velvet_1.1.03.tgz
  ~$ cd velvet_1.1.03
  ~/velvet_1.1.03$ make
  ~/velvet_1.1.03$ ./velveth output 17 -fasta -short data/test_long.fa
  ~/velvet_1.1.03$ ./velvetg output/ -exp_cov 3 -min_contig_lgth 100
  

• 使用 ​valgrind 來進行記憶體分析

  ~/velvet_1.1.03$ sudo apt-get install valgrind
  ~/velvet_1.1.03$ valgrind --tool=massif ./velveth output 17 -fasta -short data/test_long.fa
  ~/velvet_1.1.03$ ms_print massif.out.10571
  ~/velvet_1.1.03$ valgrind --tool=massif ./velvetg output/ -exp_cov 3 -min_contig_lgth 100
  ~/velvet_1.1.03$ ms_print massif.out.10561
  ~/velvet_1.1.03$ valgrind -v --leak-check=full --show-reachable=yes --read-var-info=yes ./velveth output 17 -fasta -short data/test_long.fa 2> velveth.log
  ~/velvet_1.1.03$ valgrind -v --leak-check=full --show-reachable=yes --read-var-info=yes ./velvetg output/ -exp_cov 3 -min_contig_lgth 100 2> velvetg.log
  

• 使用 gprof 進行 profile 分析

  ~/velvet_1.1.03$ cat Makefile | sed "s#^CFLAG.*#CFLAGS = -Wall -pg#" > 
  ~/velvet_1.1.03$ mv Makefile.new Makefile
  ~/velvet_1.1.03$ make
  ~/velvet_1.1.03$ ./velveth output 17 -fasta -short data/test_long.fa
  ~/velvet_1.1.03$ du -sh gmon.out 
  56K	gmon.out
  ~/velvet_1.1.03$ ./velvetg output/ -exp_cov 3 -min_contig_lgth 100
  ~/velvet_1.1.03$ du -sh gmon.out 
  152K	gmon.out
  

• 使用 ​gprof2dot 把 gprof log 轉成 graphviz 的 dot 格式，就可以畫出圖形了。

  ~/velvet_1.1.03$ sudo apt-get install graphviz python
  ~/velvet_1.1.03$ wget http://gprof2dot.jrfonseca.googlecode.com/hg/gprof2dot.py
  ~/velvet_1.1.03$ chmod a+x gprof2dot.py
  ~/velvet_1.1.03$ gprof velveth | ./gprof2dot.py | dot -Tpng -o velveth.png
  ~/velvet_1.1.03$ gprof velvetg | ./gprof2dot.py | dot -Tpng -o velvetg.png
  

• 
•