zoukankan      html  css  js  c++  java
  • Contig|scaffold|N50|L50|NG50|贪心算法|de bruiji graph|

    生物信息学

    Contigreads拼成的连续的DNA片段,连续表达一个gene。通过双端测序的contig可确定contig之间的关系得到scaffoldScaffoldreads拼成的有gapDNA片段。理想情况下,一条染色体用同一个scaffold的表达。整个genome存在很多零碎片段,可舍弃。因为duplication产生很多overlap

     

    N50L50NG50是评价genome assemblyquality的标准,评价长度时使用N50N50是一个contig的长度。不选用genome size50%是因为1.这是估计的size值不一定准;2.sequence 仅覆盖80%。评价数量使用L50L50数量越小越好。NG50表达测到genome 覆盖度。

     

     取材方法很重要,得到目标数据。

    Assembly算法有可能带来更多误差,通常二代测序和三代测序相结合。

     

     

    贪心算法原理是每一步都在找最优解,最后得到最好的结果,优点是快,缺点是不是全局最优解,出现重复序列便走不下去。

     

     

    de bruiji graph来自桥问题:

     

    比如:

    k-merk=4

    0000,成为00000000000,若下一步添0 成为0000

    0001,成为00010001001,若下一步添0 成为00100010010

    0001,成为00010001001,若下一步添1 成为00110011011

     

    在序列拼接时:

     

     

     Qualities取决于二代测序的质量;coverage&read length取决于建库方法和sequence方法。

    生物信息学处理关键是考虑研究的物种的特性,eg:某物种的duplication多,或者生物学问题的侧重点,eg:重测序。

    Genome网站:UCSC

    例子:

    植物类,希望通过“拟南芥vs抗逆抗旱植物”,清楚看到抗逆抗旱植物的相关特性,或者说希望看到的特性,做之前要估算genome size,可以看的现象是扩增基因,于是分析扩增基因的具体情况,比如对于常规部分的不同功能、通路及转录因子等,对于miRNA的探讨,分析出现的原因。

    以前的研究多重于分析生物学特征,比如某物种所有基因的罗列,现在基因组数量变大之后多研究.进化特征,比如某一个基因,eg:所有种类黄瓜的苦味,这种主观感觉定量分析的研究。

     
  • 相关阅读:
    Python文件的两种用途
    模块的搜索路径
    循环导入问题
    import和from...import
    模块的四种形式
    函数小结
    面向过程编程
    内置函数
    匿名函数
    递归
  • 原文地址:https://www.cnblogs.com/yuanjingnan/p/11725496.html
Copyright © 2011-2022 走看看