论文
Plastid Genome Assembly Using Long-read data
https://onlinelibrary.wiley.com/doi/epdf/10.1111/1755-0998.13787
还有一篇对应的论文
https://www.biorxiv.org/content/10.1101/2023.06.27.546741v1.full.pdf
Phylogenetic placement of Ceratophyllum submersum based on a complete plastome sequence derived from nanopore long read sequencing data
对应的github主页
https://github.com/Bean061/ptgaul
软件直接用conda安装就可以
使用也非常简单
ptGAUL.sh -r /path/Beta.fasta -l ont.fq -t 8 -f 3000 -o ./ptgaul/
需要指定一个参考序列，github主页上写到这个参考序列来自同属或者同科的叶绿体基因组序列都可以
软件运行非常快，基本上10分钟就能拿到组装结果，手头有数据的可以试试
这个流程首先使用minimap2将ont数据比对到参考的叶绿体基因组，然后生成paf文件，然后根据paf文件去过滤比对上的reads，然后选取50X的reads去组装，组装用到的软件是flye,我试了一下我自己手头的数据，初次比对然后过滤比对到参考基因组序列的reads大约都在1000X左右，我自己手头的ont数据文库长度应该是在20k左右，之前记得ont的8k文库好像很便宜了，不知道用8k文库的数据来组装叶绿体基因组效果怎么样
学习一下这个流程的代码，主要学习怎么根据paf格式的文件提取比对到的reads
推文记录的是自己的学习笔记，很可能存在错误，请大家批判着看
