SEQanswers

SEQanswers (http://seqanswers.com/forums/index.php)
-   Bioinformatics (http://seqanswers.com/forums/forumdisplay.php?f=18)
-   -   Error "RG ID on SAMRecord not found in header" from Picard's MarkDuplicates.jar‏ (http://seqanswers.com/forums/showthread.php?t=9564)

cliff 02-18-2011 04:44 PM

Error "RG ID on SAMRecord not found in header" from Picard's MarkDuplicates.jar‏
 
I am trying the following command:

java -Xmx4g -jar /home/picard-tools-1.35/MarkDuplicates.jar INPUT=$dir/input.bam OUTPUT=$dir/output.bam TMP_DIR=/home/temp METRICS_FILE=PCR_duplicates REMOVE_DUPLICATES=true ASSUME_SORTED=true VALIDATION_STRINGENCY=LENIENT

And getting lots of errors like "RG ID on SAMRecord not found in header".

The command is still running and I haven't seen output.bam coming out... From what I saw previously, it is most likely that I will end up getting no output when the command stops.

I used BWA for reads mapping, then SAMtools for sorting reads and merge reads from different lanes or different Illumina flowcells. For each bam file, I used the corresponding flowcell name as the prefix, e.g., FC124_1_sorted.bam.

Since I merged different flowcells from the same Illumina machine, I used "samtools merge -r". Based on SAMtools manual, "merge -r" can attach an RG tag to each alignment. The tag value is inferred from file names.

I am running commands in a cluster node which has 4GB RAM and 10Tb space and I assume the memory and space should be enough for picard.

Anyone knows how to fix the error "RG ID on SAMRecord not found in header"?


Thanks

-C

gaffa 02-19-2011 08:06 AM

Is there a line describing the read groups in your SAM header? (you can view the header on a bam file using e.g. samtools view -H)

alexbmp 11-10-2011 01:00 AM

I used picard AddOrReplaceReadGroups.jar on this problem.

I think samtools also can add (actually, replace: by samtools reheader) read groups to your bam file.

For example:
Quote:

java -Xmx4g -jar /home/lab/bin/picard/AddOrReplaceReadGroups.jar VALIDATION_STRINGENCY=LENIENT VERBOSITY=ERROR SO=unsorted INPUT=sample.cat.pic OUTPUT=sample.rga.bam RGPL=illumina RGLB=sm_1110 RGPU=barcode_1110 RGSM=sample_1110
Any other ideas?

Chipper 11-10-2011 02:47 AM

Use picard's MergeSam rather than samtools merge if the files have different headers.

alexbmp 11-10-2011 03:27 AM

@ Chipper

Yes, I also agree using picard MergeSam instead of samtools merge.

However I found a case that picard MergeSam doesn't take many inputs (500 inputs) when options are default.

So in that case I used samtools merge (the files were split; they ought to have the same header).


All times are GMT -8. The time now is 11:00 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.