My command:
$GATK HaplotypeCaller -R $ref -L $locations -I Bamfiles.txt -ploidy 1
The user error looks like this:
A USER ERROR has occurred:
Input files reference and reads have incompatible contigs: No overlapping contigs found.
reference contigs = [All of the correct contigs are listed out]
reads contigs = []
Meaning, the reads contigs are being read-in as blank
I thought maybe one of my files was aligned incorrectly. I looked at their headers
for file in $(cat Bamfiles.txt); do samtools view -H $file; done | sort | uniq -c > Out.txt
The counts for each chromosomal header were consistent and the same numbers as I have samples. The only thing that looked different: about 1/4 of my samples (from an earlier study) had VN (sam version): 1.3 while the rest had VN: 1.5. Could this be enough to cause GATK to reject my samples and if so is there an easy way to convert between samfile 1.3 and 1.5 or should I re-align them?
Hello Ace!
We believe that this post does not fit the main topic of this site.
Error resolved
For this reason we have closed your question. This allows us to keep the site focused on the topics that the community can help with.
If you disagree please tell us why in a reply below, we'll be happy to talk about it.
Cheers!
No need to close the post. That is an action used by moderators for inappropriate, incomplete, off-topic posts.
You can
accept
your own answer below to provide closure to this thread.