4.8 smORF analysis
Function
The smORF analysis module provides a complete workflow for scanning candidate small open reading frames (smORFs), filtering candidates with sequence and Kozak-context criteria, evaluating Ribo-seq translation evidence, and integrating sample-level evidence into ORF-level summary tables.
The workflow is designed for transcript-centric smORF discovery. Candidate ORFs are first generated from a genome FASTA and genePred annotation, then filtered and evaluated using P-site density tracks from Ribo-seq data.
Workflow
smorf_scanner → smorf_filter → smorf_evidence → smorf_integrate
| Step |
Command |
Main purpose |
| 1 |
smorf_scanner |
Scan transcript-centric candidate ORFs from genome FASTA and genePred annotation. |
| 2 |
smorf_filter |
Filter scanned ORFs by length, start codon, category, strand, completeness, and Kozak context. |
| 3 |
smorf_evidence |
Evaluate Ribo-seq P-site density, frame periodicity, coverage, pausing, and release evidence. |
| 4 |
smorf_integrate |
Integrate long-format sample-level evidence into ORF-level matrices and summary tables. |
| Input |
Used by |
Description |
| Genome FASTA |
smorf_scanner |
Genome sequence used to reconstruct transcript sequences. |
| genePred annotation |
smorf_scanner, optional for smorf_evidence |
Transcript annotation. Scanner uses it to define transcript structures; evidence uses it to recover ORF exon blocks when needed. |
| Scanner message table |
smorf_filter |
Candidate ORF table generated by smorf_scanner. |
| Filtered ORF table |
smorf_evidence |
Passed ORF table generated by smorf_filter. |
| P-site density files |
smorf_evidence |
Strand-specific or unstranded bedGraph/WIG files containing Ribo-seq P-site density. |
| Evidence table |
smorf_integrate |
Long-format ORF-by-sample table generated by smorf_evidence. |
Main evidence levels
| Evidence |
Meaning |
| ORF sequence |
Start codon, stop codon, ORF length, strand, and category. |
| Kozak context |
Start-codon context scored by annotated, built-in, PWM, or sequence-derived Kozak models. |
| Ribo-seq signal |
Total RPF signal, covered nucleotides/codons, and coverage ratio. |
| Periodicity |
Frame-specific signal distribution, especially frame-0 enrichment. |
| Start/stop behavior |
Start pausing, pre-stop pausing, and post-stop release signatures. |
| Coverage shape |
Uniform, disperse, or skewed RPF distribution across the ORF. |
| Multi-sample support |
Reproducibility and support level across multiple Ribo-seq samples. |
| Integration |
ORF-level summary, best sample, evidence labels, and frame-density matrix. |
Notes
smorf_scanner produces many candidate ORFs. Scanner output alone should not be treated as evidence of translation.
smorf_filter removes candidates according to sequence features and Kozak context, but it does not use Ribo-seq evidence.
smorf_evidence requires either a density list, one unstranded density file, or both plus- and minus-strand density files.
smorf_integrate is most useful when multiple Ribo-seq samples or replicates are available.
- Strand-specific P-site density files can be generated with
rpf_Bam2bw and then supplied to smorf_evidence as bedGraph or WIG tracks.