Running HaplotypeCaller in GATK

I have an issue while running HaplotypeCaller in GATK. This is the entire protocol of what i did:

  1. I have two done whole exome sequencing and i have two sequences (paired ends so forward and reverse). I indexed the human genome using bwa and also aligned the two sequences to it using bwa.
  2. I then sorted the resulting sam file using picard sortsam.
  3. I then wanted to run the haplotypeCaller of gatk on the output bam file to get the vcf file of the available SNVs. But when i run it i get an error: java.lang.IllegalArgumentException: samples cannot be empty.
  4. So i ran picards validatesamfile and this is what i got:

HISTOGRAM java.lang.String

Error Type Count
ERROR:MISSING_READ_GROUP 1
WARNING:RECORD_MISSING_READ_GROUP 131996644

What does this error mean and how do i solve this? Is there any mistake i am making in my protocol that i should i change in order to not get this error in the future?

Answers

Sign In or Register to comment.