gffread を使った transcripts fasta 転写物の配列取得

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

ゲノム配列から転写物の配列を取得

スポンサーリンク



転写物(transcript)の配列を取得したい場合、染色体上の位置(開始位置〜終了位置)をたよりに配列情報を抜き出します。位置情報を管理する方法としてよく使われるのがGFF3やGTF形式のファイルです。領域が小さくて少ない場合は自力でも切り出せますが、対象物がたくさんで構造をよくみないといけない場合は、ツールを使ったほうが便利です。cufflinksのUtilityの1つである、gffreadを使えば、転写物の配列を簡単に切り出すことができます。

gffread を使った転写物配列の切り出し

ゲノム配列(genome.fa)と位置情報(GFF3 または GTFファイル)を準備して、gffreadを実行します。オプションを豊富に持ち、以下の例では-wを指定してexonを基にした転写物の切り出し(fastaファイルの生成)を行っています。

# gffread GFF3を使った転写物(transcript)配列の切り出し
$ gffread -w transcripts.fa -g genome.fa transcripts.gff3

gffreadについて

gffreadのヘルプを参照することでオプションを確認することができます。転写物の配列切り出しの他に、gffからgtfのフォーマット変換(-Tオプション)も便利です。gffreadは、cufflinksのインストールでも使えるようになります。参考情報・ダウンロード

# gffread GFF3からGTFへの変換
$ gffread -E annotation.gff -T -o- | more
# gffread ヘルプ・オプションの参照
$ gffread --help
スポンサーリンク





ピックアップ記事

  1. R subset関数 データフレームやmatrixからの条件指定による行・列の抽出

    R の subset関数は、データフレームやマトリックスから条件にマッチした部分集合を取り出します…
  2. R 使い方 グラフの色 自作のカラーパレットとデフォルトのカラーパレット

    グラフを描画する際、色の組み合わせを1つ1つ考えるのは面倒ですよね。Rには、あらかじめ色の組み合わ…
  3. awk split サンプルでわかる列の分割とセパレータ(文字・正規表現)の指定方法

    awkのsplit関数は、セパレータを指定して文字列を分割することができます。セパレータには文字や…

人気記事

  1. IGV, 解析ツール

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

    z-score 計算方法 エクセル(Excel) 編
    統計処理で、大きく変化しているなどの判断基準にも使われる値…
  3. R データ型 - 文字列・ベクター・データフレーム・マトリックス など-, R言語, スクリプト

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

おすすめ記事

  1. awk, bash 文字列操作, シェルスクリプト

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

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

    R 使い方 軸・ラベルの調整(向き・サイズ・色など) グラフの描き方
    Rによるplot(グラフ)の描画は、手軽で大変便利です。た…