Picard AddOrReplaceReadGroups unable to read my argument?
1
0
Entering edit mode
9.4 years ago
quachtina96 ▴ 40

Hi, I am trying to run Picard's AddOrReplaceReadGroups on a bam file called exome_mtExtractremap.csort.bam (code below). the error messages states ERROR: Option 'RGPL' is required., however I have included RGPL, so I don't understand why it is occuering.

I've tried RGLB=ILLUMINA, RGLB="ILLUMINA", and RGLB=illumina. What am I doing wrong?

i=exome_mtExtractremap.csort.bam
java -Xmx2g \
  -Djava.io.tmpdir=`pwd`/tmp \
  -jar /opt/applications/picard/current/AddOrReplaceReadGroups.jar \
  INPUT="${i}" \
  OUTPUT="${i}.rg.bam" \
  SORT_ORDER=coordinate \
  RGLB="$(echo $OUTPUT | sed 's/_mtExtractremap.csort.bam.rg.bam//')" \
  RGPL="illumina" \
  RGPU="L002" \
  RGSM="${RGLB}" \
  RGID="${RGLB}_${RGPU}";
Picard Read-Group • 4.4k views
ADD COMMENT
2
Entering edit mode
9.4 years ago

I suspect there something wrong in the line just before RGPL

ADD COMMENT
0
Entering edit mode

This solved a problem for me. GATK was reporting "Invalid Argument" for my entry after --RGPU <platform.unit>. It turned out I had removed a space in the line before:

bad:

--RGPL= sample_num\
--RGPU= platform.unit \

good:

--RGPL= sample_num \
--RGPU= platform.unit \
ADD REPLY

Login before adding your answer.

Traffic: 2015 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