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
という感じ