Question: picard CalculateHSMetrics error
gravatar for Kssr
4.5 years ago by
Kssr110 wrote:

I'm trying to run picard CalculateHSMetrics step on a bam file to get hybrid selection metrics and per target coverage across a list of target and probe intervals.

I get the below error,

SequenceUtil$SequenceListsDifferException: Sequence dictionaries are not the same size (87, 85)

I checked if the @SQ header lines in my bam and interval files are the same, and they are. Is it essential that the header of bam and interval files match exactly?

metrics picard • 1.9k views
ADD COMMENTlink modified 4.5 years ago by mastal5112.0k • written 4.5 years ago by Kssr110
gravatar for mastal511
4.5 years ago by
mastal5112.0k wrote:

This is complaining about sequence dictionaries. Do your files have .dict indexes? It looks like that is what it is complaining about.

ADD COMMENTlink written 4.5 years ago by mastal5112.0k

I understand the error is about sequence dictionaries. Yes, my files -bam and target & probe intervals have .dict indexes. I'm not sure what the error message implies, @SQ lines in the header of files are same. Only difference is my bam has GO:none in the @HD line and @PG lines(which I think is okay). Is the error due to having GO:none in the @HD line of my bam? Appreciate any help.

ADD REPLYlink written 4.5 years ago by Kssr110
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 771 users visited in the last hour