error using somatic indel detector

Earendil

hello every one please i got this error while running this command

elendin@elendin-HP-Pavilion-dv6700-Notebook-PC:~/analysis of rnaseq bamfiles$ java -jar GenomeAnalysisTK.jar -R VitisVinifera.fasta -T SomaticIndelDetector -o indels.vcf -verbose indels.txt -I:normal wt.bam -I:tumor newmut.bam

MESSAGE: SAM/BAM file SAMFileReader{/home/elendin/analysis of rnaseq bamfiles/newmut.bam} is malformed: Read HWI-ST132_0461:3:2201:1211:140854#GTCCTA is either missing the read group or its read group is not defined in the BAM header, both of which are required by the GATK. Please use http://www.broadinstitute.org/gsa/wiki/index.php/ReplaceReadGroups to fix this problem

ERROR ------------------------------------------------------------------------------------------

how can i add the missing read group for HWI-ST132_0461:3:2201:1211:140854#GTCCTA or define it in the header. thanks

ebanks

This is addressed in the documentation on this forum. The error you posted specifically asks you not to post this problem to the forum because it's something you can easily fix just by reading the docs that our team has spent a lot of time preparing.

TomAmes

FYI, that link redirects to the main GATK site, and there does not seem to be any way to find the wiki entry that contains the documentation.

Geraldine_VdAuwera

We changed from the wiki to this website, but old versions of the code still produce old links.

FYI, what you're looking for is the first link you get if you use the search box with the terms "ReplaceReadGroups"...

http://gatkforums.broadinstitute.org/discussion/59/companion-utilities-replacereadgroups

TomAmes

Thanks, Geraldine. (FWIW, I'm not having the problem mentioned above--was just pointing out that the solution was not at the link given to the OP).