fastq-dump從SRA文件中提取fastq文件

fastq-dump是SRAtoolkit中使用頻率很高的命令,用于從SRA文件中拆解提取fastq文件。具體用法如下:

Usage:
  fastq-dump [options] <path> [<path>...]
  fastq-dump [options] <accession>
Use option --help for more information

一. 拆解一個sra文件

cd ~/Seqs
fastq-dump --split-files SRR6232298.sra
  • SRR6232298.sra是一個PE測序結(jié)果,所以,需要--split-files參數(shù)可以將其分解為兩個fastq文件。
  • 如果不加該參數(shù),則只有1個fastq文件(包含了兩端測序的結(jié)果)

二.批量拆解sra文件

1. 新建腳本文件

nano fqdump.sh

2. 輸入以下腳本

#!/bin/sh
for i in *sra
do
echo $i
fastq-dump --gzip --split-files $i
done

Ctrl+O; Ctrl+X保存退出。
**這里--gzip參數(shù)是為了生成壓縮的gz格式fastq文件,以節(jié)省磁盤空間

3. 運(yùn)行腳本

sh fqdump.sh
image.png

詳細(xì)說明請參加官方Documentation
https://trace.ncbi.nlm.nih.gov/Traces/sra/sra.cgi?view=toolkit_doc&f=fastq-dump
以及hiptop的簡書文章:
http://www.itdecent.cn/p/a8d70b66794c

?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
【社區(qū)內(nèi)容提示】社區(qū)部分內(nèi)容疑似由AI輔助生成,瀏覽時請結(jié)合常識與多方信息審慎甄別。
平臺聲明:文章內(nèi)容(如有圖片或視頻亦包括在內(nèi))由作者上傳并發(fā)布,文章內(nèi)容僅代表作者本人觀點(diǎn),簡書系信息發(fā)布平臺,僅提供信息存儲服務(wù)。

相關(guān)閱讀更多精彩內(nèi)容

友情鏈接更多精彩內(nèi)容