Tool:picardmetrics: Run Picard tools and collate the outputs into one file.
0
7
Entering edit mode
7.1 years ago
Kamil ★ 2.1k

picardmetrics

Picardmetrics is a Bash script that simplifies calling Picard tools and collates the different output files generated by Picard. It also has functions for generating the two input files required by CollectRnaSeqMetrics.

In order, picardmetrics run will do the following:

  1. Automatically create a sequence dictionary using your reference sequence.
  2. Create a new temporary BAM file that you can keep with option -k.
  3. Reorder the header of the BAM file to match the reference.
  4. Sort the reads in the BAM file by coordinate.
  5. Mark duplicates in the BAM file and report duplicate metrics.
  6. Run up to 8 additional Picard tools.

After running the tools, use picardmetrics collate to merge all of the generated metrics from multiple BAM files into tab-delimited files. Additionally, all of these tab-delimited files are consolidated into a single file with all metrics from all BAM files and all Picard tools.

Download

Download the latest release here: https://github.com/slowkow/picardmetrics/releases/latest

Example

Use Picard to assess the quality of your sequencing data. This example shows RNA-seq data from hundreds of glioblastoma cells and gliomasphere cell lines.

On the left, each sample is a point, and we see that samples with high mean mapping quality have the greatest number of detected genes. Further, the color of points reveals variation in the percent of reads per sample that are assigned to exons.

On the right, each sample is a bar, broken down into the percent of sequenced bases coming from different genomic regions. We see that many samples have few sequenced bases coming from coding regions relative to intergenic regions.

RNA-Seq Picard Bash R Tool • 3.3k views
ADD COMMENT
0
Entering edit mode

Thanks Kamil. This will come in handy for newbies. 

ADD REPLY
0
Entering edit mode

Looks excellent! Will come in handy for experience users as well. 

ADD REPLY

Login before adding your answer.

Traffic: 2158 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6