This is a test version of Biostars. For the public version, visit https://www.biostars.org.
run picard on all bam files in directory

I am trying to run a bash loop on all bam files in a directory and having a bit of trouble. The bash so far as well as the individual picard command is below. Thank you :).

Bash loop so far

for f in /home/cmccabe/Desktop/NGS/pool_I_090215/*.bam ; do
  bname=`basename $f`
  pref=${bname%%.bam}
  java \
    -jar /home/cmccabe/Desktop/NGS/picard-tools-1.139/picard.jar \
    CalculateHsMetrics \
    BI=/home/cmccabe/Desktop/NGS/bed/sam_sorted_unix_5column_xgen_probes.bed \
    TI=/home/cmccabe/Desktop/NGS/bed/sam_sorted_unix_5column_xgen_targets.bed \
    I=/$f \
    O=/home/cmccabe/Desktop/NGS/pool_I_090215/${pref}_all_IDT.CalculateHSmetrics
done

If the bam files were run individually here is the command

java \
  -jar /home/cmccabe/Desktop/NGS/picard-tools-1.139/picard.jar \
  CalculateHsMetrics \
  BI=/home/cmccabe/Desktop/NGS/bed/sam_sorted_unix_5column_xgen_probes.bed \
  TI=/home/cmccabe/Desktop/NGS/bed/sam_sorted_unix_5column_xgen_targets.bed \
  I=/home/cmccabe/Desktop/NGS/pool_II_090115/IonXpress_015_150901_newheader.bam \
  O=/home/cmccabe/Desktop/NGS/pool_II_090115/IonXpress_015_150901_all_IDT.CalculateHSmetrics
bam ngs picard
• 3,827 views • link
updated by Ram Ram • written by bioguy24 bioguy24

What is the problem? What is the question?

• 0 views • link
updated by Ram Ram • written by Pierre Lindenbaum Pierre Lindenbaum

What error message are you getting?

• 0 views • link
written by Alex Reynolds Alex Reynolds

Are you sharing this as a useful tip, bioguy24?

• 0 views • link
written by Ram Ram

1 answer

Some of the problems I see:

  • You're specifying parameters in an abbreviated fashion. TI should be TARGET_INTERVALS, for example.
  • The tool doesn't use BED format for the intervals files. It uses its own interval list format.
  • The output for each sample is being written to the same output file. That means you'll be overwriting it each time.

You should at minimum clearly state a question/problem and provide any error output that you obtain from the tool.

• 0 views • link
updated by Ram Ram • written by Dan D Dan D

Log in to answer this question.

More posts like this
  • GATK MarkDuplicates output and bam
    written by bioguy24 bioguy24
    The below command executes and produces the dup_metrics below, I am trying to interpret them and don't really understand the metrics and in the bam …
  • using bash variable in picard MarkDuplicates command
    written by bioguy24 bioguy24
    In the below `bash` the two variables `$bname` and `$sample` are extracted correctly however when I pass them to the `java` command `$bname` is set …
  • Bash to match and store numerical prefix to update files
    written by bioguy24 bioguy24
    In the bash below the unique headers of each vcf.gz are stored in a text file with the same name. There can be multiple vcf.gz …
  • add gene name from bed file to vcf in vcftools
    written by bioguy24 bioguy24
    I am trying to add gene information to a vcf file using vcftools: vcftools annotate \ -a genes.bed.gz \ -c CHROM,FROM,TO,GENE \ -h <(echo '##INFO=<ID=GENE,Number=1,Type=String,Description="Gene …
  • Regarding Depth coverage
    written by manojkumarbioinfo manojkumarbioinfo
    Hi. all I'm new to field of NGS. and I am trying to calculate the Dept Coverage for my alignment. For this I use picard …
  • bash loop to count variants using vcftools
    written by bioguy24 bioguy24
    I can run the below vcftools command to count individual variant types in vcf.gz files zcat /home/cmccabe/Desktop/vcf/file1.vcf.gz | vcf-annotate --fill-type | grep -oP "TYPE=\w+" | …
  • picard error version 1.139 on ubuntu 14.04
    written by bioguy24 bioguy24
    <p>I seem to have trouble with calling picard from the command line but other tools work perfectly.<br /> <br /> path:</p> <pre> /home/cmccabe/Desktop/NGS/bedops-2.14.4/bin:/home/cmccabe/Desktop/NGS/samtools-1.2:/usr/local/sbin:/usr/local/bin:/usr/sbin:/usr/bin:/sbin:/bin:/usr/games:/usr/local/games:/home/cmccabe/Desktop/NGS/bedtools2-25.0/bin:/home/cmccabe/Desktop/NGS/bbmap-35.59:/home/cmccabe/Desktop/NGS/FastQC-0.11.4:/home/cmccabe/Desktop/NGS/picard-tools-1.139</pre> <p>The code …
  • coverageBed with -g option ?
    written by bioguy24 bioguy24
    I am running a coverageBed and getting the below error. ``` cmccabe@DTV-A5211QLM:~/Desktop/NGS$ coverageBed -d -sorted -a /home/cmccabe/Desktop/NGS/bed/bedtools/xgen_targets_sorted.bed -b /home/cmccabe/Desktop/NGS/pool_I_090215/IonXpress_008_150902_newheader.bam > /home/cmccabe/Desktop/NGS/pool_I_090215/IonXpress_008_150902_output.txt ERROR: Database file /home/cmccabe/Desktop/NGS/pool_I_090215/IonXpress_008_150902_newheader.bam …
  • coverage bed error but produces output file
    written by bioguy24 bioguy24
    The below command produces a file (674 MB). ``` coverageBed \ -d \ -sorted \ -a /home/cmccabe/Desktop/NGS/bed/sorted_unix_5column_xgen_targets.bed \ -b /home/cmccabe/Desktop/NGS/pool_I_090215/IonXpress_008_150902_newheader.bam \ > /home/cmccabe/Desktop/NGS/pool_I_090215/IonXpress_008_150902_output.txt ``` Output …
  • bedtools 2.25 vs bedtools2.20 coverageBed output not consistent
    written by bioguy24 bioguy24
    Hopefully, this edit is more clear in explaining the issue and will help led to a solution. As I am lost as to why there …