美文网首页泛基因组_从头学习_
对生成的bin进行合并、注释、评估

对生成的bin进行合并、注释、评估

作者: dashan1928 | 来源:发表于2022-04-17 17:41 被阅读0次

      四个环境gtdbtk、pangenome、minigraph、drep、metawrap-env

    gtdbtk可以对所有bin进行注释,可以看到bin属于什么物种

    gtdbtk classify_wf --pplacer_cpus 30 --cpus 120 --genome_dir bins_all/ --full_tree -x fa --out_dir gtdbtk

    minigraph中利用参考基因组,加上属于同一物种bin的集合,合成一个新的泛基因组

    minigraph -xggs -t16 ref.fa sample1.fa sample2.fa > out.gfa

    gfatools gfa2fa -s graph.gfa > out.stable.fa

    pangenome对gtdbtk中注释为相同物种的bin进行合并形成泛基因组

    panaroo -i gff/*.gff -t 40 -o strict_result --clean-mode strict

    panaroo -i gff/*.gff -t 40 -o sensitive_result --clean-mode sensitive

    drep可以发现有多少不重复的基因组

    dRep dereplicate -comp 70 -con 50  <outdir>  <inportdir( including all .fa)>

    metawrap-env的checkm可以对生成的基因组进行评估

    nohup checkm lineage_wf ./ ../checkm_4mox -t 30 --pplacer_threads 10 -x fa &

    相关文章

      网友评论

        本文标题:对生成的bin进行合并、注释、评估

        本文链接:https://www.haomeiwen.com/subject/aoidertx.html