Read Duplication

determines reads duplication rate with sequence-based and mapping-based strategies (Galaxy Version 5.0.3+galaxy0)

Tool Parameters

Please provide a value for this option.
* required
(--input-file)
*
(--up-limit)
*
Minimum mapping quality for an alignment to be considered as "uniquely mapped" (--mapq)
Output the R-Script used to generate the plots

Additional Options

Send an email notification when the job completes.

Help

read_duplication.py

Two strategies were used to determine reads duplication rate:

  • Sequence based: reads with exactly the same sequence content are regarded as duplicated reads.
  • Mapping based: reads mapped to the same genomic location are regarded as duplicated reads. For splice reads, reads mapped to the same starting position and splice the same way are regarded as duplicated reads.

Inputs

Input BAM/SAM file
Alignment file in BAM/SAM format.
Upper Limit of Plotted Duplicated Times (default=500)
Only used for plotting.

Output

  1. output.dup.pos.DupRate.xls: Read duplication rate determined from mapping position of read. First column is "occurrence" or duplication times, second column is number of uniquely mapped reads.
  2. output.dup.seq.DupRate.xls: Read duplication rate determined from sequence of read. First column is "occurrence" or duplication times, second column is number of uniquely mapped reads.
  3. output.DupRate_plot.r: R script to generate pdf file
  4. output.DupRate_plot.pdf: graphical output generated from R script
shed_tool_static/toolshed.g2.bx.psu.edu/nilesh/rseqc/rseqc_read_duplication/5.0.3+galaxy0/duplicate.png

About RSeQC

The RSeQC package provides a number of useful modules that can comprehensively evaluate high throughput sequence data especially RNA-seq data. "Basic modules" quickly inspect sequence quality, nucleotide composition bias, PCR bias and GC bias, while "RNA-seq specific modules" investigate sequencing saturation status of both splicing junction detection and expression estimation, mapped reads clipping profile, mapped reads distribution, coverage uniformity over gene body, reproducibility, strand specificity and splice junction annotation.

The RSeQC package is licensed under the GNU GPL v3 license.

shed_tool_static/toolshed.g2.bx.psu.edu/nilesh/rseqc/rseqc_read_duplication/5.0.3+galaxy0/logo.png

Unnamed history

Draggable