SPAdes使用方法(基因組拼接)

第一步安裝軟件

conda install SPAdes
#通過miniconda3安裝SPAdes,可隨時調(diào)用,無需再添加環(huán)境變量扑馁。

第二步轉(zhuǎn)到工作目錄

cd /home/monkeyflower/bioworkplace

第三步盼玄,拼裝序列

spades.py --pe1-1 L4-3_1.clean.fq.gz --pe1-2 L4-3_2.clean.fq.gz -t 200 -k 21,33,55,77 -m 600 --careful --phred-offset 33 -o testfile
#--pe1-1贴彼、--pe1-2:需要拼接的前后兩個序列
#-t:線程數(shù)
#-k:kmer選擇的值,默認的值為21,33,55,77,數(shù)字必須是小于128的奇數(shù)
#-m:內(nèi)存埃儿,最大允許使用內(nèi)存(GB)
#--careful:read糾錯,減少錯誤和插入序列
#--phred-offset:堿基質(zhì)量體系,在數(shù)據(jù)糾錯中會用到器仗,illumina數(shù)據(jù)一般采用33
#-o:輸出文件目錄

第四步,導出輸出文件

#文件scaffolds.fasta和contigs.fasta便是需要的結(jié)果文件
#corrected/:包含用BayesHammer糾錯reads產(chǎn)生*.fastq.gz文件童番。
#scaffolds.fasta:包含產(chǎn)生的scaffolds
#contigs.fasta:包含產(chǎn)生的contigs
#assembly_graph.gfa:包含SPAdes組裝圖和scaffolds的路徑精钮。
#assembly_graph.fastg:包含SPAdes組裝圖。(FASTG格式)
#contigs.paths:組裝圖中包含與contigs.fasta對應的路徑剃斧。
#scaffolds.paths:組裝圖中包含與scaffolds.fasta對應的路徑
?著作權歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
  • 序言:七十年代末轨香,一起剝皮案震驚了整個濱河市,隨后出現(xiàn)的幾起案子幼东,更是在濱河造成了極大的恐慌臂容,老刑警劉巖,帶你破解...
    沈念sama閱讀 211,639評論 6 492
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件筋粗,死亡現(xiàn)場離奇詭異策橘,居然都是意外死亡炸渡,警方通過查閱死者的電腦和手機娜亿,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 90,277評論 3 385
  • 文/潘曉璐 我一進店門,熙熙樓的掌柜王于貴愁眉苦臉地迎上來蚌堵,“玉大人买决,你說我怎么就攤上這事沛婴。” “怎么了督赤?”我有些...
    開封第一講書人閱讀 157,221評論 0 348
  • 文/不壞的土叔 我叫張陵嘁灯,是天一觀的道長。 經(jīng)常有香客問我躲舌,道長丑婿,這世上最難降的妖魔是什么? 我笑而不...
    開封第一講書人閱讀 56,474評論 1 283
  • 正文 為了忘掉前任没卸,我火速辦了婚禮羹奉,結(jié)果婚禮上,老公的妹妹穿的比我還像新娘约计。我一直安慰自己诀拭,他們只是感情好,可當我...
    茶點故事閱讀 65,570評論 6 386
  • 文/花漫 我一把揭開白布煤蚌。 她就那樣靜靜地躺著耕挨,像睡著了一般。 火紅的嫁衣襯著肌膚如雪尉桩。 梳的紋絲不亂的頭發(fā)上筒占,一...
    開封第一講書人閱讀 49,816評論 1 290
  • 那天,我揣著相機與錄音魄健,去河邊找鬼赋铝。 笑死,一個胖子當著我的面吹牛沽瘦,可吹牛的內(nèi)容都是我干的革骨。 我是一名探鬼主播,決...
    沈念sama閱讀 38,957評論 3 408
  • 文/蒼蘭香墨 我猛地睜開眼析恋,長吁一口氣:“原來是場噩夢啊……” “哼良哲!你這毒婦竟也來了?” 一聲冷哼從身側(cè)響起助隧,我...
    開封第一講書人閱讀 37,718評論 0 266
  • 序言:老撾萬榮一對情侶失蹤筑凫,失蹤者是張志新(化名)和其女友劉穎,沒想到半個月后并村,有當?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體巍实,經(jīng)...
    沈念sama閱讀 44,176評論 1 303
  • 正文 獨居荒郊野嶺守林人離奇死亡,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點故事閱讀 36,511評論 2 327
  • 正文 我和宋清朗相戀三年哩牍,在試婚紗的時候發(fā)現(xiàn)自己被綠了棚潦。 大學時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片。...
    茶點故事閱讀 38,646評論 1 340
  • 序言:一個原本活蹦亂跳的男人離奇死亡膝昆,死狀恐怖丸边,靈堂內(nèi)的尸體忽然破棺而出叠必,到底是詐尸還是另有隱情,我是刑警寧澤妹窖,帶...
    沈念sama閱讀 34,322評論 4 330
  • 正文 年R本政府宣布纬朝,位于F島的核電站,受9級特大地震影響骄呼,放射性物質(zhì)發(fā)生泄漏共苛。R本人自食惡果不足惜,卻給世界環(huán)境...
    茶點故事閱讀 39,934評論 3 313
  • 文/蒙蒙 一蜓萄、第九天 我趴在偏房一處隱蔽的房頂上張望俄讹。 院中可真熱鬧,春花似錦绕德、人聲如沸患膛。這莊子的主人今日做“春日...
    開封第一講書人閱讀 30,755評論 0 21
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽踪蹬。三九已至,卻和暖如春臣咖,著一層夾襖步出監(jiān)牢的瞬間跃捣,已是汗流浹背。 一陣腳步聲響...
    開封第一講書人閱讀 31,987評論 1 266
  • 我被黑心中介騙來泰國打工夺蛇, 沒想到剛下飛機就差點兒被人妖公主榨干…… 1. 我叫王不留疚漆,地道東北人。 一個月前我還...
    沈念sama閱讀 46,358評論 2 360
  • 正文 我出身青樓刁赦,卻偏偏與公主長得像娶聘,于是被迫代替她去往敵國和親。 傳聞我的和親對象是個殘疾皇子甚脉,可洞房花燭夜當晚...
    茶點故事閱讀 43,514評論 2 348

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