MacOSXでBioinfomatics入門その６

■BLAST結果の利用

■Bioperlを使った例

perlは

シークエンスというのは所詮テキスト情報ですから、シークエンスの解析にはperlが非常に強い力を発揮します。

ちなみに、このスクリプトでは、結果をresult.outという形で一つのファイルにまとめるのでなく、test???.txtの結果はtest???.blsという形で個別のファイルになっています。batchBLAST.tcshの実行は、以下のようなコマンドをたたきます。

これで、理論上は何個のシークエンスでも上記のコマンドを１行書くだけでBLAST解析ができるようになりました。

■perlを使う。

perl -pe 's/￥r/￥n/g' （変換するファイル名）> （変換後のファイル名）

シークエンスファイルから任意の部分の配列だけ取り出すスクリプト

４日目

■Bioperlを使う。

一定のスコア以上の配列を取り出す。

５日目

■仕上げ

解析した結果を一つのテキストファイルにまとめ上げる。

ボタン一つで、これまでの内容を自動化する。

Optional

ファイルメーカーに自動的に取り込む。

ClustalWのインストール

■参考になるサイト

Installation (http://www.ncbi.nlm.nih.gov/Kuehl/prefinished/installation.html)

g

AppleのAG BLASTのページ（http://developer.apple.com/hardware/ve/acgresearch.html）

Installing NCBI's BLAST2.x Executables (http://genome.nhgri.nih.gov/blastall/blast_install/)

UNIX用に簡単に手順が書いてある。

NCBI BLASTのセットアップ手順（http://itoshi.tv/bioinformatics/blast2.html）

Linux用のものだがほぼ流用できる。

Installing BioPerl on OS X (http://www.tc.umn.edu/￣cann0010/Bioperl_OSX_install.html)

linux for molbio（http://www.geocities.co.jp/Technopolis/7075/）

linuxで各種分子生物学的ソフトウェアを動かすためのサイト

fink http://fink.sourceforge.net/

Mac OS Xにおいては大くのUnix Softwareが動作しますが、一部は変更が必要になる場合もあります。また動作するとしてもだれも試験していない場合もあります。
また./configureでインストールした場合uninstallがむずかしいsoftwareもあります。
FinkはUnixのOpen Source softwareをMac OS Xにパッケージとして移植するプロジェクト、および、移植したsoftwareを扱うためのソフトの名称です。
プロジェクトのWEBサイトはhttp://fink.sourceforge.net/です。

まず、finkをダウンロードしてインストール

echo "source /sw/bin/init.csh" >> .cshrc

nkf