AWK 数组

liaoboda

比如
B  <b>Carbohydrate Metabolism</b>
C    00010 Glycolysis / Gluconeogenesis [PATH:ko00010]
D      K00844  HK; hexokinase [EC:2.7.1.1]
D      K12407  GCK; glucokinase [EC:2.7.1.2]
D      K00845  glk; glucokinase [EC:2.7.1.2]
D      K01810  GPI, pgi; glucose-6-phosphate isomerase [EC:5.3.1.9]
D      K06859  pgi1; glucose-6-phosphate isomerase, archaeal [EC:5.3.1.9]
D      K13810  tal-pgi; transaldolase / glucose-6-phosphate isomerase [EC:2.2.1.2 5.3.1.9]
D      K00850  PFK, pfk; 6-phosphofructokinase [EC:2.7.1.11]
B  <b>HEAD Metabolism</b>
C    00020 Citrate cycle (TCA cycle) [PATH:ko00020]
D      K01647  CS, gltA; citrate synthase [EC:2.3.3.1]
D      K01648  ACLY; ATP citrate (pro-S)-lyase [EC:2.3.3.8]
D      K15230  aclA; ATP-citrate lyase alpha-subunit [EC:2.3.3.8]
D      K15231  aclB; ATP-citrate lyase beta-subunit [EC:2.3.3.8]
D      K01681  ACO, acnA; aconitate hydratase 1 [EC:4.2.1.3]
D      K01682  acnB; aconitate hydratase 2 / 2-methylisocitrate dehydratase [EC:4.2.1.3 4.2.1.99]
D是C的子集,D是C的子集. 我输入一个D中的K****的号 要找出对应的C和B.

jason680

回复 1# liaoboda


$ awk -vk="K01647" '$1~/B/{B=$0}$1~/C/{C=$0}$1~/D/{if($2==k){print B"\n"C"\n"$0;exit}}' bcd.txt
B  <b>HEAD Metabolism</b>
C    00020 Citrate cycle (TCA cycle) [PATH:ko00020]
D      K01647  CS, gltA; citrate synthase [EC:2.3.3.1]

liaoboda

多谢 多谢!