Fastq-dump:我的日常命令

原文地址:Fastq-dump: 一個神奇的軟件 - by hoptop

感謝我洲更學(xué)長~記錄一下看完學(xué)長的這篇文章之后對于我自己的fastq-dump使用建議:

默認(rèn)命令:

fastq-dump /path/to/###.sra

trinity的建議命令

在無參組裝時,使用trinity軟件在使用以上默認(rèn)命令后報錯會提示使用以下的命令:

fastq-dump --defline-seq '@$sn[_$rn]/$ri' --split-files file.sra

根據(jù)洲更學(xué)長的文章,把--split-files改成--split-3會更好一些(具體區(qū)別請移步學(xué)長的簡書文章,鏈接在上方)。

故最后得到命令是這樣的:

日常命令

fastq-dump --defline-seq '@$sn[_$rn]/$ri' --split-3 file.sra -o ../02.fastq

習(xí)慣默認(rèn)輸出到上一層目錄中的02.fastq這個文件夾中。想要輸出到別的位置記得修改-O|--outdir參數(shù)喲

選加參數(shù):

  • --gzip, --bzip2: 以壓縮文件的方式輸出結(jié)果
  • --fasta:如果下游分析只需要用到fasta文件。當(dāng)然了也有很多方法能夠把fastq轉(zhuǎn)換成fasta,比如說samtools.(嗯,偷懶必備技能~)

后記:

因為接下來會做蠻多的轉(zhuǎn)錄組或者基因組相關(guān)的工作,fastq-dump軟件簡直是大寶天天見,就像我陳飛師兄說的一樣garbage in, garbage out這第一步就出了岔子后續(xù)再怎么分析都是白搭(突然想起我那11月開始做實驗到現(xiàn)在一個多月因為一開始引物設(shè)計反了導(dǎo)致實驗一點結(jié)果都沒有的師兄。。心疼)

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

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

  • 這個不錯分享給大家,從扣上看到的,就轉(zhuǎn)過來了 《電腦專業(yè)英語》 file [fail] n. 文件;v. 保存文...
    麥子先生R閱讀 7,130評論 5 24
  • 現(xiàn)在可以用fasterq-dump, 速度更快,請閱讀都8102年了,還用fastq-dump,快換fasterq...
    xuzhougeng閱讀 69,320評論 23 100
  • Facebook中國對手是哪些 6月2日,美國社交巨頭Facebook終于在萬眾期待中遞交了IPO申請,文件除了披...
    69a9cb275028閱讀 1,051評論 0 0
  • (1) 大晚上看到九寨溝發(fā)生7級地震的消息,瞬間傻了。 萬萬沒想到,神奇美麗的九寨,被稱為“人間天堂”的九寨,竟然...
    橫渡太平洋的魚閱讀 595評論 4 7
  • 每個人都是一個圈 每個圈都生活在不同的空間 當(dāng)兩個圈自行飄落于一個平面時 必會產(chǎn)生兩種不同的結(jié)果 1.會出現(xiàn)平行移...
    沉默的泥閱讀 608評論 0 0

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