fastq-dump并行版pfastq-dump的使用

fastq-dump轉(zhuǎn)換SRA文件到fastq文件很慢,并行版本成為趨勢(shì)设联;

無論怎么換捌年,先要打好基礎(chǔ)谴轮,使用并行版本的前提是要保證NCBI的fastq-dump可以在服務(wù)器上正常運(yùn)行。

首先安裝Sratoolkit的最新版(v.2.9.2):

mkdir -p?/path-to-Sratoolkit/ && cd?/path-to-Sratoolkit/

wget? https://ftp-trace.ncbi.nlm.nih.gov/sra/sdk/2.9.2/sratoolkit.2.9.2-ubuntu64.tar.gz && \

tar zxfv sratoolkit.2.9.2-ubuntu64.tar.gz&& \

mv?sratoolkit.2.9.2-ubuntu64/* . && \ rm -rf?sratoolkit.2.9.2-ubuntu64.tar.gz?sratoolkit.2.9.2-ubuntu64


下載pfastq-dump:

git clone https://github.com/inutano/pfastq-dump && \

cd pfastq-dump && \

chmod a+x bin/pfastq-dump && \

ln -s bin/pfastq-dump? /path-to-Sratoolkit/bin


把安裝的路徑加入到賬號(hào)下的$PATH中:

echo 'PATH=/home/luna/Desktop/Software/Sratoolkit/bin:$PATH' >>?~/.bashrc && \

cp?~/.bashrc?~/.bash_profile && \

source?~/.bashrc?~/.bash_profile


使用pfastq_dump识窿,因?yàn)閜fastq_dump是基于fastq_dump寫的一個(gè)bash程序,所以參數(shù)是相同的:

對(duì)于單端數(shù)據(jù)轉(zhuǎn)換脑融,轉(zhuǎn)換后文件是fq.gz:

for id in *sra;????do pfastq-dump --threads 10 ./$id --gzip;????done?

對(duì)于雙端數(shù)據(jù)轉(zhuǎn)換喻频,轉(zhuǎn)換后文件是fq.gz:

for id in *sra;? ? ? ? do pfastq-dump? --threads 8 ./$id --split-3 --gzip;????done?


直接用sra號(hào)下載并解壓fastq文件,但是推薦下載好文件再使用fastq_dump轉(zhuǎn)換肘迎,且文件后綴是.sra(請(qǐng)注意):

單端數(shù)據(jù):

for id in SRR799545? SRR799544;????do pfastq-dump --threads 10 -s $id?--gzip;????done

雙端數(shù)據(jù):

for id in SRR799545??SRR799544;????do pfastq-dump --threads 10 -s $id?--split-3 --gzip;????done


經(jīng)過測(cè)試甥温,其實(shí)也不是那么的快啊<瞬肌姻蚓!很揪心!

最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
  • 序言:七十年代末匣沼,一起剝皮案震驚了整個(gè)濱河市狰挡,隨后出現(xiàn)的幾起案子,更是在濱河造成了極大的恐慌释涛,老刑警劉巖加叁,帶你破解...
    沈念sama閱讀 221,888評(píng)論 6 515
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件,死亡現(xiàn)場(chǎng)離奇詭異唇撬,居然都是意外死亡它匕,警方通過查閱死者的電腦和手機(jī),發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 94,677評(píng)論 3 399
  • 文/潘曉璐 我一進(jìn)店門窖认,熙熙樓的掌柜王于貴愁眉苦臉地迎上來豫柬,“玉大人,你說我怎么就攤上這事扑浸∩崭” “怎么了?”我有些...
    開封第一講書人閱讀 168,386評(píng)論 0 360
  • 文/不壞的土叔 我叫張陵首装,是天一觀的道長(zhǎng)创夜。 經(jīng)常有香客問我,道長(zhǎng)仙逻,這世上最難降的妖魔是什么驰吓? 我笑而不...
    開封第一講書人閱讀 59,726評(píng)論 1 297
  • 正文 為了忘掉前任,我火速辦了婚禮系奉,結(jié)果婚禮上檬贰,老公的妹妹穿的比我還像新娘。我一直安慰自己缺亮,他們只是感情好翁涤,可當(dāng)我...
    茶點(diǎn)故事閱讀 68,729評(píng)論 6 397
  • 文/花漫 我一把揭開白布。 她就那樣靜靜地躺著,像睡著了一般葵礼。 火紅的嫁衣襯著肌膚如雪号阿。 梳的紋絲不亂的頭發(fā)上,一...
    開封第一講書人閱讀 52,337評(píng)論 1 310
  • 那天鸳粉,我揣著相機(jī)與錄音扔涧,去河邊找鬼。 笑死届谈,一個(gè)胖子當(dāng)著我的面吹牛枯夜,可吹牛的內(nèi)容都是我干的。 我是一名探鬼主播艰山,決...
    沈念sama閱讀 40,902評(píng)論 3 421
  • 文/蒼蘭香墨 我猛地睜開眼湖雹,長(zhǎng)吁一口氣:“原來是場(chǎng)噩夢(mèng)啊……” “哼!你這毒婦竟也來了曙搬?” 一聲冷哼從身側(cè)響起摔吏,我...
    開封第一講書人閱讀 39,807評(píng)論 0 276
  • 序言:老撾萬(wàn)榮一對(duì)情侶失蹤,失蹤者是張志新(化名)和其女友劉穎织鲸,沒想到半個(gè)月后舔腾,有當(dāng)?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體,經(jīng)...
    沈念sama閱讀 46,349評(píng)論 1 318
  • 正文 獨(dú)居荒郊野嶺守林人離奇死亡搂擦,尸身上長(zhǎng)有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點(diǎn)故事閱讀 38,439評(píng)論 3 340
  • 正文 我和宋清朗相戀三年稳诚,在試婚紗的時(shí)候發(fā)現(xiàn)自己被綠了。 大學(xué)時(shí)的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片瀑踢。...
    茶點(diǎn)故事閱讀 40,567評(píng)論 1 352
  • 序言:一個(gè)原本活蹦亂跳的男人離奇死亡扳还,死狀恐怖,靈堂內(nèi)的尸體忽然破棺而出橱夭,到底是詐尸還是另有隱情氨距,我是刑警寧澤,帶...
    沈念sama閱讀 36,242評(píng)論 5 350
  • 正文 年R本政府宣布棘劣,位于F島的核電站俏让,受9級(jí)特大地震影響,放射性物質(zhì)發(fā)生泄漏茬暇。R本人自食惡果不足惜首昔,卻給世界環(huán)境...
    茶點(diǎn)故事閱讀 41,933評(píng)論 3 334
  • 文/蒙蒙 一、第九天 我趴在偏房一處隱蔽的房頂上張望糙俗。 院中可真熱鬧勒奇,春花似錦、人聲如沸巧骚。這莊子的主人今日做“春日...
    開封第一講書人閱讀 32,420評(píng)論 0 24
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽(yáng)。三九已至竣蹦,卻和暖如春顶猜,著一層夾襖步出監(jiān)牢的瞬間,已是汗流浹背痘括。 一陣腳步聲響...
    開封第一講書人閱讀 33,531評(píng)論 1 272
  • 我被黑心中介騙來泰國(guó)打工驶兜, 沒想到剛下飛機(jī)就差點(diǎn)兒被人妖公主榨干…… 1. 我叫王不留,地道東北人远寸。 一個(gè)月前我還...
    沈念sama閱讀 48,995評(píng)論 3 377
  • 正文 我出身青樓,卻偏偏與公主長(zhǎng)得像屠凶,于是被迫代替她去往敵國(guó)和親驰后。 傳聞我的和親對(duì)象是個(gè)殘疾皇子,可洞房花燭夜當(dāng)晚...
    茶點(diǎn)故事閱讀 45,585評(píng)論 2 359

推薦閱讀更多精彩內(nèi)容

  • RNAseq實(shí)際操作(實(shí)戰(zhàn)) 首先聲明矗愧,雖然是實(shí)戰(zhàn)灶芝,但是其實(shí)是學(xué)習(xí)筆記而已,初學(xué)唉韭,參考了大量大神的博客和帖子,還有...
    zd200572閱讀 4,404評(píng)論 0 29
  • 第1章 小試牛刀 $ 是普通用戶夜涕,# 表示管理員用戶 root。 shebang:#!属愤。sharp / hash ...
    巴喬書摘閱讀 6,349評(píng)論 1 4
  • Spring Cloud為開發(fā)人員提供了快速構(gòu)建分布式系統(tǒng)中一些常見模式的工具(例如配置管理女器,服務(wù)發(fā)現(xiàn),斷路器住诸,智...
    卡卡羅2017閱讀 134,702評(píng)論 18 139
  • 現(xiàn)在可以用fasterq-dump, 速度更快驾胆,請(qǐng)閱讀都8102年了,還用fastq-dump贱呐,快換fasterq...
    xuzhougeng閱讀 67,166評(píng)論 23 99
  • 系統(tǒng)管理與維護(hù)命令 date date(選項(xiàng))(參數(shù)) | 選項(xiàng) | 說明 | | :-------- | ...
    蓓蓓的萬(wàn)能男友閱讀 3,899評(píng)論 0 5