---
title: "RNA-Seq BAM File Quality Control with RSeQC (Drosophila melanogaster)"
canonical: "https://help.biobam.com/space/OBD/3055353857/RNA-Seq%20BAM%20File%20Quality%20Control%20with%20RSeQC%20(Drosophila%20melanogaster)"
format: markdown
---
# Introduction

Accession: PRJNA657899

This example dataset contains RNA-Seq reads that have been mapped (STAR) to their genome (Drosophila melanogaster). Some obtained BAM files together with the latest structural annotation (GFF) are provided here to show an exemplary usage of “BAM Quality Control” in OmicsBox.

### Dataset description

Information about the experiment.

- RNA-Seq
- 14 samples
- HiSeq2500
- Single and Paired-End

### Publication

[https://www.ncbi.nlm.nih.gov/pmc/articles/PMC7936876](https://www.ncbi.nlm.nih.gov/pmc/articles/PMC7936876/)

### Original Data

Link to the public repositories where the original data can be found:

- [https://www.ncbi.nlm.nih.gov/bioproject/PRJNA657899](https://www.ncbi.nlm.nih.gov/bioproject/PRJNA657899)
  - Convenient website to download the files from NCBI [https://sra-explorer.info/#](https://sra-explorer.info/#)
- Genome: [https://ftp.ncbi.nlm.nih.gov/genomes/all/GCF/000/001/215/GCF_000001215.4_Release_6_plus_ISO1_MT/GCF_000001215.4_Release_6_plus_ISO1_MT_genomic.fna.gz](https://ftp.ncbi.nlm.nih.gov/genomes/all/GCF/000/001/215/GCF_000001215.4_Release_6_plus_ISO1_MT/GCF_000001215.4_Release_6_plus_ISO1_MT_genomic.fna.gz)
- Annotation: [https://ftp.ncbi.nlm.nih.gov/genomes/all/GCF/000/001/215/GCF_000001215.4_Release_6_plus_ISO1_MT/GCF_000001215.4_Release_6_plus_ISO1_MT_genomic.gff.gz](https://ftp.ncbi.nlm.nih.gov/genomes/all/GCF/000/001/215/GCF_000001215.4_Release_6_plus_ISO1_MT/GCF_000001215.4_Release_6_plus_ISO1_MT_genomic.gff.gz)

# Bioinformatic Analysis

## 1- Analysis Step

### Application

[https://biobam.atlassian.net/wiki/spaces/OBD/pages/3044212737](https://biobam.atlassian.net/wiki/spaces/OBD/pages/3044212737) 

### Input

- Input BAM files
  - > Macro (lref-gdrive-file)
  - 2 samples have been selected and further reduced in size (samtools view -bo subset.bam -s 123.1 alignments.bam chr1 chr2).
- Annotation file (Gene Model)
  - [https://ftp.ncbi.nlm.nih.gov/genomes/all/GCF/000/001/215/GCF_000001215.4_Release_6_plus_ISO1_MT/GCF_000001215.4_Release_6_plus_ISO1_MT_genomic.gff.gz](https://ftp.ncbi.nlm.nih.gov/genomes/all/GCF/000/001/215/GCF_000001215.4_Release_6_plus_ISO1_MT/GCF_000001215.4_Release_6_plus_ISO1_MT_genomic.gff.gz)
  - GFF has to be extracted before providing it.

### Parameters

- Default values
- Provide the Annotations mentioned above
  - Extract the GFF before providing it

### Execution Time

- 2:20 for the 2 subsetted samples and annotations provided

### Output