美文网首页系统发育基因组分析
对BUSCO的结果进行总结

对BUSCO的结果进行总结

作者: 小黑采蘑菇 | 来源:发表于2024-01-03 16:34 被阅读0次

    平时我们对每个基因组做了BUSCO,但是官方提供的画图代码很难用,所以决定自己提取出BUSCO的结果并进行画图,这一串代码就可以批量从BUSCO的结果中提取出我们想要的结果:

    # 提取数字并写入txt文件
    echo "voucher C S D F M" > busco_summary_table.txt;
    
    for i in $(cat list.txt); do
      filename="${i}/busco/short_summary.specific.agaricales_odb10.busco.txt"
      if [ -f "$filename" ]; then
          # 提取数据 
          data=$(grep -oP 'C:([\d.]+)%\[S:([\d.]+%),D:([\d.]+)%\],F:([\d.]+)%,M:([\d.]+)%,n:\d+' "$filename" | awk -F'[:,%[]' '{print $2, $5, $8, $11, $14}')
          # 将数据写入表格
          echo "${i} $data" >> busco_summary_table; else
          echo "${i} There is no BUSCO results in ${i}!" >> busco_summary_table;
      fi
    done
    

    相关文章

      网友评论

        本文标题:对BUSCO的结果进行总结

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