泛基因集构建中JCVI做共线性分析构建完一半数据报错 core dumped

for i in `cat stepmap_ID.txt`;do

  # JCVI做共线性分析,--cpus 5 限制计算资源使用,如果不加默认--cpus 56

  python3 -m jcvi.compara.catalog ortholog \

   pangene ${i} --cpus 5  --no_strip_names --cscore 0.8;

  # 以提升后的anchors作为基因对的结果,把比对上的基因排序去重输出到mapped.txt文件

  anchor=pangene.${i}.lifted.anchors;

  awk '{ print $2}' $anchor | sort -u > ${i}.mapped.txt;

  # 未能比对上的基因合并到模板里,组成新的模板,再和下一个基因组比对

  awk '{if(ARGIND==1) {val[$0]}else{if(!($4 in val)) print $0}}' ${i}.mapped.txt ${i}.bed >> pangene.bed ;

  awk '{if(ARGIND==1) {val[$0]}else{if(!($4 in val)) print $4}}' ${i}.mapped.txt ${i}.bed |seqkit grep -f - ${i}.cds >> pangene.cds;

done

attachments-2025-03-EkTa3xhT67ea2739da1d1.png

请先 登录 后评论

2 个回答

last报错应该是因为内存不足:https://www.omicsclass.com/article/1413

请先 登录 后评论
omicsgene - 生物信息
擅长:重测序,遗传进化,转录组,GWAS

core dumped 可能是电脑系统软件编译问题,导致运行不了,建议换个电脑试试的;

或者到我们云服务器上分析数据;

请先 登录 后评论
  • 2 关注
  • 0 收藏,56 浏览
  • 落日与晚风 提出于 2天前

相似问题