#!/usr/bin/perl -w # Kevin Karplus # 24 May 2002 # make_undertaker_alignment_list < top_reported_alignments.rdb # > target-undertaker-align.script # generates a script for undertaker to use the top alignments as # sources for specific fragments and alignments. # # The script is NOT a complete undertaker script---it contains only # the ReadFragmentAlignment commands to read in the a2m files. # Setting up the target, selecting the training chains to use for # a generic fragment library, selecting a clash definition, ... # all must be done first, before doing an Include command to read # in the fragments. # BUGS: # no usage message, almost no error checking. # CHANGE LOG: # Sun Apr 9 14:02:05 PDT 2006 Kevin Karplus # Now extracing e_values as wekk $_=; die "Header missing\n" if ! /Alignment/; ; # skip the RDB header while () { @cols = split(/\t/); print "ReadFragmentAlignment $cols[0] NOFILTER SCWRL force_alignment e_value $cols[1]\n"; }