Cannot run HaplotypeCaller in GATK
2
0
Entering edit mode
20 months 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 • 493 views
ADD COMMENT
1
Entering edit mode
20 months ago
Brice Sarver ★ 3.7k

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
20 months ago
tothepoint ▴ 600
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: 2767 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