diff options
| -rwxr-xr-x | SelectiveGenomeAmplification | 4 | ||||
| -rwxr-xr-x | src/score_mers.py | 2 |
2 files changed, 3 insertions, 3 deletions
diff --git a/SelectiveGenomeAmplification b/SelectiveGenomeAmplification index 4b3c8dc..637ff2c 100755 --- a/SelectiveGenomeAmplification +++ b/SelectiveGenomeAmplification @@ -323,8 +323,8 @@ if [[ -n "$step_score" ]] || [[ -n "$all" ]]; then # output our sorted scores echo "sorting and outputting top $output_top_nb scores" echo "top scores output file: $output_directory/$current_run/top-scores" - head -n 1 $output_directory/$current_run/all-scores > $output_directory/$current_run/top-scores - tail -n +2 $output_directory/$current_run/all-scores | sort -r -t $'\t' -nk 3 | head -n $output_top_nb >> $output_directory/$current_run/top-scores + head -n 2 $output_directory/$current_run/all-scores > $output_directory/$current_run/top-scores + tail -n +3 $output_directory/$current_run/all-scores | sort -r -t $'\t' -nk 3 | head -n $output_top_nb >> $output_directory/$current_run/top-scores fi diff --git a/src/score_mers.py b/src/score_mers.py index 13050e7..fc1e81c 100755 --- a/src/score_mers.py +++ b/src/score_mers.py @@ -190,8 +190,8 @@ def percentage(part, whole, precision=2): return str(percent) + "%" def write_header(fh): + fh.write("# variables used: max_select=" + str(max_select) + " max_check=" + str(max_check) + " max_mer_distance=" + str(max_mer_distance) + " max_consecutive_binding=" + str(max_consecutive_binding) + " primer_weight=" + str(primer_weight) + "\n") fh.write("nb_primers\tCombination\tScore\tFG_mean_dist\tFG_stdev_dist\tBG_ratio\n") - def write_result(fh, score_res): combination, score_val, fg_mean_dist, fg_stddev_dist, bg_ratio = score_res fh.write(str(len(combination)) + "\t") |