海宝建议用这个拼接软件

http://www.broadinstitute.org/software/discovar/blog/?page_id=98

DISCOVAR – variant caller 适合于call variant 和拼接小基因组

DISCOVAR de novo 适合拼接大基因组

下载:

ftp://ftp.broadinstitute.org/pub/crd/DiscovarDeNovo/latest_source_code/LATEST_VERSION.tar.gz

安装:

General Instructions for Building Our Software

System Requirements

Software released from the CRD group at the Broad Institute is built and tested on a modern version of Linux for the x86_64 architecture. Our software does not run on 32-bit machines: you must have a 64-bit Linux system. Our users have successfully built and executed our software using a variety of Linux distributions including Ubuntu, RedHat, and SUSE. We expect that any flavor of x86_64 Linux will work fine, as long as it provides the necessary software prerequisites, as described below.

Basic Compiler and Library Requirements

We rely on reasonably up-to-date versions of these software packages:

  • GCC, with its associated g++ compiler for the C++ language, version 4.7 or later.  We're now using C++11 features, and require a modern GCC.
  • For ARACHNE only, the Xerces-C library. The source can be downloaded from the Xerces-c download page.  Supply the argument --with-xerces=/path/to/xercesc/installation when you run./configure.

The Build Procedure

  • Move the package that you downloaded from our FTP site to a location on you system where you'd like to build the software.
  • Extract the contents by typing: tar xzf package.tgz
  • You'll have a new subdirectory in the current directory named after the package and revision. cd package Some older releases spill all the source into the current directory, rather than creating a new subdirectory. If that's the case, ignore this step.
  • Execute the configuration script by typing: ./configure This assumes that you can copy executables to /usr/local/bin.  If you cannot, you should instead execute: ./configure --prefix=/path/to/install/directory Note: Some older packages lack a configuration script. Consider returning to the FTP site for a more recent version, or just skip this step.
  • Build the software by typing: make all
  • Install the software by typing: make install
  • The executables will be in /usr/local/bin or in /path/to/install/directory/bin.  Make sure this is on your path.

If this goes well, you're ready to go. Consult the manual for the package to learn how to set up your data and what programs to execute.

对于数据的要求:

Sequencing data requirements summary
● Illumina MiSeq or HiSeq 2500 genome sequencers
● PCR-free library preparation
● 250 base paired end reads (or longer)
● ~450 base pair fragment size
● ~60x coverage

Input files
DISCOVAR requires a BAM file containing the raw reads from the sequencer. For variant calling it also
requires a matching reference FASTA file.

call variant  命令:

DISCOVAR can currently generate variants for small regions, and not the entire genome at once. To
generate variants for a 100 kb region for example, use:
Discovar \
READS=reads.bam \
OUT_HEAD=assembly \
REGIONS=1:50000150000
\
REFERENCE=genome.fasta
The complete set of variant calls for this region is given in the text file:
assembly.final.variant

Input files

DISCOVAR requires a BAM file containing the raw reads from the sequencer. For variant calling it also
requires a matching reference FASTA file.
BAM files
The reads to assemble must be in a BAM file or files. The name of the BAM file is specified with the
required argument READS :
READS= filename
Multiple BAM files are specified using a comma separated list:
READS= filename1,filename2,...
Alternatively, the BAM files can be specified in a separate file contain a list of BAM filenames, one per
line:
READS= @listfilename

DISCOVAR calls SAMtools internally to extract reads from the BAM.

Reference file (optional)

This is only required if you are using DISCOVAR as a variant caller. The reference information is used
only for variant calling and not in the assembly process. Specifying a valid FASTA reference file is all
that is required to cause DISCOVAR to generate variants.
To specify a reference FASTA file use the optional argument REFERENCE :
REFERENCE= filename

It should be the same file that was used to generate the alignments in the input BAM file(s), or at least
should share the same coordinate system. The FASTA record names should match those in the BAM
file. Ns are allowed.In addition to the reference FASTA file, DISCOVAR also requires the associated index file ( .fai

DISCOVAR can currently de novo assemble small genomes (up to 50 Mb), with larger genome support

to come soon.

The syntax for DISCOVAR de novo assembly is:
Discovar READS= bamfilenames \
OUT_HEAD= outputfilename \
REGIONS=all

This will take as input all the reads in the BAM file reads.bam , generate an assembly, then write the
output to a set of files prefixed with assembly

安装:

安装命令(在安装包里执行):
$ CC=/opt/centos/devtoolset-1.1/root/usr/bin/gcc CXX=/opt/centos/devtoolset-1.1/root/usr/bin/c++ ./configure
$ make -j 32
$ sudo make install
 
GCC 4.7安装在/opt/centos/devtoolset-1.1,安装办法参考:
 
如果是双端测序, 运行命令:

DiscovarDeNovo

NUM_THREADS=30(线程数)

MAX_MEM_GB=300G(最大内存)

MEMORY_CHECK=False (加上这个没报错)

READS=Project_TongJi_DNAseq_THB/Sample_TongJi-DNA-1/TongJi-DNA-1_CTTGTA_L000_R1.fastq.gz,Project_TongJi_DNAseq_THB/Sample_TongJi-DNA-1/TongJi-DNA-1_CTTGTA_L000_R2.fastq.gz(要用都好隔开)

OUT_DIR=wortdic(结果所在目录)

 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 

DISCOVAR de novo的更多相关文章

  1. De novo 测序基础知识

    名词解释 De novo:拉丁文,从头开始的意思,de nove测序则是指在不需要任何参考序列的情况下对某一物种进行基因组测序,然后将测得的序列进行拼接.组装,从而绘制该物种的全基因组序列图谱. 重测 ...

  2. (转)8 reviews about de novo genome assembly

    转自:http://dskernel.blogspot.com/2012/04/8-reviews-about-de-novo-genome-assembly.html 8 reviews about ...

  3. De novo RNA-Seq Assembly Using De Bruijn Graphs

    De novo RNA-Seq Assembly Using De Bruijn Graphs  2017-06-12 09:42:47     59     0     0 在说基因组的拼接之前,可 ...

  4. 全基因组测序 从头测序(de novo sequencing) 重测序(re-sequencing)

    全基因组测序 全基因组测序分为从头测序(de novo sequencing)和重测序(re-sequencing). 从头测序(de novo)不需要任何参考基因组信息即可对某个物种的基因组进行测序 ...

  5. MCP|ZWT|Precision de novo peptide sequencing using mirror proteases of Ac-LysargiNase and trypsin for large-scale proteomics(基于Ac-LysargiNase和胰蛋白酶的蛋白组镜像de novo测序)

    一.概述 由于难以获得100%的蛋白氨基酸序列覆盖率,蛋白组de novo测序成为了蛋白测序的难点,由Ac-LysargiNase(N端蛋白酶)和胰蛋白酶构成的镜像酶组合可以解决这个问题并具有稳定性, ...

  6. chromosome interaction mapping|cis- and trans-regulation|de novo|SRS|LRS|Haplotype blocks|linkage disequilibrium

    Dissecting evolution and disease using comparative vertebrate genomics-The sequencing revolution   s ...

  7. HHP|HPLC-MS/MS|PMT|PST|de novo|

    生物医学大数据 Protein 应用 人类蛋白质组计划 Gene的存在要依靠在蛋白水平确认基因真实存在. 蛋白质组是确定时间地点的研究单元的蛋白质总体,因为时间.地点和研究单元的相互组合存在多种变化, ...

  8. Uncovering thousands of new peptides with sequence-mask-search hybrid de novo peptide sequencing framework (使用序列掩码搜索结合肽段从头测序框架发现了数千个新肽段)-解读人:刘佳维

    期刊名:Molecular & Cellular Proteomics 发表时间:(2019年12月) IF:4.828 单位: 朱拉隆功大学 费城威斯塔研究所 物种:人 技术:de novo ...

  9. 基因组 de novo 组装原理

    Falcon软件的组装流程 为了错误校正,将原始子reads进行overlap 预组装和错误校正 错误校正后reads的overlap检测 overlap的过滤 从overlap构建图 从图构建con ...

随机推荐

  1. linux percona-toolkit的安装

    percona-toolkit是一套linux下的工具集,我们需要使用到 pt-query-digest 对mysql慢日志查询来做分析. 1.下载 http://www.percona.com/do ...

  2. jQuery EasyUI Combobox无法检索中文输入的问题

    在项目里使用了EasyUI的Combobox,当ComboBox的item是英文时,都能正常检索出对应项,但是如果使用中文输入法输入几个字母然后通过按shift键输入时,奇怪的事情发生了,combob ...

  3. POJ 1719 二分图最大匹配(记录路径)

    Shooting Contest Time Limit: 1000MS   Memory Limit: 10000K Total Submissions: 4097   Accepted: 1499 ...

  4. UVa 10870 - Recurrences

    http://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&page=show_problem&p ...

  5. mouseOver与rollOver

    区别: 当父容器监听这两个事件,鼠标从父容器移到子容器再移回父容器时,会触发mouseOver.mouseout事件,但是不会触发rollover.rollout事件.

  6. Android的R.java文件

    1.Android资源管理简介: Android应用程序资源可以分为两大类,分别放在assets和res文件夹下.assets目录下保存的是一些原始的文件,可以以任何方式来进行组织.这些文件最终会被原 ...

  7. [pjsip]Pjlib中的链表结构

    Pjlib的链表结构跟常见的链表结构有所区别,如下图所示: 图1:一般链表结构 图2:pjlib中的链表结构 可以看到一般的双向链表是链表节点包含数据域,而pjlib中是数据域包含链表节点.一般的链表 ...

  8. Windows Azure 实操 —— 迁移本地SharePoint服务器到Azure

    博客地址 http://blog.csdn.net/foxdave 注意:如果你是第二代虚拟机,那就别看这个了,老老实实在Azure上重新创建吧,Azure不支持第二代虚拟机. 写在之前,对Azure ...

  9. 使用WebClient Post方式模拟上传文件和数据

    假如某网站有个表单,例如(url: http://localhost/login.aspx):帐号  密码 我们需要在程序中提交数据到这个表单,对于这种表单,我们可以使用 WebClient.Uplo ...

  10. Centos6升级内核2.6到3.x过程

    最近公司有一个应用,安装需要内核版本3.1以后,不得已,需要升级下内核版本: 1. 安装必要依赖 # yum groupinstall "Development Tools" #y ...