Skip to content

Commit

Permalink
updating the nucleotide call to fasttree to include -nt and -gtr
Browse files Browse the repository at this point in the history
  • Loading branch information
Mike Lee authored and Mike Lee committed Oct 7, 2024
1 parent 824be8e commit 00d4d13
Showing 1 changed file with 10 additions and 2 deletions.
12 changes: 10 additions & 2 deletions bin/GToTree
Original file line number Diff line number Diff line change
Expand Up @@ -3436,7 +3436,11 @@ if [ $align_only == 'false' ]; then
# setting env variable for FastTreeMP
export OMP_NUM_THREADS=${num_jobs}

FastTreeMP ${alignment_file} > ${tree_file} 2> >(tee -a $gtotree_log >&2)
if [ $nucleotide == 'false' ]; then
FastTreeMP ${alignment_file} > ${tree_file} 2> >(tee -a $gtotree_log >&2)
else
FastTreeMP -nt -gtr ${alignment_file} > ${tree_file} 2> >(tee -a $gtotree_log >&2)
fi

fasttree_used="true" # setting to report citations at end

Expand All @@ -3454,7 +3458,11 @@ if [ $align_only == 'false' ]; then
printf " It is currently $curr_time; the process started at $start_time.\n" | tee >( sed 's/\x1b\[[0-9;]*m//g' >> ${gtotree_log} )
printf " Current process runtime: $(($duration / 60 / 60)) hours and $((($duration / 60) % 60)) minutes.\n\n" | tee >( sed 's/\x1b\[[0-9;]*m//g' >> ${gtotree_log} )

FastTree ${alignment_file} > ${tree_file} 2> >(tee -a $gtotree_log >&2)
if [ $nucleotide == 'false' ]; then
FastTree ${alignment_file} > ${tree_file} 2> >(tee -a $gtotree_log >&2)
else
FastTree -nt -gtr ${alignment_file} > ${tree_file} 2> >(tee -a $gtotree_log >&2)
fi

fasttree_used="true" # setting to report citations at end

Expand Down

0 comments on commit 00d4d13

Please sign in to comment.