清华大佬耗费三个月吐血整理的几百G的资源,免费分享!....>>>
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 "$out1\t$out2\t$out3\t$out4\t$out5\t$out6\t$out7\t$out8\t$out9\t$out10\t$out11\t$out12\n";
}
}