Hello! I am trying to conduct a SNP calling with HaplotypeCaller but i got an error message that theres something wrong with my *.BAM files. After checking what is wrong i got a message that my *BAM files do not have read groups. My question is how do i deal with that error? Can i add 'dummy read groups' so i could force running HaplotypeCaller? How the results would change if i would add read groups correctly? I dont quite understand the required parameters for AddOrReplaceReadGroups tool.
Here are headers of my *.BAM files (used: samtools viev -H ...):
P.S. I mapped whole human genome with 3 different tools. Raw NGS reads were pair-end reads.
Thanks in advance!
•
written
by
ImBioMario •
2 answers
AddOrReplaceReadGroups in the Picard manual should tell you everything that you need to know. There's also a link to the GATK's definition list for more detailed explanations. You should also consider searching on Biostars - this has come up many times before.
java -jar /pathtopicard/AddOrReplaceReadGroups.jar I= input.bam O=output.bam RGID=4 RGLB=libx RGPL=illumina RGSM=namesample
Check this
Log in to answer this question.
-
Error with Gatk Haplotypecaller
written by Sowmya Pulapet •
-
Error while using GATK Haplotypecaller
written by Sowmya Pulapet •
-
Malformed BAM while running GATK
written by shivangi.agarwal800
-
Using GATK AddOrReplaceReadGroups
written by vctrm67 •
-
Error doing VCF files with GATK
written by theo.stefan.1 •
-
Adding Read Groups to Merged Bams
written by dthorbur
-
GATK: SAM file doesn't have any read groups defined in the header.
written by mahnoornaseer97 •
-
a VERY BASIC question about "add or replace groups" for BAM files
written by CrazyB
-
Redundant @Sq Lines In Bam File
written by deepthithomaskannan
-
Picard Addreplacereadgroup Error: Rg Id Not Found In Header
written by William
