blast结果的m7格式转换为m8格式

以下代码由PHP站长网 52php.cn收集自互联网

现在PHP站长网小编把它分享给大家，仅供参考

while (<>)
{
   if (/Iteration_query-def/)
   {
   />(.*?)</;
   $out1=$1;
   }
   if (/Hit_def/)
   {
    />(.*?)</;
   $out2=$1;
   }
   if (/Hsp_bit-score/)
   {
    />(.*?)</;
   $out12=$1;
   }
   if (/Hsp_evalue/)
   {
   />(.*?)</;
   $out11=$1;
   }
   if (/Hsp_query-from/)
   {
     />(.*?)</;
   $out7=$1;
   }
   if (/Hsp_query-to/)
   {
   />(.*?)</;
   $out8=$1;
   }
   if (/Hsp_hit-from/)
   {
    />(.*?)</;
   $out9=$1;
   }
   if (/Hsp_hit-to/)
   {
    />(.*?)</;
   $out10=$1;
   }
   if (/Hsp_align-len/)
   {
   />(.*?)</;
   $out4=$1;
   }
   if (/Hsp_gaps/)
   {
    />(.*?)</;
   $out6=$1;
   }
     if (/Hsp_identity/)
   {
    />(.*?)</;
   $tmp=$1;
   }

   if (/Hsp_qseq/)
   {
    $out3=100*$tmp/$out4;
    $out3=sprintf"%0.2f",$out3;
    $out5=$out4-$tmp-$out6;
    print "$out1t$out2t$out3t$out4t$out5t$out6t$out7t$out8t$out9t$out10t$out11t$out12n";
   }

}

以上内容由PHP站长网【52php.cn】收集整理供大家参考研究

如果以上内容对您有帮助，欢迎收藏、点赞、推荐、分享。

（编辑：李大同）

【声明】本站内容均来自网络，其相关言论仅代表作者个人观点，不代表本站立场。若无意侵犯到您的权利，请及时与联系站长删除相关内容!