RSeQC is a comprehensive quality control tool for RNA-seq data analysis. It provides a set of modules for assessing various aspects of RNA-seq experiments.
mamba install -c bioconda rseqc
Before running RSeQC modules, ensure you have:
- A sorted and indexed BAM file
- A reference gene model in BED format
-
Basic Alignment Statistics:
bam_stat.py -i sample_sorted.bam > bam_stat_output.txt
-
Read Distribution:
read_distribution.py -i sample_sorted.bam -r reference.bed > read_distribution_output.txt
-
Gene Body Coverage:
geneBody_coverage.py -r reference.bed -i sample_sorted.bam -o gene_body_coverage_output
-
Infer Experiment:
infer_experiment.py -i sample_sorted.bam -r reference.bed > infer_experiment_output.txt
-
Inner Distance:
inner_distance.py -i sample_sorted.bam -o inner_distance_output -r reference.bed
-
Junction Saturation:
junction_saturation.py -i sample_sorted.bam -r reference.bed -o junction_saturation_output
-
Read Duplication:
read_duplication.py -i sample_sorted.bam -o read_duplication_output
-
GC Content:
read_GC.py -i sample_sorted.bam -o read_GC_output
-
Read Quality:
read_quality.py -i sample_sorted.bam -o read_quality_output
-
Nucleotide Composition:
read_NVC.py -i sample_sorted.bam -o read_NVC_output