blast 使い方 エラーと解決 tblastx

バイオインフォ道場、くまぞうです。

バイオインフォマティクスの解析では、様々なツールを使います。たまにはエラーに遭遇してしまい、解決するのに時間がかかったりします。BLASTは、大変使用頻度の高いツールの1つですが、見慣れないエラーに遭遇したので、エラーと対応方法をまとめます。

スポンサーリンク



BLASTについて

BLASTとは、塩基配列やアミノ酸配列のアライメントと比較を高速に行うためのアルゴリズム・プログラムです。BLASTは、解析の用途に応じて下記のプログラムを提供しています。手持ちの情報・比較したい内容などに応じてこれらを選択します。今回使用したのは、ncbi-blast-2.5.0tblastxおよびmakeblastdbです。

プログラムQueryDB説明
megablast塩基配列塩基配列塩基配列の比較。デフォルト設定でblastnより高速。
blastn塩基配列塩基配列塩基配列の比較。
tblastnアミノ酸配列塩基配列DBを6通りの読み枠で翻訳・比較。
tblastx塩基配列塩基配列Query・DBを6通りの読み枠で翻訳・比較。
blastpアミノ酸配列アミノ酸配列アミノ酸配列の比較。
blastx塩基配列アミノ酸配列Queryを6通りの読み枠で翻訳・比較。

遭遇したエラーと対応方法

遭遇したエラーと原因

インストール後、最初にmakeblastdbを実行しました(正常終了)。その後、tblastxを実行したら、下記のようなエラーが発生しました。原因は、FASTAファイルのフォーマットが崩れていたことでした。FASTAファイルのヘッダー行に>がありませんでした。

Error: NCBI C++ Exception:
T0 "/home/PATH.../src/objtools/readers/fasta.cpp", line 2204: Error: ncbi::objects::CFastaReader::PostWarning() 
- CFastaReader: Near line 1, there's a line that doesn't look like plausible data, but it's not marked as defline or comment. 

対応方法

FASTAファイルのヘッダー行の先頭に>を付与しました。

# ヘッダーを数字を含む行と定義し、>&で先頭に>を挿入。
cat fasta | sed -e '/[0-9]/s/.*/>&/' > fasta2  

フォーマットの崩れを見直し1(エラー発生)

$ cat my.fasta
ENTRY1 | Symbol: ...  # ヘッダー行の先頭に>が無い
AGCAT...

$ tblastx -query my.fasta -db db/mydb  # 異常終了
Error: NCBI C++ Exception:
T0 "/home/PATH.../src/objtools/readers/fasta.cpp", line 2204: Error: ncbi::objects::CFastaReader::PostWarning() 
- CFastaReader: Near line 1, there's a line that doesn't look like plausible data, but it's not marked as defline or comment. 
$ cat my.fasta | sed -e '/[0-9]/s/.*/>&/' | tee my2.fasta  # ヘッダーに">"追加
>ENTRY1 | Symbol: ...
AGCAT...

$ tblastx -query my2.fasta -db db/mydb  # 正常終了

その他のエラーと対応方法

今回のエラーに対する対策を調べているとき、その他のエラーで困っている人の意見が参考になりました。確認できたものを載せておきます。

フォーマットの崩れを見直し2(不要なスペース):エラーなし

FASTAのフォーマットでは、「ヘッダー行の>の後に不要なスペースを入れない」とあります。フォーマット不正かもしれませんが、エラーにはなりませんでした。但し、同様にFASTA形式を要求する他のプログラムでは、エラーになることがあります。

$ cat my.fasta
>   ENTRY1 | Symbol: ...  # >の後に、不要なスペースが挿入されている
AGCAT...

$ tblastx -query my.fasta -db db/mydb  # 正常終了

データベースに関するエラー1(makeblastなし)(エラー発生)

BLASTのデータベースを作成するには、makeblastdbプログラムを実行します。データベースを作成しないで、もととなるFASTAファイルをdbオプションに指定したら、「BLAST Database error」が出力されました。makeblastdbを使って、データベースを作成してからblastプログラムを動かせば良いです。

$ tblastx -query my.fasta -db db/mydb.fasta  # makeblastdb無し
BLAST Database error: No alias or index file found for nucleotide database [...] in search path

$ makeblastdb -in db/mydb.fasta -dbtype nucl -out mydb  # makeblastdb実行
$ tblastx -query my.fasta -db db/mydb  # 正常終了

データベースに関するエラー2(dbtype違い)(エラー発生)

BLASTのデータベースを作成するには、makeblastdbプログラムを実行します。その時、dbtypeオプションでnucl(塩基配列の場合)またはprot(アミノ酸配列の場合)を指定します。tblastxなどの指定した解析プログラムに対応したデータベース(tblastxは塩基データベースを要求)が正しく参照できないとエラーになり、「BLAST Database error」が出力されます。makeblastdbを使って、対応するデータベースを作成してからblastプログラムを動かせば良いです。

$ makeblastdb -in db/mydb.prot -dbtype prot -out myprot  # makeblastdb(prot)
$ tblastx -query my.fasta -db db/myprot  # アミノ酸データベースを参照するように指示
BLAST Database error: No alias or index file found for nucleotide database [...] in search path
スポンサーリンク





カテゴリー

今週よく読まれている記事

  1. 学会・会議で英語が絶望的にできなくて困るケース | オンライン開催で「英語力のなさ」を痛感した場合の対処法

    学会・会議で英語ができなくてに困るケース学会やイベント・会議などが、オンラインで開催される…
  2. r tidyverse 使い方 | 列 filter 絞り込み select関数 – dplyrパッケージ

    tidyverseで1つのデータフレームの列の絞り込みは、dplyrパッケージのselect関数を…
  3. プログラミングで疲れた脳をリフレッシュ 〜 鬼滅の刃「感動」と「やる気アップ」でストレス発散!

    ストレス発散は鬼滅で。「50%OFF」で読む!脳のパフォーマンスを上げるには、適度な休憩と…
  4. AWS ディスク容量不足 新しいボリュームを追加する

    バイオインフォマティクスでは大きなファイルを扱うことがあるので、ディスク不足に陥ることがあります。…
  5. 「知っている」と「知らない」とでは、もしものとき、大違いになる – コロナうつ対策

    コロナうつなどという言葉を聞くようになりましたが、派遣切り、解雇、リストラは、これから本格化します…

人気記事

  1. R言語

    R subset関数 データフレームやmatrixからの条件指定による行・列の抽出
    R の subset関数は、データフレームやマトリックスか…
  2. IGV, 解析ツール

    IGV 使い方 インストール〜便利な使い方まで | リファレンス・マッピングデータ・アノテーションを読み込んで表示しよう
    IGV(Integrative Genomics View…
  3. Excel, 統計

    z-score 計算方法 エクセル(Excel) 編
    統計処理で、大きく変化しているなどの判断基準にも使われる値…

おすすめ記事

  1. シェルスクリプト

    シェルスクリプト | ファイル存在チェック・空ファイルチェック
    bashでスクリプトを作成するときに、よく使うのがファイル…
  2. R言語, グラフ

    R 使い方 軸・ラベルの調整(向き・サイズ・色など) グラフの描き方
    Rによるplot(グラフ)の描画は、手軽で大変便利です。た…
  3. awk, シェルスクリプト

    bash 部分文字列・置換・長さ・連結・分割の文字列処理
    bashのよく使う文字列処理、部分文字列・置換・連結・長さ…