Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • Marius
    8armed
    • Dec 2010
    • 30

    #1

    Problem using BASH to pipe Piccard command

    I'm new to Bash, but would like to use it.
    I would like to pipe two Picard (GATK pipeline) commands using Bash. Each of these commands works fine when they are run as individual Bash scripts (see below). However, when I try to combine them within a single bash script, it no longer works. I really cannot find the bug though. Help is much appreciated.


    BASH Script 1
    #!/bin/bash

    # Set paths
    picardtools='PATH/picard.jar'
    bams="PATH/BAM"
    outSorted='PATH/outSorted'

    for i in BamFile1 BamFile2
    do
    java -jar $picardtools SortSam INPUT=$bams/"$i" OUTPUT=$outSorted/"$i" SORT_ORDER=coordinate CREATE_INDEX=TRUE VALIDATION_STRINGENCY=LENIENT
    done


    BASH Script 2
    #!/bin/bash

    # Set paths
    picardtools='PATH/picard.jar'
    outSorted='PATH/outSorted'
    outSorted2='PATH/outSorted2'

    for i in BamFile1 BamFile2
    do
    java -jar $picardtools AddOrReplaceReadGroups INPUT=$outSorted/$i OUTPUT=$outSorted2/"$i" SORT_ORDER=coordinate RGID="$i" RGLB="$i" RGPL=illumina RGPU="$i" RGSM="$i" CREATE_INDEX=True
    done


    ######################

    Here my attempts to pipe the two Scripts into one. But neither of these attempts worked. Can you help me find out why?

    Attempt 1
    #!/bin/bash

    # Set paths
    picardtools='PATH/picard.jar'
    bams="PATH/BAM"
    outSorted='PATH/outSorted'
    outSorted2='PATH/outSorted2'

    for i in BamFile1 BamFile2
    do
    java -jar $picardtools SortSam INPUT=$bams/"$i" OUTPUT=$outSorted/"$i" SORT_ORDER=coordinate CREATE_INDEX=TRUE VALIDATION_STRINGENCY=LENIENT
    java -jar $picardtools AddOrReplaceReadGroups INPUT=$outSorted/$i OUTPUT=$outSorted2/"$i" SORT_ORDER=coordinate RGID="$i" RGLB="$i" RGPL=illumina RGPU="$i" RGSM="$i" CREATE_INDEX=True
    done


    Attempt 2
    #!/bin/bash

    # Set paths
    picardtools='PATH/picard.jar'
    bams="PATH/BAM"
    outSorted2='PATH/outSorted2'

    for i in BamFile1 BamFile2
    do
    java -jar $picardtools SortSam INPUT=$bams/"$i" OUTPUT=/dev/stdout SORT_ORDER=coordinate CREATE_INDEX=TRUE VALIDATION_STRINGENCY=LENIENT | \
    java -jar $picardtools AddOrReplaceReadGroups INPUT=/dev/stdin OUTPUT=$outSorted2/"$i" SORT_ORDER=coordinate RGID="$i" RGLB="$i" RGPL=illumina RGPU="$i" RGSM="$i" CREATE_INDEX=True
    done

Latest Articles

Collapse

  • SEQadmin2
    Beyond CRISPR/Cas9: Understand, Choose, and Use the Right Genome Editing Tool
    by SEQadmin2



    CRISPR/Cas9 sparked the gene editing revolution for both research and therapeutics.1 But this system still showed severe issues that limited its applications. The most prominent were the heavy reliance on PAM sequences, delivery limitations, double-stranded breaks that prompt unintended edits and cell death, and editing inefficiency (both in targeting and in knock-in reliability).

    Despite this, “CRISPR helped turn genome editing from a specialized technique into
    ...
    07-31-2026, 11:01 AM
  • SEQadmin2
    Proteomic Platforms: How to Choose the Right Analytical Strategy to Improve Detection and Clinical Applications
    by SEQadmin2


    Proteomics platforms are evolving rapidly, with advances in mass spectrometry and affinity-based approaches expanding what researchers can detect and at what scale. As the field moves toward deeper proteome coverage and clinical applications, scientists face an increasingly complex landscape of tools. This article will explore how researchers are navigating these choices to find the right platform for their work.

    The systematic characterization of the human proteome has
    ...
    07-20-2026, 11:48 AM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by SEQadmin2, Yesterday, 12:22 PM
0 responses
14 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 08-11-2026, 10:35 AM
0 responses
16 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 08-06-2026, 07:41 AM
0 responses
31 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 08-03-2026, 10:13 AM
0 responses
50 views
0 reactions
Last Post SEQadmin2  
Working...