@@ -94,25 +94,25 @@ gmsc-mapper createdb -i examples/target.faa -o examples/ -m mmseqs
94
94
- Input is genome contig sequences.
95
95
96
96
``` bash
97
- gmsc-mapper -i examples/example.fa -o examples_output/ --db examples/targetdb.dmnd --habitat examples/ref_habitat.txt --quality examples/ref_quality.txt --taxonomy examples/ref_taxonomy.txt
97
+ gmsc-mapper -i ./ examples/example.fa -o ./ examples_output/ --db ./ examples/targetdb.dmnd --habitat ./ examples/ref_habitat.npy --habitat-index ./examples/ref_habitat_index.tsv -- quality ./ examples/ref_quality.txt --taxonomy ./ examples/ref_taxonomy.npy --taxonomy-index ./examples/ref_taxonomy_index.tsv
98
98
```
99
99
100
100
- Input is amino acid sequences.
101
101
102
102
``` bash
103
- gmsc-mapper --aa-genes examples/example.faa -o examples_output/ --db examples/targetdb.dmnd --habitat examples/ref_habitat.txt --quality examples/ref_quality.txt --taxonomy examples/ref_taxonomy.txt
103
+ gmsc-mapper --aa-genes ./ examples/example.faa -o ./ examples_output/ --db ./ examples/targetdb.dmnd --habitat ./ examples/ref_habitat.npy --habitat-index ./examples/ref_habitat_index.tsv -- quality ./ examples/ref_quality.txt --taxonomy ./ examples/ref_taxonomy.npy --taxonomy-index ./examples/ref_taxonomy_index.tsv
104
104
```
105
105
106
106
- Input is nucleotide gene sequences.
107
107
108
108
``` bash
109
- gmsc-mapper --nt-genes examples/example.fna -o examples_output/ --db examples/targetdb.dmnd --habitat examples/ref_habitat.txt --quality examples/ref_quality.txt --taxonomy examples/ref_taxonomy.txt
109
+ gmsc-mapper --nt-genes ./ examples/example.fna -o ./ examples_output/ --db ./ examples/targetdb.dmnd --habitat ./ examples/ref_habitat.npy --habitat-index ./examples/ref_habitat_index.tsv -- quality ./ examples/ref_quality.txt --taxonomy ./ examples/ref_taxonomy.npy --taxonomy-index ./examples/ref_taxonomy_index.tsv
110
110
```
111
111
112
112
- Check another alignment tool: MMseqs2
113
113
114
114
``` bash
115
- gmsc-mapper -i examples/example.fa -o examples_output/ --db examples/targetdb --habitat examples/ref_habitat.txt --quality examples/ref_quality.txt --taxonomy examples/ref_taxonomy.txt --tool mmseqs
115
+ gmsc-mapper -i ./ examples/example.fa -o ./ examples_output/ --db ./ examples/targetdb --habitat ./ examples/ref_habitat.npy --habitat-index ./examples/ref_habitat_index.tsv -- quality ./ examples/ref_quality.txt --taxonomy ./ examples/ref_taxonomy.npy --taxonomy-index ./examples/ref_taxonomy_index.tsv --tool mmseqs
116
116
```
117
117
118
118
## Usage
@@ -280,14 +280,6 @@ The output folder will contain
280
280
281
281
* ` --quiet ` : Disable alignment console output. (default: False )
282
282
283
- * ` --db ` : Path to the GMSC database file. (default: ../db/targetdb.dmnd)
284
-
285
- * ` --habitat ` : Path to the habitat file. (default: ../db/ref_habitat.tsv.xz)
286
-
287
- * ` --taxonomy ` : Path to the taxonomy file. (default: ../db/ref_taxonomy.tsv.xz)
288
-
289
- * ` --quality ` : Path to the quality file. (default: ../db/ref_quality.tsv.xz)
290
-
291
283
### Subcommands and Parameters
292
284
Subcommands: ` gmsc-mapper createdb `
293
285
@@ -298,6 +290,3 @@ Subcommands: `gmsc-mapper createdb`
298
290
* ` -m/--mode ` : Alignment tool (Diamond / MMseqs2).
299
291
300
292
* ` --quiet ` : Disable alignment console output. (default: False )
301
-
302
- ## Sensitivity choices considering time and memory usage
303
- To be done
0 commit comments