Skip to content

Commit

Permalink
fixed bug for hallmark taxonomy file
Browse files Browse the repository at this point in the history
  • Loading branch information
mtisza1 committed May 9, 2022
1 parent 83aa1f8 commit d9e6a82
Show file tree
Hide file tree
Showing 4 changed files with 11 additions and 11 deletions.
20 changes: 10 additions & 10 deletions README.md
Original file line number Diff line number Diff line change
Expand Up @@ -139,7 +139,7 @@ conda activate cenote-taker2_env
cd Cenote-Taker2
# run the update script:
`python update_ct2_databases.py --hmm True --protein True`
python update_ct2_databases.py --hmm True --protein True
```

## Schematic
Expand All @@ -152,15 +152,15 @@ Cenote-Taker 2 currently runs in a python wrapper.
Check environments:

`conda info --envs`
```
#Default:
Default:

`conda activate cenote-taker2_env`

Or if you've put your conda environment in a custom location:
conda activate cenote-taker2_env
`conda activate /path/to/better/directory/cenote-taker2_env`
#Or if you've put your conda environment in a custom location:
conda activate /path/to/better/directory/cenote-taker2_env
```
2. Run the python script to get the help menu (see options below).

```
Expand All @@ -175,7 +175,7 @@ python /path/to/Cenote-Taker2/run_cenote-taker2.py -h
```
python /path/to/Cenote-Taker2/run_cenote-taker2.py -c MY_CONTIGS.fasta -r my_contigs1_ct -m 32 -t 32 -p true -db virion
#Or, if you want to save a log of the run, add "2>&1 | tee output.log" to the end of the command):
#Or, if you want to save a log of the run, add "2>&1 | tee output.log" to the end of the command:
python /path/to/Cenote-Taker2/run_cenote-taker2.py -c MY_CONTIGS.fasta -r my_contigs1_ct -m 32 -t 32 -p true -db virion 2>&1 | tee output.log
```
Expand All @@ -197,7 +197,7 @@ python /path/to/Cenote-Taker2/run_cenote-taker2.py -c MY_VIRUSES.fasta -r viruse

For very divergent genomes, setting `-hh hhsearch` will marginally improve number of genes that are annotated. This setting increasese the run time quite a bit. On the other hand, setting `-hh none` will skip the time consuming hhblits step. With this, you'll still get pretty good genome maps, and might be most appropriate for very large virus genome databases, or for runs where you just want to do a quick check.

#### *Discovery*
#### *Discovery + Annotation*

**Virus-like particle (VLP) prep assembly:**

Expand All @@ -222,7 +222,7 @@ python /path/to/Cenote-Taker2/run_cenote-taker2.py -c MY_WGS_ASSEMBLY.fasta -r m
```


**Bacterial reference genome or MAG**
**Bacterial isolate genome or MAG**

`-p True -db virion --minimum_length_linear 3000 --lin_minimum_hallmark_genes 2`

Expand Down
2 changes: 1 addition & 1 deletion cenote-taker2.1.5.sh
Original file line number Diff line number Diff line change
Expand Up @@ -1601,7 +1601,7 @@ fi
cd ${base_directory}/${run_title}

if [ "$HALLMARK_TAX" == "True" ] && [ -d DTR_contigs_with_viral_domain/ ] ;then
HALLMARK_FILES=$( find DTR_contigs_with_viral_domain/ -maxdepth 1 -type f -name ".*rotate.AA.hmmscan.sort.out" | sed 's/\.\///g' )
HALLMARK_FILES=$( find DTR_contigs_with_viral_domain/ -maxdepth 1 -type f -name "*rotate.AA.hmmscan.sort.out" | sed 's/\.\///g' )
if [ -n "${HALLMARK_FILES}" ] ; then
MDYT=$( date +"%m-%d-%y---%T" )
echo "time update: reporting taxonomy for each hallmark gene, circular contigs " $MDYT
Expand Down
Binary file modified cenote_DNA_test_contigs1.png
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Binary file modified cenote_RNA_test_contigs1.png
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.

0 comments on commit d9e6a82

Please sign in to comment.