HETATMのresidue nameがない。。
あるソフトで分子構造をPDBに書き出すと、
HETATMレコードのresidue名がomitされてた。
別のプログラムでエラーになる。。
スペース決め打ちで分割して、XXXを挿入したら、
とりあえず読めた。
#!/usr/bin/perl -w if($#ARGV == -1){ printf(" Usage: %s <filename> \n", $0); exit -1; } open(INFILE, $ARGV[0]) || die " Can\'t open INFILE: $!\n"; while(<INFILE>) { if(/^HETATM?/) { @block = split(/\s{7}?/,$_); print $block[0], " XXX ", $block[1], $block[2]; }else{ print $_; } } close(INFILE); exit 0;