Aide-memoire

Bio-infomaticsで困った事の備忘録的なサムシング

自分が困った事をまとめた備忘録的なサムシング

mclblastlineを使いたい

一つのFastaファイルの中で,配列類似性に基づいてクラスタリングしたい

→mcl (http://micans.org/mcl/)を使う

 

いろんなやりかたあるけど,mclblastlineを使ったら楽かもしれない?

## protein配列

formatdb -i input.fa -p T

 

blastall -p blastp -i input.fa -d input.fa -o input.fa.blast -m 8 -e 1e-10

mclblastline --mcl-I=1.0 --blast-m9  input.fa.blast

  # --mcl-te=20 でコア数指定

 

cytoscape用のsifファイルを作りたい

mcxdump -imx out.input.fa.blast.I10 -dump-sif A -tabr input.fa.blast.ta

 

という感じ