https://www.genome.jp/kegg-bin/get_htext?hsa00001+3101
C開(kāi)頭的就是kegg的pathway的ID所在行，D開(kāi)頭的就是屬于它的kegg的所有的基因
perl -alne '{if(/^C/){/PATH:hsa(\d+)/;$kegg=$1}else{print "$kegg\t$F[1]" if /^D/ and $kegg;}}' hsa00001.keg >kegg2gene.txt
++++++++++++++++++++++++++++++++
#!usr/bin/perluse warnings;use strict; my ($path, $num);open IN, 'hsa00001.keg';open OUT, '>kegg_sorting'; while (){
? chomp;
? ? if (/^C/){
? ? ? ? ($num)=$_=~/C\s*(\d+).*/;
? ? ? ? #print OUT"$num\n";
? ? }
? ? elsif(/^D/){
? ? ? ? ($path)=$_=~/D\s+(\d+).*/;
? ? ? ? print OUT "$num\t$path\n";
? ? }
}
