Nextflow Modules
Showing module(s) with keyword "smorf"
| Module | Keywords | Description |
|---|---|---|
| nf-core/custom/orfcollapse | orf ribo-seq catalogue smorf deduplication | Collapse small ORFs that share an amino-acid sequence cluster into a single catalogue entry. Pair with `custom/orfmerge` (coordinate-based catalogue), `bedtools/getfasta` + `seqkit/translate` (AA FASTA keyed by orf_id), and `mmseqs/easycluster` (AA clusters) upstream. The coordinate-based merge in `custom/orfmerge` only groups ORFs that overlap on the genome, so the same micropeptide encoded at several distinct, non-overlapping loci (typically repetitive regions) survives as separate rows. This adopts the peptide-level deduplication and 0.9 amino-acid-similarity threshold of the GENCODE Ribo-seq ORF consolidation (Mudge et al. 2022, Nat Biotechnol, doi:10.1038/s41587-022-01369-0; gencode-riboseqORFs collapse_cutoff 0.9), implemented here with MMseqs2 sequence-identity clustering rather than that tool's longest-shared-string / P-site-overlap metric. Small ORFs (`aa_length` <= `--smorf-max-aa`, default 100) are clustered by amino-acid identity upstream and this module folds each multi-member cluster down to one representative. Only small ORFs are collapsed; larger ORFs are passed through untouched. Eligibility is the catalogue's `is_smorf` flag, independent of `orf_class`, so a short uORF and a short novel ORF are both candidates; `--smorf-max-aa` re-derives the flag and aborts on disagreement. Among the members of a cluster the representative is chosen by class specificity, then longest aa_length, then orf_id, so the result does not depend on which sequence MMseqs2 labelled the cluster representative. Catalogue row order is preserved; dropped members fold their `called_by_<caller>` / `score_<caller>` evidence, `n_samples` / `samples` recurrence and gene mappings into the survivor. |