Skip to content

Commit 7aaf5c6

Browse files
committed
Help update
1 parent 36557be commit 7aaf5c6

File tree

3 files changed

+10
-2
lines changed

3 files changed

+10
-2
lines changed

README.md

Lines changed: 4 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -54,6 +54,9 @@ cd FAMSA && make
5454

5555
# export a pairwise identity (PID) matrix to the CSV format (square)
5656
./famsa -dist_export -pid -square_matrix ./test/adeno_fiber/adeno_fiber pid.csv
57+
58+
# profile-profile alignment without refining output
59+
./famsa -refine_mode off ./test/adeno_fiber/upgma.no_refine.part1.fasta ./test/adeno_fiber/upgma.no_refine.part2.fasta pp.fasta
5760
```
5861

5962

@@ -120,6 +123,7 @@ Options:
120123
* `-keep-duplicates` - keep duplicated sequences during alignment (default: disabled - duplicates are removed prior and restored after the alignment)
121124
* `-gz` - enable gzipped output (default: disabled)
122125
* `-gz-lev <value>` - gzip compression level [0-9] (default: 7)
126+
* `-refine_mode <on | off | auto>` - refinement mode (default: `auto` - the refinement is enabled for sets <= 1000 seq.)
123127

124128

125129
### Guide tree import and export

run.sh

Lines changed: 3 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -20,3 +20,6 @@
2020

2121
# export pairwise identity (PID) matrix to CSV format (square)
2222
./famsa -dist_export -pid -square_matrix ./test/adeno_fiber/adeno_fiber pid.csv
23+
24+
# profile-profile alignment without refining output
25+
./famsa -refine_mode off ./test/adeno_fiber/upgma.no_refine.part1.fasta ./test/adeno_fiber/upgma.no_refine.part2.fasta pp.fasta

src/core/params.cpp

Lines changed: 3 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -89,12 +89,13 @@ void CParams::show_usage(bool expert)
8989
<< " (default: disabled - duplicates are removed prior and restored after the alignment).\n\n"
9090

9191
<< " -gz - enable gzipped output (default: " << bool2str[gzippd_output] << ")\n"
92-
<< " -gz-lev <value> - gzip compression level [0-9] (default: " << gzip_level << ")\n\n";
92+
<< " -gz-lev <value> - gzip compression level [0-9] (default: " << gzip_level << ")\n"
93+
<< " -refine_mode <on | off | auto> - refinement mode (default: auto - the refinement is enabled for sets <= " << thr_refinement << " seq.)\n\n";
94+
9395

9496
if (expert) {
9597
LOG_NORMAL << "Advanced options:\n"
9698
<< " -r <value> - no. of refinement iterations (default: " << n_refinements << ")\n"
97-
<< " -refine_mode <on|off|auto> - refinement mode (default: auto - the refinement is disabled for sets larger than " << thr_refinement << " seq.)\n"
9899
<< " -go <value> - gap open cost (default: " << gap_open << ")\n"
99100
<< " -ge <value> - gap extension cost (default: " << gap_ext << ")\n"
100101
<< " -tgo <value> - terminal gap open cost (default: " << gap_term_open << ")\n"

0 commit comments

Comments
 (0)