T2T基因组组装,去除自冗余,发现是minimap2 -D -cx asm5 -t 10 tmp1.fa tmp1.fa > tmp1.self.paf 命令产生的tmp1.self.pa文件中有各别行输出结果有问题(17[M::worker_pipeline::1636.730*6.86] mapped 12 sequences),后续awk命令报错。如下图

attachments-2025-08-HD20m9vX68ac0c8425ecd.jpgattachments-2025-08-NxyF6Ym468ac0cac67dfd.jpg


minimap2 -D -cx asm5  -t 10 tmp1.fa tmp1.fa  > tmp1.self.paf

attachments-2025-08-9GPYB02q68ac10a039e68.jpg

请先 登录 后评论

最佳答案 1天前

因为你nohup这个命令,但是没有把这个命令写成sh的脚本去运行,这样log文件会输出到你的结果文件里面。命令写成sh文件,再提交运行。得到新得结果。

请先 登录 后评论

其它 0 个回答

  • 0 关注
  • 0 收藏,75 浏览
  • 青春无悔 提出于 2天前

相似问题