Skip to content

Commit

Permalink
RLS Version 1.5.0
Browse files Browse the repository at this point in the history
Add support for local querying of the AMPSphere

Full `ChangeLog`:

- Add support for local searching
- Slightly change output format for AMPSphere matching
- Add --verbose and --quiet flags
  • Loading branch information
luispedro committed Sep 20, 2024
1 parent 79481be commit 83fe40b
Show file tree
Hide file tree
Showing 14 changed files with 38 additions and 15 deletions.
2 changes: 1 addition & 1 deletion ChangeLog
Original file line number Diff line number Diff line change
@@ -1,4 +1,4 @@
Unreleased
Version 1.5.0 2024-09-20
* Add support for local searching
* Slightly change output format for AMPSphere matching
* Add --verbose and --quiet flags
Expand Down
21 changes: 18 additions & 3 deletions README.md
Original file line number Diff line number Diff line change
Expand Up @@ -128,17 +128,32 @@ at `test/abundances/expected.abundance.txt`.

### AMPSphere Querying

Macrel also supports querying the [AMPSphere database](https://ampsphere.big-data-biology.org/) (described in [Santos-Júnior et al., 2024](https://doi.org/10.1016/j.cell.2024.05.013)). To do so, use the `query-ampsphere` subcommand:
Macrel also supports querying the [AMPSphere
database](https://ampsphere.big-data-biology.org/) (described in [Santos-Júnior
et al., 2024](https://doi.org/10.1016/j.cell.2024.05.013)). To do so, use the
`query-ampsphere` subcommand:

```bash
macrel query-ampsphere \
--fasta example_seqs/pep8.faa \
--output out_ampsphere
```

Note that this command requires internet access as it uses the AMPSphere API.
Note that, by default, this command requires internet access as it uses the
AMPSphere API. Alternatively, you can use the `--local` flag to download a copy
of the database and run the queries locally. This only requires the network the
first time.

By default it performs exact matching, but you can also use MMSeqs2 to perform approximate matching by using the `--query-mode=mmseqs` (or `--query-mode=hmm` for HMMER).
```bash
macrel query-ampsphere \
--local \
--fasta example_seqs/pep8.faa \
--output out_ampsphere
```

By default it performs exact matching, but you can also use MMSeqs2 to perform
approximate matching by using the `--query-mode=mmseqs` (or `--query-mode=hmm`
for HMMER).

### Community

Expand Down
8 changes: 8 additions & 0 deletions docs/whatsnew.md
Original file line number Diff line number Diff line change
@@ -1,5 +1,13 @@
# What's new? (History)

## Version 1.5.0

*Released 20 September 2024*

- Add support for local searching when querying the [AMPSphere database](https://ampsphere.big-data-biology.org/)
- Slightly change output format for AMPSphere matching
- Add --verbose and --quiet flags

## Version 1.4.0

*Released 27 June 2024*
Expand Down
2 changes: 1 addition & 1 deletion macrel/macrel_version.py
Original file line number Diff line number Diff line change
@@ -1 +1 @@
__version__ = '1.5.0.dev0'
__version__ = '1.5.0'
2 changes: 1 addition & 1 deletion tests/contigs.cluster/expected.prediction
Original file line number Diff line number Diff line change
@@ -1,4 +1,4 @@
# Prediction from macrel v1.5.0.dev0
# Prediction from macrel v1.5.0
Access Sequence AMP_family AMP_probability Hemolytic Hemolytic_probability
smORF_2 RFLIKMVKVNLMNGKLIRKISLM CLP 0.634 Hemo 0.871
smORF_19 FFNDGKGTIYYGIKKYFRIYF CLP 0.673 Hemo 0.822
2 changes: 1 addition & 1 deletion tests/contigs.nosmorfs/expected.percontigs
Original file line number Diff line number Diff line change
@@ -1,4 +1,4 @@
# Prediction from macrel v1.5.0.dev0
# Prediction from macrel v1.5.0
# Macrel calculated for the sample a density of 0.000 AMPs / Mbp.
contig length ORFs smORFs AMPs
scaffold2530_2_MH0058 1324 1 0 0
2 changes: 1 addition & 1 deletion tests/contigs.nosmorfs/expected.prediction
Original file line number Diff line number Diff line change
@@ -1,2 +1,2 @@
# Prediction from macrel v1.5.0.dev0
# Prediction from macrel v1.5.0
Access Sequence AMP_family AMP_probability Hemolytic Hemolytic_probability
2 changes: 1 addition & 1 deletion tests/contigs/expected.percontigs
Original file line number Diff line number Diff line change
@@ -1,4 +1,4 @@
# Prediction from macrel v1.5.0.dev0
# Prediction from macrel v1.5.0
# Macrel calculated for the sample a density of 45.062 AMPs / Mbp.
contig length ORFs smORFs AMPs
C4060843_1_MH0058 518 1 1 0
Expand Down
2 changes: 1 addition & 1 deletion tests/contigs/expected.prediction
Original file line number Diff line number Diff line change
@@ -1,4 +1,4 @@
# Prediction from macrel v1.5.0.dev0
# Prediction from macrel v1.5.0
Access Sequence AMP_family AMP_probability Hemolytic Hemolytic_probability
scaffold75334_1_MH0058_1 RFLIKMVKVNLMNGKLIRKISLM CLP 0.634 Hemo 0.871
scaffold33693_17_MH0058_2 FFNDGKGTIYYGIKKYFRIYF CLP 0.673 Hemo 0.822
2 changes: 1 addition & 1 deletion tests/peptides/expected.prediction
Original file line number Diff line number Diff line change
@@ -1,4 +1,4 @@
# Prediction from macrel v1.5.0.dev0
# Prediction from macrel v1.5.0
Access Sequence AMP_family AMP_probability Hemolytic Hemolytic_probability
AP00002|AMP YVPLPNVPQPGRRPFPTFPGQGPFNPKIKWPQGY CLP 0.861 Hemo 0.663
AP00007|AMP GNNRPVYIPQPRPPHPRL CLP 0.970 Hemo 0.515
Expand Down
2 changes: 1 addition & 1 deletion tests/reads.se/expected.percontigs
Original file line number Diff line number Diff line change
@@ -1,4 +1,4 @@
# Prediction from macrel v1.5.0.dev0
# Prediction from macrel v1.5.0
# Macrel calculated for the sample a density of 59.743 AMPs / Mbp.
contig length ORFs smORFs AMPs
k47_0 3379 4 2 0
Expand Down
2 changes: 1 addition & 1 deletion tests/reads.se/expected.prediction
Original file line number Diff line number Diff line change
@@ -1,4 +1,4 @@
# Prediction from macrel v1.5.0.dev0
# Prediction from macrel v1.5.0
Access Sequence AMP_family AMP_probability Hemolytic Hemolytic_probability
k47_10_1 RFLIKMVKVNLMNGKLIRKISLM CLP 0.634 Hemo 0.871
k47_11_1 FFNDGKGTIYYGIKKYFRIYF CLP 0.673 Hemo 0.822
2 changes: 1 addition & 1 deletion tests/reads/expected.percontigs
Original file line number Diff line number Diff line change
@@ -1,4 +1,4 @@
# Prediction from macrel v1.5.0.dev0
# Prediction from macrel v1.5.0
# Macrel calculated for the sample a density of 57.627 AMPs / Mbp.
contig length ORFs smORFs AMPs
k77_11 1303 2 1 0
Expand Down
2 changes: 1 addition & 1 deletion tests/reads/expected.prediction
Original file line number Diff line number Diff line change
@@ -1,4 +1,4 @@
# Prediction from macrel v1.5.0.dev0
# Prediction from macrel v1.5.0
Access Sequence AMP_family AMP_probability Hemolytic Hemolytic_probability
k77_12_1 RFLIKMVKVNLMNGKLIRKISLM CLP 0.634 Hemo 0.871
k77_15_1 FFNDGKGTIYYGIKKYFRIYF CLP 0.673 Hemo 0.822

0 comments on commit 83fe40b

Please sign in to comment.