• 别人处理二代测序的流程


     做了基因组研究的summary,请各位大老帮忙看看,GCE评估180M(二倍体),杂合率为0.0119/21kmer,用ALLPATH-LG的ErrorCorrectReads.pl进行评估83M(单倍体),SNP率1/48,接着用canu进行组装,组装出165M的基因组,(1)直接进行FinisherSC和pilon进行升级,得150M基因组,用二代reads进行回帖,99.2%的回帖率,BUSCO 96.9%的完整率。
    (2)得到canu的结果后,使用HaploMerger2获得单倍体93M,接着用FinisherSC和pilon进行升级,得到91M的基因组~二代回帖92%,BUSCO评估95.8%

    我要先做比较基因组的东西,NCBI中已有上传的基因组为子实体测序(单倍体),我打算用(2)所得的单倍体基因组进行比较
    我们还测菌索,和菌丝的RNA-seq,为了兼顾杂合率,我们用的(1)所得基因组进行基因预测,然后使用预测的GTF文件进行定量差异表达分析,然后GO富集,KEGG富集等分析

    我有的疑问在于我用得到的单倍体基因组进行比较基因组学可以吗?
    我用二倍体的基因组进行RNA-seq,因为二倍体组装多少会把一些纯合位点合并,那么并不是所有基因都是二倍的,会对差异基因寻找,GO富集,KEGG,等有影响吗?
  • 相关阅读:
    .gitignore语法
    每日阅读
    css摘要
    ubuntu安装qq、微信
    django中views中方法的request参数
    js html标签select 中option 删除除了第一行外的其他行
    js 新增标签、标签属性
    python中None与0、Null、false区别
    python class中__init__函数、self
    for foreach循环
  • 原文地址:https://www.cnblogs.com/modaidai/p/9579637.html
Copyright © 2020-2023  润新知