#!/usr/bin/perl # convert start into .prot format for cyana $x=999.00; $error=0.010; $cnt=1; printf"#Atom_num Shift_ppm error_ppm Res_name Res_atom Res_num\n"; while(<>) { @wd = split ' '; if($wd[0] =~ /^ATOM/ && $wd[2]=~/^H/ ){ printf("%s\t %3.2f\t %1.2f\t %s\t %s\t %s\t %s\t",$cnt,$x,$error, $wd[3], $wd[2], $wd[5]); printf"\n"; $cnt++; } } close();