Cannot run HaplotypeCaller in GATK
2
0
Entering edit mode
4.8 years ago
ImBioMario • 0

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

https://pastebin.com/fdSyyfPW

https://pastebin.com/dwK9h14F

https://pastebin.com/QGwkd05X

P.S. I mapped whole human genome with 3 different tools. Raw NGS reads were pair-end reads.

Thanks in advance!

snp gatk haplotypecaller read groups • 1.2k views
ADD COMMENT
1
Entering edit mode
4.8 years ago
Brice Sarver ★ 3.8k

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.

ADD COMMENT
1
Entering edit mode
4.8 years ago
tothepoint ▴ 940
java -jar /pathtopicard/AddOrReplaceReadGroups.jar I= input.bam O=output.bam RGID=4 RGLB=libx RGPL=illumina RGSM=namesample

Check this

ADD COMMENT

Login before adding your answer.

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