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

  四个环境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 &

©著作权归作者所有,转载或内容合作请联系作者
【社区内容提示】社区部分内容疑似由AI辅助生成,浏览时请结合常识与多方信息审慎甄别。
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

相关阅读更多精彩内容

友情链接更多精彩内容