利用masurca进行混合二、三代测序数据基因组组装

现在组装新基因组,一般用pacbio的hifi测序模式进行olc组装+Hi-c辅助组装染色体就行了,但是我们实验室以前测了一些基因组是低覆盖度的pacbio的clr模式测序+50x以上二代,没有HI-C数据做的泛基因组,还是有需求对这一部分组装进行改进。

该软件结合了debreijn 图和 Overlap-Layout-Consensus的方法,3.2.1版本以来masurca可以混合二代、三代数据进行基因组组装。

(按我的理解,2代组装用了debreijn 图,组装成super reads.然后以3代为骨架进行map比对,比对上的super reads 用OLC进行组装)

对计算资源的要求是:

1 10Mb : 16Gb 内存,8 + 核,10Gb 磁盘空间 少于1小时

2 500Mb: 128Gb 内存,16个内核,1Tb 磁盘空间 1-2天

3 1Gb : 256Gb 内存,32 + 核,2Tb 磁盘空间 4-5天

4 3Gb: 512Gb 内存,32 + 核心,5Tb 磁盘空间 15-20天

5 30Gb : 1Tb 内存,64 + 核,10Tb + 磁盘空间 60-90天

所以各位自己看一下怎么搞。

安装

有人制作了bioconda的包,但是官方明确说明不能用bioconda安装……

1 安装numactl(必须要有,否则flye组装部分无法执行,然后程序弹出)

sudo apt-get install numactl

我搞了3-4天没组装成功,一直显示3代测序数据有问题,要我检查,但是也没查出毛病来,后面在某个论坛的旮旯里找到了一个回答,让安装这个,然后就没问题了。

2 下载boost最新版本,进行安装,

Boost C++ Libraries

下载完了,解压,进入文件夹

sudo ./boostrap.sh

sudo ./b2 install

然后就安装完毕了

3  安装masurca

GitHub - alekseyzimin/masurca 下载最新安装包

解压,进入文件夹

sh ./install.sh 会自动安装依赖软件包

以上步骤就完成了全部的安装,反正我这里这样处理完就可以组装了

使用

这个软件进化到4.03版本的时候可以使用一条简单的命令行进行组装

a 没有三代测序

/path_to_MaSuRCA/bin/masurca -t 32 -i /path_to/pe_R1.fa,/path_to/pe_R2.fa

b 有三代测序

/path_to_MaSuRCA/bin/masurca -t 32 -i /path_to/pe_R1.fa,/path_to/pe_R2.fa -r /path_to/nanopore.fastq.gz

一般情况下还是使用安装文件夹里面内置的config文件进行配置

主要就是2代、3代测序文件的位置,预期基因组大小、能用多少个核。具体就不细说了,文档写的很详细。

配置文件写好以后,$/install _ path/MaSuRCA-X.X. x/bin/masurca config.txt用这个命令能够生成一个assemble.sh的文件,然后运行就可以了。

需要注意的是ubuntu内置的sh是dash,不能执行这个shell,必须要bash  ./assemble.sh才行。

另外 flye组装步骤只支持最高128线程,但是前序步骤我这里160都没问题,可能可以支持更好,所以个人经验是config里第一次线程数写你的机器能用的最高线程,等跑到flye步骤出错后,再吧线程改为128以下,重新接着跑。

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 159,015评论 4 362
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 67,262评论 1 292
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 108,727评论 0 243
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 43,986评论 0 205
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 52,363评论 3 287
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 40,610评论 1 219
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 31,871评论 2 312
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 30,582评论 0 198
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 34,297评论 1 242
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 30,551评论 2 246
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 32,053评论 1 260
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 28,385评论 2 253
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 33,035评论 3 236
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 26,079评论 0 8
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 26,841评论 0 195
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 35,648评论 2 274
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 35,550评论 2 270

推荐阅读更多精彩内容