Skip to content
Merged
Show file tree
Hide file tree
Changes from all commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
Original file line number Diff line number Diff line change
Expand Up @@ -3,9 +3,16 @@
Minimal crafted ORF catalogue used as input to the `custom/orfcollapse` module
test. The real chr20 data does not contain the case orfcollapse handles (the
same micropeptide encoded at distinct, non-overlapping loci), so this small
fixture is hand-built to exercise it: two identical-peptide smORFs on opposite
strands/loci, one canonical CDS and one unique smORF.
fixture is hand-built to exercise it: two identical-peptide small ORFs on
opposite strands/loci, one canonical CDS and one unique small ORF.

`orf_class` is positional only — it records where an ORF sits relative to the
annotated CDS and never encodes length. The `is_smorf` column carries the
length flag (`aa_length` <= `--smorf-max-aa`), and `orf_type_native` preserves
each caller's own ORF-type label. The two identical-peptide rows share a class
deliberately: orfcollapse prefers the more specific class when a peptide cluster
spans several, so differing classes would change which row survives the fold.

- `cohort.catalogue.{bed12,tsv}`, `cohort.orf_to_gene.tsv`, `cohort.catalogue.mqc.tsv`: orfmerge-shaped catalogue.
- `cohort.catalogue.aa.fasta`: catalogue peptides (headers carry the `(+)/(-)` suffix bedtools getfasta -nameOnly -s emits).
- `cohort_cluster.tsv`: MMseqs2 cluster TSV grouping the two identical smORFs.
- `cohort_cluster.tsv`: MMseqs2 cluster TSV grouping the two identical small ORFs.
Original file line number Diff line number Diff line change
@@ -1,3 +1,9 @@
Class Count
smORF 3
canonical_cds 1
uORF 2
uoORF 0
dORF 1
doORF 0
intORF 0
novel_u 0
other 0
Original file line number Diff line number Diff line change
@@ -1,5 +1,5 @@
orf_id chrom start end strand gene_id transcript_id orf_class aa_length called_by_ribotish called_by_ribocode called_by_ribotricer called_by_rpbp called_by_price score_ribotish score_ribocode score_ribotricer score_rpbp score_price n_samples samples
orf_00000001 chr20 100 160 + g1 t1 smORF 20 1 0 0 0 0 0.01 1 sampleA
orf_00000002 chr20 9000 9060 - g5 t5 smORF 20 0 1 0 0 0 0.02 1 sampleB
orf_00000003 chr20 2000 2900 + g2 t2 canonical_cds 300 1 0 0 0 0 0.5 1 sampleA
orf_00000004 chr20 50 125 + g3 t3 smORF 25 1 0 0 0 0 0.3 1 sampleA
orf_id chrom start end strand gene_id transcript_id orf_class aa_length called_by_ribotish called_by_ribocode called_by_ribotricer called_by_rpbp called_by_price score_ribotish score_ribocode score_ribotricer score_rpbp score_price n_samples samples orf_type_native is_smorf
orf_00000001 chr20 100 160 + g1 t1 uORF 20 1 0 0 0 0 0.01 1 sampleA 5'UTR 1
orf_00000002 chr20 9000 9060 - g5 t5 uORF 20 0 1 0 0 0 0.02 1 sampleB uORF 1
orf_00000003 chr20 2000 2900 + g2 t2 canonical_cds 300 1 0 0 0 0 0.5 1 sampleA Annotated 0
orf_00000004 chr20 50 125 + g3 t3 dORF 25 1 0 0 0 0 0.3 1 sampleA 3'UTR 1