この記事では、特定のSRA IDを使用してFASTQ形式のHi-Cデータをダウンロードするスクリプトの作成と実行方法について説明します。スクリプトは、失敗した場合1分後にリトライする機能を備えています。
conda仮想環境構築
まず、必要なパッケージをインストールするためのconda仮想環境を作成します。
conda create -n sra_env_x86_64 python=3.8
conda activate sra_env_x86_64
SRAtoolsのインストール
次に、SRAtoolkitを仮想環境にインストールします。
conda install -c bioconda sra-tools
インストールの確認
インストールが正しく行われたか確認するために、以下のコマンドを実行します。
fastq-dump --version
スクリプトの作成
以下のスクリプトを任意のエディタで作成し、download_fastq.shという名前で保存します。
#!/bin/bash
SRA_ID=$1
OUTPUT_DIR="/path/to/your/output_directory"
LOG_FILE="${OUTPUT_DIR}/${SRA_ID}.log"
# conda環境のアクティベート
source $(conda info --base)/etc/profile.d/conda.sh
conda activate sra_env_x86_64
while true; do
nohup fastq-dump --split-files --outdir ${OUTPUT_DIR} ${SRA_ID} > ${LOG_FILE} 2>&1
if [ $? -eq 0 ]; then
echo "Download completed successfully."
break
else
echo "Download failed, retrying in 1 minute..."
sleep 60
fi
done
実行方法
スクリプトを実行するには、ターミナルを開き、スクリプトが保存されているディレクトリに移動して以下のコマンドを実行します。
bash download_fastq.sh <SRA_ID>
<SRA_ID>の部分は、ダウンロードしたいSRAのIDに置き換えてください。
実行中の動作
スクリプトを実行すると、以下のようなメッセージがターミナルに表示されます。
- ダウンロードが成功した場合:
Download completed successfully. - ダウンロードが失敗し、リトライする場合:
Download failed, retrying in 1 minute...
これにより、Hi-Cデータのダウンロードプロセスが自動化され、失敗した場合も再試行することで、データを確実に取得することができます。
ログファイルについて
スクリプトは、ダウンロードプロセスの詳細なログを指定されたログファイルに保存します。ログファイルのパスは、LOG_FILE変数で指定されており、各SRA IDに対応したファイル名が生成されます。これにより、ダウンロードプロセス中に発生したエラーや進行状況を確認することができます。
例えば、SRA IDがSRR123456の場合、ログファイルは/path/to/your/output_directory/SRR123456.logという名前で保存されます。ログファイルには、fastq-dumpコマンドの出力が全て記録されます。