The current GATK version is 3.6-0
Examples: Monday, today, last week, Mar 26, 3/26/04

Howdy, Stranger!

It looks like you're new here. If you want to get involved, click one of these buttons!

Powered by Vanilla. Made with Bootstrap.

error using somatic indel detector

EarendilEarendil Member Posts: 2
edited January 2013 in Ask the GATK team

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

Post edited by Geraldine_VdAuwera on

Answers

  • ebanksebanks Broad InstituteMember, Administrator, Broadie, Moderator, Dev Posts: 698 admin

    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.

    Eric Banks, PhD -- Director, Data Sciences and Data Engineering, Broad Institute of Harvard and MIT

  • TomAmesTomAmes Member Posts: 6

    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_VdAuweraGeraldine_VdAuwera Administrator, Dev Posts: 10,677 admin

    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

    Geraldine Van der Auwera, PhD

  • TomAmesTomAmes Member Posts: 6

    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).

Sign In or Register to comment.