bowtie 使い方 マッピングとオプション

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

bowtie は、代表的なマッピングソフトです。高速にマッピング処理を行うことが出来ます。ただし、ギャップは考慮しません。

スポンサーリンク



Bowtie – 使い方

  1. まずは、インデックスファイルを作る

    マッピングのリファレンスとなるFastaファイルからインデックスファイルを作成します。インデックスファイルの作成には時間がかかります。bowtieのサイトに予め作成されたインデックスファイルが準備されているので、ダウンロードして使うことができます。Pre-built indexes

    bowtie-build -f ref_sequence index_name
  2. マッピングを行う

    bowtie-buildで作成したインデックスと、リードファイルを指定してマッピングを行います。

    bowtie -S -p 4 index_name single-read.fastq mapped.sam
    bowtie -S -p 4 index_name -1 paired-read1.fastq -2 paired-read2.fastq mapped.sam

Bowtie – よく用いられるオプション

アライメント

  • -v int
    リード全体で許容される最大ミスマッチ数。-v 0は、完全一致のみ。
  • -n int
    -lで指定された長さについて許容されるミスマッチ数
  • -l int
    シード配列の長さ

レポート

  • -k int
    リード当たりのレポートするヒット数の最大数。
  • -a
    全てのヒットをレポートする。
  • -m int
    指定数以上にヒットが得られる場合はレポートしない。-m 1ならば、マッピングがユニークなリードのみ。
  • –best
    n モードでは、backtrackを125->800に増やす。
    v モードでは、良いアライメントから順にレポートする。
  • –strata
    アライメント候補グループ(stratum)が複数存在する場合、より良いグループ(stratum)をレポートする。–best も一緒に指定する。

その他

  • -p int
    スレッド数
  • -S
    SAM出力

よく使われるパラメータ

  • デフォルト
    “-n 2 -l 28 -e 70”
  • 最も少ないミスマッチでマップされるリードのみ出力
    “-a –best –strata”

参考情報

bowtie マニュアル
bowtie インストール

スポンサーリンク





カテゴリー

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

  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. Excel, 統計

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

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

おすすめ記事

  1. R言語, グラフ

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

    bash 部分文字列・置換・長さ・連結・分割の文字列処理
    bashのよく使う文字列処理、部分文字列・置換・連結・長さ…
  3. シェルスクリプト

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