Trusted, Bureaucrats, Administrators
1,660
edits
Changes
Jump to navigation
Jump to search
protein 16 Q 16 N 3.410860e+06 5.823456e+06 3.427248e+06
protein 22 E 22 N 3.681818e+06 6.378086e+06 3.663392e+06
</source>
=== Select residues ===
<source lang="text">
set IN=table_ser_files_model.txt ;
set SELRESIS=sel_resi.txt ;
awk '{print $2}' $IN > $SELRESIS ;
gedit $SELRESIS
</source>
After this step, you should remove those lines with residues you dont want to have in the new file.
<source lang="text">
set IN=table_ser_files_model.txt ;
set SELRESIS=sel_resi.txt ;
set RESIS=`cat $SELRESIS` ;
set DATE=`date '+%Y%m%d_%H%M'` ;
set TMP=${IN}_${DATE} ;
cp $IN $TMP ;
rm $IN ;
set T=$ ;
echo $RESIS ;
foreach RESI ($RESIS) ;
cat $TMP | awk "{if (${T}2 == $RESI ) { print ${T}0; } }" >> $IN
end ;
cat $IN ;
relax_disp -g -l LOGFILE.txt table_ser_files_relax.py
</source>
[[Category:Peak_lists]]
