Scavenger
Scavenger recovers unaligned RNA-seq reads by identifying potential alignment locations through sequence similarity to already aligned reads, improving the sensitivity of downstream RNA-seq analyses.
Key Features:
- Recovery of Unaligned Reads: Identifies potential alignment locations for previously unmapped reads by comparing them to aligned reads and proposing new alignments based on sequence similarity.
- Python Implementation: Implemented in Python.
- Applicability to RNA-seq Types: Evaluated on simulated and real datasets, including single-cell RNA-seq, demonstrating applicability across diverse RNA-seq data types.
- Detection of Genetic Variants: Recovered reads often exhibit more genetic variants relative to the reference genome, implicating personal genomic variation in alignment failures.
- Impact on Downstream Analyses: Integration of recovered reads alters gene expression estimation and differential expression results, notably affecting lowly expressed genes such as pseudogenes.
Scientific Applications:
- Improving Read Alignment Sensitivity: Addresses false-negative non-alignment to increase the comprehensiveness of RNA-seq alignments.
- Gene Expression Profiling and Differential Expression: Enhances estimation of gene expression levels and differential expression analysis by incorporating recovered reads.
- Variant Detection and Personalized Genomics: Reveals additional genetic variants relative to the reference genome that may be missed by initial alignments.
Methodology:
Identifies potential alignment locations for unaligned reads by detecting sequence similarity with aligned reads and integrates the recovered reads into downstream analyses.
Topics
Details
- License:
- MIT
- Programming Languages:
- Shell, Python
- Added:
- 11/14/2019
- Last Updated:
- 12/17/2020
Operations
Publications
Yang A, Tang JYS, Troup M, Ho JWK. Scavenger: A pipeline for recovery of unaligned reads utilising similarity with aligned reads. F1000Research. 2019;8:1587. doi:10.12688/f1000research.19426.1.
Funding: - National Health and Medical Research Council: 1105271
- National Heart Foundation of Australia: 100848