home
MENUPREVNEXT

MacOSXでBioinfomatics入門その7

そのほかのソフトウェア

 

■BLAST結果の利用

 

 

 

■Bioperlを使った例

perlは

 

 

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

batchBLAST.tcsh

#!/bin/tcsh
set seqdir=$1
cd $seqdir
foreach file (*.Seq)
blastall -p blastn -d /bio/blast/db/nt -i $file > ${file:r}.bls
end

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

% tcsh /bio/scripts/batchBLAST.tcsh /bio/test

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

 

■perlを使う。

perl -pe 's/¥r/¥n/g' (変換するファイル名)> (変換後のファイル名)

 

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

4日目

■Bioperlを使う。

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

5日目

■仕上げ

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

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

 

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


 
     
更新記録

●2004年4月1日:新規掲載

home
MENUPREVNEXT