Picard-tools

From wiki
Revision as of 18:22, 15 February 2017 by Rf (talk | contribs) (Created page with "= Introduction = Started as a companion package to GATK, so that GATK could concentrate on its core functions, but has blossomed beyond that to be quite substnatial in its ow...")
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to: navigation, search

Introduction

Started as a companion package to GATK, so that GATK could concentrate on its core functions, but has blossomed beyond that to be quite substnatial in its own right, and is sometimes better than samtools.

Its toolbox takes the form of subcommand which are fed to the chief picard-tools program. As it is java, this will be a jar file.

Note that Java 1.8 is a prerequisite for picard-tools.

Usage

picard-tools follows conventions set by GATK. This means that it is a bit awkward clumsy for one-liner command-lines and sits better in a script.

Typical script

Here we use the MergeSamFiles tool.

#!/bin/bash 
A=( $* )
NA=${#A[@]}
INPS=""
for i in ${A[@]}; do
    INPS="$INPS -I $i"
done
INPS="$INPS -O allmerged.bam"
java -jar ~/swmake/jars/picard.jar MergeSamFiles SORT_ORDER=coordinate${INPS}