マッピング
使用するツール Cat, Trimmomatic, Picard FastqToSam, Bedtools bamtofastq 今回は何をする? 最近、カバレッジの不足を補うために、別々に実施したWGSデータを合体させて解析に使う機会がありました。 「単純にcatコマンドで合成するだけでできるよ」と共…
FastQC, trimmomatic, bwa, samtools 導入難易度★☆☆☆☆ 使用難易度★★★☆☆ 今回は何をする? 前回の記事で収集したマラリア原虫のWGSデータを使ってマッピングを行い、BAM形式のファイルを生成します。内容は以前に投稿した下記の記事と重複しますが、これから…
SRA形式からBAMファイル出力までのロードマップ 記事にした作業内容について、自分のメモとして作成している図表を時折アップしていきます。 覚書程度の内容です。 Twitterで記事の更新をお知らせしているので、興味を持たれた方は是非フォローをお願いしま…
BWAの使い方。 bwa mem, samtools sort, samtools index これまでに準備したファイルを使って、BWAを用いたリファレンスゲノムへのWGSデータのマッピングを行います。 今回使う主なコマンドは 1. bwa mem 2. samtools sort 3. samtools index の3つ。 Twitt…
リファレンスゲノムのダウンロード、インデックスファイルの作成、SRAからfastqおよびfastq.gzへの変換する。 BWA index, fastq-dump, pigzの使い方。 マッピングに必要な準備としてリファレンスゲノムのダウンロード、インデックスファイルの作成、SRA形式…
prefetchの使い方 複数のSRAファイルをprefetchを使って一括ダウンロードする方法を紹介します。 個々のSRAファイルは、accession numberを検索してブラウザから一つずつダウンロードすることも可能ですが、多数のデータを用いる場合に非効率的です。 そこで…