Skip to content

Commit

Permalink
✨ Added a folder to put the genome by genome analysis to be more clear
Browse files Browse the repository at this point in the history
  • Loading branch information
rdenise committed Nov 28, 2023
1 parent 9cdf2eb commit f83bd66
Show file tree
Hide file tree
Showing 3 changed files with 27 additions and 26 deletions.
4 changes: 2 additions & 2 deletions Dockerfile
Original file line number Diff line number Diff line change
@@ -1,8 +1,8 @@
FROM mambaorg/micromamba:1.4.9
LABEL maintainer="Remi Denise (rdenise@ucc.ie)"
LABEL version="0.2.4"
LABEL version="0.2.5"

RUN micromamba install -y -n base -c conda-forge -c bioconda taxmyphage==0.2.4 && \
RUN micromamba install -y -n base -c conda-forge -c bioconda taxmyphage==0.2.5 && \
micromamba clean --all --yes
WORKDIR /app
ENTRYPOINT ["/usr/local/bin/_entrypoint.sh", "taxmyphage"]
47 changes: 24 additions & 23 deletions README.md
Original file line number Diff line number Diff line change
Expand Up @@ -267,29 +267,30 @@ The time to classify a phage will depend on the number of hits and number of pha
```
[output_folder] <- General output folder
├── Summary_taxonomy.tsv <- Summary of the analysis for all the genomes (summarises what was printed to screen)
└── [genome query_id] <- Results output for the genome query_id
├── Output_of_taxonomy.csv <- Output of the taxonomy classification
├── Summary_file.txt <- Summary of the analysis (summarises what was printed to screen)
├── heatmap.pdf <- Heatmap of the similarity to the closest relatives (pdf)
├── heatmap.png <- Heatmap of the similarity to the closest relatives (png)
├── heatmap.svg <- Heatmap of the similarity to the closest relatives (svg)
├── known_taxa.fa <- Fasta file of the closest relatives
├── mash.txt <- Output of the MASH analysis
├── query.fasta <- Input fasta file
├── similarities.tsv <- Similarities to the closest relatives
├── top_right_matrix.tsv <- Top right matrix of similarity to closest relatives (same as heatmap)
└── pmv <- VIRIDIC-like analysis
├── pmv_in.fa <- Input fasta file
├── pmv_in.fa.blastn_vs2_self.tab.gz <- Blastn output of the input fasta file against itself
├── pmv_in.fa.genus_species_clusters.tsv <- Clusters of the closest relatives
├── pmv_in.fa.ndb <- Blastn database of the closest relatives
├── pmv_in.fa.nhr
├── pmv_in.fa.nin
├── pmv_in.fa.njs
├── pmv_in.fa.not
├── pmv_in.fa.nsq
├── pmv_in.fa.ntf
└── pmv_in.fa.nto
└── Results_per_genome <- Folder containing the results for each genome
└── [genome query_id] <- Results output for the genome query_id
├── Output_of_taxonomy.csv <- Output of the taxonomy classification
├── Summary_file.txt <- Summary of the analysis (summarises what was printed to screen)
├── heatmap.pdf <- Heatmap of the similarity to the closest relatives (pdf)
├── heatmap.png <- Heatmap of the similarity to the closest relatives (png)
├── heatmap.svg <- Heatmap of the similarity to the closest relatives (svg)
├── known_taxa.fa <- Fasta file of the closest relatives
├── mash.txt <- Output of the MASH analysis
├── query.fasta <- Input fasta file
├── similarities.tsv <- Similarities to the closest relatives
├── top_right_matrix.tsv <- Top right matrix of similarity to closest relatives (same as heatmap)
└── pmv <- VIRIDIC-like analysis
├── pmv_in.fa <- Input fasta file
├── pmv_in.fa.blastn_vs2_self.tab.gz <- Blastn output of the input fasta file against itself
├── pmv_in.fa.genus_species_clusters.tsv <- Clusters of the closest relatives
├── pmv_in.fa.ndb <- Blastn database of the closest relatives
├── pmv_in.fa.nhr
├── pmv_in.fa.nin
├── pmv_in.fa.njs
├── pmv_in.fa.not
├── pmv_in.fa.nsq
├── pmv_in.fa.ntf
└── pmv_in.fa.nto
```

##### Output files explained
Expand Down
2 changes: 1 addition & 1 deletion taxmyphage/actions.py
Original file line number Diff line number Diff line change
Expand Up @@ -97,7 +97,7 @@ def all_classification(
if char in genome_id:
genome_id = genome_id.replace(char, "_")

results_path = os.path.join(args.output, genome_id)
results_path = os.path.join(args.output, "Results_per_genome", genome_id)
print_ok(f"\nClassifying {genome.id} in result folder {results_path}...")

timer_start = time.time()
Expand Down

0 comments on commit f83bd66

Please sign in to comment.