云之南

风声，雨声，读书声，声声入耳；家事，国事，天下事，事事关心

日志

关于我

云之南

专业背景：计算机科学研究方向与兴趣: JavaEE-Web软件开发, 生物信息学, 数据挖掘与机器学习, 智能信息系统目前工作: 基因组, 转录组, NGS高通量数据分析, 生物数据挖掘, 植物系统发育和比较进化基因组学

文章分类

LOFTER精选

八招诀窍，教你实力撩妹 >

网易考拉推荐

新一代测序技术（NGS）组装简介三

2010-01-01 11:17:48| 分类：生物信息学 | 标签： |举报 |字号大中小订阅

下载LOFTER 我的照片书 |

MIRA-related(An Automated Genome and EST Assembler, 作者：Bastien Chevreux 最新版本：V2.9.46, 22/06 2009）

MIRA是一个多通道（multipass）组装软件，主要用于基因组或EST（表达序列标签）数据. 它能够处理各种混合型数据，如Solexa, 454, 3730等，并且它能用于检测重复序列和SNPs；

但MIRA需要比较长的时间运行，尤其你在处理mate pair（双末端测序序列）信息时；

MIRA采取逐步递进式精化组装的contigs。在每条consensus序列中的每个碱基，它都会产生一个repeat ration的柱状图，并给予打分和标记；

MIRA采用了一个跟DNPTrapper处理SNP的相似策略；

MIRA会变得混乱，如果你的输入的DNA序列是来自不同个体的 - 他们会有着很大不同的SNPs位点。此时你则必须要告诉MIRA你的DNA序列是来自不同个体的基因组序列。

一些意见：

用AMOS检验

预先筛选基因组的n-mer重复序列. 这将给你的组装带来期望的结果。

预先筛选您的序列为你预料到目前和屏蔽这些重复序列。

用Blast比对你的contigs - 看结果是否有意义的?

用Blast比对你的singlets - 看看是不是你所期望看到的(e.g. sequences that translate into missing elements of pathways?)

Assembly and analysis组装软件

gsAssembler - aka Newbler, from Roche. Version 2.3 was just released before the workshop. GenePool users can download the software and manuals.

Velvet

Curtain

MIRA

Abyss

The R libraries ShortRead and Rolexa

Minimus2

Phrap

Celera assembler and CABOG

Bambus - plus a document on how to use Bambus with MIRA

Visualising assemblies组装可视化软件

gsAssembler - aka Newbler, has a graphical interface.

Hawkeye

Tablet

Consed

Gap5

CLCBio

EagleView

Other software其他相关应用软件

simhtsd - Given a reference sequence, simhtsd creates a large set of short nucleotide reads, simulating the output from high throughput DNA sequencers such as the Illumina Genome Analyzer II.

MAQ simulate - another way to simulate short nucleotide reads. Part of the MAQ suite.

Make - a utility program, usually used for specifying the building of executable programs. Also suggested as a useful

scripting tool for running pipelines such as assembly pipelines.

Maker -a genome annotation pipeline for smaller eukaryotic and prokaryotic genome projects to annotate their

genomes and to create genome databases.

RAST annotation server

我自己也写了一个模拟序列程序，可以模拟DNA序列，并同时考虑有variation(SNPS, Indels, SVs)，甲基化(Methylation)，测序错误，各种NGS形式的序列，单末端或双末端，BAC序列，加adpater，primier，vector等--simulateSeq (v0.1.9，但尚未发布）

后续工作：

Outcomes

Script everything. This helps cut down unnecessary work, and aids in reproducibility and accountability. One tool mentioned as particularly useful in this regard is Make, as it understands dependencies, and can be used to avoid re-running things that have already been done, as well as keeping a record of steps taken.

Share scripts. A discussion was held about how and where a useful script repository could be set up. Tools such as trac and subversion were mentioned. This is probably a topic to be followed up on. In the meantime, the NERC Environmental Bioinformatics Centre has a webpage where utility scripts for handling new sequencing data types are being recorded. People are invited to submit links to their scripts, or the scripts themselves if they do not have capacity to host them locally, to be added to this webpage, until an alternative script repository is established.

Use metrics. Mark Blaxter prepared an overview of the types of metrics discussed and trialled during the workshop. These are a place to start when assembling sequence data:

Base level metrics

high length of contigs (and full length relative to expected genome size)

high N50 few contigs in N50

many contigs over 1 kb over 10kb

longest contig

All of these need to be balanced with quality of the assembly of course.

Comparative assembly metrics

eveness of coverage of contigs

same assembly achieved with several parameter sets

same assembly achieved with different programs/algorithms

same assembly achieved with different data sets.

Read distribution metrics

low propoation of reads rejected as singletons (genome only)

even coverages of reads over assembly (congurent with expected fold coverage)

correct spacing of paired reads in assembly

Biological affirmation

synteny with related genome

breaks/contig ends map to likely repetitive elements (Tn)

congruent with transcriptome data (exon mapping)

congruent with restriction map (optical) map or genetic map

源自：

Summary report on the Next Generation Sequcen Assembly workshop and NextGenBug meeting held from

Nov 30 - Dec 2, 2009 in Edinburgh.

Next Generation Sequencing Assembly Workshop

eScience Institute, Edinburgh, November 30 - December 2, 2009

Full program and list of presenters. All of the presentations will be available on the eScience Institute website.

This workshop was funded by the Scottish Bioinformatics Forum, the National eScience Institute and the GenePool

评论这张

转发至微博

阅读(1754)| 评论(0)

推荐转载

历史上的今天

在LOFTER的更多文章

this.p={  m:2,
              b:2,
              loftPermalink:'',
              id:'fks_083070082086089068081085084095085094087074093087083071',
              blogTitle:'新一代测序技术（NGS）组装简介 三',
              blogAbstract:'<P style=\"TEXT-INDENT: 2em;\"  \></P\>  <P style=\"TEXT-INDENT: 2em;\"  \>MIRA-related(An Automated Genome and EST Assembler, 作者：<A rel=\"nofollow\" href=\"http://mailto:%EF%BB%BFbach@chevreux.org/\"  \>Bastien Chevreux</A\> 最新版本：V2.9.46, 22/06 2009）</P\>  <P style=\"TEXT-INDENT: 2em;\"  \>MIRA是一个多通道（multipass）组装软件，主要用于基因组或EST（表达序列标签）数据. 它能够处理各种混合型数据，如Solexa, 454, 3730等，并且它能用于检测重复序列和SNPs；</P\>  <P style=\"TEXT-INDENT: 2em;\"  \>但MIRA需要比较长的时间运行，尤其你在处理mate pair（双末端测序序列）信息时；</P\>',
              blogTag:'',
              blogUrl:'blog/static/18699154201001111748535',
              isPublished:1,
              istop:false,
              type:0,
              modifyTime:1316912665685,
              publishTime:1262315868535,
              permalink:'blog/static/18699154201001111748535',
              commentCount:0,
              mainCommentCount:0,
              recommendCount:0,
              bsrk:-100,
              publisherId:0,
              recomBlogHome:false,
              currentRecomBlog:false,
              attachmentsFileIds:[],
              vote:{},
              groupInfo:{},
              friendstatus:'none',
              followstatus:'unFollow',
              pubSucc:'',
              visitorProvince:'',
              visitorCity:'',
              visitorNewUser:false,
              postAddInfo:{},
              mset:'000',
              mcon:'',
              srk:-100,
              remindgoodnightblog:false,
              isBlackVisitor:false,
              isShowYodaoAd:true,
              hostIntro:'专业背景：计算机科学                           \n\n研究方向与兴趣: JavaEE-Web软件开发,\n生物信息学, 数据挖掘与机器学习, 智能信息系统                                                                \n目前工作: 基因组, 转录组, NGS高通量数据分析, 生物数据挖掘, 植物系统发育和比较进化基因组学',
              hmcon:'1',
              selfRecomBlogCount:'0',
              lofter_single:'<a hidefocus="true" target="_blank" href="http://inpin.blog.163.com/blog/static/2633040402016620112520219/"><img onload="var nimg = new Image(); nimg.src=\'http://blog.163.com/newpage/images/analyse.png?lofter_single&t=\'+new Date().getTime()" border="0" src="http://imglf.nosdn.127.net/img/TlZCWHdFWFpidkQ1dGtJK1N5cFRBNzBLUFFPYnJEVi9LWk1KTzJlMTU2Yz0.jpg"/></a>'
            }

{list a as x}
    {if !!x}
    <div class="iblock nbw-fce nbw-f40">
      <a class="fc03 noul" target="_blank" hidefocus="true" href="http://blog.163.com/${x.visitorName}/">
      {if x.visitorName==visitor.userName}
      <img alt="${x.visitorNickname|escape}" onerror="this.src=location.f40" class="cwd bdwa bdc0" src="${fn1(x.visitorName)}&r=${visitor.imageUpdateTime}"/>
      {else}
      <img alt="${x.visitorNickname|escape}" onerror="this.src=location.f40" class="cwd bdwa bdc0" src="${fn1(x.visitorName)}"/>
      {/if}
      </a>
      <div class="cwd vname thide">
        {if x.moveFrom=='wap'}
          <a class="noul pnt" target="_blank" href="http://blog.163.com/services/wapblog.html?frompersonalbloghome"><span title="来自网易手机博客" class="iblock wapIcon"> </span></a>
        {elseif x.moveFrom=='iphone'}
          <a class="noul pnt" target="_blank"><span title="来自iPhone客户端" class="iblock iphoneIcon"> </span></a>
        {elseif x.moveFrom=='android'}
          <a class="noul pnt" target="_blank"><span title="来自Android客户端" class="iblock androidIcon"> </span></a>
        {elseif x.moveFrom=='mobile'}
          <a class="noul pnt" target="_blank" href="http://blog.163.com/services/emsblog.html?frompersonalbloghome"><span title="来自网易短信写博" class="iblock wapIcon"> </span></a>
        {/if}
        <a class="fc03 m2a"  target="_blank" hidefocus="true" href="http://blog.163.com/${x.visitorName}/">
          ${fn(x.visitorNickname,8)|escape}
        </a>
      </div>
    </div>
    {/if}
    {/list}

<#--最新日志，群博日志--> <#--推荐日志-->

<p class="fc06">推荐过这篇日志的人：</p>
    <div>
      {list a as x}
      {if !!x}
      <div class="iblock nbw-fce nbw-f40">
        <a class="fc03 noul" target="_blank" hidefocus="true" href="http://blog.163.com/${x.recommenderName}/">
        <img alt="${x.recommenderNickname|escape}" onerror="this.src=location.f40" class="cwd bdwa bdc0" src="${fn1(x.recommenderName)}"/>
        </a>
        <div class="cwd thide">
          <a class="fc03 m2a" target="_blank" hidefocus="true" href="http://blog.163.com/${x.recommenderName}/">
            ${fn(x.recommenderNickname,6)|escape}
          </a>
        </div>
      </div>
      {/if}
      {/list}
    </div>
    {if !!b&&b.length>0}
    <p  class="fc06">他们还推荐了：</p>
    <ul>
    {list b as y}
      {if !!y}
        <li class="rrb"><span class="iblock">·</span><a class="fc03 m2a" target="_blank" href="http://blog.163.com/${y.recommendBlogPermalink}/?from=blog/static/18699154201001111748535">${y.recommendBlogTitle|escape}</a></li>
      {/if}
    {/list}
    </ul>
    {/if}

<#--引用记录--> <#--博主推荐--> <#--随机阅读--> <#--首页推荐--> <#--历史上的今天--> <#--被推荐日志--> <#--上一篇，下一篇--> <#-- 热度 -->

{list a as x}
    {if !!x}
    <div class="hotItem iblock nbw-fce nbw-f40">
      <a class="fc03 noul" target="_blank" hidefocus="true" href="http://blog.163.com/${x.publisherUsername}/">
      {if x.publisherUsername==visitor.userName}
      <img alt="${x.publisherNickname|escape}" onerror="this.src=location.f40" class="cwd bdwa bdc0" src="${fn1(x.publisherUsername)}&r=${visitor.imageUpdateTime}"/>
      {else}
      <img alt="${x.publisherNickname|escape}" onerror="this.src=location.f40" class="cwd bdwa bdc0" src="${fn1(x.publisherUsername)}"/>
      {/if}
      </a>
      <div class="cwd vname thide">
        <a class="fc03 m2a"  target="_blank" hidefocus="true" href="http://blog.163.com/${x.publisherUsername}/">
          ${fn(x.publisherNickname,8)|escape}
        </a>
      </div>
      <a class="f-myLikeIcons hottype {if x.type==1} js-liketype{elseif x.type==2} js-reblogtype{elseif x.type==3} js-sharetype{else}{/if}" target="_blank" hidefocus="true" href="http://blog.163.com/${x.publisherUsername}/"> </a>
    </div>
    {/if}
    {/list}

<#-- 网易新闻广告 --> <#--右边模块结构--> <#--评论模块结构--> <#--引用模块结构--> <#--博主发起的投票-->

页脚

我的照片书 - 博客风格 - 手机博客 - 下载LOFTER APP - 订阅此博客

云之南

导航