Skip to content

Commit

Permalink
bootstraping & DNA-mode
Browse files Browse the repository at this point in the history
  • Loading branch information
sinamajidian authored Dec 20, 2023
1 parent ebd11ad commit 644a4a4
Showing 1 changed file with 17 additions and 0 deletions.
17 changes: 17 additions & 0 deletions README.md
Original file line number Diff line number Diff line change
Expand Up @@ -112,6 +112,23 @@ read2tree --standalone_path marker_genes/ --output_path output --reads species3_
read2tree --standalone_path marker_genes/ --output_path output --merge_all_mappings --tree
```

### bootstraping

To have bootstrap values a metric for quality of internal nodes, you can run the following
```
thread=20
iqtree -T ${thread} -s output/concat_*_aa.phy -bb 1000
```
The `.phy` file is either `concat_sample_aa.phy` or `concat_merge_aa.phy` corresponding to single- or multi-species mode.

It is also possible to us [trimal](http://trimal.cgenomics.org/use_of_the_command_line_trimal_v1.2) for trimming msa `trimal -in <inputfile> -out <outputfile> -automated1`

### DNA-mode tree inferrence
For closely related species, the user can infer tree using MSA of nucleotide sequences.
```
thread=20
iqtree -T ${thread} -s output/concat_*_dna.phy
```

## Test example

Expand Down

0 comments on commit 644a4a4

Please sign in to comment.